Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional diffusion models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Measurements of the electron energy distribution function in partially magnetized low temperature plasmas

While Langmuir probes (LPs) are relatively simple and inexpensive plasma diagnostics for the electron density, temperature, and the electron energy distribution function (EEDF), the interpretation of the measured current–voltage (I–V) characteristic is complicated considerably by the presence of a magnetic field. In regimes where the electron gyroradius is comparable to the probe radius, the electron flux to the probe surface is retarded by reduced mobility across field lines and can no longer be described by a thermal model. Predicting the current collected by the probe in these regimes requires accurate estimates of the plasma diffusion coefficients, which are usually difficult to obtain. In this work, we measure electron energy distribution functions in E×B plasmas with magnetized electrons and non-magnetized ions in argon and krypton gases, using both a LP and laser Thomson scattering (LTS) at various magnetic fields. Using the LTS measurements to provide a robust benchmark for comparison, we compare existing theories describing the flux to the probe under magnetized conditions. We find that even when the electron gyroradius associated with the effective electron temperature is small compared to the probe radius, the EEDF computed using classical probe theory is still robust at energies higher than the energy at which the gyroradius becomes larger than the probe size. For plasmas that are not strongly non-Maxwellian, we formulate a method to extract robust density and temperature measurements using physics-informed fitting techniques to analyze EEDFs computed using classical theory.

Devin, E. G. [Princeton Plasma Physics Laboratory ↗
DOE OSTI2016

THE EFFECTS OF GRAIN SIZE AND TEMPERATURE DISTRIBUTIONS ON THE FORMATION OF INTERSTELLAR ICE MANTLES

Computational models of interstellar gas-grain chemistry have historically adopted a single dust-grain size of 0.1 micron, assumed to be representative of the size distribution present in the interstellar medium. Here, we investigate the effects of a broad grain-size distribution on the chemistry of dust-grain surfaces and the subsequent build-up of molecular ices on the grains, using a three-phase gas-grain chemical model of a quiescent dark cloud. We include an explicit treatment of the grain temperatures, governed both by the visual extinction of the cloud and the size of each individual grain-size population. We find that the temperature difference plays a significant role in determining the total bulk ice composition across the grain-size distribution, while the effects of geometrical differences between size populations appear marginal. We also consider collapse from a diffuse to a dark cloud, allowing dust temperatures to fall. Under the initial diffuse conditions, small grains are too warm to promote grain-mantle build-up, with most ices forming on the mid-sized grains. As collapse proceeds, the more abundant, smallest grains cool and become the dominant ice carriers; the large population of small grains means that this ice is distributed across many grains, with perhaps no more than 40 monolayers of ice each (versus several hundred assuming a single grain size). This effect may be important for the subsequent processing and desorption of the ice during the hot-core phase of star formation, exposing a significant proportion of the ice to the gas phase, increasing the importance of ice-surface chemistry and surface–gas interactions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

The COS-Halos survey: origins of the highly ionized circumgalactic medium of star-forming galaxies

The total contribution of diffuse halo gas to the galaxy baryon budget strongly depends on its dominant ionization state. In this paper, we address the physical conditions in the highly ionized circumgalactic medium (CGM) traced by O VI absorption lines observed in COS-Halos spectra. We analyze the observed ionic column densities, absorption-line widths and relative velocities, along with the ratios of N V/O VI for 39 fitted Voigt profile components of O vi. We compare these quantities with the predictions given by a wide range of ionization models. Photoionization models that include only extragalactic UV background radiation are ruled out; conservatively, the upper limits to N V/O VI and measurements of N{sub OVI} imply unphysically large path lengths ≳100 kpc. Furthermore, very broad O VI absorption (b > 40 km s{sup −1}) is a defining characteristic of the CGM of star-forming L* galaxies. We highlight two possible origins for the bulk of the observed O VI: (1) highly structured gas clouds photoionized primarily by local high-energy sources or (2) gas radiatively cooling on large scales behind a supersonic wind. Approximately 20% of circumgalactic O vi does not align with any low-ionization state gas within ±50 km s{sup −1} and is found only in halos with M{sub halo} < 10{sup 12} M{sub ⊙}. We suggest that this type of unmatched O vi absorption traces the hot corona itself at a characteristic temperature of 10{sup 5.5} K. We discuss the implications of these very distinct physical origins for the dynamical state, gas cooling rates, and total baryonic content of L* gaseous halos.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Effect of Radiation on Chromospheric Magnetic Reconnection: Reactive and Collisional Multi-fluid Simulations

We study magnetic reconnection under chromospheric conditions in five different ionization levels from 0.5% to 50% using a self-consistent two-fluid (ions + neutrals) model that accounts for compressibility, collisional effects, chemical inequilibrium, and anisotropic heat conduction. Results with and without radiation are compared, using two models for the radiative losses: an optically thin radiation loss function, and an approximation of the radiative losses of a plasma with photospheric abundances. The results without radiation show that reconnection occurs faster for the weakly ionized cases as a result of the effect of ambipolar diffusion and fast recombination. The tearing mode instability appears earlier in the low ionized cases and grows rapidly. We find that radiative losses have a stronger effect than was found in previous results as the cooling changes the plasma pressure and the concentration of ions inside the current sheet. This affects the ambipolar diffusion and the chemical equilibrium, resulting in thin current sheets and enhanced reconnection. The results quantify this complex nonlinear interaction by showing that a strong cooling produces faster reconnections than have been found in models without radiation. The results accounting for radiation show timescales and outflows comparable to spicules and chromospheric jets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

On the importance of the heat and mass transfer resistances in internally-cooled liquid desiccant dehumidifiers and regenerators

Liquid desiccant heat and mass exchangers are a promising technology for efficient humidity control in buildings. Many researchers have investigated these exchangers, often using numerical models to predict their performance. However, there is a lack of information in the literature on the magnitude of the heat and mass transfer resistances, both for the dehumidifier (which absorbs moisture from the air) and the regenerator (which heats the liquid desiccant to re-concentrate it). This article focuses on internally-cooled, 3-fluid exchangers in a parallel plate geometry. Water heats or cools a desiccant across a plate, and the desiccant absorbs or releases water into an airstream through a membrane. A sensitivity analysis was used to estimate the importance of each of the heat and mass transfer resistances (air, membrane, desiccant, plate, water), and how it changes with different design geometries. The results show that, for most designs, the latent and sensible heat transfer of the dehumidifier is dominated by the air mass transfer resistance and air heat transfer resistance, respectively. The air mass transfer resistance is also important for the regenerator, but much less so; the change in the desiccant equilibrium humidity ratio due to a change in either temperature or desiccant mass fraction is much higher at the regenerator's higher temperatures. This increases the importance of (1) getting heat from the water to the desiccant/membrane interface, and (2) diffusing salt ions quickly away from the desiccant/membrane interface. The membrane heat transfer and water heat transfer resistances were found to be the least important. These results can help inform decisions about what simplifying assumptions to make in numerical models, and can also help in designing these exchangers by understanding which resistances are most important.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Effects of elastic softening and helium accumulation kinetics on surface morphological evolution of plasma-facing tungsten

Based on a continuous-domain model, capable of accessing the spatiotemporal scales relevant to fuzz formation on the surface of plasma-facing component (PFC) tungsten, we report self-consistent simulation results that elucidate the effects of elastic softening and helium (He) accumulation kinetics on the surface morphological response of PFC tungsten. The model accounts for the softening of the elastic moduli in the near-surface region of PFC tungsten, including both thermal softening at high temperature and softening due to He accumulation upon He implantation. The dependence of the elastic moduli on the He content follows an exponential scaling relation predicted by molecular-dynamics simulations, while the He content in the near-surface region of PFC tungsten evolves according to a first-order saturation kinetics, consistent with experimental and simulation results reported in the literature. We establish that this elastic softening accelerates both nanotendril growth on the PFC surface and the onset of fuzz formation. We also explore the role of the rate of He accumulation to a saturation level in the near-surface region of irradiated tungsten in the onset of fuzz formation. For PFC tungsten surfaces such as W(110) where, under typical irradiation conditions, the characteristic time scale for stress-driven surface diffusion is comparable to the characteristic time scale for He accumulation, we find that accelerating the rate of He accumulation accelerates the growth rate of nanotendrils emanating from the surface. Additionally, we present a systematic parametric study of the PFC surface morphological response to explore its dependence on the He accumulation kinetics that is controlled by the irradiation conditions for low-energy implantation. Finally, we introduce an incubation time for nanotendril growth on the PFC surface, a concept equivalent to that of incubation fluence discussed in the literature, to predict and explain the minimum exposure time required to observe fuzz formation on PFC tungsten surfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Topology optimization for the full-cell design of porous electrodes in electrochemical energy storage devices

In this paper, we introduce a density-based topology optimization framework to design porous electrodes for maximum energy storage. We simulate the full cell with a model that incorporates electronic potential, ionic potential, and electrolyte concentration. The system consists of three materials, namely pure liquid electrolyte and the porous solids of the anode and cathode, for which we determine the optimal placement. We use separate electronic potentials to model each electrode, which allows interdigitated designs. As a result, a penalization is required to ensure that the anode and cathode do not touch, i.e., causing a short circuit. We compare multiple 2D designs generated for different fixed conditions, e.g. material properties. A 3D design with complex channel and interlocked structure is also created. All optimized designs are far superior to the traditional monolithic electrode design with respect to energy storage metrics. We observe up to a 750% increase in energy storage for cases with slow effective ionic diffusion within the porous electrode.

25 ENERGY STORAGE↗
DOE OSTI2021

Modeling of High-Temperature Corrosion of Zirconium Alloys Using the eXtended Finite Element Method (X-FEM)

Oxidation modeling in modern nuclear fuel performance codes is currently limited by the lack of coupling with mechanics, thus preventing proper description of how high-temperature oxidation impacts mechanical properties. This is mostly due to the fact that the finite difference formalism adopted in corrosion models is incompatible with the direct coupling with mechanics in the finite element modeling employed in modern nuclear fuel performance codes. In this study, a physically based zirconium alloy corrosion model called the Coupled-Current Charge Compensation (C4) model, which was initially developed for operating temperature conditions, has been updated to include high-temperature corrosion in order to provide additional critical information (e.g., oxygen concentration profile) under loss-of-coolant accident (LOCA) conditions—information lacking in existing empirical models. The C4 model was implemented in the MOOSE finite-element framework developed at Idaho National Laboratory, enabling it to be used in the BISON nuclear fuel performance code based on the MOOSE framework. To precisely track the different interfaces at a relatively low computational cost, the eXtended Finite Element Method (X-FEM) was applied in MOOSE. The model’s results were compared to those of existing empirical models as well as metallographic analysis of high-temperature oxidized Zircaloy-4 coupons. Oxygen diffusivities in the a and ß phases resulting from this comparison closely agree with those found in the literature. The C4 model implemented with X-FEM in MOOSE now has the capability to accurately predict oxide, oxygen-stabilized a, and prior ß phase layer growth kinetics under isothermal exposure at high temperature (1000–1500°C). Furthermore, in contrast with the empirical models, the C4 model accounts for the finite thickness of the fuel cladding. It can predict the oxygen concentration profile evolution through the whole cladding, enabling evaluation of the remaining ductile thickness—a crucial variable for modeling the mechanical behavior of the fuel cladding under LOCA. Furthermore, this implementation allows direct coupling with mechanics, at a low computing cost, using finite-element-based nuclear fuel performance codes such as BISON.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Computational framework for modeling membrane processes without process and solution property simplifications

Accurately modeling membrane processes is critical to evaluating novel process configurations, designing scalable membrane systems, informing process cost estimates, and directing future research. Most membrane process models trade accuracy for computational efficiency by employing simplified approximations of the process (i.e. no salt flux, no pressure drop) and solution properties (i.e. ideal solution, and constant density, viscosity, and diffusivity). Our report introduces a detailed one-dimensional finite difference model for evaluating membrane processes that avoids these common simplifications. We apply this model to quantify the error introduced by these simplifications for case studies of reverse osmosis, osmotically assisted reverse osmosis, forward osmosis, and pressure retarded osmosis. While the magnitude of error introduced by these simplifications is reliant on the case study parameters and specifications, we find that existing model formulations can underestimate or overestimate average water flux by nearly 50% for some membrane processes operating under standard conditions. Lastly we investigate the error introduced by simplified inlet-outlet models that do not solve the governing system of differential equations, and we assess the accuracy of novel inlet-outlet formulations that use a log and geometric mean, instead of the typical arithmetic mean, to represent non-linear water flux profiles.

42 ENGINEERING↗
DOE OSTI2016

IceCube constraints on fast-spinning pulsars as high-energy neutrino sources

Relativistic winds of fast-spinning pulsars have been proposed as a potential site for cosmic-ray acceleration from very high energies (VHE) to ultrahigh energies (UHE). We re-examine conditions for high-energy neutrino production, considering the interaction of accelerated particles with baryons of the expanding supernova ejecta and the radiation fields in the wind nebula. We make use of the current IceCube sensitivity in diffusive high-energy neutrino background, in order to constrain the parameter space of the most extreme neutron stars as sources of VHE and UHE cosmic rays. We demonstrate that the current non-observation of 10{sup 18} eV neutrinos put stringent constraints on the pulsar scenario. For a given model, birthrates, ejecta mass and acceleration efficiency of the magnetar sources can be constrained. When we assume a proton cosmic ray composition and spherical supernovae ejecta, we find that the IceCube limits almost exclude their significant contribution to the observed UHE cosmic-ray flux. Furthermore, we consider scenarios where a fraction of cosmic rays can escape from jet-like structures piercing the ejecta, without significant interactions. Such scenarios would enable the production of UHE cosmic rays and help remove the tension between their EeV neutrino production and the observational data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Extension of the PINN diffusion model to k-eigenvalue problems

This paper extends our recent work on the Physics-Informed Neural Networks (PINN) approach for the fixed source diffusion models and applies it to the diffusion theory based k-eigenvalue problems. To make the PINN equitable for the eigenvalue problems, we introduce a novel integral regularization term to the loss function in the framework, and allow the direct inference of the principal eigenvalue and the associated eigenfunction. The regularization term enforces a pre-defined value on the integration of the model predictions, and this value can be directly related to a physical property of the system. We also introduce an additional learnable parameter to approximate the principal eigenvalue. As a proof of principle, we solve the one-group two-dimensional k-eigenvalue neutron diffusion equation in this work. We then provide two numerical examples to demonstrate the applicability of the PINN approach. In each example, we solve the k-eigenvalue diffusion equation in a multi-region configuration constrained with a set of Robin boundary conditions for generality. We use a FEM solution based on the power-iteration method to verify the results of the PINN solution. The results showed relative percentage error in the predicted eigenvalue of about 0.77% and about 1.2% for example 1 and example 2, respectively. The mean absolute error in the predicted flux for example 1 is ∼ 0.002 and for example 2 is ∼ 0.0024. These results indicate some preliminary successes of the PINN application to k-eigenvalue problems. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2013

The eigenmode perspective of NMR spin relaxation in proteins

We developed in recent years the two-body (protein and probe) coupled-rotator slowly relaxing local structure (SRLS) approach for elucidating protein dynamics from NMR spin relaxation. So far we used as descriptors the set of physical parameters that enter the SRLS model. They include the global (protein-related) diffusion tensor, D{sub 1}, the local (probe-related) diffusion tensor, D{sub 2}, and the local coupling/ordering potential, u. As common in analyzes based on mesoscopic dynamic models, these parameters have been determined with data-fitting techniques. In this study, we describe structural dynamics in terms of the eigenmodes comprising the SRLS time correlation functions (TCFs) generated by using the best-fit parameters as input to the Smoluchowski equation. An eigenmode is a weighted exponential with decay constant given by an eigenvalue of the Smoluchowski operator, and weighting factor determined by the corresponding eigenvector. Obviously, both quantities depend on the SRLS parameters as determined by the SRLS model. Unlike the set of best-fit parameters, the eigenmodes represent patterns of motion of the probe-protein system. The following new information is obtained for the typical probe, the {sup 15}N−{sup 1}H bond. Two eigenmodes, associated with the protein and the probe, dominate when the time scale separation is large (i.e., D{sub 2} ≫ D{sub 1}), the tensorial properties are simple, and the local potential is either very strong or very weak. When the potential exceeds these limits while the remaining conditions are preserved, new eigenmodes arise. The multi-exponentiality of the TCFs is associated in this case with the restricted nature of the local motion. When the time scale separation is no longer large, the rotational degrees of freedom of the protein and the probe become statistically dependent (coupled dynamically). The multi-exponentiality of the TCFs is associated in this case with the restricted nature of both the local and the global motion. The effects of local diffusion axiality, potential strength, and extent of mode-coupling on the eigenmode setup are investigated. We detect largely global motional or largely local motional eigenmodes. In addition, we detect mixed eigenmodes associated with correlated/prograde or anti-correlated/retrograde rotations of the global (D{sub 1}) and local (D{sub 2}) motional modes. The eigenmode paradigm is applied to N−H bond dynamics in the β-sheet residue K19, and the α-helix residue A34, of the third immunoglobulin-binding domain of streptococcal protein G. The largest contribution to the SRLS TCFs is made by mixed anti-correlated D{sub 1} and D{sub 2} eigenmodes. The next largest contribution is made by D{sub 1}-dominated eigenmodes. Eigenmodes dominated by the local motion contribute appreciably to A34 and marginally to K19. Correlated D{sub 1} and D{sub 2} eigenmodes contribute exclusively to K19 and do not contribute above 1% to A34. The differences between K19 and A34 are delineated and rationalized in terms of the best-fit SRLS parameters and mode-mixing. It may be concluded that eigenmode analysis is complementary and supplementary to data-fitting-based analysis.

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

Consideration of Grain Size Distribution in the Diffusion of Fission Gas to Grain Boundaries

We analyze the accumulation of fission gas on grain boundaries in a polycrystalline microstructure with a distribution of grain sizes. The diffusion equation is solved throughout the microstructure to evolve the gas concentration in space and time. Grain boundaries are treated as infinite sinks for the gas concentration, and we monitor the cumulative gas inventory on each grain boundary throughout time. We consider two important cases: first, a uniform initial distribution of gas concentration without gas production (correlating with post-irradiation annealing), and second, a constant gas production rate with no initial gas concentration (correlating with in-reactor conditions). The results show that a single-grain-size model, such as the Booth model, over predicts the gas accumulation on grain boundaries compared with a polycrystal with a grain size distribution. Also, a considerable degree of scatter, or variability, exists in the grain boundary gas accumulation when comparing all of the grain boundaries in the microstructure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2013

Exchange and spin states in quantum dots under strong spatial correlations. Computer simulation by the Feynman path integral method

The fundamental laws in the behavior of electrons in model quantum dots that are caused by exchange and strong Coulomb correlations are studied. The ab initio path integral method is used to numerically simulate systems of two, three, four, and six interacting identical electrons confined in a three-dimensional spherical potential well with a parabolic confining potential against the background of thermal fluctuations. The temperature dependences of spin and collective spin magnetic susceptibility are calculated for model quantum dots of various spatial sizes. A basically exact procedure is proposed for taking into account the permutation symmetry and the spin state of electrons, which makes it possible to perform numerical calculations using modern computer facilities. The conditions of applicability of a virial energy estimator and its optimum form in exchange systems are determined. A correlation estimator of kinetic energy, which is an alternative to a basic estimator, is suggested. A fundamental relation between the kinetic energy of a quantum particle and the character of its virtual diffusion in imaginary time is demonstrated. The process of natural 'pairing' of electron spins during the compression of a quantum dot and cooling of a system is numerically reproduced in terms of path integrals. The temperature dependences of the spin magnetic susceptibility of electron pairs with a characteristic maximum caused by spin pairing are obtained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Automated Parameterized CFD Simulations of Phase-Change Material Embedded Heat Exchangers

Thermal storage can be implemented using Phase-Change Materials (PCM), which absorb significant latent heat with a relatively small temperature change. PCM phase-change processes are transient and are driven by thermal diffusion and natural convection – the latter, especially for melting process. Modeling and simulation of PCM heat exchangers is typically computationally intensive due to the relative complex time-dependent underlying physics. This paper presents a study investigating a single straight tube with circular transverse fins in the vertical orientation, using PCM’s with 37oC nominal melting temperature, leveraging automated parametric CFD simulations. The computational domain is reduced to 2D axisymmetric for a segment of the HX comprising a region in between two fins with constant wall temperature as boundary condition. This approach allows for faster evaluations of the transient thermal characteristics for a given design space of interest. Correlations and reduced order models (ROM’s) are developed. Results indicate that the ROM’s reproduce the CFD simulations with high accuracy for designs from the training data set, however, designs from cross-validation and test sets do not always have good agreement. The limited number of designs in the training set is a possible cause for the lack of robustness of the ROM’s and thus additional training designs will be added in future work. The ROM’s solve the simulations on average 500,000 times faster than CFD.

25 ENERGY STORAGE↗
DOE OSTI2013

An ab initio model for the modulation of galactic cosmic-ray electrons

The modulation of galactic cosmic-ray electrons is studied using an ab initio three-dimensional steady state cosmic-ray modulation code in which the effects of turbulence on both the diffusion and drift of these cosmic-rays are treated as self-consistently as possible. A significant refinement is that a recent two-component turbulence transport model is used. This model yields results in reasonable agreement with observations of turbulence quantities throughout the heliosphere. The sensitivity of computed galactic electron intensities to choices of various turbulence parameters pertaining to the dissipation range of the slab turbulence spectrum, and to the choice of model of dynamical turbulence, is demonstrated using diffusion coefficients derived from the quasi-linear and extended nonlinear guiding center theories. Computed electron intensities and latitude gradients are also compared with spacecraft observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Dynamics of Nanoparticles in Entangled Polymer Solutions

The mean square displacement < r 2 > of nanoparticle probes dispersed in simple isotropic liquids and in polymer solutions is interrogated using fluorescence correlation spectroscopy and single-particle tracking (SPT) experiments. Probe dynamics in different regimes of particle diameter (d), relative to characteristic polymer length scales, including the correlation length (ξ), the entanglement mesh size (a), and the radius of gyration (R g ), are investigated. In simple fluids and for polymer solutions in which d >> R g , long-time particle dynamics obey random-walk statistics < r 2 >:t, with the bulk zero-shear viscosity of the polymer solution determining the frictional resistance to particle motion. In contrast, in polymer solutions with d < R g , polymer molecules in solution exert noncontinuum resistances to particle motion and nanoparticle probes appear to interact hydrodynamically only with a local fluid medium with effective drag comparable to that of a solution of polymer chain segments with sizes similar to those of the nanoparticle probes. Under these conditions, the nanoparticles exhibit orders of magnitude faster dynamics than those expected from continuum predictions based on the Stokes–Einstein relation. SPT measurements further show that when d > a, nanoparticle dynamics transition from diffusive to subdiffusive on long timescales, reminiscent of particle transport in a field with obstructions. This last finding is in stark contrast to the nanoparticle dynamics observed in entangled polymer melts, where X-ray photon correlation spectroscopy measurements reveal faster but hyperdiffusive dynamics. As a result, we analyze these results with the help of the hopping model for particle dynamics in polymers proposed by Cai et al. and, on that basis, discuss the physical origins of the local drag experienced by the nanoparticles in entangled polymer solutions.

77 NANOSCIENCE AND NANOTECHNOLOGY↗