Search NASA⌕ Search

SEARCH · Search NASA

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Trends in runaway electron plateau partial recombination by massive H 2 or D 2 injection in DIII-D and JET and first extrapolations to ITER and SPARC

Abstract Experimental trends in thermal plasma partial recombination resulting from massive D 2 injection into high-Z (Ar) containing runaway electron (RE) plateaus in DIII-D and JET are studied for the purpose of achieving sufficiently low electron density ( n e ≈ 10 18 m − 3 ) to increase RE final loss MHD levels. In both DIII-D and JET, thermal electron density n e is found to drop by ∼100 × when the thermal plasma partially recombines, with a minimum at a vacuum vessel-averaged D 2 density in the range 10 20 − 10 21 m − 3 . RE effective resistivity also drops after partial recombination, indicating expulsion of the Ar content. The n e level after partial recombination is found to increase as RE current is increased. The amount of initial Ar in the RE plateau is not observed to have a strong effect on partial recombination. Partial recombination timescales of order 5 ms in DIII-D and 15 ms in JET are observed. These basic trends and timescales are matched with a 1D diffusion model, which is then used to extrapolate to ITER and SPARC tokamaks. Within the approximations of this model, it is predicted that ITER will be able to achieve sufficiently low n e values on time scales faster than expected RE plateau vertical drift timescales (of order 100 ms), provided sufficient D 2 or H 2 is injected. In SPARC, it is predicted that achieving significant n e recombination will be challenging, due to the very high RE current density. In both ITER and SPARC, it is predicted that achieving low n e will be easier with Ar as a background impurity (rather than Ne).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Direct link between disorder and magnetoresistance in topological semimetals

The extent to which disorder influences the properties of topological semimetals is relevant to the understanding of topological states and their use in practical applications. Using molecular beam epitaxy, we achieve systematic control of point defect concentrations in the prototypical Dirac semimetal Cd 3 As 2 to gain insight into the role of disorder on electron transport behavior. Using the guiding center diffusion model for linear magnetoresistance, we extract point defect densities as a function of deposition conditions. We find that reducing cadmium defect concentrations by an order of magnitude results in an 2x increase in the magnetoresistance from 450% to 900%. This finding yields important information in the quest to identify the origin of linear magnetoresistance in a wider range of materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2011

Onset of nondiffusive phonon transport in transient thermal grating decay

The relaxation of a spatially sinusoidal temperature perturbation in a dielectric crystal at a temperature comparable to or higher than the Debye temperature is investigated theoretically. We assume that most phonons contributing to the specific heat have a mean free path (MFP) much shorter than the thermal transport distance and can be described by the thermal diffusion model. Low-frequency phonons that may have MFP comparable to or longer than the grating period are described by the Boltzmann transport equation. These low-frequency phonons are assumed to interact with the thermal reservoir of high-frequency phonons but not with each other. Within the single mode relaxation time approximation, an analytical expression for the thermal grating relaxation rate is obtained. We show that the contribution of “ballistic” phonons with long MFP to the effective thermal conductivity governing the grating decay is suppressed compared to their contribution to thermal transport at long distances. In conclusion, the reduction in the effective thermal conductivity in Si at room temperature is found to be significant at grating periods as large as 10 μm.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Localization landscape theory of disorder in semiconductors. III. Application to carrier transport and recombination in light emitting diodes

In this paper we introduce a novel method to account for quantum disorder effects into the classical drift-diffusion model of semiconductor transport through the localization landscape theory. Quantum confinement and quantum tunneling in the disordered system change dramatically the energy barriers acting on the perpendicular transport of heterostructures. In addition, they lead to percolative transport through paths of minimal energy in the two-dimensional (2D) landscape of disordered energies of multiple 2D quantum wells. This model solves the carrier dynamics with quantum effects self-consistently and provides a computationally much faster solver when compared with the Schrödinger equation resolution. The theory also provides a good approximation to the density of states for the disordered system over the full range of energies required to account for transport at room temperature. The current-voltage characteristics modeled by three-dimensional simulation of a full nitride-based light emitting diode (LED) structure with compositional material fluctuations closely match the experimental behavior of high-quality blue LEDs. The model allows also a fine analysis of the quantum effects involved in carrier transport through such complex heterostructures. Finally, details of carrier population and recombination in the different quantum wells are given.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Ion Thermal Decoupling and Species Separation in Shock-Driven Implosions

Anomalous reduction of the fusion yields by 50% and anomalous scaling of the burn-averaged ion temperatures with the ion-species fraction has been observed for the first time in DHe 3 -filled shock-driven inertial confinement fusion implosions. Two ion kinetic mechanisms are used to explain the anomalous observations: thermal decoupling of the D and He 3 populations and diffusive species separation. The observed insensitivity of ion temperature to a varying deuterium fraction is shown to be a signature of ion thermal decoupling in shock-heated plasmas. The burn-averaged deuterium fraction calculated from the experimental data demonstrates a reduction in the average core deuterium density, as predicted by simulations that use a diffusion model. Accounting for each of these effects in simulations reproduces the observed yield trends.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Beam loading compensation of traveling wave linacs through the time dependence of the rf drive

Beam loading in traveling-wave linear accelerating structures leads to unacceptable spread of particle energies across an extended train of bunched particles due to beam-induced field and dispersion. Methods for modulating the rf power driving linacs are effective at reducing energy spread, but for general linacs do not have a clear analytic foundation. We report here methods for calculating how to modulate the rf drive in arbitrarily nonuniform traveling-wave linacs within the convective-transport (power-diffusion) model that results in no additional energy spread due to beam loading (but not dispersion). Varying group velocity, loss factor, and cell quality factor within a structure, and nonzero particle velocity, are handled.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Latent diffusion can map beam loss to two-dimensional phase-space projections

Beam loss monitors (BLMs) and beam current monitors (BCMs) are ubiquitous at particle accelerators around the world. These simple devices provide noninvasive high-level beam measurements but give no insight into the detailed 6D (𝑥,𝑦,𝑧,𝑝 𝑥 ,𝑝 𝑦 ,𝑝 𝑧 ) beam phase-space distributions or dynamics. We show that generative conditional latent diffusion models can learn intricate patterns to solve the extreme inverse problem of mapping waveforms of tens of BLMs or BCMs along an accelerator to detailed 2D projections of a charged particle beam’s 6D phase-space density. This transformational method can be used at any particle accelerator to transform simple noninvasive devices into detailed beam phase-space diagnostics. We demonstrate this concept via multiparticle simulations of the high-intensity beam in the kilometer-long Los Alamos Neutron Science Center linear proton accelerator.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Localization of fluorescent targets in deep tissue with expanded beam illumination for studies of cancer and the brain

Imaging fluorescence through millimeters or centimeters of tissue has important in vivo applications, such as guiding surgery and studying the brain. Often, the important information is the location of one of more optical reporters, rather than the specifics of the local geometry, motivating the need for a localization method that provides this information. We present an optimization approach based on a diffusion model for the fast localization of fluorescent inhomogeneities in deep tissue with expanded beam illumination that simplifies the experiment and the reconstruction. We show that the position of a fluorescent inhomogeneity can be estimated while assuming homogeneous tissue parameters and without having to model the excitation profile, reducing the computational burden and improving the utility of the method. We perform two experiments as a demonstration. First, a tumor in a mouse is localized using a near infrared folate-targeted fluorescent agent (OTL38). This result shows that localization can quickly provide tumor depth information, which could reduce damage to healthy tissue during fluorescence-guided surgery. Second, another near infrared fluorescent agent (ATTO647N) is injected into the brain of a rat, and localized through the intact skull and surface tissue. This result will enable studies of protein aggregation and neuron signaling.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Numerical Improvements in Magnetohydrodynamic, Pulsed Power Simulations of Near-Target Plasmas

Magnetohydrodynamic (MHD) simulations of pulsed power experiments frequently result in unphysical runaway heating when extended to low values of mass density. Traditionally, this has been addressed by the use of conductivity floors below which the plasma is given an arbitrarily low conductivity value. Here, a low-density treatment is presented that allows for low-density material to carry current while maintaining a stable temperature. This treatment is implemented in the Ares multiphysics code. Numerical modifications to the conventional vacuum treatment include an energy-conserving density floor, a modified averaging procedure to determine the thermal conductivity at the edge of adjacent zones, and a modified averaging procedure for determining the thermal conductivity in mixed zones that contain material from differing regions. Additionally, anomalous resistivity and Bohm diffusion models are implemented as simplified models for microphysics-induced enhancement of collisional transport. Importantly, the advantage of these various improvements are illustrated in a simple 1-D pulsed power target where the combined changes result in stable temperatures within the lower density regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Observed volatilization behavior of silicon carbide in flowing hydrogen above 2000 K

In this work, the intrinsic compatibility of silicon carbide (SiC) and hydrogen (H 2 ) at high temperatures (2000-2473 K) and pressure near one atmosphere was evaluated through a combination of thermodynamic calculations and hot hydrogen exposure testing. Thermodynamic calculations predict the decomposition of SiC in a hydrogen environment to form free silicon (Si) and free carbon (C). Free Si is predicted to vaporize from the surface as a volatile species, while free C may interact with H 2 to form the hydrocarbons CH 4 (T < 2100 K) or C 2 H 2 (T > 2100 K). Coupons of high purity chemical vapor deposition (CVD) β-SiC were exposed to slowly flowing hydrogen at temperatures ranging between 2000 and 2473 K. SiC experienced active attack as the result of H 2 exposure, exhibiting linear weight loss kinetics and an Arrhenius dependence of weight loss on exposure temperature. The linear volatilization constant was experimentally evaluated to correspond with an activation energy of 370 ± 18 kJ/mol. Due to the dependence of observed corrosion rates on gas velocity, corrosion of SiC in flowing H 2 was determined to be governed by external mass transfer of volatile Si species through the boundary layer. Experimentally derived mass losses were in good agreement with mass losses predicted by a boundary layer limited gas diffusion model.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Oxygen isotope thermometry using quartz inclusions in garnet

Oxygen isotope ratios of quartz inclusions (QI) within garnet from granulite and amphibolite facies gneisses in the Adirondack Mountains, NY were analysed and used to determine metamorphic temperatures. Primary QI for eight of 12 samples have δ 18 O values significantly lower than matrix quartz (MQ). The primary QI retain δ 18 O values representative of thermal conditions during garnet crystallization, whereas the δ 18 O values of MQ were raised by diffusive exchange with other matrix minerals (e.g. mica and feldspar) during cooling. The δ 18 O differences between QI and MQ show that garnet (a mineral with slow diffusion of oxygen) can armour QI from isotopic exchange with surrounding matrix, even during slow cooling. These differences between δ 18 O in MQ and QI can further be used to test cooling rates by Fast Grain Boundary diffusion modelling. Criteria for identifying QI that preserve primary compositions and are suitable for thermometry were developed based on comparative tests. Relations between δ 18 O and inclusion size, distance of inclusion to host–garnet rim, core–rim zonation of individual inclusions, and presence or absence of petrological features (healed cracks in QI, inclusions in contact with garnet cracks lined by secondary minerals, and secondary minerals along the inclusion grain boundary) were investigated. In this study, 61% of QI preserve primary δ 18 O and 39% were associated with features that were linked to reset δ 18 O values. If δ 18 O in garnet is homogeneous and inclusions are removed, laser-fluorination δ 18 O values of bulk garnet are more precise, more accurate, and best for thermometry. Intragrain δ 18 O(Grt) profiles measured in situ by ion microprobe show no δ 18 O zonation. Almandine–rich garnet (Alm 60–75 ) from each sample was measured by laser-fluorination mass-spectrometry (LF-MS) for δ 18 O and compared with ion microprobe measurements of δ 18 O in QI for thermometry. The Δ 18 OO(Qz–Grt) values for Adirondack samples range from 2.66 to 3.24‰, corresponding to temperatures of 640–740 °C (A[Qz–Alm] = 2.71). Out of 12 samples that were used for thermometry, nine are consistent with previous estimates of peak temperature (625–800 °C) based on petrological and carbon–isotope thermometry for regional granulite and upper amphibolite facies metamorphism. The three samples that disagree with independent thermometry for peak metamorphism are from the anorthosite–mangerite–charnockite–granite suite in the central Adirondacks and yield temperatures of 640–665 °C, ~100 °C lower than previous estimates. Finally, these low temperatures could be interpreted as thermal conditions during late (post-peak) crystallization of garnet on the retrograde path.

58 GEOSCIENCES↗
DOE OSTI2021

Application of Scalar Filtered Density Function to Turbulent Flows Under Supercritical Condition

Abstract The scalar filtered density function (FDF) methodology is extended and employed for large eddy simulation (LES) of turbulent flows under supercritical condition. To describe real fluid behavior, the extended methodology incorporates the generalized heat and mass diffusion models along with real fluid thermodynamic relations which are derived using the cubic Peng–Robinson equation of state. These models are implemented within the stochastic differential equations comprising the scalar FDF transport. Simulations are conducted of a temporally developing mixing layer under supercritical condition and the results are assessed by comparing with data generated by direct numerical simulation (DNS) of the same layer. The consistency of the proposed FDF methodology is assessed. The LES-FDF predictions are shown to agree favorably with the DNS data and exhibit several key features pertaining to supercritical turbulent flows.

Energy & Fuels↗
DOE OSTI2019

The role of permafrost on the morphology of an MIS 3 moraine from the southern Laurentide Ice Sheet

The Laurentide Ice Sheet (LIS) represents the single largest contributor to global sea level during the late Pleistocene glacial-interglacial cycle. Fluctuations of global sea level prior to the Last Glacial Maximum (LGM; 26–19 ka) are well documented, but the terrestrial extent of the LIS prior to the LGM is uncertain. In central Wisconsin, United States, three north-south–trending end moraines of the LIS are preserved. The two eastward moraines are constrained to the LGM through radiocarbon and surface exposure dating methods, while the westward moraine is constrained only by relative dating techniques and is estimated to be up to >125 ka. Here, we report surface exposure ages from this westward moraine (the Arnott moraine) and demonstrate that the LIS reached an equivalent extent to its LGM position at ca. 35 ka, during Marine Isotope Stage (MIS) 3. These new age constraints for the Arnott moraine are significantly younger than prior estimates, which were based on sediment weathering and low-relief moraine morphology. To address this dichotomy between the relatively young exposure dates and the low relief, we applied a landscape diffusion model to the Arnott moraine and found that during permafrost conditions, intense erosion can act to efficiently and quickly smooth the moraine surface. Here, this study connects the LIS as a potential source for sea-level lowering during MIS 3 while simultaneously demonstrating the inherent uncertainty associated with relative dating techniques if nonlinear erosional conditions are not accounted for in periglacial terrains.

58 GEOSCIENCES↗
DOE OSTI2024

Rehydrated glass embayments record the cooling of a Yellowstone ignimbrite

Abstract Hydration fronts penetrate 50–135 μm into glassy rhyolite embayments hosted in quartz crystals from the Mesa Falls Tuff in the Yellowstone Plateau volcanic field. The hydration fronts occur as steep enrichments that reach 2.4 ± 0.6 wt% H2O at the embayment opening, representing much higher values than interior concentrations of 0.9 ± 0.2 wt% H2O. Molecular water accounts for most of the water enrichment. Water speciation indicates the hydration fronts comprise absorbed meteoric water that modified the original magmatic composition of the rhyolitic glass. We used finite difference diffusion models to demonstrate that glass rehydration was likely produced over a few decades as the ignimbrite cooled. Such temperatures and time scales are consistent with rare firsthand observations of decadal hydrothermal systems associated with cooling ignimbrites at Mount Pinatubo (Philippines) and the Valley of Ten Thousand Smokes (Alaska).

Geology↗
DOE OSTI2017

Mathematical simulation of heat and mass transfer in a cylindrical channel versus the tangential momentum accommodation coefficient

The process of heat and mass transfer in a long cylindrical channel has been considered in terms of the mirror–diffuse model of the Maxwell boundary condition. The Hazen–Williams equation is used as a basic equation of the process kinetics. A constant temperature gradient is maintained in the channel. The heat and mass fluxes through the cross section of the channel versus the tangential momentum accommodation coefficient have been calculated in a wide range of the Knudsen number. The heat flux profiles have been constructed. A comparison with relevant published data has been carried out.

42 ENGINEERING↗
DOE OSTI2023

High-Level Synthesis of Irregular Applications: A Case Study on Influence Maximization

The Influence Maximization problem is the problem of identifying a small cohort of actors from a broader population that, when initially activated in a diffusion process, are expected to result in a large number of activations in the population. While the problem is known to be NP-hard, several approximation algorithms have been devised by leveraging its submodular structure. While these algorithms are theoretically efficient, they are computationally very expensive in practice. This work advances the current state-of-the-art parallelization scheme for the IMM algorithm by devising the adoption of custom hardware accelerators implemented on FPGAs by leveraging High Level Synthesis from OpenCL. We study the performance of our proposed approach by exploring optimizations tailored at improving the parallel efficiency of the accelerators and highlight their effects and limitations in accelerating complex graph analytic applications. Our experimental evaluation shows that FPGA acceleration can improve the performance of the LT diffusion model up to 1.72x for the entire application and up to 2.90x for its most important kernel with respect to a CPU only parallel execution. The FPGA acceleration of the LT model shows also a 1.54x reduction in energy consumption when compared to a parallel CPU only run.

Neff, Reece W.↗
DOE OSTI2016

Tritium migration to the surfaces of Type 316 stainless steel; aluminum 6061; and oxygen-free, high-conductivity copper

The migration of tritium to the surfaces of aluminum 6061, oxygen-free, high-conductivity copper (OFHC), and stainless-steel 316 from the bulk metal was studied using low-pressure Tonks–Langmuir argon plasma. The plasma is shown to be effective at removing tritium from metal surfaces in a controlled manner. Tritium is removed in decreasing quantities with successive plasma exposures, which suggests a depletion of the surface and near-surface tritium inventories. A diffusion model was developed to predict tritium migration from the bulk and its accumulation in the water layers present on the metal surface. The model reproduces the rate of tritium re-growth on the surface for all three metals and can be used to calculate the triton solubility in the water layers present on metal surfaces. The ratio of surface-to-bulk solubilities at the water-layer/bulk-metal interface uniquely determines the concentration ratio between these two media. Removing the tritium-rich water layers induces tritium to migrate from the bulk to the surface. Furthermore, this process is driven by a concentration gradient that develops in the bulk because of the perturbation on the surface.

36 MATERIALS SCIENCE↗