Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Analysis of Anisotropic Nonlocal Diffusion Models: Well-posedness of Fractional Problems for Anomalous Transport

We analyze the well-posedness of an anisotropic, nonlocal diffusion equation. Establishing an equivalence between weighted and unweighted anisotropic nonlocal diffusion operators in the vein of unified nonlocal vector calculus, we apply our analysis to a class of fractional-order operators and present rigorous estimates for the solution of the corresponding anisotropic anomalous diffusion equation. Furthermore, we extend our analysis to the anisotropic diffusion-advection equation and prove well-posedness for fractional orders s ∊ [0.5, 1). We also present an application of the advection-diffusion equation to anomalous transport of solutes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Parameters estimation using the first passage times method in a jump-diffusion model

The main purposes of this paper are two contributions: (1) it presents a new method, which is the first passage time (FPT method) generalized for all passage times (GPT method), in order to estimate the parameters of stochastic Jump-Diffusion process. (2) it compares in a time series model, share price of gold, the empirical results of the estimation and forecasts obtained with the GPT method and those obtained by the moments method and the FPT method applied to the Merton Jump-Diffusion (MJD) model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A baseline for ensemble-based, time-resolved inflow reconstruction for a single turbine using large-eddy simulations and latent diffusion models

We are interested in reconstructing winds flowing through a turbine on a second-by-second basis over a 10 min window. Previously, we developed a machine learning algorithm that takes in a snapshot of wind speed measurements and generates ensembles of three-dimensional wind field estimates. Here, we use these estimates as initial conditions in large-eddy simulations and reconstruct atmospheric and turbine response dynamic quantities in a synthetic field campaign. In doing so, we establish a baseline for model validation that future time-aware data assimilation techniques will be compared to. In turbine-free case studies, ground truth wind speeds consistently fall within our estimated wind speed distribution for the first 100 s after the simulation start. In simulations with turbines, the wind estimates show a small bias of 0.10 m s 1 and good correlation of 0.80 during the first 100 s. During this window, our estimates of the Blade 1 bending moment and generator power typically span the ground truth, with the estimate of the former performing better overall. In summary, this approach shows promise as a stand-alone technique for reconstructing real-world inflow and turbine dynamics in 1-2 min windows and as a foundation for future time-aware data assimilation techniques.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A generative artificial intelligence framework based on a molecular diffusion model for the design of metal-organic frameworks for carbon capture

Metal-organic frameworks (MOFs) exhibit great promise for CO 2 capture. However, finding the best performing materials poses computational and experimental grand challenges in view of the vast chemical space of potential building blocks. Here, we introduce GHP-MOFassemble, a generative artificial intelligence (AI), high performance framework for the rational and accelerated design of MOFs with high CO 2 adsorption capacity and synthesizable linkers. GHP-MOFassemble generates novel linkers, assembled with one of three pre-selected metal nodes (Cu paddlewheel, Zn paddlewheel, Zn tetramer) into MOFs in a primitive cubic topology. GHP-MOFassemble screens and validates AI-generated MOFs for uniqueness, synthesizability, structural validity, uses molecular dynamics simulations to study their stability and chemical consistency, and crystal graph neural networks and Grand Canonical Monte Carlo simulations to quantify their CO 2 adsorption capacities. We present the top six AI-generated MOFs with CO 2 capacities greater than 2m mol g -1 , i.e., higher than 96.9% of structures in the hypothetical MOF dataset.

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

Moisture ingress in commercial steel drums: Water content determination, diffusion modelling and predicted permeation rates

Commercial steel drums underpin the global economy, playing a pivotal role in the storage and transportation of critical materials. Transported and stored materials, such as food, chemical and nuclear waste, can be sensitive to ambient conditions, particularly moisture that can enhance negative effects such as corrosion and material degradation. Although international standards and regulations are in place for the qualification of steel drums, there are no current testing requirements, established limits or boundaries for the permeation of moisture into the drums during transportation or storage. This work aims to provide insights into the moisture ingress over time into properly sealed steel drums and provides estimated moisture ingress rates over time through extrapolation. Water vapour transmission rate (WVTR) measurements through the gasket material at 10–40°C were 0.11–2.1 g/m 2 /day resulting in a permeation activation energy of 30.2 kJ/mol. Water sorption measurements and Karl Fischer titration (KFT) on ethylene propylene diene monomer (EPDM) gasket material revealed a decrease in equilibrium moisture saturation with increasing temperature. KFT measurements also revealed the presence of moisture within the adhesive and drum wall after exposure to ambient conditions. KFT and Fourier transform infrared spectroscopy (FTIR) show that moisture will desorb from the EPDM and drum wall after exposure to desiccating conditions, although a minimal amount of moisture will remain present. When sealed to the manufacturer's recommendations, the steel drums are effective in minimizing moisture ingress. In conclusion, in sealed empty drums, moisture ingress rates for 19-L drums were 0.4–1.5 mg/day at 25°C 15% relative humidity (RH) and increased to 7.1–8.8 mg/day at 40°C 90% RH, and moisture ingress rates for 210-L drums were 2.5 and 3.5 mg/day at field deployment conditions of 15.5°C 51.5% RH and 23°C 40% RH, respectively.

42 ENGINEERING↗
DOE OSTI2014

SU-E-T-196: Heat Diffusion Modeling for Digital Holographic Interferometry Dosimetry

Purpose: We have previously demonstrated that with Digital Holographic Interferometry (DHI) 2D spatial calorimetric measurements of high dose rate radiation sources can be obtained. The impact of heat transfer must be considered when undertaking any form of calorimetric measurement, as the radiation induced temperature distributions are subject to degradation due to heat diffusion. Unaccounted for, this limits the accuracy of the approach especially for long delivery times. Methods: 3D modelling of the heat diffusion in water was undertaken, and two different approaches developed to account for this effect. The mathematical framework to describe heat diffusion in 3D was applied, with the differential equations solved numerically using an implicit method. The first approach involved the comparison of the DHI measurements to an independent dose model of the source. The model was forward modeled to account for the heat diffusion during irradiation, allowing a direct comparison to validate the measured results. The second approach involved the correction of the measured data directly, by comparing the temperature distribution of two instances and subtracting the effects of heat diffusion of the first distribution from the second instance. This required the use of the Abel transform to approximate the 3D dose distribution from the 2D DHI results, thus limiting the approach to radiation applications possessing cylindrical symmetry. Results: The first approach resulted in higher accuracy and was more straightforward, but has a major limitation in that the measured results are only able to be utilized in comparison with an independent dose model. The applicability of the second approach is affected by noise in the measurement data and introduces higher uncertainties, but results in higher usability of the final data. Conclusion: Both approaches were implemented, and if used in conjunction would provide the most utility for the interpretation and use of DHI measurements.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2015

A Process for Modelling Diffuse Scattering from Disordered Molecular Crystals, Illustrated by Application to Monoclinic 9-Chloro-10-methylanthracene

Diffuse scattering from a crystal contains valuable information about the two-body correlations (related to the nanoscale order) in the material. Despite years of development, the detailed analysis of single crystal diffuse scattering (SCDS) has yet to become part of the everyday toolbox of the structural scientist. Recent decades have seen the pair distribution function approach to diffuse scattering (in fact, total scattering) from powders become a relatively routine tool. However, analysing the detailed, complex, and often highly anisotropic three-dimensional distribution of SCDS remains valuable yet rare because there is no routine method for undertaking the analysis. At present, analysis requires significant investment of time to develop specialist expertise, which means that the analysis of diffuse scattering, which has much to offer, is not incorporated thorough studies of many compounds even though it has the potential to be a very useful adjunct to existing techniques. This article endeavours to outline in some detail how the diffuse scattering from a molecular crystal can be modelled relatively quickly and largely using existing software tools. It is hoped this will provide a template for other studies. To enable this, the entire simulation is included as deposited material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Discovering a reaction–diffusion model for Alzheimer’s disease by combining PINNs with symbolic regression

Misfolded tau proteins play a critical role in the progression and pathology of Alzheimer's disease. Recent studies suggest that the spatio-temporal pattern of misfolded tau follows a reaction-diffusion type equation. However, the precise mathematical model and parameters that characterize the progression of misfolded protein across the brain remain incompletely understood. Here, we use deep learning and artificial intelligence to discover a mathematical model for the progression of Alzheimer's disease using longitudinal tau positron emission tomography from the Alzheimer's Disease Neuroimaging Initiative database. Specifically, we integrate physics informed neural networks (PINNs) and symbolic regression to discover a reaction-diffusion type partial differential equation for tau protein misfolding and spreading. First, we demonstrate the potential of our model and parameter discovery on synthetic data. Then, we apply our method to discover the best model and parameters to explain tau imaging data from 46 individuals who are likely to develop Alzheimer's disease and 30 healthy controls. Our symbolic regression discovers different misfolding models f(c) for two groups, with a faster misfolding for the Alzheimer's group, f(c) = 0.23c 3 – 1.34c 2 + 1.11c, than for the healthy control group, f(c) = –c 3 + 0.62c 2 + 0.39c. Our results suggest that PINNs, supplemented by symbolic regression, can discover a reaction-diffusion type model to explain misfolded tau protein concentrations in Alzheimer's disease. Furthermore, we expect our study to be the starting point for a more holistic analysis to provide image-based technologies for early diagnosis, and ideally early treatment of neurodegeneration in Alzheimer's disease and possibly other misfolding-protein based neurodegenerative disorders.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Reaction–drift–diffusion models from master equations: application to material defects

We present a general method to produce well-conditioned continuum reaction–drift–diffusion equations directly from master equations on a discrete, periodic state space. We assume the underlying data to be kinetic Monte Carlo models (i.e. continuous-time Markov chains) produced from atomic sampling of point defects in locally periodic environments, such as perfect lattices, ordered surface structures or dislocation cores, possibly under the influence of a slowly varying external field. Our approach also applies to any discrete, periodic Markov chain. Here, the analysis identifies a previously omitted non-equilibrium drift term, present even in the absence of external forces, which can compete in magnitude with the reaction rates, thus being essential to correctly capture the kinetics. To remove fast modes which hinder time integration, we use a generalized Bloch relation to efficiently calculate the eigenspectrum of the master equation. A well conditioned continuum equation then emerges by searching for spectral gaps in the long wavelength limit, using an established kinetic clustering algorithm to define a proper reduced, Markovian state space.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Enabling fast charging of lithium-ion batteries through secondary- /dual- pore network: Part I - Analytical diffusion model

Battery performance is strongly correlated with electrode microstructural properties. Enabling fast charging of lithium-ion batteries requires improved through-plane ionic diffusion that can be achieved through, among other strategies, structured electrodes with a secondary- or dual-pore network (SPN). In this work, an analytical model investigates the impact of such an SPN on ionic diffusion with a composite electrode, considering various pore-channel geometries and comparing to standard electrodes with identical gravimetric- and volumetric-specific theoretical capacities. Relevant SPN design parameters and tortuosity coefficients are identified according to three optimization objectives that aim to balance the improved overall through-plane diffusion, thanks to the coarse aligned channels, and degraded in-plane diffusion because of the porous matrix densification required to maintain gravimetric- and volumetric-specific theoretical capacities. The model indicates that a relatively low amount of SPN is required and that electrodes with high through-plane tortuosity and low in-plane tortuosity benefit most from such architecture.

25 ENERGY STORAGE↗
DOE OSTI2014

Perona Malik anisotropic diffusion model using Peaceman Rachford scheme on digital radiographic image

In image processing, it is important to remove noise without affecting the image structure as well as preserving all the edges. Perona Malik Anisotropic Diffusion (PMAD) is a PDE-based model which is suitable for image denoising and edge detection problems. In this paper, the Peaceman Rachford scheme is applied on PMAD to remove unwanted noise as the scheme is efficient and unconditionally stable. The capability of the scheme to remove noise is evaluated on several digital radiography weld defect images computed using MATLAB R2009a. Experimental results obtained show that the Peaceman Rachford scheme improves the image quality substantially well based on the Peak Signal to Noise Ratio (PSNR). The Peaceman Rachford scheme used in solving the PMAD model successfully removes unwanted noise in digital radiographic image.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Reproducing the observed energy-dependent structure of Earth's electron radiation belts during storm recovery with an event-specific diffusion model

Here, we present dynamic simulations of energy-dependent losses in the radiation belt “slot region” and the formation of the two-belt structure for the quiet days after the 1 March storm. The simulations combine radial diffusion with a realistic scattering model, based data-driven spatially and temporally resolved whistler-mode hiss wave observations from the Van Allen Probes satellites. The simulations reproduce Van Allen Probes observations for all energies and L shells (2–6) including (a) the strong energy dependence to the radiation belt dynamics (b) an energy-dependent outer boundary to the inner zone that extends to higher L shells at lower energies and (c) an “S-shaped” energy-dependent inner boundary to the outer zone that results from the competition between diffusive radial transport and losses. We find that the characteristic energy-dependent structure of the radiation belts and slot region is dynamic and can be formed gradually in ~15 days, although the “S shape” can also be reproduced by assuming equilibrium conditions. The highest-energy electrons (E > 300 keV) of the inner region of the outer belt (L ~ 4–5) also constantly decay, demonstrating that hiss wave scattering affects the outer belt during times of extended plasmasphere. Through these simulations, we explain the full structure in energy and L shell of the belts and the slot formation by hiss scattering during storm recovery. We show the power and complexity of looking dynamically at the effects over all energies and L shells and the need for using data-driven and event-specific conditions.

58 GEOSCIENCES↗
DOE OSTI2021

A multifidelity method for a nonlocal diffusion model

Nonlocal models feature a finite length scale, referred to as the horizon, such that points separated by a distance smaller than the horizon interact with each other. Such models have proven to be useful in a variety of settings. However, due to the reduced sparsity of discretizations, they are also generally computationally more expensive compared to their local differential equation counterparts. In this work, we introduce a multifidelity Monte Carlo method that combines the high-fidelity nonlocal model of interest with surrogate models that use coarser grids and/or smaller horizons and thus have lower fidelities and lower costs. Using the multifidelity method, the overall computational cost of uncertainty quantification is reduced without compromising accuracy. It is shown for a one-dimensional nonlocal diffusion example that speedups of up to two orders of magnitude can be achieved using the multifidelity method to estimate the expectation of an output of interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Phase transitions in a reaction-diffusion model on a line with boundaries

A one-dimensional model on a line of length L is investigated, which involves particle diffusion as well as single particle annihilation. There are also creation and annihilation at the boundaries. The static and dynamical behaviors of the system are studied. It is seen that the system could exhibit a dynamical phase transition. For small drift velocities, the relaxation time does not depend on the absorption rates at the boundaries. This is the fast phase. For large velocities, the smaller of the absorption rates at boundaries enter the relaxation rate and makes it longer. This is the slow phase. Finally, the effect of a random particle creation in the bulk is also investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Secondary Pore Network Design Optimization Analytical Diffusion Model for Lithium Ion Battery

Fast ion transport along the electrode thickness is required for Lithium ion battery (LIB) to sustain continuous high-rate discharge or charge currents, which is a desired feature for the fast charging of electric vehicle batteries. Advance in LIB electrode manufacturing allows to control, to a certain extent, the electrode microstructure architecture in order to improve the macroscopic performances of these devices. For instance, several groups have experimentally tested architectures tailored specifically to increase the through plane ionic diffusion in regards for fast charge application by adding channels or groove-like large pores oriented along the electrode thickness. These voids do not replace the existing pores of the electrode microstructure but consist in a second, coarser, pore network, called the secondary pore network (SPN). The model quantifies the impact of SPN on ionic diffusion with a composite electrode through a simple analytical approach, considering various pore channel geometries and comparing to standard electrodes with identical gravimetric and identical volumetric specific theoretical capacity. The goal is to find optimal SPN design parameters according to three optimization objectives aiming to balance the improved overall through-plane, thanks to the coarse aligned channels, and the degraded in-plane diffusion, due to the porous matrix densification required to maintain gravimetric and volumetric specific theoretical capacity. The analytical model must be considered as a pre-screening tool to quickly investigate large design space, that can be refined later with more complex electrochemical CPU-expensive numerical models.

Usseglio Viretta, Francois↗
DOE OSTI2024

Score-based diffusion models for generating liquid argon time projection chamber images

For the first time, we show high-fidelity generation of Liquid Argon Time Projection Chamber (LArTPC-like) data using a generative neural network. This demonstrates that methods developed for natural images do transfer to LArTPC-produced images, which, in contrast to natural images, are globally sparse but locally dense. We present the score-based diffusion method employed. We evaluate the fidelity of the generated images using several quality metrics, including modified measures used to evaluate natural images, comparisons between high-dimensional distributions, and comparisons relevant to LArTPC experiments. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2021

Incorporating Electrochemical Halide Oxidation into Drift-Diffusion Models to Explain Performance Losses in Perovskite Solar Cells under Prolonged Reverse Bias

Partial shading of a solar module can induce a set of cells within the module to operate under reverse bias. Studies have shown that metal halide perovskite solar cells with a wide variety of compositions and contacts exhibit interesting behavior in reverse bias that includes both reversible performance loss and non-reversible degradation. In this paper, an advanced drift-diffusion approach incorporating an electrochemical term to explain the short-circuit, open circuit and fill factor losses that are experimentally measured after prolonged reverse bias is used. It is shown that holes can tunnel into the perovskite due to sharp band bending near the contact, accumulate within the bulk of the perovskite absorber, and trigger the oxidation of halides to form neutral halogens. The density of neutral halogens is much higher in reverse bias because there are hardly any electrons available to reduce the iodine. The resulting halogens act as bulk recombination centers. While the interstitial halogen density does decay when the cell is operated in forward bias, permanent degradation can occur if the iodine diffuses out of the perovskite layer. Finally, the ways in which changing parameters such as the mobile ion density or the series resistance at the contact can influence device performance and stability are discussed.

degradation↗