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

Evaluation of the Effect of Burnup Acceleration on UO 2 Microstructure Evolution

Accelerated fuel qualification (AFQ) is a methodology by which new nuclear fuels are developed in an accelerated time frame compared with historical fuel qualification approaches. AFQ generally relies on high-fidelity physics-based modeling and simulation tools to adequately describe fuel performance as well as on revolutionary methods to accelerate burnup accumulation and collect relevant data more quickly. This report summarizes the use of advanced fuel modeling and simulation tools to evaluate microstructures from commercially irradiated fuel and microstructures from proposed MiniFuel irradiations, in which burnup accumulation is accelerated while prototypic temperature conditions are maintained. In this milestone, we used the mesoscale fuel performance code MARMOT to model the evolution of irradiated UO 2 microstructures and their potential restructuring at high burnup. The simulation conditions were informed by BISON models of both commercially irradiated fuel and MiniFuel. A first set of simulations investigated the recrystallization behavior of fully dense microstructures and showed full recrystallization at burnups as low as 52 MWd/kgU at 950°C. However, these simulations did not account for the presence of fission gas bubbles (predicted by BISON). Therefore, a second set of simulations including fission gas bubbles was performed and indicated that at the lowest temperature considered (650°C), the porous UO 2 microstructures have the highest total Gibbs free energies and are likely to recrystallize earlier than higher temperature cases (800 and 950°C), which agrees with high-burnup fuel characterization data. The results also showed that at lower temperature (650°C), the total free energies of the PWR fuel and MiniFuel microstructures are not significantly different. However, at the highest temperature (950°C), MiniFuel microstructures have a lower free energy than that of the PWR fuel microstructure. The competing effects between the temperature-dependent grain nucleation rate and the reduction of the free energy of the microstructure at higher temperature as a result of diffusion indicated that restructuring may occur at even higher temperatures than those considered in this study. In addition to the microstructure evolution modeling efforts, the burnup gradient across a single fuel specimen was also considered. This evaluation was for the VXF-15 position of the High Flux Isotope Reactor (HFIR) using the code suite HFIRCON, which was developed to automate the workflow for evaluating targets and fuel as they are irradiated in HFIR. The burnup gradient evaluation showed a dependence on both the axial and radial locations within the specimen, with a maximum difference of 1.7 between the inner and outermost radial layers. This relationship was further supported by considering the fission product speciation with respect to location within the specimen, which showed a higher concentration of 239 Pu, 240 Pu, and 241 Pu on the outer radial locations of the specimen than the center. The findings of the burnup and speciation evaluation show that some amount of self-shielding is occurring in the specimen when irradiated in the high-flux environment of HFIR; however, this impact is more pronounced for natural uranium when compared to 6% enrichment due to the higher ratio of 238 U in the specimen. Further analyses are required to understand the sensitivity of this gradient to spatial mesh and enrichment of the specimen.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Numerical predictions of cryogenic hydrogen vertical jets

Here, a comparison of Computational Fluid Dynamics (CFD) predictions with measurements is presented for cryo-compressed hydrogen vertical jets. The stagnation conditions of the experiments are characteristic of unintended leaks from pipe systems that connect cryogenic hydrogen storage tanks and could be encountered at a fuel cell refueling station. Jets with pressure up to 5 bar and temperatures just above the saturation liquid temperature were examined. Comparisons are made to the centerline mass fraction and temperature decay rates, the radial profiles of mass fraction and the contours of volume fraction. Two notional nozzle approaches are tested to model the under-expanded jet that was formed in the tests with pressures above 2 bar. In both approaches the mass and momentum balance from the throat to the notional nozzle are solved, while the temperature at the notional nozzle was assumed equal to the nozzle temperature in the first approach and was calculated by an energy balance in the second approach. The two approaches gave identical results. Satisfactory agreement with the measurements was found in terms of centerline mass fraction and temperature. However, for test with 3 and 4 bar release the concentration was overpredicted. Furthermore, a wider radial spread was observed in the predictions possibly revealing higher degree of diffusion using the k-ε turbulence model. An integral model for cryogenic jets was also developed and provided good results. Finally, a test simulation was performed with an ambient temperature jet and compared to the cold jet showing that warm jets decay faster than cold jets.

08 HYDROGEN↗
DOE OSTI2016

Multi-scale diffuse interface modeling of multi-component two-phase flow with partial miscibility

In this paper, we introduce a diffuse interface model to simulate multi-component two-phase flow with partial miscibility based on a realistic equation of state (e.g. Peng–Robinson equation of state). Because of partial miscibility, thermodynamic relations are used to model not only interfacial properties but also bulk properties, including density, composition, pressure, and realistic viscosity. As far as we know, this effort is the first time to use diffuse interface modeling based on equation of state for modeling of multi-component two-phase flow with partial miscibility. In numerical simulation, the key issue is to resolve the high contrast of scales from the microscopic interface composition to macroscale bulk fluid motion since the interface has a nanoscale thickness only. To efficiently solve this challenging problem, we develop a multi-scale simulation method. At the microscopic scale, we deduce a reduced interfacial equation under reasonable assumptions, and then we propose a formulation of capillary pressure, which is consistent with macroscale flow equations. Moreover, we show that Young–Laplace equation is an approximation of this capillarity formulation, and this formulation is also consistent with the concept of Tolman length, which is a correction of Young–Laplace equation. At the macroscopical scale, the interfaces are treated as discontinuous surfaces separating two phases of fluids. Our approach differs from conventional sharp-interface two-phase flow model in that we use the capillary pressure directly instead of a combination of surface tension and Young–Laplace equation because capillarity can be calculated from our proposed capillarity formulation. A compatible condition is also derived for the pressure in flow equations. Furthermore, based on the proposed capillarity formulation, we design an efficient numerical method for directly computing the capillary pressure between two fluids composed of multiple components. Finally, numerical tests are carried out to verify the effectiveness of the proposed multi-scale method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Plasma edge and plasma-wall interaction modelling: Lessons learned from metallic devices

Robust power exhaust schemes employing impurity seeding are needed for target operational scenarios in present day tokamak devices with metallic plasma-facing components (PFCs). For an electricity-producing fusion power plant at power density P sep /R > 15 MW/m divertor detachment is a requirement for heat load mitigation. 2D plasma edge transport codes like the SOLPS code as well as plasma-wall interaction (PWI) codes are key to disentangle relevant physical processes in power and particle exhaust. With increased quantitative credibility in such codes more realistic and physically sound estimates of the life-time expectations and performance of metallic PFCs can be accomplished for divertor conditions relevant for ITER and DEMO. An overview is given on the recent progress of plasma edge and PWI modelling activities for (carbon-free) metallic devices, that include results from JET with the ITER-like wall, ASDEX Upgrade and Alcator C-mod. It is observed that metallic devices offer an opportunity to progress the understanding of underlying plasma physics processes in the edge. The validation of models can be substantially improved by eliminating carbon from the experiment as well as from the numerical system with reduced degrees of freedom as no chemical sputtering from amorphous carbon layers and no carbon or hydro-carbon transport are present. With the absence of carbon as the primary plasma impurity and given the fact that the physics of the PWI at metallic walls is less complex it is possible to isolate the crucial plasma physics processes relevant for particle and power exhaust. Finally, for a reliable 2D dissipative plasma exhaust model these are: cross-field drifts, complete kinetic neutral physics, geometry effects (including main-chamber, divertor and sub-divertor structures), SOL transport reflecting also the non-diffusive nature of anomalous transport, as well as transport within the pedestal region in case of significant edge impurity radiation affecting pedestal pressure and hence P sep .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Modelling of noble anaesthetic gases and high hydrostatic pressure effects in lipid bilayers

Our objective was to study molecular processes that might be responsible for inert gas narcosis and high-pressure nervous syndrome. The classical molecular dynamics trajectories (200 ns-long) of dioleoylphosphatidylcholine (DOPC) bilayers simulated by the Berger force field were evaluated for water and the atomic distribution of noble gases around DOPC molecules at a pressure range of 1 - 1000 bar and temperature of 310 Kelvin. Xenon and argon have been tested as model gases for general anesthetics, and neon has been investigated for distortions that are potentially responsible for neurological tremor at hyperbaric conditions. The analysis of stacked radial pair distribution functions of DOPC headgroup atoms revealed the explicit solvation potential of gas molecules, which correlates with their dimensions. The orientational dynamics of water molecules at the biomolecular interface should be considered as an influential factor; while excessive solvation effects appearing in the lumen of membrane-embedded ion channels could be a possible cause of inert gas narcosis. All the noble gases tested exhibit similar patterns of the order parameter for both DOPC acyl chains, which is opposite to the patterns found for the order parameter curve at high hydrostatic pressures in intact bilayers. This finding supports the ‘critical volume’ hypothesis of anesthesia pressure reversal. The irregular lipid headgroup-water boundary observed in DOPC bilayers saturated with neon in the pressure range of 1 - 100 bar could be associated with the possible manifestation of neurological tremor at the atomic scale. The non-immobilizer neon also demonstrated the highest momentum impact on the normal component of the DOPC diffusion coefficient representing monolayers undulations rate, which indicates enhanced diffusivity, rather than atom size, as the key factor.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Kelvin Probe Force Microscopy in liquid using Electrochemical Force Microscopy

Conventional closed loop-Kelvin probe force microscopy (KPFM) has emerged as a powerful technique for probing electric and transport phenomena at the solid-gas interface. The extension of KPFM capabilities to probe electrostatic and electrochemical phenomena at the solid–liquid interface is of interest for a broad range of applications from energy storage to biological systems. However, the operation of KPFM implicitly relies on the presence of a linear lossless dielectric in the probe-sample gap, a condition which is violated for ionically-active liquids (e.g., when diffuse charge dynamics are present). Here, electrostatic and electrochemical measurements are demonstrated in ionically-active (polar isopropanol, milli-Q water and aqueous NaCl) and ionically-inactive (non-polar decane) liquids by electrochemical force microscopy (EcFM), a multidimensional (i.e., bias- and time-resolved) spectroscopy method. In the absence of mobile charges (ambient and non-polar liquids), KPFM and EcFM are both feasible, yielding comparable contact potential difference (CPD) values. In ionically-active liquids, KPFM is not possible and EcFM can be used to measure the dynamic CPD and a rich spectrum of information pertaining to charge screening, ion diffusion, and electrochemical processes (e.g., Faradaic reactions). EcFM measurements conducted in isopropanol and milli-Q water over Au and highly ordered pyrolytic graphite electrodes demonstrate both sample- and solvent-dependent features. Finally, the feasibility of using EcFM as a local force-based mapping technique of material-dependent electrostatic and electrochemical response is investigated. The resultant high dimensional dataset is visualized using a purely statistical approach that does not require a priori physical models, allowing for qualitative mapping of electrostatic and electrochemical material properties at the solid–liquid interface.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

BISON Simulated and Experimental Fission Product Release Comparisons from Reradiated AGR-3/4 Compacts During High Temperature Heating Tests

The fuel performance modeling code BISON was used to predict the release of fission products iodine-131 (131I), xenon-133 (133Xe), and krypton-85 (85Kr) from four re-irradiated AGR-3/4 fuel compacts containing tristructural isotropic (TRISO) coated particles during high-temperature isothermal heating tests. The AGR-3/4 fuel compacts were irradiated in the Advanced Test Reactor (ATR) as part of the third and fourth series of planned experiments to support the Advanced Gas Reactor (AGR) Program. They were subsequently stored and re-irradiated in the Neutron Radiography (NRAD) reactor for approximately five days and then stored for another five to eight days before being subjected to isothermal heating tests in the Fuel Accident Condition Simulation (FACS furnace) for 200 to 300 hours at temperatures between 1000°C and 1600°C to evaluate fission product release at elevated temperatures. New nuclide-specific fission product source term models for the three nuclides of interest were developed using the reactor multiphysics code Griffin and implemented into BISON to support this work. The new source term models were incorporated into coupled compact- and particle-scale BISON simulations, which predict spatially- and temporally-resolved radionuclide generation, radioactive decay, transport, and release throughout the entire irradiation history, including the initial ATR irradiation, NRAD re-irradiations, FACS heating tests, and intermediate periods spent in storage. The experimentally measured fission product release from the heating tests were compared to modeling release predictions calculated by BISON to evaluate how well the code compares to experimental results. Overall, the experimental measured and BISON predicted comparative release results varied but generally agreed to within 5 particle equivalents. Comparative release results identified general observations to take into consideration to help refine future models and reduce uncertainties associated with both the measurement results and predictive results. This includes developing new uranium oxycarbide (UCO) specific kernel diffusivities for the three isotopes examined to more accurately reflect the material properties of the fuel form. Deriving new diffusivities will aid in producing a more informed BISON model

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Analytical models of hydrogen transport in graphite

The importance of graphite-hydrogen chemical reactions to fusion, fission, and hydrogen storage applications, combined with the rapidly evolving knowledge on the underlying mechanisms, has led to the development of multiple models to describe hydrogen transport in graphite. Significant differences exist among these models, resulting from discrepancies in the modeling assumptions, intended degree of fidelity, and conditions of applicability. This paper attempts at reconciling these apparent differences by providing a comprehensive description of the constitutive equations governing hydrogen transport in graphite at high-temperature, identifying outstanding gaps in knowledge, illustrating how these different models approach them, and proposing alternative analytical formulations grounded on experimental results from hydrogen-graphite studies. Governing equations, closing relations, and simplifying assumptions are discussed for hydrogen transport at the inter-granular and intra-granular level, accompanied by compiled experimental data and illustrated energy diagrams associated to the proposed transport mechanisms. Analytical formulations are provided to reproduce competing hypotheses on the mechanisms, supporting the development of a range of computational models that can enable resolution of outstanding knowledge gaps through comparative testing against experimental data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Computational Modeling of Fluctuations in Energy and Metabolic Pathways of Methanogenic Archaea

The methanogenic archaea, anaerobic microbes that convert CO2 and H2 and/or other small organic fermentation products into methane, play an unusually large role in the global carbon cycle. As they perform the final step in the anaerobic breakdown of biomass, methanogens are a biogenic source of an estimated one billion tons methane each year. Depending on the location, produced methane can be considered as either a greenhouse gas (agricultural byproduct), sequestered carbon storage (methane hydrate deposits), or a potential energy source (organic wastewater treatment). These microbes therefore represent an important target for biotechnology applications. Computational models of methanogens with predictive power are useful aids in the adaptation of methanogenic systems, but need to connect processes of wide-ranging time and length scales. In this project, we developed several computational methodologies for modeling the dynamic behavior of entire cells that connects stochastic reaction-diffusion dynamics of individual biochemical pathways with genome-scale modeling of metabolic networks. While each of these techniques were in the realm of well-defined computational methods, here we integrated them to develop several entirely new approaches to systems biology. The first scientific aim of the project was to model how noise in a biochemical pathway propagates into cellular phenotypes. Genetic circuits have been optimized by evolution to regulate molecular processes despite stochastic noise, but the effect of such noise on a cellular biochemical networks is currently unknown. An integrated stochastic/systems model of Escherichia coli species was created to analyze how noise in protein expression gives—and therefore noise in metabolic fluxes—gives rise to multiple cellular phenotype in isogenic population. After the initial work developing and validating methods that allow characterization of the heterogeneity in the model organism E. coli, the project shifted toward investigations of the methanogen Methanosarcina acetivorans. By integrating an unprecedented transcriptomics dataset for growth of the methanogen on many substrates with an in silico model, heterogeneity in metabolic pathway usage and methane production were examined. This lent insight into the physiological requirements of the organism under different environmental conditions and uncovered the unique regulatory role that mRNA half-life has in shaping metabolic flux distributions in this organism.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

EGR Dilution and Fuel Property Effects on High-Efficiency Spark-Ignition Flames

Modern spark ignition internal combustion engines rely on fast combustion rates and high dilution to achieve high brake thermal efficiencies. To accomplish this, new engine designs have moved towards increased tumble ratios and stroke-to-bore ratios. Increased tumble ratios correlate positively with increases in turbulent kinetic energy and improved fuel and residual gas mixing, all of which favor faster and more efficient combustion. Longer stroke-to-bore ratios allow higher geometric compression ratios and use of late intake valve closing to control peak compression pressures and temperatures. The addition of dilution to improve efficiency is limited by the resulting increase in combustion instabilities manifested by cycle-to-cycle variability. A number of effects - preferential diffusion, turbulence-combustion interactions, stochastic flow patterns, laminar-turbulent flame kernel transitions, and relative length and velocity scales between flame and turbulence - are believed to be responsible for the increase in cycle-to-cycle variations, where their contributions are likely interlinked. Several studies have shown the influence of stochastic flow characteristics on the nature of combustion instabilities, such as velocity patterns on flame kernel formation and cycle-to-cycle variations in residual gas. However, few have focused on the specific effects of fuel properties. The objective of this work is to contrast the effects of dilution on propane stoichiometric combustion against gasoline. Dilution tolerance experiments were conducted in a purpose-built high stroke-to-bore ratio single cylinder engine with both gasoline and LPG. Three-dimensional full cycle computational fluid dynamics (CFD) simulations employing a level-set combustion approach and Reynolds averaged Navier-Stokes (RANS) turbulence modeling was used to qualitatively assess the changes in length and velocity scales for turbulence and the flame. The experimental results showed that LPG can tolerate higher exhaust gas recirculation (EGR) dilution under a variety of conditions. Analysis of CFD simulations showed that propane flames are likely less sensitive to influences from the flow field due to less thickening of the flame and higher effective flame speeds.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Applications of Nickelate perovskites for neuromorphic computing from electronic structure and Machine Learning

While the limit of Moore's law is presently being reached with current microelectronic technologies, we need to develop new paradigms that overcome this limitation. In that respect, neuromorphic computing is a concept that emulates the neural behavior and response of the human brain, and it has been recognized as a promising alternative approach. In this research project, we will perform multi-fidelity scale bridging to explore the potential use of materials with metal to insulator transition for neuromorphic applications. In particular, rare earth nickelates are promising for such purposes, as the transition in these materials is quite sensitive to a broad set of different external stimuli. Our multi-fidelity approach will bridge the high-fidelity electronic structure calculations with classical potentials. We will bridge dynamical mean field theory with a classical atomistic representation via a deep learning force field. The neural network is trained with energies, charges, and forces obtained by accurate electronic structure theories based on Dynamical Mean Field Theory. The configurational space is generated from known crystal phases, ab initio molecular dynamics with exchange-correlation functionals corrected with the Hubbard model, disordered phases with different concentrations of oxygen vacancies, and nonsymmetrical positions and induced strain by grain interfaces or contact with a substrate. Strategies to train the model with a reduced number of training examples are obtained from active learning methods, and new structures for improving the learning process are generated by using machine learning autoencoders. This classical potential will be validated through a diversity of electronic structure methods and represents an important step to combine the flexibility and accuracy of first-principles with the speed of classical potentials. The generated multi-fidelity surrogate model will be used to understand the role of strain, oxygen vacancies, proton doping, the variation of the crystal phase, substrate effects, vibrational effects as the octahedral rotation, grain boundaries and defect effects on the response of a Metal to Insulator Transition (MIT) in correlated materials. Long time and large-scale simulations will help understand the role of different stimuli to control the hysteresis of the MIT, as it has been experimentally suggested. Selected configurations will be analyzed with higher-level theories to provide an accurate electronic description and to study how the orbitals and charges are rearranged under different conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Characterizing the Subsurface Using Deformation from Pumping and Surface Loading Tests

Changes in pore pressure or applied load cause deformation that can be used to characterize the subsurface. A series of field tests were conducted by measuring the deformation and pressure response from a variety of stresses at our test site in South Carolina. The site is underlain by saprolite to a depth of 30 m, and the water table in the unconfined aquifer occurs at a depth of 9 m. The types of tests performed included constant-rate pumping tests and surface point loading tests using the weight of a person or vehicle. Vertical displacement of the saprolite in response to these stresses was measured at depths of 3, 6, 9 and 11 meters in both saturated and unsaturated conditions using an extensometer with a resolution of +/- 10 nm. Results from 2.5-hr-long constant-rate pumping tests indicate a repeatable, compressive vertical strain of roughly 0.1 microstrain in the vadose zone and 1-10 microstrain in the saturated zone. Drawdown was 5 m in the pumping well and 0.1 m in a nearby (r = 9m) observation well. Under saturated conditions, the strain rate gradually decreased throughout the test in a pattern that resembles the drawdown curve in the monitoring well. In contrast, the strain in the vadose zone increased quickly early in the test, but the strain rate diminished and was negligible after approximately 30 minutes of pumping. Interestingly, drawdown in the pumping well also increased quickly over 30 minutes and remained roughly constant after that. Additional tests were conducted by applying a load (weight of a person or vehicle) at the ground surface above the extensometer and allowing the signal to stabilize before removing the load. In the vadose zone, the soil compressed rapidly and then stabilized. In the saturated zone, sudden compression was followed by additional compression at a much slower rate. A similar response was observed upon removal of the load. This transient signal in the saturated zone is interpreted to be a poroelastic effect. When a load is applied, the solid skeleton compresses and pressurizes the pore fluid, causing flow. The strain rate is then a function of pressure diffusion. Simulation of the deformation from both pumping and loading tests using an unsaturated, poroelastic finite-element model indicates that strains are sensitive to the Young’s modulus and hydraulic conductivity as functions of saturation.

well tests, strain, surface loading, hydromechanic↗
DOE OSTI2021

Reactive force fields for aqueous and interfacial magnesium carbonate formation

Here we develop Mg/C/O/H ReaxFF parameter sets for two environments: an aqueous force field for magnesium ions in solution and an interfacial force field for minerals and mineral–water interfaces. Since magnesium is highly ionic, we choose to fix the magnesium charge and model its interaction with C/O/H through Coulomb, Lennard-Jones, and Buckingham potentials. We parameterize the forcefields against several crystal structures, including brucite, magnesite, magnesia, magnesium hydride, and magnesium carbide, as well as Mg 2+ water binding energies for the aqueous forcefield. Then, we test the forcefield for other magnesium-containing crystals, solvent separated and contact ion-pairs and single-molecule/multilayer water adsorption energies on mineral surfaces. We also apply the forcefield to the forsterite–water and brucite–water interface that contains a bicarbonate ion. We observe that a long-range proton transfer mechanism deprotonates the bicarbonate ion to carbonate at the interface. Free energy calculations show that carbonate can attach to the magnesium surface with an energy barrier of about 0.22 eV, consistent with the free energy required for aqueous Mg–CO 3 ion pairing. Also, the diffusion constant of the hydroxide ions in the water layers formed on the forsterite surface are shown to be anisotropic and heterogeneous. These findings can help explain the experimentally observed fast nucleation and growth of magnesite at low temperature at the mineral–water–CO 2 interface in water-poor conditions.

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

A force-level theory of the rheology of entangled rod and chain polymer liquids. I. Tube deformation, microscopic yielding, and the nonlinear elastic limit

In this paper, we employ a first-principles-based, force-level approach to construct the anharmonic tube confinement field for entangled fluids of rigid needles, and also for chains described at the primitive-path (PP) level in two limiting situations where chain stretch is assumed to either be completely equilibrated or unrelaxed. The influence of shear and extensional deformation and polymer orientation is determined in a nonlinear elastic limit where dissipative relaxation processes are intentionally neglected. For needles and PP-level chains, a self-consistent analysis of transverse polymer harmonic dynamical fluctuations predicts that deformation-induced orientation leads to tube weakening or widening. In contrast, for deformed polymers in which chain stretch does not relax, we find tube strengthening or compression. For all three systems, a finite maximum transverse entanglement force localizing the polymers in effective tubes is predicted. The conditions when this entanglement force can be overcome by an externally applied force associated with macroscopic deformation can be crisply defined in the nonlinear elastic limit, and the possibility of a “microscopic absolute yielding” event destroying the tube confinement can be analyzed. For needles and contour-relaxed PP chains, this force imbalance occurs at a stress of order the equilibrium shear modulus and a strain of order unity, corresponding to a mechanically fragile entanglement tube field. However, for unrelaxed stretched chains, tube compression stabilizes transverse polymer confinement, and there appears to be no force imbalance. These results collectively suggest that the crossover from elastic to irreversible viscous response requires chain retraction to initiate disentanglement. Finally, we qualitatively discuss comparisons with existing phenomenological models for nonlinear startup shear, step strain, and creep rheology experiments.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2014

RADIAL TRANSPORT OF LARGE-SCALE MAGNETIC FIELDS IN ACCRETION DISKS. II. RELAXATION TO STEADY STATES

We study the time evolution of a large-scale magnetic flux threading an accretion disk. The induction equation of the mean poloidal field is solved under the standard viscous disk model. Magnetic flux evolution is controlled by two timescales: one is the timescale of the inward advection of the magnetic flux, τ{sub adv}. This is induced by the dragging of the flux by the accreting gas. The other is the outward diffusion timescale of the magnetic flux τ{sub dif}. We consider diffusion due to the Ohmic resistivity. These timescales can be significantly different from the disk viscous timescale τ{sub disk}. The behaviors of the magnetic flux evolution are quite different depending on the magnitude relationship of the timescales τ{sub adv}, τ{sub dif}, and τ{sub disk}. The most interesting phenomena occur when τ{sub adv} << τ{sub dif}, τ{sub disk}. In such a case, the magnetic flux distribution approaches a quasi-steady profile much faster than the viscous evolution of the gas disk, and the magnetic flux has also been tightly bundled to the inner part of the disk. In the inner part, although the poloidal magnetic field becomes much stronger than the interstellar magnetic field, the field strength is limited to the maximum value that is analytically given by our previous work. We also find a condition for the initial large magnetic flux, which is a fossil of the magnetic field dragging during the early phase of star formation that survives for a duration in which significant gas disk evolution proceeds.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The response of the Amazon ecosystem to the photosynthetically active radiation fields: integrating impacts of biomass burning aerosol and clouds in the NASA GEOS Earth system model

The Amazon experiences fires every year, and the resulting biomass burning aerosols, together with cloud particles, influence the penetration of sunlight through the atmosphere, increasing the ratio of diffuse to direct photosynthetically active radiation (PAR) reaching the vegetation canopy and thereby potentially increasing ecosystem productivity. In this study, we use the NASA Goddard Earth Observing System (GEOS) model with coupled aerosol, cloud, radiation, and ecosystem modules to investigate the impact of Amazon biomass burning aerosols on ecosystem productivity, as well as the role of the Amazon's clouds in tempering this impact. The study focuses on a 7-year period (2010–2016) during which the Amazon experienced a variety of dynamic environments (e.g., La Niña, normal years, and El Niño). The direct radiative impact of biomass burning aerosols on ecosystem productivity – called here the aerosol diffuse radiation fertilization effect – is found to increase Amazonian gross primary production (GPP) by 2.6 % via a 3.8 % increase in diffuse PAR (DFPAR) despite a 5.4 % decrease in direct PAR (DRPAR) on multiyear average during burning seasons. On a monthly basis, this increase in GPP can be as large as 9.9 % (occurring in August 2010). Consequently, the net primary production (NPP) in the Amazon is increased by 1.5 %, or ~92 Tg C yr -1 – equivalent to ~37 % of the average carbon lost due to Amazon fires over the 7 years considered. Clouds, however, strongly regulate the effectiveness of the aerosol diffuse radiation fertilization effect. The efficiency of this fertilization effect is the highest in cloud-free conditions and linearly decreases with increasing cloud amount until the cloud fraction reaches ~0.8, at which point the aerosol-influenced light changes from being a stimulator to an inhibitor of plant growth. Nevertheless, interannual changes in the overall strength of the aerosol diffuse radiation fertilization effect are primarily controlled by the large interannual changes in biomass burning aerosols rather than by changes in cloudiness during the studied period.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Probing the Structure and Dynamics of Fluid Mixtures in Porous Materials Through Ultrafast Vibrational Spectro-Microscopy and Many-Body Molecular Dynamics

This research program focused on the characterization of both molecular structure and dynamics of aqueous solutions in various metal-organic frameworks (MOFs) and organic polymers that have recently been proposed for applications in water treatment technologies. This was accomplished by integrating ultrafast vibrational spectroscopy and microscopy, which are sensitive to local environments, molecular orientation, and dynamical couplings, with many-body molecular dynamics (MB-MD) simulations, which enable realistic modeling of aqueous systems and porous materials under different thermodynamic conditions. Our studies allowed for identifying the physical mechanisms and characterizing the underlying molecular interactions between water and host materials that determine the adsorption and transport processes of water in various MOFs and β-cyclodextrin polymers – prototypical examples of porous materials and polymers, that have recently been proposed for technological applications in water purification. By combining bulk- and surface-sensitive, spatially resolved ultrafast vibrational spectroscopy with MB-MD simulations, it was possible to gain broad insight into the structure and mobility of water under heterogeneous confinement at various length scales. In particular, measurements of surface- and bulk-sensitive vibrational spectra at single crystal level, integrated with MB-MD simulations were used to characterize domain-specific, structure–binding affinity relationships. Besides providing specific information about the behavior of aqueous solutions in prototypical MOFs and organic polymers for potential applications in water treatment technologies, our studies also provided fundamental insights into the properties of aqueous solutions in confined heterogeneous environments, which have broad implications in many areas, ranging from heterogeneous catalysis to oil recovery, ion transport processes, and reaction-diffusion processes in crowded systems.

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

Radionuclide Waste Disposal: Development of Multi-scale Experimental and Modeling Capabilities (Final Report)

The DOE EPSCoR Implementation project “Radioactive waste management: Development of multi-scale experimental and modeling capabilities” helped to develop a team of scientists and engineers from Clemson University, South Carolina State University, and the University of South Carolina to address the disposition of nuclear wastes and study the transport of radioisotopes from a waste repository in the near and far field. The project involved 20 faculty from the three institutions as well as 13 postdoctoral fellows, 32 graduate students, and 32 undergraduate students and was active from 2014-2019. Additionally, we forged new collaborations with nine researchers from DOE laboratories SRNL, LLNL, and ANL as well as the University of Manchester and the China Academy of Engineering Physics. The overarching goal of the project was to understand the conditions under which important classes of co-reactants, ranging from counter ions in crystal lattices to dissolved oxygen in pores, control the chemistry and transport characteristics of radionuclides in engineered waste forms and natural soils. Our approach was to characterize the time and length scales over which non-equilibrium states are maintained by rate-limiting, or rate-enhancing, reactions between radionuclides and co-reactants due to interactions between physical mass-transfer processes (i.e., advection, diffusion) and (biogeo) chemical reactions. We have focused our project on three specific classes of reactions relevant to radionuclide transport at DOE legacy sites: ion exchange/substitution, ligand complexation, and redox-mediated reactions. Understanding radionuclide migration requires detailed knowledge of how changes to a system – whether engineered or natural – drive the behavior of co-reactants, which in turn provide the geochemical context controlling radionuclide transport. Student engagement and training were a primary focus of the project in order to create a pipeline of researchers who could work in the area of nuclear waste disposition to support the state and the nation. Over the duration of the project we worked with 32 undergraduate, graduated 18 M.S. students and 14 Ph.D. students, and advised 13 postdoctoral fellows. The Ph.D students and postdocs have primarily taken positions at DOE laboratories, academia, and industry. Through our collaborative team, numerous follow on projects have been started with over $5M in sponsored research. Additionally, thus far the team has published 53 peer reviewed papers and given over 75 technical presentations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗