Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion model”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2023

Outlook on texture evolution in additively manufactured stainless steels: Prospects for hydrogen embrittlement resistance, overview of mechanical, and solidification behavior

Abstract Realizing application specific manufacture with fusion-based additive manufacturing (F-BAM) processes requires understanding of the physical phenomena that drive evolution of microstructural attributes, such as texture. Current approaches for understanding texture evolution in F-BAM are majorly considerate of the phenomena occurring only during solidification. This hinders the comprehensive understanding and control of texture during F-BAM. In this perspective article, we discuss several physical phenomena occurring during and after solidification that can determine texture in F-BAM processed stainless steels (SS). A crystal plasticity-coupled hydrogen adsorption-diffusion modeling framework is also leveraged to demonstrate the prospects of grain boundary engineering with F-BAM for enhanced hydrogen embrittlement resistance of SS. Implications of varying thermokinetics in F-BAM for solidification behavior of SS are discussed. Additionally, microstructural attributes that are key to high temperature mechanical performance of SS are highlighted. Considerations as outlined in this perspective article will enable grain boundary engineering and application specific microstructural design of SS with F-BAM. Graphical abstract

36 MATERIALS SCIENCE↗
DOE OSTI2024

High-dimensional and permutation invariant anomaly detection

Methods for anomaly detection of new physics processes are often limited to low-dimensional spaces due to the difficulty of learning high-dimensional probability densities. Particularly at the constituent level, incorporating desirable properties such as permutation invariance and variable-length inputs becomes difficult within popular density estimation methods. In this work, we introduce a permutation-invariant density estimator for particle physics data based on diffusion models, specifically designed to handle variable-length inputs. We demonstrate the efficacy of our methodology by utilizing the learned density as a permutation-invariant anomaly detection score, effectively identifying jets with low likelihood under the background-only hypothesis. To validate our density estimation method, we investigate the ratio of learned densities and compare to those obtained by a supervised classification algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Natural Tracers and Multi-Scale Assessment of Caprock Sealing Behavior: A Case Study of the Kirtland Formation, San Juan Basin

The assessment of caprocks for geologic CO{sub 2} storage is a multi-scale endeavor. Investigation of a regional caprock - the Kirtland Formation, San Juan Basin, USA - at the pore-network scale indicates high capillary sealing capacity and low permeabilities. Core and wellscale data, however, indicate a potential seal bypass system as evidenced by multiple mineralized fractures and methane gas saturations within the caprock. Our interpretation of {sup 4}He concentrations, measured at the top and bottom of the caprock, suggests low fluid fluxes through the caprock: (1) Of the total {sup 4}He produced in situ (i.e., at the locations of sampling) by uranium and thorium decay since deposition of the Kirtland Formation, a large portion still resides in the pore fluids. (2) Simple advection-only and advection-diffusion models, using the measured {sup 4}He concentrations, indicate low permeability ({approx}10-20 m{sup 2} or lower) for the thickness of the Kirtland Formation. These findings, however, do not guarantee the lack of a large-scale bypass system. The measured data, located near the boundary conditions of the models (i.e., the overlying and underlying aquifers), limit our testing of conceptual models and the sensitivity of model parameterization. Thus, we suggest approaches for future studies to better assess the presence or lack of a seal bypass system at this particular site and for other sites in general.

03 NATURAL GAS↗
DOE OSTI2011

On the Origin of Radiation Growth of HCP Crystals

The aim of the present work is to study theoretically the radiation growth (RG) of hcp-type materials with a particular focus on the effect of one-dimensionally (1-D) migrating clusters of self-interstitial atoms (SIAs), which are steadily produced in displacement cascades under neutron or heavy-ion irradiation. A reaction-diffusion model is developed for the description of RG in single hcp-type metallic crystals. The model reproduces all RG stages observed in neutron-irradiated annealed samples of pure Zr and Zr alloys, such as high strain rate at low, strain saturation at intermediate and breakaway growth at relatively high irradiation doses. In addition, it accounts for the striking observations of negative strains in prismatic directions and coexistence of vacancy- and SIA-type prismatic loops. The role of cold work in RG behavior and alignment of the vacancy-type loops along basal planes are revealed and the maximum strain rate is estimated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Partial Support for the Federal Committee for Meteorological Services and Supporting Research (Final Report)

The OFCM sponsored a session within the 16th annual George Mason University Atmospheric Transport and Diffusion Modeling Conference. The theme of the session was Identifying the Research Needs and Current/Planned Research Projects within the Operational and Research Planetary Boundary Layer Interagency Community. Subject matter experts (SME) from the NWS, U.S. Department of Energy, and the National Center for Atmospheric Research participated in the session’s panel discussion. The OFCM-sponsored special session SMEs identified the following areas as operational and research priorities/requirements: Transitioning from planetary boundary layer (PBL) schemes to large-scale eddy simulation; Separating out waves from turbulence and understanding when and how to simulate internal waves degenerating into turbulence; Treating stable boundary layers and unstable boundary layers as separate entities; Increasing the capacity to observe heterogenous boundary layers and to use these observations to improve data assimilation; Dealing with terrain (in all its variations) in a more appropriate way; Conducting a coordinated, joint field program to address the priorities/requirements listed above. The Federal Coordinator brought the results of the successful forum to the attention of the broader interagency PBL community with the goal of eliciting from that community how best to proceed in meeting the operational and research requirements articulated at the special session. On December 16, 2011, the Federal Coordinator also convened a meeting of the Interdepartmental Committee for Meteorological Services and Supporting Research (ICMSSR), a committee with members from 12 federal agencies, including DOE. One of the topics addressed at the meeting was the status of coupling ATD models and the impact across the federal meteorological community. A status report was provided to all ICMSSR members.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Entropy Stable Discretization of Compressible Flows in Thermochemical Nonequilibrium

Entropy stable numerical methods for compressible flow have been demonstrated to exhibit better robustness than purely linearly stable methods and need less overall artificial dissipation for long simulations in subsonic and transonic flows. In this work we seek to extend these benefits to multicomponent, multitemperature flows in thermochemical nonequilibrium such as combustion and hypersonic flight. We first derive entropy functions that symmetrize the governing equations and allow stability proofs for such systems. The impact of diffusion model selection on provable entropy stability is considered in detail, including both rigorous models of irreversible thermodynamics and simplified models of greater practical interest. Based on the proven entropy functions we develop affordable, entropy conservative two-point flux functions for solution in conservation form. We derive entropy conservative fluxes for calorically and thermally perfect mixtures, with heat capacities described by either polynomials of the temperature or formulas from statistical thermodynamics.

42 ENGINEERING↗
DOE OSTI2013

Probing Hohlraum Plasma Conditions by Neutronics From Fusion in the Hohlraum Volume

We propose a novel technique to measure the ion temperature in ignition scale hohlraums, namely by filling the hohlraum with DT gas. Quantifying the neutron production within the hohlraum volume (in space and time) will help measure the ion temperature. Varying the amount of cross beam energy transfer (CBET) is predicted to affect the ion temperature due to saturation effects. By varying CBET (via wavelength choices in the beams) we could monitor the DT production as a measure of this ion heating effect. Ion diffusion models in the simulations can also be tested in this manner.

Rosen, M. D.↗
DOE OSTI2025

Developing the Science Basis for Understanding Polymer Encapsulant Degradation Mechanisms: DuraMAT 2.0 Final Project Report

Polymeric encapsulants are essential materials in photovoltaic modules, protecting sensitive electronics from the environment while providing mechanical integrity to the multilayered assembly. However, these polymeric materials are susceptible to degradation processes driven by the ingress of environmental species, ultraviolet radiation, thermal stresses, and mechanical loading. In this study, we employ a combined atomistic simulation and accelerated aging experimental approach to study the molecular-scale mechanisms of encapsulant degradation. Classical molecular dynamics simulations quantify the diffusion of environmental and degradation species through the polymer matrix, producing composition-specific diffusion coefficients. Reactive simulations characterize activation energy barriers and reaction rate constants for key chemical pathways. In parallel, thermal-desorption analyses coupled with mass spectrometry monitor the emergence and concentration profiles of degradation products under controlled stressor conditions. By integrating simulation and experiment, we establish quantitative correlations between polymer composition, species diffusivity, and chemical reactivity. We anticipate that these relations and quantitative values could serve as high-fidelity inputs to reaction-diffusion models, enabling physics-informed lifetime predictions and guiding the design of more durable encapsulant materials for solar energy applications.

36 MATERIALS SCIENCE↗