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,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Fictitious phase separation in Li layered oxides driven by electro-autocatalysis

Layered oxides widely used as lithium-ion battery electrodes are designed to be cycled under conditions that avoid phase transitions. Although the desired single-phase composition ranges are well established near equilibrium, operando diffraction studies on many-particle porous electrodes have suggested phase separation during delithiation. Notably, the separation is not always observed, and never during lithiation. These anomalies have been attributed to irreversible processes during the first delithiation or reversible concentration-dependent diffusion. However, these explanations are not consistent with all experimental observations such as rate and path dependencies and particle-by-particle lithium concentration changes. Here, we show that the apparent phase separation is a dynamical artefact occurring in a many-particle system driven by autocatalytic electrochemical reactions, that is, an interfacial exchange current that increases with the extent of delithiation. We experimentally validate this population-dynamics model using the single-phase material L ix (Ni 1/3 Mn 1/3 Co 1/3 )O 2 (0.5 < x < 1) and demonstrate generality with other transition-metal compositions. Operando diffraction and nanoscale oxidation-state mapping unambiguously prove that this fictitious phase separation is a repeatable non-equilibrium effect. We quantitatively confirm the theory with multiple-datastream-driven model extraction. More generally, our study experimentally demonstrates the control of ensemble stability by electro-autocatalysis, highlighting the importance of population dynamics in battery electrodes (even non-phase-separating ones).

36 MATERIALS SCIENCE↗
DOE OSTI2021

Generation of mean flows in rotating anisotropic turbulence: The case of solar near-surface shear layer

Results from helioseismology indicate that the radial gradient of the rotation rate in the near-surface shear layer (NSSL) of the Sun is independent of latitude and radius. Theoretical models using the mean-field approach have been successful in explaining this property of the NSSL, while global direct or large-eddy magnetoconvection models have so far been unable to reproduce this. We investigate the reason for this discrepancy by measuring the mean flows, Reynolds stress, and turbulent transport coefficients under conditions mimicking those in the solar NSSL. Simulations with as few ingredients as possible to generate mean flows were studied. These ingredients are inhomogeneity due to boundaries, anisotropic turbulence, and rotation. The parameters of the simulations were chosen such that they matched the weakly rotationally constrained NSSL. The simulations probe locally Cartesian patches of the star at a given depth and latitude. The depth of the patch was varied by changing the rotation rate such that the resulting Coriolis numbers covered the same range as in the NSSL. We measured the turbulent transport coefficient relevant for the nondiffusive (Λ-effect) and diffusive (turbulent viscosity) parts of the Reynolds stress and compared them with predictions of current mean-field theories. A negative radial gradient of the mean flow is generated only at the equator where meridional flows are absent. At other latitudes, the meridional flow is comparable to the mean flow corresponding to differential rotation. We also find that the meridional components of the Reynolds stress cannot be ignored. Additionally, we find that the turbulent viscosity is quenched by rotation by about 50% from the surface to the bottom of the NSSL. Our local simulations do not validate the explanation for the generation of the NSSL from mean-field theory where meridional flows and stresses are neglected. However, the rotational dependence of the turbulent viscosity in our simulations agrees well with theoretical predictions. Moreover, our results agree qualitatively with global convection simulations in that an NSSL can only be obtained near the equator.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

The appearance, motion, and disappearance of three-dimensional magnetic null points

While theoretical models and simulations of magnetic reconnection often assume symmetry such that the magnetic null point when present is co-located with a flow stagnation point, the introduction of asymmetry typically leads to non-ideal flows across the null point. To understand this behavior, we present exact expressions for the motion of three-dimensional linear null points. The most general expression shows that linear null points move in the direction along which the magnetic field and its time derivative are antiparallel. Null point motion in resistive magnetohydrodynamics results from advection by the bulk plasma flow and resistive diffusion of the magnetic field, which allows non-ideal flows across topological boundaries. Null point motion is described intrinsically by parameters evaluated locally; however, global dynamics help set the local conditions at the null point. During a bifurcation of a degenerate null point into a null-null pair or the reverse, the instantaneous velocity of separation or convergence of the null-null pair will typically be infinite along the null space of the Jacobian matrix of the magnetic field, but with finite components in the directions orthogonal to the null space. Not all bifurcating null-null pairs are connected by a separator. Furthermore, except under special circumstances, there will not exist a straight line separator connecting a bifurcating null-null pair. The motion of separators cannot be described using solely local parameters because the identification of a particular field line as a separator may change as a result of non-ideal behavior elsewhere along the field line.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Limiting Surface Density Method Adapted to Large Arrays of Heterogeneous Shipping Packages with Nonlinear Responses

The classic Limiting Surface Density (LSD) method is an empirical calculation technique for analyzing and setting mass limits for fissile items in storage arrays. LSD is a desirable method because it can reduce or eliminate the need for lengthy detailed Monte Carlo models of storage arrays. The original (or classic) method was developed based on idealized arrays of bare spherical metal items in air-spaced cubic units in a water-reflected cubic array. In this case, the geometric and material-based surface densities were acceptably correlated by linear functions. Later updates to the method were made to allow for concrete reflection rather than water, cylindrical masses rather than spheres, different material forms, and noncubic arrays. However, in the intervening four decades since those updates, little work has been done to update the method, especially for use with contemporary highly heterogeneous shipping packages that are noncubic and stored in noncubic arrays. In this work, the LSD method is reevaluated for application to highly heterogeneous shipping packages for fissile material. The package modeled is the 9975 shipping package, currently the primary package used to store fissile material at Savannah River Site’s K-Area Complex. The package is neither cubic nor rectangular but resembles nested cylinders of stainless steel, lead, aluminum, and Celotex. The fissile content is assumed to be a cylinder of plutonium metal. The packages may be arranged in arrays with both an equal number of packages per side (package cubic) and an unequal number of packages per side (noncubic). The cubic arrangements are used to derive the 9975-specific material and geometry constants for the classic linear form LSD method. The linear form of the LSD, with noncubic array adjustment, is applied and evaluated against computational models for these packages to determine the critical unit fissile mass. Sensitivity equations are derived from the classic method, and these are also used to make projections of the critical unit fissile mass. It was discovered that the heterogeneous packages have a nonlinear surface density versus critical mass relationship compared to the acceptably linear response of bare spherical fissile masses. Methodology is developed to address the nonlinear response. In so doing, the solution to the nonlinear LSD method becomes decoupled from the critical mass of a single unit, adding to its flexibility. The ability of the method to predict changes in neutron multiplication due to perturbations in a parameter is examined to provide a basis for analyzing upset conditions. In conclusion, a full rederivation of the classic LSD method from diffusion theory is also included as this was found to be lacking in the available literature.

42 ENGINEERING↗
DOE OSTI2013

Zeolite Membrane Reactor for Water Gas Shift Reaction for Hydrogen Production

Gasification of biomass or heavy feedstock to produce hydrogen fuel gas using current technology is costly and energy-intensive. The technology includes water gas shift reaction in two or more reactor stages with inter-cooling to maximize conversion for a given catalyst volume. This project is focused on developing a membrane reactor for efficient conversion of water gas shift reaction to produce a hydrogen stream as a fuel and a carbon dioxide stream suitable for sequestration. The project was focused on synthesizing stable, hydrogen perm-selective MFI zeolite membranes for high temperature hydrogen separation; fabricating tubular MFI zeolite membrane reactor and stable water gas shift catalyst for membrane reactor applications, and identifying experimental conditions for water gas shift reaction in the zeolite membrane reactor that will produce a high purity hydrogen stream. The project has improved understanding of zeolite membrane synthesis, high temperature gas diffusion and separation mechanisms for zeolite membranes, synthesis and properties of sulfur resistant catalysts, fabrication and structure optimization of membrane supports, and fundamentals of coupling reaction with separation in zeolite membrane reactor for water gas shift reaction. Through the fundamental study, the research teams have developed MFI zeolite membranes with good perm-selectivity for hydrogen over carbon dioxide, carbon monoxide and water vapor, and high stability for operation in syngas mixture containing 500 part per million hydrogen sulfide at high temperatures around 500°C. The research teams also developed a sulfur resistant catalyst for water gas shift reaction. Modeling and experimental studies on the zeolite membrane reactor for water gas shift reaction have demonstrated the effective use of the zeolite membrane reactor for production of high purity hydrogen stream.

08 HYDROGEN↗
DOE OSTI2026

WellPINN: Accurate Well Representation for Transient Fluid Pressure Diffusion in Subsurface Reservoirs With Physics‐Informed Neural Networks

Accurate representation of pumping wells is essential for reliable reservoir characterization and simulation of operational scenarios in subsurface flow models. Physics-informed neural networks (PINNs) are emerging as a promising alternative to numerical models for reservoir modeling, offering seamless integration of monitoring data and governing physical equations. However, existing PINN-based studies face major challenges in capturing fluid pressure near wells when using a source/sink term, particularly during the early stages after pumping begins. We address this problem by introducing WellPINN, a workflow in which an initially trained PINN infers fluid pressure across the entire reservoir domain using a large equivalent well radius. This initial PINN solution is then locally refined around the well by a set of subdomain PINNs that are trained for smaller equivalent well radii. Continuity across these subdomain interfaces as well as at the initial condition is ensured by hard-constraining each PINN on its subdomain boundary. Our results demonstrate WellPINN as the first workflow of its kind to focus on accurate inference of fluid pressure from pumping rates throughout the entire injection period, significantly advancing the potential of PINNs for inverse modeling and operational scenario simulations. All data and code for this paper are openly available at https://doi.org/10.20350/DIGITALCSIC/17260.

58 GEOSCIENCES↗
DOE OSTI2019

Brownmillerite Phase Formation and Evolution in Epitaxial Strontium Ferrite Heterostructures

By controlling the synthesis and process conditions, we show that the oxygen vacancy channels (OVCs) of Brownmillerite-structured SrFeO2.5+? (0=??0.5) thin films can be oriented differently on the same LaAlO3(001) substrate due to either a bottom up or top down synthesis route. The brownmillerite structure formation and evolution processes yield distinctively different optical and electronic properties, owing to the difference in the oxygen content and OVC configuration in the resultant films. Moreover, we also show that in situ transmission electron microscopy (TEM) can induce oxygen loss and drive a phase transition from perovskite SrFeO3 to brownmillerite SrFeO2.5+? with mixed horizontally and vertically aligned OVCs. The dimensionality change during the TEM sample preparation makes the bottom up and top down model not applicable. Our results offer further insight into the phase stability and oxygen-diffusion mechanisms, which is important for the predictive synthesis of novel functional materials.

SrFeO3-d, compressive strain, oxygen vacancy chann↗
DOE OSTI2025

Automated Segmentation of Twin Boundaries in TRISO Silicon Carbide Using Deep Neural Networks

Coated particle fuels, such as the tristructural isotropic (TRISO) fuel particle, are essential for high-temperature gas reactor (HTGR) applications due to their efficiency and stability under normal and off-normal conditions. However, widespread commercialization and deployment of this technology for next-generation nuclear applications require robust quality assurance and quality control (QA/QC) methods linking fabrication, properties, and performance. Of the many important metrics for TRISO QA/QC, quantification of the silicon carbide (SiC) microstructure is critical because it correlates with fission product retention during irradiation. Previous work has shown extensive twinning of the SiC microstructure, which strongly affects microstructural metrics; however, twin grain boundaries are not expected play a significant role in fission product diffusion. This report summarizes the initial development, training, and testing of a machine learning image processing algorithm to detect twin grain boundaries in a backscattered electron image, which can be removed so that microstructural metrics can be recalculated for legacy data. Further development and deployment of this model will provide automated, scalable improvement of potential QA/QC methods for the SiC layer of TRISO particles.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

SHOCK CONNECTIVITY IN THE 2010 AUGUST AND 2012 JULY SOLAR ENERGETIC PARTICLE EVENTS INFERRED FROM OBSERVATIONS AND ENLIL MODELING

During periods of increased solar activity, coronal mass ejections (CMEs) can occur in close succession and proximity to one another. This can lead to the interaction and merger of CME ejecta as they propagate in the heliosphere. The particles accelerated in these shocks can result in complex solar energetic particle (SEP) events, as observing spacecraft form both remote and local shock connections. It can be challenging to understand these complex SEP events from in situ profiles alone. Multipoint observations of CMEs in the near-Sun environment, from the Solar Terrestrial Relations Observatory –Sun Earth Connection Coronal and Heliospheric Investigation and the Solar and Heliospheric Observatory Large Angle and Spectrometric Coronagraph, greatly improve our chances of identifying the origin of these accelerated particles. However, contextual information on conditions in the heliosphere, including the background solar wind conditions and shock structures, is essential for understanding SEP properties well enough to forecast their characteristics. Wang–Sheeley–Arge WSA-ENLIL + Cone modeling provides a tool to interpret major SEP event periods in the context of a realistic heliospheric model and to determine how much of what is observed in large SEP events depends on nonlocal magnetic connections to shock sources. We discuss observations of the SEP-rich periods of 2010 August and 2012 July in conjunction with ENLIL modeling. We find that much SEP activity can only be understood in the light of such models, and in particular from knowing about both remote and local shock source connections. These results must be folded into the investigations of the physics underlying the longitudinal extent of SEP events, and the source connection versus diffusion pictures of interpretations of SEP events.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Predictions of diffusion rates of large organic molecules in secondary organic aerosols using the Stokes–Einstein and fractional Stokes–Einstein relations

Information on the rate of diffusion of organic molecules within secondary organic aerosol (SOA) is needed to accurately predict the effects of SOA on climate and air quality. Diffusion can be important for predicting the growth, evaporation, and reaction rates of SOA under certain atmospheric conditions. Often, researchers have predicted diffusion rates of organic molecules within SOA using measurements of viscosity and the Stokes–Einstein relation ($D∝1/η$, where D is the diffusion coefficient and η is viscosity). However, the accuracy of this relation for predicting diffusion in SOA remains uncertain. Using rectangular area fluorescence recovery after photobleaching (rFRAP), we determined diffusion coefficients of fluorescent organic molecules over 8 orders in magnitude in proxies of SOA including citric acid, sorbitol,and a sucrose–citric acid mixture. These results were combined with literature data to evaluate the Stokes–Einstein relation for predicting the diffusion of organic molecules in SOA. Although almost all the data agree with the Stokes–Einstein relation within a factor of 10, a fractional Stokes–Einstein relation ($D∝1/η^ξ$) with ξ =0.93 is a better model for predicting the diffusion of organic molecules in the SOA proxies studied. In addition, based on the output from a chemical transport model, the Stokes–Einstein relation can overpredict mixing times of organic molecules within SOA by as much as 1 order of magnitude at an altitude of ~3 km compared to the fractional Stokes–Einstein relation with ξ =0.93. These results also have implications for other areas such as in food sciences and the preservation of biomolecules.

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

Interplay between Ion Transport, Applied Bias, and Degradation under Illumination in Hybrid Perovskite p-i-n Devices

In this work, we studied ion transport in hybrid organic–inorganic perovskite p-i-n devices as a function of applied bias under device operating conditions. Using electrochemical impedance spectroscopy (EIS) and equivalent circuit modeling, we elucidated various resistive and capacitive elements in the device. We show that ion migration is predictably influenced by a low applied forward bias, characterized by an increased capacitance at the hole-transporting (HTM) and electron-transporting material (ETM) interfaces, as well as in bulk. However, unlike observations in n-i-p devices, we found that there is a capacitive discharge leading to ion redistribution in the bulk at high forward biases. Furthermore, we show that a chemical double-layer capacitance buildup as a result of ion accumulation impacts the electronic properties of the device, likely by inducing either charge pinning or charge screening, depending on the direction of the ion-induced field. Finally, we extrapolate ion diffusion coefficients (~10 –7 cm 2 s –1 ) and ionic conductivities (~10 –7 S cm –1 ) from the Warburg mass (ion) diffusion response and show that, as the device degrades, there is an overall depletion of capacitive effects coupled with increased ion mobility.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Dissolution Behavior and Phase Evolution During Aluminum Oxide Dissolution in BOF Slag

Basic oxygen furnace (BOF) slag is a byproduct in the steel making process and is produced in considerable amounts. Thus, its valorization is important. Additives such as Al{sub 2}O{sub 3} can be utilized to modify the composition and mineralogy of the slag provided that a fast and effective dissolution of these additives in molten slags can be ensured. With this in mind, the dissolution of Al{sub 2}O{sub 3} particles in molten BOF slag has been investigated between 1500 and 1600 °C. In all cases, a dissolution zone forms around the dissolving particle which limits the dissolution rate in the melt under natural convection conditions. The dissolution in the current system is found to be an indirect process, wherein crystalline reaction products like CaAl{sub 4}O{sub 7} (CA{sub 2}) and CaAl{sub 2}O{sub 4} (CA) form on the surface of the particles. CA forms at 1500 °C and CA{sub 2} at 1600 °C. These reaction products form a continuous layer on the surface of the dissolving particles but the thickness of this reaction layer remains limited. The experimental data have been utilized to evaluate the dissolution paths. Near the Al{sub 2}O{sub 3}/slag interface, the inter-diffusion of Al{sub 2}O{sub 3} and CaO occurs due to which the dissolution path shifts towards the Al{sub 2}O{sub 3}-CaO binary system. A phenomenological model is presented to explain the evolution of the dissolution process in the current system, wherein the kinetics of the process is used to explain the reaction products formed at 1500 °C and 1600 °C.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Evaluating particle-suspension reactor designs for Z-scheme solar water splitting via transport and kinetic modeling

Sunlight-driven water splitting to produce hydrogen and oxygen provides a pathway to store available solar energy in the form of stable, energy-dense chemical bonds. Here we investigate a tandem particle-suspension reactor design for solar water splitting comprising micron-scale photocatalyst particles suspended in an aqueous solution with soluble redox shuttles. A porous separator facilitates redox species transport between the hydrogen and oxygen evolution reaction compartments while averting gas crossover. A two-dimensional, transient model of the reactor is presented to illustrate the coupling between light absorption, interfacial electron-transfer kinetics and species transport, and their combined impacts on overall solar-to-hydrogen conversion efficiency. The volumetric reactivity of the suspended semiconductor particles is dictated by combining the (photo)current-voltage behavior of a photodiode with Butler-Volmer electron-transfer kinetics. For the first time, a quantitative approach to determine the impacts of surface-dependent redox shuttle kinetic parameters on reaction selectivity in a Z-scheme system is established. Model results provide insights on the effects of optical, transport and kinetic properties of the semiconductor particles and the redox shuttles on the overall reactor performance. Solar-to-hydrogen reactor efficiencies predicted with BiVO 4 particles for oxygen evolution are at least two times larger than efficiencies achieved with wider band-gap TiO 2 particles due to enhanced visible light absorption; hydrogen evolution with SrTiO 3 :Rh particles was considered for both cases. Superior performance is predicted with proton-coupled electron transfer redox shuttles (para-benzoquinone/hydroquinone and iodide/iodate) that absorb little-to-no visible light, while also facilitating operation at near-neutral pH conditions, as compared to the non-proton-coupled triiodide/iodide and iron(iii)/iron(ii) redox shuttles. For 1 cm tall reaction compartments, diffusive species transport is fast enough to sustain reactor operation at a 1% solar-to-hydrogen conversion efficiency for both para-benzoquinone/hydroquinone and iodate/iodide redox shuttles with less than 2.2 mg L -1 of each of BiVO 4 and SrTiO 3 :Rh particles in the solution.

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

Dehydration of multilamellar fatty acid membranes: Towards a computational model of the stratum corneum

The level of hydration controls the cohesion between apposed lamellae of saturated free fatty acids found in the lipid matrix of stratum corneum, the outermost layer of mammalian skin. This multilamellar lipid matrix is highly impermeable to water and ions, so that the local hydration shell of its fatty acids may not always be in equilibrium with the acidity and relative humidity, which significantly change over a course of days during skin growth. The homeostasis of the stratum corneum at each moment of its growth likely requires a balance between two factors, which affect in opposite ways the diffusion of hydrophilic species through the stratum corneum: (i) an increase in water order as the lipid lamellae come in closer contact, and (ii) a decrease in water order as the fraction of charged fatty acids is lowered by pH. Herein molecular dynamics simulations are employed to estimate the impact of both effects on water molecules confined between lamellae of fatty acids. Under conditions where membrane undulations are energetically favorable, the charged fatty acids are able to sequester cations around points of contact between lamellae that are fully dehydrated, while essentially maintaining a multilamellar structure for the entire system. This observation suggests that the undulations of the fatty acid lamellae control the diffusion of hydrophilic species through the water phase by altering the positional and rotational order of water molecules in the embedded/occluded “droplets.”

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

Stochastic and mixed density functional theory within the projector augmented wave formalism for simulation of warm dense matter

Stochastic density functional theory (DFT) and mixed stochastic-deterministic DFT are burgeoning approaches for the calculation of the equation of state and transport properties in materials under extreme conditions. In the intermediate warm dense matter regime, a state between correlated condensed matter and kinetic plasma, electrons can range from being highly localized around nuclei to delocalized over the whole simulation cell. The plane-wave basis pseudopotential approach is thus the typical tool of choice for modeling such systems at the DFT level. Unfortunately, stochastic DFT methods scale as the square of the maximum plane-wave energy in this basis. To reduce the effect of this scaling and improve the overall description of the electrons within the pseudopotential approximation, we present stochastic and mixed DFT approaches developed and implemented within the projector augmented wave formalism. In conclusion, we compare results between the different DFT approaches for both single-point and molecular dynamics trajectories and present calculations of self-diffusion coefficients of solid density carbon from 1 to 50 eV.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Numerical Modeling and Dynamic Analysis of a Wave-Powered Reverse-Osmosis System

A wave energy converter (WEC) system has the potential to convert the wave energy resource directly into the high-pressure flow that is needed by the desalination system to pump saltwater to the reverse-osmosis membrane and provide the required pressure level to generate freshwater. In this study, a wave-to-water numerical model was developed to investigate the potential use of a wave-powered desalination system (WPDS) for water production. The model was developed by coupling a time-domain radiation-and-diffraction method-based numerical tool (WEC-Sim) for predicting the hydrodynamic performance of WECs with a solution-diffusion model that was used to simulate the reverse-osmosis (RO) process. The objective of this research is to evaluate the WPDS dynamics and the overall efficiency of the system. To evaluate the feasibility of the WPDS, the wave-to-water numerical model was applied to simulate a desalination system that used an oscillating surge WEC device to pump seawater through the system. The hydrodynamics WEC-Sim simulation results for the oscillating surge WEC device were validated against existing experimental data. The RO simulation was verified by comparing the results to those from the Dow Chemical Company’s reverse osmosis system analysis (ROSA) model, which has been widely used to design and simulate RO systems. The wave-to-water model was then used to analyze the WPDS under a range of wave conditions and for a two-WECs-coupled RO system to evaluate the influence of pressure and flow rate fluctuation on the WPDS performance. The results show that the instantaneous energy fluctuation from waves has a significant influence on the responding hydraulic pressure and flow rate, as well as the recovery ratio and, ultimately, the water-production quality. Nevertheless, it is possible to reduce the hydraulic fluctuation for different sea states while maintaining a certain level of freshwater production, and a WEC array that produces water can be a viable, near-term solution to the nation’s water supply. A discussion on the dynamic impact of hydraulic fluctuation on the WPDS performance and potential options to reduce the fluctuation and their trade-offs is also presented.

16 TIDAL AND WAVE POWER↗
DOE OSTI2013

Comparison of actinide production in traveling wave and pressurized water reactors

The geopolitical problems associated with civilian nuclear energy production arise in part from the accumulation of transuranics in spent nuclear fuel. A traveling wave reactor is a type of breed-burn reactor that could, if feasible, reduce the overall production of transuranics. In one possible configuration, a cylinder of natural or depleted uranium would be subjected to a fast neutron flux at one end. The neutrons would transmute the uranium, producing plutonium and higher actinides. Under the right conditions, the reactor could become critical, at which point a self-stabilizing fission wave would form and propagate down the length of the reactor cylinder. The neutrons from the fission wave would burn the fissile nuclides and transmute uranium ahead of the wave to produce additional fuel. Fission waves in uranium are driven largely by the production and fission of {sup 239}Pu. Simulations have shown that the fuel burnup can reach values greater than 400 MWd/kgIHM, before fission products poison the reaction. In this work we compare the production of plutonium and minor actinides produced in a fission wave to that of a UOX fueled light water reactor, both on an energy normalized basis. The nuclide concentrations in the spent traveling wave reactor fuel are computed using a one-group diffusion model and are verified using Monte Carlo simulations. In the case of the pressurized water reactor, a multi-group collision probability model is used to generate the nuclide quantities. We find that the traveling wave reactor produces about 0.187 g/MWd/kgIHM of transuranics compared to 0.413 g/MWd/kgIHM for a pressurized water reactor running fuel enriched to 4.95 % and burned to 50 MWd/kgIHM. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2014

Diffusion-controlled reactions modeling in Geant4-DNA

Context Under irradiation, a biological system undergoes a cascade of chemical reactions that can lead to an alteration of its normal operation. There are different types of radiation and many competing reactions. As a result the kinetics of chemical species is extremely complex. The simulation becomes then a powerful tool which, by describing the basic principles of chemical reactions, can reveal the dynamics of the macroscopic system. To understand the dynamics of biological systems under radiation, since the 80s there have been on-going efforts carried out by several research groups to establish a mechanistic model that consists in describing all the physical, chemical and biological phenomena following the irradiation of single cells. This approach is generally divided into a succession of stages that follow each other in time: (1) the physical stage, where the ionizing particles interact directly with the biological material; (2) the physico-chemical stage, where the targeted molecules release their energy by dissociating, creating new chemical species; (3) the chemical stage, where the new chemical species interact with each other or with the biomolecules; (4) the biological stage, where the repairing mechanisms of the cell come into play. This article focuses on the modeling of the chemical stage. Method This article presents a general method of speeding-up chemical reaction simulations in fluids based on the Smoluchowski equation and Monte-Carlo methods, where all molecules are explicitly simulated and the solvent is treated as a continuum. The model describes diffusion-controlled reactions. This method has been implemented in Geant4-DNA. The keys to the new algorithm include: (1) the combination of a method to compute time steps dynamically with a Brownian bridge process to account for chemical reactions, which avoids costly fixed time step simulations; (2) a k–d tree data structure for quickly locating, for a given molecule, its closest reactants. The performance advantage is presented in terms of complexity, and the accuracy of the new algorithm is demonstrated by simulating radiation chemistry in the context of the Geant4-DNA project. Application The time-dependent radiolytic yields of the main chemical species formed after irradiation are computed for incident protons at different energies (from 50 MeV to 500 keV). Both the time-evolution and energy dependency of the yields are discussed. The evolution, at one microsecond, of the yields of hydroxyls and solvated electrons with respect to the linear energy transfer is compared to theoretical and experimental data. According to our results, at high linear energy transfer, modeling radiation chemistry in the trading compartment representation might be adopted.

60 APPLIED LIFE SCIENCES↗