Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion problems”

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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Floating Island International PHASE I FINAL TECHNICAL REPORT (Advancing the Development of Floating Solar-Powered Nanobubble Aeration Systems for Use with Floating Treatment Wetlands in Natural and Man-Made Waterbodies)

More and more freshwater lakes are suffering from algae blooms, which deprive water of oxygen and lead to widespread loss of aquatic life and production of dangerous toxins. When oxygen is lacking, methane is generated in the sediments as algae is decomposed; it has been calculated that more than half of global methane emissions come from nutrient-impaired freshwater systems (Beaulieu, 2021). Nutrients, mainly from agricultural run-off, feed algae blooms and are exacerbated by warming related to climate change. Artificial aeration is frequently used to restore oxygen to a waterbody. But diffuser aeration is inefficient and expensive, needs on-site grid power, and is failing to keep pace with the demands of water as it gets warmer with climate change. In the last five years, nanobubble aeration has been introduced into freshwater applications and shows promise for rapidly increasing dissolved oxygen effectively throughout the water column. Nanobubbles deliver oxygen in bubbles that have no buoyancy, so they stay in water longer and release their oxygen more fully, compared to large bubbles that rise and burst at the surface. The prospect of this new technology answering the deficiencies of diffuser aeration drove us to initiate the current project. During our Phase I period, we have successfully tested a nanobubble aeration system that can super-saturate oxygen levels. It is operated on solar power, and both the nanobubbler and solar array are mounted on a proprietary floating island platform, that also performs biological nutrient recycling. A rudimentary system was assembled and tested on a 6.5-acre research lake at our headquarters in Montana, where it was subjected to a summer drought that reduced water levels significantly, periods of severe cold in winter, and spring conditions that included heavy rainfall and violent thunderstorms. Our twice-weekly sampling throughout the project showed that dissolved oxygen levels rose rapidly and were maintained throughout winter and spring, well into June. The nanobubbler was able to run on solar power for long periods. Methane levels were tested periodically and found to decrease as oxygen increased. The nanobubbles appeared to have no adverse impact on fish or other aquatic life exposed to them. The stability of the installation survived the weather conditions. The many problems we encountered taught us what we need to improve in the next iteration. Towards the end of our Phase I, we were awarded a supplementary state grant that enabled us to acquire a new nanobubble system for testing that runs on DC, is extremely efficient and has few moving parts. We plan to team this with low-profile solar panels, lithium batteries and a controller that provides real-time data. The ultimate goal of this project is to commercialize an affordable and scalable lake management system that will oxygenate water from the surface down to the sludge. It will enliven fisheries, prevent toxic algae blooms and – most importantly for the fate of our planet - inhibit the production and release of methane from the sediments. We anticipate that as carbon credits expand to include methane, the cost of remediating waterbodies that produce methane will be offset. We firmly believe that this is a practical technology that works and will make a big difference when fully implemented.

14 SOLAR ENERGY↗
DOE OSTI2017

Helium interactions with alumina formed by atomic layer deposition show potential for mitigating problems with excess helium in spent nuclear fuel

Helium gas accumulation from alpha decay during extended storage of spent fuel has potential to compromise the structural integrity the fuel. We report results obtained with surrogate nickel particles which suggest that alumina formed by atomic layer deposition can serve as a low volume-fraction, uniformly-distributed phase for retention of helium generated in fuel particles such as uranium oxide. Thin alumina layers may also form transport paths for helium in the fuel rod, which would otherwise be impermeable. Micron-scale nickel particles, representative of uranium oxide particles in their low helium solubility and compatibility with the alumina synthesis process, were homogeneously coated with alumina approximately 3–20 nm by particle atomic layer deposition (ALD) using a fluidized bed reactor. Particles were then loaded with helium at 800 °C in a tube furnace. Subsequent helium spectroscopy measurements showed that the alumina phase, or more likely a related nickel/alumina interface structure, retains helium at a density of at least 10 17 atoms/cm 3 . High resolution transmission electron microscopy revealed that the thermal treatment increased the alumina thickness and generated additional porosity. Results from Monte Carlo simulations on amorphous alumina predict the helium retention concentration at room temperature could reach 10 21 atoms/cm 3 at 400 MPa, a pressure predicted by others to be developed in uranium oxide without an alumina secondary phase. This concentration is sufficient to eliminate bubble formation in the nuclear fuel for long-term storage scenarios, for example. Finally, measurements by others of the diffusion coefficient in polycrystalline alumina indicate values several orders of magnitude higher than in uranium oxide, which then can also allow for helium transport out of the spent fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Extended hierarchical solvent perturbations from curved surfaces of mesoporous silica particles in a deep eutectic solvent

Here, many applications of deep eutectic solvents (DES) rely on exploitation of their unique yet complex liquid structures. Due to the ionic nature of the DES components, their diffuse structures are perturbed in the presence of a charged surface. We hypothesize that it is possible to perturb the bulk DES structure far (> 100 nm) from a curved, charged surface with mesoscopic dimensions. We performed in situ, synchrotron-based ultra-small angle X-ray scattering (USAXS) experiments to study the solvent distribution near the surface of charged mesoporous silica particles (MPS) (≈ 0.5 µm in diameter) suspended in both water and a common type of DES (1:2 choline Cl-:ethylene glycol). A careful USAXS analysis reveals that the perturbation of electron density distribution within the DES extends ≈ 1 μm beyond the particle surface, and that this perturbation can be manipulated by the addition of salt ions (AgCl). The concentration of the pore-filling fluid is greatly reduced in the DES. Notably, we extracted the real-space structures of these fluctuations from the USAXS data using a simulated annealing approach that does not require a priori knowledge about the scattering form factor, and can be generalized to a wide range of complex small-angle scattering problems.

36 MATERIALS SCIENCE↗
DOE OSTI2017

grim: A Flexible, Conservative Scheme for Relativistic Fluid Theories

Hot, diffuse, relativistic plasmas such as sub-Eddington black-hole accretion flows are expected to be collisionless, yet are commonly modeled as a fluid using ideal general relativistic magnetohydrodynamics (GRMHD). Dissipative effects such as heat conduction and viscosity can be important in a collisionless plasma and will potentially alter the dynamics and radiative properties of the flow from that in ideal fluid models; we refer to models that include these processes as Extended GRMHD. Here we describe a new conservative code, grim, that enables all of the above and additional physics to be efficiently incorporated. grim combines time evolution and primitive variable inversion needed for conservative schemes into a single step using an algorithm that only requires the residuals of the governing equations as inputs. This algorithm enables the code to be physics agnostic as well as flexibility regarding time-stepping schemes. grim runs on CPUs, as well as on GPUs, using the same code. We formulate a performance model and use it to show that our implementation runs optimally on both architectures. grim correctly captures classical GRMHD test problems as well as a new suite of linear and nonlinear test problems with anisotropic conduction and viscosity in special and general relativity. As tests and example applications, we resolve the shock substructure due to the presence of dissipation, and report on relativistic versions of the magneto-thermal instability and heat flux driven buoyancy instability, which arise due to anisotropic heat conduction, and of the firehose instability, which occurs due to anisotropic pressure (i.e., viscosity). Finally, we show an example integration of an accretion flow around a Kerr black hole, using Extended GRMHD.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Optimization of electrical treatment strategy for surface roughness reduction in conducting thin films

The surface roughness of deposited conducting thin films is responsible for various materials reliability problems in nanoelectronics and nanofabrication technologies. In this work, we report a modeling and simulation study that aims at optimizing the electrical surface treatment of deposited conducting thin films as a physical processing strategy for their surface roughness reduction. Our research is based on a continuum model of film surface morphological evolution that accounts for the residual stress in the deposited conducting thin film, the film’s wetting of the substrate layer that it is deposited on, film texture and surface diffusional anisotropy, and surface electromigration. Through systematic linear stability analysis and dynamical simulation protocols, we examine in detail the effects of film surface crystallographic orientation and applied electric field direction toward minimizing the electric field strength required for film surface smoothening. We believe that the critical electric field strength requirement for surface roughness reduction on {110}, {100}, and {111} surfaces of face-centered cubic crystalline conducting thin films exhibits a very strong dependence on the applied electric field direction, expressed as the electric field misalignment with respect to the principal residual stress directions in the film and the fast surface diffusion directions. Due to these findings, we optimize the electrical treatment strategy for surface roughness reduction of conducting thin films with respect to all relevant processing and material parameters.

42 ENGINEERING↗
DOE OSTI2020

Choosing the Best Modeling Platform for Radiological Risk Assessment Models - 20468

Radiological risk assessments, in the form of performance or safety assessments, are often required under regulations or guidance for remediation of contaminated land, decommissioning of contaminated buildings or structures, and radioactive waste disposal. These risk assessments are usually supported by fate and transport models that address decay and ingrowth of radionuclides, as well as their movement through engineered systems and the natural environment. These models are often projected thousands, or more, years into the future, largely because the radioactive species change through decay and ingrowth, and hence the magnitude of the radioactive effect changes with time. There are many computer codes that are available to address this type of modeling. They range from addressing specific pathways or processes such as infiltration of water, groundwater, surface water, air, biota, diffusion and advection of water and gases, to those that try to couple all processes together to evaluate the impact of fate and transport through the entire system to places in space and time to which access is assumed. These different types of codes are sometimes separated with the monikers process-level and systems-level codes, although it is often not clear that this separation does justice to the capabilities of the many codes that are available to evaluate fate and transport of radionuclides. The focus of this paper is the latter group of modeling codes. Several systems level modeling codes exist and are used. There are differences between these codes in terms of utility, flexibility, complexity and cost. The purpose of this paper is to compare a few of these codes in the context of work currently being performed by the International Atomic Energy Agency (IAEA) Modeling and Data for Radiological Impact Assessments (MODARIA) II Working Group 1 (WG1). The MODARIA II WG1's main focus is how stakeholder engaged decision analysis can, or should, be applied to radiological contamination problems so that better, longstanding, sustainable, solutions are reached. However, the WG1 also recognizes the potential impact of the modeling tools that are chosen to address radiological risk, which is often a primary objective of decision making for radiological problems. Other objectives might also be important, such as constraining costs, obtaining financing, minimizing impact on ecosystems, saving cultural resources, saving jobs, farmland, environmental justice, etc., in a full decision analysis for a given radiological contamination problem, but none of these other objectives have the same types of complex modeling needs as minimize radiological dose. Consequently, a further focus of WG1 is to evaluate the potential impacts on decision making of the choice of fate and transport, and risk assessment, modeling codes that are used to support decision making. The WG1 will produce a report at the end of 2020 that will focus on an approach to effective decision making and stakeholder engagement. The report will also consider the role that performance assessment modeling should play in the decision-making process, including the impact of the choice of modeling tools or computer codes on risk-informed decision making. Several sites around the World have been made available by Member States for these model comparisons, and several modeling tools have been considered. However, the focus of this paper is on two of the sites, one in Belgium and one in Ukraine, and on three of the tools: NORMALYSA (NORM And Legacy Site Assessment); GoldSim{sup C}, and AMBER{sup C}. The final report from this working group will also cover other modeling tools, including RESRAD, and PC-Cream{sup R}. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Comparison of momentum transport models for numerical relativity

The main problems of nonvacuum numerical relativity, compact binary mergers and stellar collapse, involve hydromagnetic instabilities and turbulent flows, so that kinetic energy at small scales leads to mean effects at large scale that drive the secular evolution. Notable among these effects is momentum transport. We investigate two models of this transport effect, a relativistic Navier-Stokes system and a turbulent mean stress model, that are similar to all of the prescriptions that have been attempted to date for treating subgrid effects on binary neutron star mergers and their aftermath. Our investigation involves both stability analysis and numerical experimentation on star and disk systems. We also begin the investigation of the effects of particle and heat transport on postmerger simulations. Here, we find that correct handling of turbulent heating is crucial for avoiding unphysical instabilities. Given such appropriate handling, the evolution of a differentially rotating star and the accretion rate of a disk are reassuringly insensitive to the choice of prescription. However, disk outflows can be sensitive to the choice of method, even for the same effective viscous strength. We also consider the effects of eddy diffusion in the evolution of an accretion disk and show that it can interestingly affect the composition of outflows.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Bayesian learning with Gaussian processes for low-dimensional representations of time-dependent nonlinear systems

This work presents a data-driven method for learning low-dimensional time-dependent physics-based surrogate models whose predictions are endowed with uncertainty estimates. We use the operator inference approach to model reduction that poses the problem of learning low-dimensional model terms as a regression of state space data and corresponding time derivatives by minimizing the residual of reduced system equations. Standard operator inference models perform well with accurate training data that are dense in time, but producing stable and accurate models when the state data are noisy and/or sparse in time remains a challenge. Another challenge is the lack of uncertainty estimation for the predictions from the operator inference models. Our approach addresses these challenges by incorporating Gaussian process surrogates into the operator inference framework to (1) probabilistically describe uncertainties in the state predictions and (2) procure analytical time derivative estimates with quantified uncertainties. The formulation leads to a generalized least-squares regression and, ultimately, reduced-order models that are described probabilistically with a closed-form expression for the posterior distribution of the operators. The resulting probabilistic surrogate model propagates uncertainties from the observed state data to reduced-order predictions. Furthermore, we demonstrate the method is effective for constructing low-dimensional models of two nonlinear partial differential equations representing a compressible flow and a nonlinear diffusion–reaction process, as well as for estimating the parameters of a low-dimensional system of nonlinear ordinary differential equations representing compartmental models in epidemiology.

Data-driven model reduction↗
DOE OSTI2025

Nonlinear Ensemble Filtering with Diffusion Models: Application to the Surface Quasigeostrophic Dynamics

The intersection between classical data assimilation methods and novel machine learning techniques has attracted significant interest in recent years. Here, we explore another promising solution in which diffusion models are used to formulate a robust nonlinear ensemble filter for sequential data assimilation. Unlike standard machine learning methods, the proposed ensemble score filter (EnSF) is completely training free and can efficiently generate a set of analysis ensemble members. Here, in this study, we apply the EnSF to a surface quasigeostrophic model and compare its performance against the popular local ensemble transform Kalman filter (LETKF), which makes Gaussian assumptions in the analysis step. Numerical tests demonstrate that EnSF maintains stable performance in the absence of localization and for a variety of experimental settings. We find that while LETKF maintains optimal performance in the case of linear observations of the entire state and a perfect model, EnSF shows improvements over LETKF when nonlinear observations are assimilated and the system is subject to unexpected model errors. A spectral decomposition of the analysis results in this nonlinear observation regime shows that the largest improvements over LETKF occur at large scales (small wavenumbers), where LETKF lacks sufficient ensemble spread. Overall, this initial application of EnSF to a geophysical model of intermediate complexity motivates further development of the algorithm for more realistic problems.

Artificial intelligence↗
DOE OSTI2014

Rigidity and resistance of larval- and adult schistosomes-medium interface

Graphical abstract: - Highlights: • Schistosoma larvae and worms are studied by neutron scattering. • Measurements on larvae were repeated after one day and by increasing temperature. • The flexibility properties of larvae and adult parasites are compared. • The parasite rigidity is related to their resistance to the hostile environment. • Insight into the parasite defense mechanisms to the immune system attack is achieved. - Abstract: Schistosomiasis is second only to malaria in prevalence and severity, and is still a major health problem in many tropical countries worldwide with about 200–300 million cases and with more than 800 million people at risk of infection. Based on these data, the World Health Organization recommends fostering research efforts for understanding at any level the mechanisms of the infection and then decreasing the social and economical impact of schistosomiasis. A key role is played by the parasite apical lipid membrane, which is entirely impervious to the surrounding elements of the immune system. We have previously demonstrated that the interaction between schistosomes and surrounding medium is governed by a parasite surface membrane sphingomyelin-based hydrogen barrier. In the present article, the elastic contribution to the total motion as a function of the exchanged wave-vector Q and the mean square displacement values for Schistosoma mansoni larvae and worms and Schistosomahaematobium worms have been evaluated by quasi elastic neutron scattering (QENS). The results point out that S. mansoni larvae show a smaller mean square displacement in comparison to S. mansoni and S. haematobium worms. These values increased by repeating the measurements after one day. These differences, which are analogous to those observed for the diffusion coefficient we previously evaluated, are interpreted in terms of rigidity of the parasite-medium interaction. S. mansoni larvae are the most rigid systems, while S. haematobium worms are the most flexible. In addition, temperature and hypoxia induce a weakening of the schistosome-medium interaction. These evidences are related to the strength of the hydrogen-bonded interaction between parasites and environment that we previously determined. We have shown that S. mansoni worms are characterized by a weakened interaction in respect to the larvae, while the S. haematobium worms more weakly interact with the surrounding medium than S. mansoni. The present QENS analysis allowed us to characterize the rigidity of larval- and adult S. mansoni and S. haematobium-host interface and to relate it to the parasite resistance to the hostile elements of the surrounding medium and to the immune effectors attack.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Absorption of ac fields in amorphous indium-oxide films

Absorption data from applied ac fields in Anderson-localized amorphous indium-oxide (In{sub x}O) films are shown to be frequency and disorder dependent. The absorption shows a roll-off at a frequency which is much lower than the electron-electron scattering rate of the material when it is in the diffusive regime. This is interpreted as evidence for discreteness of the energy spectrum of the deeply localized regime. This is consistent with recent many-body localization scenarios. As the metal-insulator transition is approached, the absorption shifts to higher frequencies. Comparing with the previously obtained results on the crystalline version of indium-oxide (In{sub 2}O{sub 3−x}) implies a considerably higher inelastic electron-phonon scattering rate in the amorphous material. The range over which the absorption versus frequency decreases may indicate that a wide distribution of localization length is a common feature in these systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Local nonlinear rf forces in inhomogeneous magnetized plasmas

The local nonlinear forces induced by radio frequency (rf) waves are derived in inhomogeneous magnetized plasmas, where the inhomogeneity exists in the rf fields, in the static magnetic field as well as in the equilibrium density and temperature. The local parallel force is completely resonant, but a novel component dependent on those inhomogeneities is obtained as the result of the inhomogeneous transport of parallel resonant-absorbed momentum by the nonlinear perpendicular drift flux. In the local poloidal force, the component induced by the inhomogeneity of rf power absorption is also confirmed and it can be recognized as the residual effect from the incomplete cancellation between the rate of the diamagnetic poloidal momentum gain and the Lorentz force due to the radial diffusion-like flux. The compact expression for radial force is also obtained for the first time, whose nonresonant component is expressed as the sum of the ponderomotive force on particles and the gradients of the nonresonant perpendicular pressure and of the nonresonant momentum flux due to the finite temperature effect. Numerical calculations in a 1-D slab model show that the resonant component dependent on the inhomogeneities may be significant when the ion absorption dominates the resonant wave-particle interaction. A quantitative estimation shows that the novel component in the parallel force is important to understand the experiments of the ion-cyclotron-frequency mode-conversion flow drive.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Three species one-dimensional kinetic model for weakly ionized plasmas

A three species one-dimensional kinetic model is presented for a spatially homogeneous weakly ionized plasma subjected to the action of a time varying electric field. Planar geometry is assumed, which means that the plasma evolves in the privileged direction of the field. The energy transmitted to the electric charges is channelized to the neutrals thanks to collisions, a mechanism that influences the plasma dynamics. Charge-charge interactions have been designed as a one-dimensional collision term equivalent to the Landau operator used for fully ionized plasmas. Charge-neutral collisions are modelled by a conservative drift-diffusion operator in the Dougherty's form. The resulting set of coupled integro-differential equations is solved with the stable and robust propagator integral method. This semi–analytical method feasibility accounts for non–linear effects without appealing to linearisation or simplifications, providing conservative physically meaningful solutions even for initial or emerging sharp velocity distribution function profiles. It is found that charge-neutral collisions exert a significant effect since a quite different plasma evolution arises if compared to the collisionless limit. In addition, substantial differences in the system motion are found for constant and temperature dependent collision frequencies cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

AGR-3/4 Irradiation Test Train Disassembly and Component Metrology First Look Report

The AGR-3/4 experiment was designed to study fission product transport within graphitic matrix material and nuclear-grade graphite. To this end, this experiment consisted of 12 capsules, each fueled with 4 compacts containing UCO TRISO particles as driver fuel and 20 UCO designed-to-fail (DTF) fuel particles in each compact. The DTF fuel was fabricated with a thin pyrocarbon layer which was intended to fail during irradiation and provide a source of fission products. These fission products could then migrate through the compact and into the surrounding concentric rings of graphitic matrix material and/or nuclear graphite. Through post-irradiation examination (PIE) of the rings (including physical sampling and gamma scanning) fission product concentration profiles within the rings can be determined. These data can be used to elucidate fission product transport parameters (e.g. diffusion coefficients within the test materials) which will be used to inform and refine models of fission product transport. After irradiation in the Advanced Test Reactor (ATR) had been completed in April 2014, the AGR-3/4 experiment was shipped to the Hot Fuel Examination Facility (HFEF) at the Materials and Fuels Complex (MFC) for inspection, disassembly, and metrology. The AGR-3/4 test train was received at MFC in two separate shipments between February and April 2015. Visual examinations of the test train exterior did not indicate dimensional distortion, and only two small discolored areas were observed at the bottom of Capsules 8 and 9. No corresponding discoloration was found on the inside of these capsules, however. Prior to disassembly, the two test train sections were subject to analysis via the Precision Gamma Scanner (PGS), which did not indicate that any gross fuel relocation had occurred. A series of specialized tools (including clamps, cutters, and drills) had been designed and fabricated in order to carry out test train disassembly and recovery of capsule components (graphite rings and fuel compacts). This equipment performed well for separating each capsule in the test train and extracting the capsule components. Only a few problems were encountered. In one case, the outermost ring (the sink ring) was cracked during removal of the capsule through tubes. Although the sink ring will be analyzed in order to obtain a mass balance of fission products in the experiment, these cracks do not pose a major concern because the sink ring will not be analyzed in detail to obtain the spatial distribution of fission products. In Capsules 4 and 5, the compacts could not be removed from the inner rings. Strategies for removing the compacts are being evaluated. Sampling the inner rings with the compacts in-place is also an option. Dimensional measurements were made on the compacts, inner rings, outer rings, and sink rings. The diameters of all compacts decreased by 0.5 to 2.0 %. Generally, the extent of diametric shrinkage increased linearly with increasing neutron fluence. Most compact lengths also decreased. Compact lengths decreased with increasing fluence, reaching maximum shrinkage of about 0.9 % at a fast fluence of 4.0x10 25 n/m 2 E > 0.18 MeV. Above this fluence, the extent of length shrinkage appeared to decrease with fluence, and two compacts from Capsule 7 were found to have slightly increased in length (< 0.1 %) after a fluence of 5.2x10 25 n/m 2 .

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

AGR-3/4 Irradiation Test Train Disassembly and Component Metrology First Look Report

The AGR-3/4 experiment was designed to study fission product transport within graphitic matrix material and nuclear-grade graphite. To this end, this experiment consisted of 12 capsules, each fueled with 4 compacts containing UCO TRISO particles as driver fuel and 20 UCO designed-to-fail (DTF) fuel particles in each compact. The DTF fuel was fabricated with a thin pyrocarbon layer which was intended to fail during irradiation and provide a source of fission products. These fission products could then migrate through the compact and into the surrounding concentric rings of graphitic matrix material and/or nuclear graphite. Through post-irradiation examination (PIE) of the rings (including physical sampling and gamma scanning) fission product concentration profiles within the rings can be determined. These data can be used to elucidate fission product transport parameters (e.g. diffusion coefficients within the test materials) which will be used to inform and refine models of fission product transport. After irradiation in the Advanced Test Reactor (ATR) had been completed in April 2014, the AGR-3/4 experiment was shipped to the Hot Fuel Examination Facility (HFEF) at the Materials and Fuels Complex (MFC) for inspection, disassembly, and metrology. The AGR-3/4 test train was received at MFC in two separate shipments between February and April 2015. Visual examinations of the test train exterior did not indicate dimensional distortion, and only two small discolored areas were observed at the bottom of Capsules 8 and 9. No corresponding discoloration was found on the inside of these capsules, however. Prior to disassembly, the two test train sections were subject to analysis via the Precision Gamma Scanner (PGS), which did not indicate that any gross fuel relocation had occurred. A series of specialized tools (including clamps, cutters, and drills) had been designed and fabricated in order to carry out test train disassembly and recovery of capsule components (graphite rings and fuel compacts). This equipment performed well for separating each capsule in the test train and extracting the capsule components. Only a few problems were encountered. In one case, the outermost ring (the sink ring) was cracked during removal of the capsule through tubes. Although the sink ring will be analyzed in order to obtain a mass balance of fission products in the experiment, these cracks do not pose a major concern because the sink ring will not be analyzed in detail to obtain the spatial distribution of fission products. In Capsules 4 and 5, the compacts could not be removed from the inner rings. Strategies for removing the compacts are being evaluated. Sampling the inner rings with the compacts in-place is also an option. Dimensional measurements were made on the compacts, inner rings, outer rings, and sink rings. The diameters of all compacts decreased by 0.5 to 2.0 %. Generally, the extent of diametric shrinkage increased linearly with increasing neutron fluence. Most compact lengths also decreased. Compact lengths decreased with increasing fluence, reaching maximum shrinkage of about 0.9 % at a fast fluence of 4.0x1025 n/m2 E > 0.18 MeV. Above this fluence, the extent of length shrinkage appeared to decrease with fluence, and two compacts from Capsule 7 were found to have slightly increased in length (< 0.1 %) after a fluence of 5.2x1025 n/m2.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Neutronics Calculation Advances at Los Alamos: Manhattan Project to Monte Carlo

The history and advances of neutronics calculations at Los Alamos during the Manhattan Project through the present are reviewed. Substantial improvements to neutron diffusion methods and the invention of both the Monte Carlo neutron transport methods in 1947 and deterministic discrete ordinates Sn in 1953 were all made at Los Alamos just after the Manhattan Project. We briefly summarize early simpler and more approximate neutronics methods and then describe the need to better predict neutronics behavior through consideration of theoretical equations, models and algorithms, experimental measurements, and available computing capabilities and their limitations. This paper briefly covers key advances in deterministic methods during the Manhattan Project. These capabilities, coupled with increasing postwar defense needs and the invention of electronic computing with the Electronic Numeric Integrator and Computer, known as ENIAC, and the Mathematical Analyzer Numerical Integrator and Automatic Computer Model, known as MANIAC, led to the creation of Monte Carlo and deterministic discrete ordinates neutronics transport methods. We note the important role that the scientific comradery between the Los Alamos scientists played in the process. This paper briefly covers the early methods, algorithms, computers, and electronic and women pioneers that enabled Monte Carlo to spread to all areas of science. We focus heavily on these early developments and the subsequent creation of the MCNP® code, advances in its associated nuclear data, and its applications to problems of national defense at Los Alamos.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Collapse of neutrino wave functions under Penrose gravitational reduction

Models of spontaneous wave function collapse have been postulated to address the measurement problem in quantum mechanics. Their primary function is to convert coherent quantum superpositions into incoherent ones, with the result that macroscopic objects cannot be placed into widely separated superpositions for observably prolonged times. Many of these processes will also lead to loss of coherence in neutrino oscillations, producing observable signatures in the flavor profile of neutrinos at long travel distances. The majority of studies of neutrino oscillation coherence to date have focused on variants of the continuous state localization model, whereby an effective decoherence strength parameter is used to model the rate of coherence loss with an assumed energy dependence. Another class of collapse models that have been proposed posit connections to the configuration of gravitational field accompanying the mass distribution associated with each wave function that is in the superposition. A particularly interesting and prescriptive model is Penrose’s description of gravitational collapse which proposes a decoherence time τ determined through E g τ ∼ ℏ , where E g is a calculable function of the Newtonian gravitational potential. Here we explore application of the Penrose collapse model to neutrino oscillations, reinterpreting previous experimental limits on neutrino decoherence in terms of this model. We identify effects associated with both spatial collapse and momentum diffusion, finding that the latter is ruled out in data from the IceCube South Pole Neutrino Observatory so long as the neutrino wave packet width at production is σ ν , x ≤ 2 × 10 − 12 m . Published by the American Physical Society 2024

Astronomy & Astrophysics↗