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 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Helium bubble formation in aged plutonium metal at ambient conditions

A model of helium bubble formation is developed that assumes the initial nascent bubbles are generated when two diffusing helium atoms come within a distance of about the interatomic spacing of the plutonium atoms. These nascent bubbles are assumed to be locked in place and grow as additional helium atoms diffuse across their surfaces. Calculations of the helium-bubble number density and sizes as a function of age are compared to experiment. This leads to an inferred helium diffusion coefficient of-the-order of 10 -20 m 2 /s, supports the idea that most of the helium from alpha-decay end up in bubbles, and that previous transmission electron microscopy size estimates are biased low. The energy that helium can transfer to 40-year-old plutonium, upon expansion, is determined to be <0.1 Joules per gram of metal. The model predicts differences between the bubble properties in naturally and accelerated aged plutonium.

36 MATERIALS SCIENCE↗

DOE OSTI2015

Chaotic dynamics and diffusion in a piecewise linear equation

Genetic interactions are often modeled by logical networks in which time is discrete and all gene activity states update simultaneously. However, there is no synchronizing clock in organisms. An alternative model assumes that the logical network is preserved and plays a key role in driving the dynamics in piecewise nonlinear differential equations. We examine dynamics in a particular 4-dimensional equation of this class. In the equation, two of the variables form a negative feedback loop that drives a second negative feedback loop. By modifying the original equations by eliminating exponential decay, we generate a modified system that is amenable to detailed analysis. In the modified system, we can determine in detail the Poincaré (return) map on a cross section to the flow. By analyzing the eigenvalues of the map for the different trajectories, we are able to show that except for a set of measure 0, the flow must necessarily have an eigenvalue greater than 1 and hence there is sensitive dependence on initial conditions. Further, there is an irregular oscillation whose amplitude is described by a diffusive process that is well-modeled by the Irwin-Hall distribution. There is a large class of other piecewise-linear networks that might be analyzed using similar methods. The analysis gives insight into possible origins of chaotic dynamics in periodically forced dynamical systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

Modeling, control, and dynamic performance analysis of a reverse osmosis desalination plant integrated within hybrid energy systems

An RO (reverse osmosis) desalination plant is proposed as an effective, FLR (flexible load resource) to be integrated into HES (hybrid energy systems) to support various types of ancillary services to the electric grid, under variable operating conditions. To study the dynamic (transient) analysis of such system, among the various unit operations within HES, special attention is given here to the detailed dynamic modeling and control design of RO desalination process with a spiral-wound membrane module. The model incorporates key physical phenomena that have been investigated individually into a dynamic integrated model framework. In particular, the solution-diffusion model modified with the concentration polarization theory is applied to predict RO performance over a large range of operating conditions. Simulation results involving several case studies suggest that an RO desalination plant, acting as a FLR, can provide operational flexibility to participate in energy management at the utility scale by dynamically optimizing the use of excess electrical energy. Here, the incorporation of additional commodity (fresh water) produced from a FLR allows a broader range of HES operations for maximizing overall system performance and profitability. For the purpose of assessing the incorporation of health assessment into process operations, an online condition monitoring approach for RO membrane fouling supervision is addressed in the case study presented.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2017

Heterogeneous Neutron Leakage Model for PWR Pin-by-pin Calculation

In the last few decades, the nodal diffusion calculation accompanied with the assembly homogenization has been widely employed. The single assembly calculation with the reflective boundary condition (BC) would be carried out to determine the assembly-homogenized parameters. Because of the reflective BC, no axial and/or radial leakage rates have been taken into accounted for the lattice calculation. To represent the global influence of the reactor core on an assembly calculation, one mostly frequently uses a procedure called leakage model to approximately take into account the absorption due to the macroscopic structure of the physical flux in a heterogeneous situation. For the assembly homogenization, this absorption is introduced as the product of the leakage coefficient and the buckling obtained for an equivalent homogeneous medium using the well-known homogeneous B1 theory. Clearly, this is an imperfect representation of leakages in a heterogeneous medium. However, the errors introduced by the foregoing procedure are small for the whole? -assembly homogenization in a PWR. When it comes to the pin-by-pin calculation, it is not the case. The leakage rate? will be affected by the neutron streaming effect caused by the heterogeneity of the assembly. Especially when poison pins or empty channels are present, the streaming effect becomes important. The errors introduced appear particularly visible. The heterogeneous leakage model? which assumes that the neutron flux is variational and the leakage rate is not uniform over the lattice pitch is a more appropriate choice. When the heterogeneous leakage model is applied in the lattice calculation for the pin-cell homogenization, the pin-homogenized problem would be changed compared with the original assembly problem without any homogenization technique. So, the Cooperation of the heterogeneous leakage model with the pin-cell homogenization theory is unsettled. In this paper, both homogenous and heterogeneous leakage models are evaluated in pin-by-pin calculation. Their implements used in pin-cell homogenization of the lattice calculation are studied. The diffusion coefficient obtained from transport cross-section and the space-dependent leakage coefficient together with the two leakage models are analyzed for the practical reactor-core applications.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Re-evaluating cloud chamber constraints on depositional ice growth in cirrus clouds – Part 1: Model description and sensitivity tests

Abstract. Ice growth from vapor deposition is an important process for the evolution of cirrus clouds, but the physics of depositional ice growth at the low temperatures (<235 K) characteristic of the upper troposphere and lower stratosphere is not well understood. Surface attachment kinetics, generally parameterized as a deposition coefficient αD, control ice crystal habit and also may limit growth rates in certain cases, but significant discrepancies between experimental measurements have not been satisfactorily explained. Experiments on single ice crystals have previously indicated the deposition coefficient is a function of temperature and supersaturation, consistent with growth mechanisms controlled by the crystal's surface characteristics. Here we use observations from cloud chamber experiments in the Aerosol Interactions and Dynamics in the Atmosphere (AIDA) aerosol and cloud chamber to evaluate surface kinetic models in realistic cirrus conditions. These experiments have rapidly changing temperature, pressure, and ice supersaturation such that depositional ice growth may evolve from diffusion limited to surface kinetics limited over the course of a single experiment. In Part 1, we describe the adaptation of a Lagrangian parcel model with the Diffusion Surface Kinetics Ice Crystal Evolution (DiSKICE) model (Zhang and Harrington, 2014) to the AIDA chamber experiments. We compare the observed ice water content and saturation ratios to that derived under varying assumptions for ice surface growth mechanisms for experiments simulating ice clouds between 180 and 235 K and pressures between 150 and 300 hPa. We found that both heterogeneous and homogeneous nucleation experiments at higher temperatures (>205 K) could generally be modeled consistently with either a constant deposition coefficient or the DiSKICE model assuming growth on isometric crystals via abundant surface dislocations. Lower-temperature experiments showed more significant deviations from any depositional growth model, with different ice growth rates for heterogeneous and homogeneous nucleation experiments.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

A fractional model for anomalous diffusion with increased variability: Analysis, algorithms and applications to interface problems

Fractional equations have become the model of choice in several applications where heterogeneities at the microstructure result in anomalous diffusive behavior at the macroscale. Here, we introduce a new fractional operator characterized by a doubly-variable fractional order and possibly truncated interactions. Under certain conditions on the model parameters and on the regularity of the fractional order we show that the corresponding Poisson problem is well-posed. Additionally, we introduce a finite element discretization and describe an efficient implementation of the finite-element matrix assembly in the case of piecewise constant fractional order. Through several numerical tests, we illustrate the improved descriptive power of this new operator across media interfaces. Furthermore, we present one-dimensional and two-dimensional h-convergence results that show that the variable-order model has the same convergence behavior as the constant-order model.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Generative design of de novo proteins based on secondary-structure constraints using an attention-based diffusion model

We report two generative deep-learning models that predict amino acid sequences and 3D protein structures on the basis of secondary-structure design objectives via either the overall content or the per-residue structure. Both models are robust regarding imperfect inputs and offer de novo design capacity because they can discover new protein sequences not yet discovered from natural mechanisms or systems. The residue-level secondary-structure design model generally yields higher accuracy and more diverse sequences. These findings suggest unexplored opportunities for protein designs and functional outcomes within the vast amino acid sequences beyond known proteins. Our models, based on an attention-based diffusion model and trained on a dataset extracted from experimentally known 3D protein structures, offer numerous downstream applications in the conditional generative design of various biological or engineering systems. Future work could include additional conditioning and an exploration of other functional properties of the generated proteins for various properties beyond structural objectives.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2014

Numerical study of transition to supersonic flows in the edge plasma

The plasma scrape-off layer (SOL) in a tokamak is characterized by ion flow down a long narrow flux tube terminating on a solid surface. The ion flow velocity along a magnetic field line can be equal to or greater than sonic at the entrance of a Debye sheath or upstream in the presheath. This paper presents a numerical study of the transition between subsonic and supersonics flows. A quasineutral one-dimensional (1D) fluid code has been used for modeling of plasma transport in the SOL along magnetic field lines, both in steady state and under transient conditions. The model uses coupled equations for continuity, momentum, and energy balance with ionization, radiation, charge exchange, and recombination processes. The recycled neutrals are described in the diffusion approximation. Standard Bohm sheath criterion is used as boundary conditions at the material surface. Three conditions conducive for the generation of supersonic flows in SOL plasmas have been explored. It is found that in steady state high (attached) and low (detached) divertor temperatures cases, the role of particle, momentum, and energy loss is critical. For attached case, the appearance of shock waves in the divertor region if the incoming plasma flow is supersonic and its effect on impurity retention is presented. In the third case, plasma expansion along the magnetic field can yield time-dependent supersonic solutions in the quasineutral rarefaction wave. Such situations can arise in the parallel transport of intermittent structures such as blobs and edge localized mode filaments along field lines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Cooperative diffusion in body-centered cubic iron in Earth and super-Earths’ inner core conditions

Abstract The physical chemistry of iron at the inner-core conditions is key to understanding the evolution and habitability of Earth and super-Earth planets. Based on full first-principles simulations, we report cooperative diffusion along the longitudinally fast ⟨ 111 ⟩ directions of body-centered cubic (bcc) iron in temperature ranges of up to 2000–4000 K below melting and pressures of ∼300–4000 GPa. The diffusion is due to the low energy barrier in the corresponding direction and is accompanied by mechanical and dynamical stability, as well as strong elastic anisotropy of bcc iron. These findings provide a possible explanation for seismological signatures of the Earth’s inner core, particularly the positive correlation between P wave velocity and attenuation. The diffusion can also change the detailed mechanism of core convection by increasing the diffusivity and electrical conductivity and lowering the viscosity. The results need to be considered in future geophysical and planetary models and should motivate future studies of materials under extreme conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2016

Broadband light sources based on InAs/InGaAs metamorphic quantum dots

We propose a design for a semiconductor structure emitting broadband light in the infrared, based on InAs quantum dots (QDs) embedded into a metamorphic step-graded In{sub x}Ga{sub 1−x}As buffer. We developed a model to calculate the metamorphic QD energy levels based on the realistic QD parameters and on the strain-dependent material properties; we validated the results of simulations by comparison with the experimental values. On this basis, we designed a p-i-n heterostructure with a graded index profile toward the realization of an electrically pumped guided wave device. This has been done by adding layers where QDs are embedded in In{sub x}Al{sub y}Ga{sub 1−x−y}As layers, to obtain a symmetric structure from a band profile point of view. To assess the room temperature electro-luminescence emission spectrum under realistic electrical injection conditions, we performed device-level simulations based on a coupled drift-diffusion and QD rate equation model. On the basis of the device simulation results, we conclude that the present proposal is a viable option to realize broadband light-emitting devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

Report on simulation of fission gas and fission product diffusion in UO 2

In UO 2 nuclear fuel, the retention and release of fission gas atoms such as xenon (Xe) are important for nuclear fuel performance by, for example, reducing the fuel thermal conductivity, causing fuel swelling that leads to mechanical interaction with the clad, increasing the plenum pressure and reducing the fuel–clad gap thermal conductivity. We use multi-scale simulations to determine fission gas diffusion mechanisms as well as the corresponding rates in UO 2 under both intrinsic and irradiation conditions. In addition to Xe and Kr, the fission products Zr, Ru, Ce, Y, La, Sr and Ba have been investigated. Density functional theory (DFT) calculations are used to study formation, binding and migration energies of small clusters of Xe atoms and vacancies. Empirical potential calculations enable us to determine the corresponding entropies and attempt frequencies for migration as well as investigate the properties of large clusters or small fission gas bubbles. A continuum reaction-diffusion model is developed for Xe and point defects based on the mechanisms and rates obtained from atomistic simulations. Effective fission gas diffusivities are then obtained by solving this set of equations for different chemical and irradiation conditions using the MARMOT phase field code. The predictions are compared to available experimental data. The importance of the large Xe U3O cluster (a Xe atom in a uranium + oxygen vacancy trap site with two bound uranium vacancies) is emphasized, which is a consequence of its high mobility and high binding energy. We find that the Xe U3O cluster gives Xe diffusion coefficients that are higher for intrinsic conditions than under irradiation over a wide range of temperatures. Under irradiation the fast-moving Xe U3O cluster recombines quickly with irradiation-induced interstitial U ions, while this mechanism is less important for intrinsic conditions. The net result is higher concentration of the Xe U3O cluster for intrinsic conditions than under irradiation. We speculate that differences in the irradiation conditions and their impact on the Xe U3O cluster can explain the wide range of diffusivities reported in experimental studies. However, all vacancy-mediated mechanisms underestimate the Xe diffusivity compared to the empirical radiation-enhanced rate used in most fission gas release models. We investigate the possibility that diffusion of small fission gas bubbles or extended Xe-vacancy clusters may give rise to the observed radiation-enhanced diffusion coefficient. These studies highlight the importance of U divacancies and an octahedron coordination of uranium vacancies encompassing a Xe fission gas atom. The latter cluster can migrate via a multistep mechanism with a rather low effective barrier, which together with irradiation-induced clusters of uranium vacancies, gives rise to the irradiation-enhanced diffusion coefficient observed in experiments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Ohm’s Law, the Reconnection Rate, and Energy Conversion in Collisionless Magnetic Reconnection

Magnetic reconnection is a ubiquitous plasma process that transforms magnetic energy into particle energy during eruptive events throughout the universe. Reconnection not only converts energy during solar flares and geomagnetic substorms that drive space weather near Earth, but it may also play critical roles in the high energy emissions from the magnetospheres of neutron stars and black holes. In this review article, we focus on collisionless plasmas that are most relevant to reconnection in many space and astrophysical plasmas. Guided by first-principles kinetic simulations and spaceborne in-situ observations, we highlight the most recent progress in understanding this fundamental plasma process. We start by discussing the non-ideal electric field in the generalized Ohm’s law that breaks the frozen-in flux condition in ideal magnetohydrodynamics and allows magnetic reconnection to occur. We point out that this same reconnection electric field also plays an important role in sustaining the current and pressure in the current sheet and then discuss the determination of its magnitude (i.e., the reconnection rate), based on force balance and energy conservation. This approach to determining the reconnection rate is applied to kinetic current sheets with a wide variety of magnetic geometries, parameters, and background conditions. We also briefly review the key diagnostics and modeling of energy conversion around the reconnection diffusion region, seeking insights from recently developed theories. Finally, future prospects and open questions are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Phase-field model of alloy solidification far from chemical equilibrium at the solid-liquid interface

We further develop a recently introduced phase-field model of far-from-equilibrium alloy solidification under additive manufacturing conditions [K. Ji et al., Phys. Rev. Lett. 130, 026203 (2023)]. This model utilizes enhanced solute diffusivity within the spatially diffuse interface region to quantitatively capture solute trapping with a larger interface width, thereby making simulations on experimentally relevant length and timescales computationally feasible. The main developments presented here include testing the robustness of different variational formulations, extending the model to concentrated alloys by incorporating solid and liquid free energies from thermodynamic databases, as illustrated for hypoeutectic Al-Ag alloys with CALPHAD, extending convergence tests as a function of interface width to 3D, and carrying out simulations in both 2D and 3D to examine existing theories of microstructure development. Our results indicate that the simplest variational formulation that interpolates the bulk free-energy density between its solid and liquid forms is the most robust. Remarkably, for hypoeutectic Al-Ag alloys, this formulation yields a high-velocity nonequilibrium phase diagram that is independent of interface width, thereby demonstrating that the framework of enhanced solute diffusivity can be nontrivially extended to concentrated alloys. Other variational formulations have restricted ranges of materials or processing parameters that can be reliably modeled. We use 2D simulations to construct high-velocity microstructure selection maps for dilute Al-Cu alloys. The results validate the important role of latent heat rejection at the interface and extend the limited predictions of linear stability analysis [A. Karma and A. Sarkissian, Phys. Rev. E 47, 513 (1993)] and sharp-interface 1D simulations to fully nonlinear regimes. Furthermore, 3D simulations, carried out using a computationally tractable axisymmetric cellular/dendritic interface shape, demonstrate a good convergence similar to that observed in 2D as a function of interface width. Full 3D simulations, in turn, reveal that the standard theory of absolute stability is a good predictor of the upper critical velocity beyond which steady-state growth becomes unstable, despite the different morphological manifestations of this instability in 2D and 3D.

36 MATERIALS SCIENCE↗

DOE OSTI2017

Thermo-mechanical performance assessment of selected plates from MP-1 high power experiments

Monolithic plate-type fuel is a fuel form that is being developed for the conversion of high performance research and test reactors to low-enrichment uranium fuels. These fuel-plates are comprised of a high density, low enrichment, U-Mo alloy based fuel foil encapsulated in an aluminum cladding. To benchmark this new design, number of plates has been irradiated with satisfactory performance. As a part of continuing evaluation efforts, a set of plates covering range of operational parameters is scheduled to be tested during MP-1 irradiation experiments. It is necessary to evaluate the thermo-mechanical performance of plates during irradiation. For this, selected plates with distinct operational histories; covering low power, high power and high fission density were simulated. Fully coupled three-dimensional models of plates with a capability to evolve mechanical and thermal properties of constituent materials with irradiation time and burn-up were developed. The models input used projected parameters, including plate geometry, irradiation history and coolant conditions as input. The model output included temperature, displacement and stresses in the fuel, cladding and diffusion barrier. The fuel behavioral model considered inelastic behavior including volumetric swelling due to solid and gaseous products, irradiation induced creep, thermal expansion, conductivity degradation and plasticity. A visco-plastic behavioral model was used for the cladding that included thermal creep, irradiation hardening, growth due to fast neutrons and Mises plasticity. The plates were then simulated by using projected irradiation parameters. The resulting temperature, displacement and stress-strains were comparatively evaluated on the selected paths. The results were then compared with those of plates from previous RERTR experiments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2018

Thermo-mechanical performance assessment of selected plates from MP-1 low power experiments

42 ENGINEERING↗

DOE OSTI2022

Advection versus diffusion in Richtmyer-Meshkov mixing

The Richtmyer-Meshkov (RM) instability is one of the most severe degradation mechanisms for inertial confinement fusion (ICF), and mitigating it has been a priority for the global ICF effort. In this Letter, the instability's ability to atomically mix is linked to its background decay of residual turbulent energy. Here, we show how recently derived inequalities from the mathematical theory of PDEs constrain the evolution. A model RM process at leading order may diffusively mix or retain imprints of its initial structures indefinitely, depending on initial conditions, and there exists a theoretical range of zero-mixing for certain values of parameters. The results may apply to other systems resembling scalar transport in decaying turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Non-equilibrium simulations of hydraulic permeation: Role of mechanical boundary conditions in dense membranes

Understanding the mechanisms of water transport in reverse osmosis membranes is critical for improving membrane performance and guiding material design. While classical models describe transport as either solution-diffusion (SD)—involving concentration-driven diffusion through a homogeneous medium—or pore-flow (PF)—involving pressure-driven convection through percolated water channels—their applicability to crosslinked polyamide membranes remains debated. Here, using non-equilibrium molecular dynamics simulations, we investigate the impact of mechanical support conditions on pressure-driven water transport in polyamide membranes across varying crosslink densities and pressure differentials (1000–5000 bar). Two support conditions are considered: graphene-restrained, representing experimentally relevant supported membranes, and freeze-restrained, mimicking a self-supported structure. In graphene-restrained systems, water concentration gradients and constant pressure profiles emerge, consistent with SD theory and incompatible with PF assumptions due to the absence of percolated pores and sub-nanometer voids. In contrast, freeze-restrained systems display uniform water concentration and linearly decreasing pressure at 1000 bar, and exhibit compressibility-induced water gradients and partial percolation at 5000 bar, resembling PF-like behavior. However, the underlying assumptions of PF theory—continuous solvent pathways and pressure transmission through water-filled pores—are not met under most conditions. Our results demonstrate that accurate modeling of reverse osmosis membranes must incorporate realistic mechanical boundary conditions to distinguish between transport mechanisms. For dense polyamide membranes supported by porous substrates, graphene-restrained simulations best reflect experimental setups and support the SD model as the dominant mechanism of water permeation.

molecular dynamics↗

DOE OSTI2024

Vacancy-mediated transport and segregation tendencies of solutes in fcc nickel under diffusional creep: A density functional theory study

The Nabarro-Herring (NH) diffusional creep theory postulates the vacancy-mediated transport of atoms under a stress gradient as the creep mechanism under low-stress and high-temperature conditions. In multicomponent alloys, we premise that this stress-assisted flow of vacancies to and from grain boundaries will produce elemental segregation. An observation of such segregation, validated with theoretical predictions, can provide the necessary experimental evidence for the occurrence of NH creep. Theoretical calculations of the segregation tendencies via analyzing the dominant solute diffusion mechanisms and the difference in diffusivities of the elements are therefore essential. To this end, this study applies density functional theory calculations of migration barriers and solute-vacancy binding energies as input to the self-consistent mean-field theory to assess the vacancy-mediated diffusion mechanisms, transport coefficients, and segregation tendencies of Co, Cr, Mo, Re, Ta, and W solutes in face-centered-cubic Ni. We find Co, Re, and W to be slow diffusers at high temperatures and Cr, Mo, and Ta to be fast diffusers. Further analysis shows that the slow diffusers tend to always enrich at vacancy sinks over a wide range of temperatures. In contrast, the fast diffusers show a transition from depletion to enrichment as the temperature lowers. Furthermore, our analysis of the segregation tendencies under tensile hydrostatic strains shows that slow diffusers are largely unaffected by the strain and favor enrichment. On the other hand, the fast diffusers exhibit high sensitivity to strain and their segregation tendency can transition from depletion to enrichment at a given temperature. Finally, the transport coefficients calculated in this work are expected to serve as input to mesoscale microstructure models to provide a more rigorous assessment of solute segregation under NH creep conditions.

36 MATERIALS SCIENCE↗