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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Catalytic resonance of ammonia synthesis by simulated dynamic ruthenium crystal strain

Ammonia affords dense storage for renewable energy as a fungible liquid fuel, provided it can be efficiently synthesized from hydrogen and nitrogen. In this work, the catalysis of ammonia synthesis was computationally explored beyond the Sabatier limit by dynamically straining a ruthenium crystal (±4%) at the resonant frequencies (10 2 to 10 5+ Hz) of N 2 surface dissociation and hydrogenation. Density functional theory calculations at different strain conditions indicated that the energies of NH x surface intermediates and transition states scale linearly, allowing the description of ammonia synthesis at a continuum of strain conditions. A microkinetic model including multiple sites and surface diffusion between step and Ru(0001) terrace sites of varying ratios for nanoparticles of differing size revealed that dynamic strain yields catalytic ammonia synthesis conversion and turnover frequency comparable to industrial reactors (400°C, 200 atm) but at lower temperature (320°C) and an order of magnitude lower pressure (20 atm).

10 SYNTHETIC FUELS↗

DOE OSTI2018

Calculation of the Influence of the Ion Current Density and Temperature on the Accumulation Kinetics of Point Defects under the Irradiation of Si with Light Ions

Numerical calculations of the accumulation kinetics of point defects—vacancies and divacancies—under the irradiation of Si with ions having masses of M{sub 1} ≤ 31 amu and energies of E ≤ 100 keV under various irradiation conditions are performed. The previously proposed diffusion-coagulation model is used without application of the “weak diffusion” approximation, which was performed during its analytical implementation. The main peculiarities of the dependences of the concentrations of vacancies and divacancies on the dose, ion-current density, and temperature under irradiation are analyzed. A physical interpretation of these results is given. The developed computing complex is rather flexible and makes it possible to analyze the influence of model input parameters by means of their variation and include additional processes into consideration if necessary.

36 MATERIALS SCIENCE↗

NASA NTRS1994

Comparison of PDF and Moment Closure Methods in the Modeling of Turbulent Reacting Flows

In modeling turbulent reactive flows, Probability Density Function (PDF) methods have an advantage over the more traditional moment closure schemes in that the PDF formulation treats the chemical reaction source terms exactly, while moment closure methods are required to model the mean reaction rate. The common model used is the laminar chemistry approximation, where the effects of turbulence on the reaction are assumed negligible. For flows with low turbulence levels and fast chemistry, the difference between the two methods can be expected to be small. However for flows with finite rate chemistry and high turbulence levels, significant errors can be expected in the moment closure method. In this paper, the ability of the PDF method and the moment closure scheme to accurately model a turbulent reacting flow is tested. To accomplish this, both schemes were used to model a CO/H2/N2- air piloted diffusion flame near extinction. Identical thermochemistry, turbulence models, initial conditions and boundary conditions are employed to ensure a consistent comparison can be made. The results of the two methods are compared to experimental data as well as to each other. The comparison reveals that the PDF method provides good agreement with the experimental data, while the moment closure scheme incorrectly shows a broad, laminar-like flame structure.

Norris, Andrew T.↗

NASA NTRS1995

Theoretical Study of Response of InGaAs Metal-Semiconductor-Metal Photodetectors

We present a theoretical investigation of the response of metal-semiconductor-metal (MSM) photodetectors made of InGaAs lattice-matched to InP using a two-dimensional drift-diffusion model with a thermionic-field emission boundary condition for the heterojunctions. The effect of including a top InAlAs layer to increase the effective barrier height of the metal fingers on the InGaAs active layer is thoroughly examined and found to limit the collection of the photocurrent signal due to the electron and hole barriers that it forms with InGaAs. Due to the thickness and height of the InAlAs barrier layer in existing designs, the tunneling current obtained from the model is found to be negligibly small to significantly affect the output signal current. In an attempt to obtain a better response, different design structures including one where a quasi-Schottky contact is utilized are studied and their speed of response, breakdown voltage, and dark current are compared to that of the usual InGaAs device.

Salem, Ali F.↗

NASA NTRS2023

Numerical Investigation of Film Coefficient Approximation for Chemically Reacting Boundary-Layer Flows

Aerothermal analysis of spacecraft planetary entry is heavily dependent on heritage engineering models. The film coefficient heat transfer model examined in this paper estimates the convective heating to the vehicle for a laminar, dissociated, chemically reacting boundary layer for an Earth atmosphere. This model requires information about the vehicle and flowfield for a given trajectory point and estimates a proportional relationship between enthalpy potential and convective heat flux. In practice it is the aerothermal engineer who must decide which assumptions are appropriate for his/her application. This work looks at numerous CFD simulations for an arbitrary, axisymmetric flight vehicle to analyze the relative importance of both the mass and energy constraints imposed at the wall boundary, as well as the effect of various diffusion models. Within the subset of tested energy boundary conditions, it is found that the most desirable energy boundary condition is the radiative equilibrium boundary condition, which permits conservative estimates of convective heat flux, but also generates flowfield-dependent spatial thermal distributions along the surface. Other key findings are presented in an effort to make the film coefficient engineering model readily available to design engineers across industry.

Film coefficient↗

DOE OSTI2019

Transport and Loss of Ring Current Electrons Inside Geosynchronous Orbit During the 17 March 2013 Storm

Ring current electrons (1–100 keV) have received significant attention in recent decades, but many questions regarding their major transport and loss mechanisms remain open. In this study, we use the four-dimensional Versatile Electron Radiation Belt code to model the enhancement of phase space density that occurred during the 17 March 2013 storm. Our model includes global convection, radial diffusion, and scattering into the Earth's atmosphere driven by whistler-mode hiss and chorus waves. We study the sensitivity of the model to the boundary conditions, global electric field, the electric field associated with subauroral polarization streams, electron loss rates, and radial diffusion coefficients. The results of the code are almost insensitive to the model parameters above 4.5 R E R E , which indicates that the general dynamics of the electrons between 4.5 R E and the geostationary orbit can be explained by global convection. We found that the major discrepancies between the model and data can stem from the inaccurate electric field model and uncertainties in lifetimes. We show that additional mechanisms that are responsible for radial transport are required to explain the dynamics of ≥40-keV electrons, and the inclusion of the radial diffusion rates that are typically assumed in radiation belt studies leads to a better agreement with the data. The overall effect of subauroral polarization streams on the electron phase space density profiles seems to be smaller than the uncertainties in other input parameters. This study is an initial step toward understanding the dynamics of these particles inside the geostationary orbit.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

The accuracy of multi-group models for nonlocal electron transport in magnetized plasmas

In the extreme conditions of inertial confinement fusion experiments, heat flow plays a vital role, but local diffusive models frequently break down and overestimate the heat flow. The situation becomes more complicated again in the significant magnetic fields generated during laser–plasma interactions or in magnetized fusion schemes. Accurate non-local and magnetized heat flow computations can be carried out using Vlasov–Fokker–Planck (VFP) simulations, but these are computationally expensive. There is, therefore, significant interest in using faster multi-group models to accurately calculate the non-local heat flow in magnetized plasmas. We benchmark two such multi-group models for calculating the heat flow, M1 and hybrid-AWBS-BGK, against diffusive models and full VFP simulations, before applying the models to realistic example test cases, both magnetized and unmagnetized. We find that the multi-group models generally perform very well for moderate non-localities up to kλmfp∼0.01, but the computational cost increases dramatically. hybrid-AWBS-BGK performs more effectively than M1 at high non-localities, up to kλmfp∼1, due to its adaptive solver and robust P1 closure, but tends to fail in very strong magnetic fields. Both codes are much faster than VFP simulations but are still slow in steep temperature gradients.

Arran, C. (ORCID:0000000286448118)↗

DOE OSTI2016

Experimental Studies and Model Analysis of Noble Gas Fractionation in Porous Media

Core Ideas Ambient noble gases in porous media fractionate from through‐flow of other gases. CO 2 fluxes and noble‐gas fractionation are measured in simple sandbox experiments. Modeled fractionation provides indirect proxy, e.g., of steady CO 2 leakage fluxes. The noble gases, which are chemically inert under normal terrestrial conditions but vary systematically across a wide range of atomic mass and diffusivity, offer a multicomponent approach to investigating gas dynamics in unsaturated soil horizons, including transfer of gas between saturated zones, unsaturated zones, and the atmosphere. To evaluate the degree to which fractionation of noble gases in the presence of an advective–diffusive flux agrees with existing theory, a simple laboratory sand column experiment was conducted. Pure CO 2 was injected at the base of the column, providing a series of constant CO 2 fluxes through the column. At five fixed sampling depths within the system, samples were collected for CO 2 and noble gas analyses, and ambient pressures were measured. Both the advection–diffusion and dusty gas models were used to simulate the behavior of CO 2 and noble gases under the experimental conditions, and the simulations were compared with the measured depth‐dependent concentration profiles of the gases. Given the relatively high permeability of the sand column (5 × 10 −11 m 2 ), Knudsen diffusion terms were small, and both the dusty gas model and the advection–diffusion model accurately predicted the concentration profiles of the CO 2 and atmospheric noble gases across a range of CO 2 flux from ∼700 to 10,000 g m −2 d −1 . The agreement between predicted and measured gas concentrations demonstrated that, when applied to natural systems, the multi‐component capability provided by the noble gases can be exploited to constrain component and total gas fluxes of non‐conserved (CO 2 ) and conserved (noble gas) species or attributes of the soil column relevant to gas transport, such as porosity, tortuosity, and gas saturation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2018

Quantum Transmission Conditions for Diffusive Transport in Graphene with Steep Potentials

We present a formal derivation of a drift-diffusion model for stationary electron transport in graphene, in presence of sharp potential profiles, such as barriers and steps. Assuming the electric potential to have steep variations within a strip of vanishing width on a macroscopic scale, such strip is viewed as a quantum interface that couples the classical regions at its left and right sides. In the two classical regions, where the potential is assumed to be smooth, electron and hole transport is described in terms of semiclassical kinetic equations. The diffusive limit of the kinetic model is derived by means of a Hilbert expansion and a boundary layer analysis, and consists of drift-diffusion equations in the classical regions, coupled by quantum diffusive transmission conditions through the interface. The boundary layer analysis leads to the discussion of a four-fold Milne (half-space, half-range) transport problem.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2021

Tools for Assessing Performance Project: FY2021 Quarter 4 Report

The rotatable building located at Texas Tech’s Reese Technology Center is approximately 14 m wide (width being defined as more normal to the wind than parallel), 9 m long (aligned more with the wind than perpendicular), and 4 m tall. By placing 29 sonic anemometers downwind of the building, see Fig. 1, this facility provided a chance to collect data concerning both the wake velocity deficit distribution in the downstream and crosssteam directions (relative to the mean wind) behind an isolated building. A preliminary comparison between this data and the recently- developed fast-running diffusive wake model, which was developed based on wind tunnel and LES (performed with JOULES) simulations, for the purpose of either validating this model of understanding potential persistent differences between the idealized wind tunnel or LES conditions and those of full-scale real world phenomena. For this preliminary exploration, we utilized the diffusive wake model implemented in the QUIC model.

42 ENGINEERING↗

NASA NTRS1991

Analysis of redox additive-based overcharge protection for rechargeable lithium batteries

The overcharge condition in secondary lithium batteries employing redox additives for overcharge protection, has been theoretically analyzed in terms of a finite linear diffusion model. The analysis leads to expressions relating the steady-state overcharge current density and cell voltage to the concentration, diffusion coefficient, standard reduction potential of the redox couple, and interelectrode distance. The model permits the estimation of the maximum permissible overcharge rate for any chosen set of system conditions. Digital simulation of the overcharge experiment leads to numerical representation of the potential transients, and estimate of the influence of diffusion coefficient and interelectrode distance on the transient attainment of the steady state during overcharge. The model has been experimentally verified using 1,1-prime-dimethyl ferrocene as a redox additive. The analysis of the experimental results in terms of the theory allows the calculation of the diffusion coefficient and the formal potential of the redox couple. The model and the theoretical results may be exploited in the design and optimization of overcharge protection by the redox additive approach.

Narayanan, S. R.↗

DOE OSTI2017

Diffusive mass transport in agglomerated glassy fallout from a near-surface nuclear test

Aerodynamically-shaped glassy fallout is formed when vapor phase constituents from the nuclear device are incorporated into molten carriers (i.e. fallout precursor materials derived from soil or other near-field environmental debris). The effects of speciation and diffusive transport of condensing constituents are not well defined in models of fallout formation. Previously we reported observations of diffuse micrometer scale layers enriched in Na, Fe, Ca, and 235 U, and depleted in Al and Ti, at the interfaces of agglomerated fallout objects. Here in this paper, we derive the timescales of uranium mass transport in such fallout as it cools from 2500 K to 1500 K by applying a 1-dimensional planar diffusion model to the observed 235 U/ 30 Si variation at the interfaces. By modeling the thermal transport between the fireball and the carrier materials, the time of mass transport is calculated to be <0.6 s, <1 s, <2 s, and <3.5 s for fireball yields of 0.1 kt, 1 kt, 10 kt, and 100 kt respectively. Based on the calculated times of mass transport, a maximum temperature of deposition of uranium onto the carrier material of ~2200 K is inferred (1σ uncertainty of ~200 K). We also determine that the occurrence of micrometer scale layers of material enriched in relatively volatile Na-species as well as more refractory Ca-species provides evidence for an oxygen-rich fireball based on the vapor pressure of the two species under oxidizing conditions. These results represent the first application of diffusion-based modeling to derive material transport, thermal environments, and oxidation-speciation in near-surface nuclear detonation environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2014

Simulation of xenon, uranium vacancy and interstitial diffusion and grain boundary segregation in UO 2

In light water reactor fuel, gaseous fission products segregate to grain boundaries, resulting in the nucleation and growth of large intergranular fission gas bubbles. Based on the mechanisms established from density functional theory (DFT) and empirical potential calculations 1, continuum models for diffusion of xenon (Xe), uranium (U) vacancies and U interstitials in UO 2 have been derived for both intrinsic conditions and under irradiation. Segregation of Xe to grain boundaries is described by combining the bulk diffusion model with a model for the interaction between Xe atoms and three different grain boundaries in UO 2 ( Σ5 tilt, Σ5 twist and a high angle random boundary),as derived from atomistic calculations. All models are implemented in the MARMOT phase field code, which is used to calculate effective Xe and U diffusivities as well as redistribution for a few simple microstructures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS2001

Flight Experiment to Study Double-Diffusive Instabilities in Silver-Dopped Lead Bromide Crystals

The main objective of the present program is to understand thermosolutal convection during crystal growth of PbBr2-AgBr alloys. This involves identification of the growth conditions for microgravity experiments delineating the microsegregation, observation of convecto-diffusive instabilities and comparison with theoretical models. The overall objectives can be summarized as follows: 1. Observe and study the double diffusive and morphological instabilities in controlled conditions and to compare with theoretically predicted convective and morphological instability curves. 2. Study the three-dimensional morphological instabilities and resulting cellular growth that occur near the onset of morphological instability in the bulk samples under purely diffusive conditions. 3. Understand the micro-and macro-segregation of silver dopant in lead bromide crystals in microgravity. 4. Provide basic data on convective behavior in alloy crystals grown by the commercially important Bridgman crystal growth process.

Singh, N. B.↗

DOE OSTI2013

Entropy-based separation of yeast cells using a microfluidic system of conjoined spheres

A physical model is derived to create a biological cell separator that is based on controlling the entropy in a microfluidic system having conjoined spherical structures. A one-dimensional simplified model of this three-dimensional problem in terms of the corresponding effects of entropy on the Brownian motion of particles is presented. This dynamic mechanism is based on the Langevin equation from statistical thermodynamics and takes advantage of the characteristics of the Fokker-Planck equation. This mechanism can be applied to manipulate biological particles inside a microfluidic system with identical, conjoined, spherical compartments. This theoretical analysis is verified by performing a rapid and a simple technique for separating yeast cells in these conjoined, spherical microfluidic structures. The experimental results basically match with our theoretical model and we further analyze the parameters which can be used to control this separation mechanism. Both numerical simulations and experimental results show that the motion of the particles depends on the geometrical boundary conditions of the microfluidic system and the initial concentration of the diffusing material. This theoretical model can be implemented in future biophysics devices for the optimized design of passive cell sorters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2017

Turbulent mixing and removal of ozone within an Amazon rainforest canopy

Simultaneous profiles of turbulence statistics and mean ozone mixing ratio are used to establish a relation between eddy diffusivity and ozone mixing within the Amazon forest. A one–dimensional diffusion model is proposed and used to infer mixing time scales from the eddy diffusivity profiles. Data and model results indicate that during daytime conditions, the upper (lower) half of the canopy is well (partially) mixed most of the time and that most of the vertical extent of the forest can be mixed in less than an hour. During nighttime, most of the canopy is predominantly poorly mixed, except for periods with bursts of intermittent turbulence. Even though turbulence is faster than chemistry during daytime, both processes have comparable time scales in the lower canopy layers during nighttime conditions. Nonchemical loss time scales (associated with stomatal uptake and dry deposition) for the entire forest are comparable to turbulent mixing time scale in the lower canopy during the day and in the entire canopy during the night, indicating a tight coupling between turbulent transport and dry deposition and stomatal uptake processes. Because of the significant time of day and height variability of the turbulent mixing time scale inside the canopy, it is important to take it into account when studying chemical and biophysical processes happening in the forest environment. In conclusion, the method proposed here to estimate turbulent mixing time scales is a reliable alternative to currently used models, especially for situations in which the vertical distribution of the time scale is relevant.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Understanding the E balance for water management in hydroxide exchange membrane fuel cells

As a key challenge for hydroxide exchange membrane fuel cells (HEMFCs), understanding water management deeply is critical for performance and durability improvement. Here, in this work, we show that the balance between water generation and phase-change-induced (PCI) flow is the main factor controlling cell voltage over a wide range of conditions where the cells are vulnerable to severe flooding. The regulation of flooding by heat generation can be reflected by a stable voltage with increase of current density in the polarization curve. This PCI balance voltage (E balance ), which determines the shape of polarization curve, is modeled as a function of cell operating conditions and gas diffusion layer (GDL) properties. Temperature, pressure, and thermal conductivity of GDL were found important for tuning the level of E balance , both in HEMFCs and proton exchange membrane fuel cells. The understanding of E balance will help to generate strategies for resolving flooding in HEMFCs, and it offers a new perspective for polarization curve interpretation.

25 ENERGY STORAGE↗

DOE OSTI2022

Mason–Weaver theory: Revised and extended for a semi-infinite domain

Mason and Weaver developed equations to describe small particles settling under the influence of gravity and Brownian motion, including the limiting case for an infinitely deep suspension. We encountered this common convection–diffusion equation and no-flux boundary conditions in a model for dynamics of adsorbed polymers in dead end pores of a depolymerization catalyst. Close examination reveals that the Mason–Weaver solution is not correct for the infinite domain with a non-uniform initial condition. In this paper, we obtain the time dependent Green’s function for the no flux boundary condition and also for a more general reactive boundary condition. We demonstrate how the results provide solutions, via superposition, which provide solutions for several boundary conditions and all initial conditions.

Chen, Ziqiu↗