Search NASA⌕ Search

SEARCH · Search NASA

Results for “solute transport”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

"Source Term Modeling for Advanced Gas Micro-Reactors"

Maintaining the safety of the public, environment, and operating personnel is the most important factor in designing, operating, maintaining, and decommissioning nuclear reactors. In recent years, there has been a growing interest in the development of micro-reactors employing TRi-structural ISOtropic (TRISO)-coated particle fuel. In gas reactors, TRISO fuel plays an important role in the safety case for high temperature reactors because of the fission product retention properties of the fuel. This ability enables the use of a functional containment strategy for the reactor where multiple barriers are used to prevent fission product release to the environment. Part of the safety analysis of these advanced reactors is the assessment of radionuclide releases under normal and accident conditions through the multiple credited safety barriers. Using conservative assumptions, a mechanistic analysis can be performed to quantify these releases that combines the probabilistic assessment of failure with analytic solutions to radionuclide transport equations. Source term modeling for TRISO fuel has been performed for previous reactor designs; however, these models are outdated, in many cases proprietary, and need updates to be applied to the current state of TRISO fuel technology and alternative gas reactor core configurations [1]. Currently, the only publicly available source term assessment for gas reactors is an expert-based Monte Carlo simulation based on the effectiveness of the fuel kernel, coating layers, and graphite block in a modular high temperature gas reactor [2]. Thus, there is a need to develop a simple, versatile, and mechanistic model of fission product release and transport in gas reactor cores that could be applied to a variety of reactors through user inputs and reactor-specific radionuclide inventories. The release is calculated by the diffusion of the key safety important fission products through the kernel, silicon carbide (SiC), graphite for both intact and defective TRISO particles based on fuel and graphite temperatures in the reactor under normal operation. These releases from the fuel enter the coolant where they can plate-out on cooler surfaces. A clean-up model is included for designs with a coolant purification system to remove fission gases. This initial distribution of fission products in the reactor serves as an initial condition for potential releases under postulated accident conditions. The model then can calculate the fission product release for any transient temperature profile and fission product releases can then be used to assess radiological dose to the workers and the public using conventional dose tools. Data on the diffusion of fission products is based on historic German TRISO experiments and the more current Department of Energy (DOE) Advanced Gas Reactor (AGR) TRISO fuel development program. The model is coded in python with inputs and outputs in excel spreadsheets, as well as python plotting utilities to aid in the interpretation of the results. References: [1] INL, NGNP Mechanistic Source Term White Paper, INL-10-17997, July 2010. [2] David A. Petti, Richard R. Hobbins, Peter Lowry, Hans Gougar, “Representative Source Terms and The Influence of Reactor Attributes on Functional Containment in Modular High Temperature Gas-cooled Reactors,” Nuclear Technology, Vol. 184, p. 181-197, Nov. 2013.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2023

Water transport in reverse osmosis membranes is governed by pore flow, not a solution-diffusion mechanism

We performed nonequilibrium molecular dynamics (NEMD) simulations and solvent permeation experiments to unravel the mechanism of water transport in reverse osmosis (RO) membranes. The NEMD simulations reveal that water transport is driven by a pressure gradient within the membranes, not by a water concentration gradient, in marked contrast to the classic solution-diffusion model. We further show that water molecules travel as clusters through a network of pores that are transiently connected. Permeation experiments with water and organic solvents using polyamide and cellulose triacetate RO membranes showed that solvent permeance depends on the membrane pore size, kinetic diameter of solvent molecules, and solvent viscosity. This observation is not consistent with the solution-diffusion model, where permeance depends on the solvent solubility. Motivated by these observations, we demonstrate that the solution-friction model, in which transport is driven by a pressure gradient, can describe water and solvent transport in RO membranes.

42 ENGINEERING↗
DOE OSTI2020

Mobile Sorption-based Thermal Battery for Harvesting Low-Temperature Geothermal Energy

Around 20% of the total primary energy in the United States is consumed for thermal demands of buildings such as space cooling, dehumidification, and space heating (EIA 2018). Low-temperature geothermal energy is abundant and can effectively satisfy buildings’ thermal demands. However, low-temperature geothermal energy is underutilized because the energy density of geothermal fluid is too low to justify the costs associated with transporting it between existing geothermal resources and buildings. The mobile sorption-based thermal battery (MSTB) system has been developed using three-phase (i.e., vapor–liquid, solution–solid, crystal) sorption technology to harvest low-temperature heat and store it with a much higher energy density than the geothermal fluid. The energy density of salt crystals is over six times higher than geothermal fluid, which makes long-distance transportation of salt crystals economically feasible. Salt crystals can be used to dehumidify air or provide space cooling in buildings, which alleviates peak demand on the electricity grid by offsetting electricity use for these end uses. This helps improve the grid’s stability and resilience. High-energy storage density, fast crystallization, and dissolution of salt crystals are all critical to the viability and performance of the MSTB system. Therefore, the design and operation of MSTB systems need to ensure effective generation and dissolution of salt crystals inside the MSTB. To achieve this target, this seedling project developed an experimental apparatus for characterizing the crystallization and dissolution processes. The energy density and potential latent cooling capacity of the MSTB are also evaluated based on lab test results. The crystallization results showed that the generated lithium chloride hydrate crystals are fluffy, the crystallization process lasts about 50 min, and the maximum crystal fraction (i.e., the ratio of crystal mass to the mass in the MSTB) can be up to 51.1% of the total mass in the MSTB at a solution flow rate of 1.58 g/s. The dissolution results show that the salt crystals in the MSTB can be fully dissolved within 15–28 min, based on different test conditions. Reducing solution flow rate and cooling water temperature can achieve increased energy storage density and crystal fraction. While the increase in the discharge rate (i.e., latent cooling capacity for dehumidifying air) is achieved by increasing flow rate and temperature of inlet diluted solution, as well as by using a pump for internal solution circulation, the discharge rate increases by 38%, from 0.95 kW to 1.31 kW. Compared with increasing the inlet solution flow rate, power consumption of salt solution transportation can be reduced by using a pump for internal solution circulation. The crystallization test results also showed that the maximum energy storage density is 981.8 kJ/kg, and the maximum discharge rate of the dissolution tests is ≤1.79 kW. Both are above the target values of 900 kJ/kg and 1.75 kW) for this project. The work reported here proves the feasibility and advancement of the MSTB system, which is helpful to the further study and improvement of the MSTB system.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Single Grid Error Estimation for Neutron Transport Solvers

The method of nearby problems (MNP) is a solution verification technique that does not require the use of multiple spatial grids. To estimate spatial discretization error without requiring a high-fidelity spatial grid, an analytical curve fit is interpolated from the numerical solution. The residual between the curve fit solution and numerical solution is calculated and added as an additional source term to the governing equation. The nearby solution is estimated using the updated source term and boundary conditions to remain consistent with the curve fit interpolation. The nearby solution can be compared to the curve fit solution as a discretization error estimation while using a single spatial grid. Without the use of higher fidelity spatial grids, the MNP is able to approximate the spatial discretization error, a facet of solution verification. The application of the method of nearby problems is presented for one- and two-dimensional neutron transport problems for both fixed source and criticality problems on the spatial variable. The fixed source results demonstrate the effectiveness of nearby problems for spatial error identification using the discrete ordinates method. Criticality results are shown to identify area of high spatial error for the C5G7 problem as well as for the discrete ordinates solver. A novel approach of combining the capabilities of Monte Carlo with the discrete ordinates nearby problems is presented for one- and two-dimensional fixed source problems. In conclusion, the MNP demonstrates its effectiveness at identifying spatial error on a single structured grid with a wide variety of neutron transport problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Transport error estimation using residual Monte Carlo

The residual Monte Carlo (RMC) method is also known in the literature as sequential Monte Carlo and reduced-source Monte Carlo. Given a Monte Carlo method for solving a linear equation and an approximate solution to that system, the residual method enables use of essentially the same Monte Carlo algorithm to directly compute the additive error or “defect” associated with the approximate solution. As the size of the defect decreases relative to the size of the solution, the residual Monte Carlo method becomes increasingly efficient relative to the standard Monte Carlo (SMC) method. Here we present a new RMC algorithm for evaluating the space-angle error in S n radiation transport solutions, and provide computational examples demonstrating that it can be far more efficient than SMC for this purpose. Herein we also describe a particular pitfall that must be avoided if RMC is to be efficient, and explain why the performance of RMC can significantly differ between different transport problems and different quantities of interest for the same problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

OECD/NEA COUPLED NEUTRONIC/THERMAL-FLUIDS BENCHMARK OF THE MHTGR-350 MW CORE DESIGN: RESULTS FOR PHASE-I EXERCISE 1

This report presents the results from code-to-code comparisons for a prismatic modular high-temperature gas-cooled reactor 350-MW design for steady-state neutron transport calculations. Integral parameters as well as power and flux distributions are provided. Differences between solutions from diffusion and transport operators are analysed. Effects of control rod homogenization between full block and one sixth block are discussed.

42 ENGINEERING↗
DOE OSTI2023

Probability of Initiation in Neutron Transport

We discuss the numerical solution of the nonlinear integro-differential equation for the probability of a divergent neutron chain in a stationary system (i.e., the probability of initiation (POI)). We follow the development described in Bell’s classic paper on the stochastic theory of neutron transport. As noted by Bell, the linearized form of this equation resembles the linear adjoint neutron transport equation. A matrix formalism for the discretized steady state (or forward) neutron equation in slab geometry is first developed and is then used to derive the discrete adjoint equation. A main advantage of this discrete development is that the resulting discrete adjoint equation does not depend upon how the multigroup cross sections for the forward problem are obtained. That is, we derive the discrete adjoint directly from the discrete forward equations rather than discretizing directly the adjoint equation. This also guarantees that the discrete adjoint operator is consistent with the inner product used to define the adjoint operator. We discuss three approaches for the numerical solution of the POI equations, and present numerical results on several test problems. The three solution methods are a simple fixed-point iteration, a second approach that is akin to a nonlinear Power iteration, and a third approach which uses a Newton-Krylov nonlinear solver. We also give sufficient conditions to guarantee the existence and uniqueness of nontrivial solutions to our discrete POI equations when the discrete system is supercritical, and that only the trivial solution exists when the discrete system is subcritical. Our approach is modeled after the analysis presented for the continuous POI equations by Mokhtar-Kharroubi and Jarmouni-Idrissi, and by Pazy and Rabinowitz.

42 ENGINEERING↗
DOE OSTI2024

A Stochastic Calculus Approach to Boltzmann Transport

Traditional Monte Carlo methods for particle transport utilize source iteration to express the solution, the flux density, of the transport equation as a Neumann series. Our contribution is to show that the particle paths simulated within source iteration are associated with the adjoint flux density and the adjoint particle paths are associated with the flux density. Here, we make our assertion rigorous through the use of stochastic calculus by representing the particle path used in source iteration as a solution to a stochastic differential equation (SDE). The solution to the adjoint Boltzmann equation is then expressed in terms of the same SDE, and the solution to the Boltzmann equation is expressed in terms of the SDE associated with the adjoint particle process. An important consequence is that the particle paths used within source iteration simultaneously provide Monte Carlo samples of the flux density and adjoint flux density in the detector and source regions, respectively. The significant practical implication is that particle trajectories can be reused to obtain both forward and adjoint quantities of interest. To the best our knowledge, the reuse of entire particles paths has not appeared in the literature. Monte Carlo simulations are presented to support the reuse of the particle paths.

Boltzmann transport↗
DOE OSTI2023

Quantifying selective solvent transport under an electric field in mixed-solvent electrolytes

Electrolytes in lithium-ion batteries comprise solvent mixtures, but analysis of ion transport is always based on treating the solvents as a single-entity. We combine electrophoretic NMR (eNMR) measurements and molecular dynamics (MD) simulations to quantify electric-field-induced transport in a concentrated solution containing LiPF 6 salt dissolved in an ethylene carbonate/ethyl methyl carbonate (EC/EMC) mixture. The selective transport of EC relative to EMC is reflected in the difference between two transference numbers, defined as the fraction of current carried by cations relative to the velocity of each solvent species. This difference arises from the preferential solvation of cations by EC and its dynamic consequences. The simulations reveal the presence of a large variety of transient solvent-containing clusters which migrate at different velocities. Rigorous averaging over different solvation environments is essential for comparing simulated and measured transference numbers. Our study emphasizes the necessity of acknowledging the presence of four species in mixed-solvent electrolytes.

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

Incorporating a naphthalene diimide polymer into a fullerene electron-transport layer to improve the fracture energy of perovskite solar cells

By blending a naphthalene diimide polymer into C 60 , we made a solution-processed electron-transport layer (ETL) for perovskite solar cells with fracture energies of 1.25 J m −2 , over 3× higher than that of thermally evaporated C 60 . Fracture energies were measured in a double cantilever beam configuration, and fracture surface images showed a fracture location near the ETL/perovskite interface, indicating a toughening of the interface between the ETL and Ag. We show that this modification to the ETL has no adverse effect on solar cell performance, and highlight the additional benefit of reduced parasitic absorption; a finding relevant for tandem solar cells.

electron transport layer↗
DOE OSTI2022

New Approaches to Improved PEM Fuel Cell Catalyst Layers

Polymer-electrolyte membrane (PEM) fuel-cells are one of the most promising energy conversion technologies for renewable clean energy applications. A major challenge preventing the widespread use and commercialization of PEM fuel cells is achieving high performance with low-loadings of platinum group metal (PGM) catalysts. One of the factors driving performance limitations in the cell is the mass transport losses within the cathode catalyst layers (CCL) due to sluggish oxygen-reduction reactions occurring at the platinum-ionomer interface. Any resistance to transport of these ionic and gaseous species within the CCL results in mass-transport limitations and performance losses, especially at high current densities. It is known that mass-transport losses increase with reduced platinum loading, thereby creating a performance-cost tradeoff for fuel cells. A viable solution to reduce the transport resistances in the catalyst layers is to create new ionomers that can provide good ion and oxygen transport needed to accomplish high-performing fuel cell catalysts. Using this approach Tetramer Technologies and LBNL propose improved fuel-cell catalyst ionomers based on Tetramers proprietary polymer chemistry.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

Analysis of Anisotropic Nonlocal Diffusion Models: Well-posedness of Fractional Problems for Anomalous Transport

We analyze the well-posedness of an anisotropic, nonlocal diffusion equation. Establishing an equivalence between weighted and unweighted anisotropic nonlocal diffusion operators in the vein of unified nonlocal vector calculus, we apply our analysis to a class of fractional-order operators and present rigorous estimates for the solution of the corresponding anisotropic anomalous diffusion equation. Furthermore, we extend our analysis to the anisotropic diffusion-advection equation and prove well-posedness for fractional orders s ∊ [0.5, 1). We also present an application of the advection-diffusion equation to anomalous transport of solutes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Innovative rail transport of a supersized land-based wind turbine blade

Wind turbine blade logistic providers are being challenged with escalating costs and routing complexities as one-piece blade approach lengths of 75 m in various regions of the U.S. land-based market. New lower cost solutions are needed to enable further reductions in the levelized cost of energy (LCOE) and continued market expansion. In this paper, a novel method of using existing U.S. rail infrastructure to deploy 100-m, one-piece blades to U.S. land-based wind sites is numerically investigated. The study removes the constraint that blades must be kept rigid during transport, and it allows bending to keep blades within a clearance profile while navigating horizontal and vertical curvatures. Novel system optimization and blade design processes consider blade structural constraints and rail logistic constraints in parallel to develop a highly flexible, rail-transportable blade. Results indicate maximum deployment potential in the Interior region of the United States and limited deployment potential in other regions. The study concludes that innovative rail transportation solutions combined with advanced rotor technologies can provide a feasible alternative to segmentation and support continued LCOE reductions in the U.S. land-based wind energy market.

17 WIND ENERGY↗
DOE OSTI2022

Plutonium reactive transport in fractured granite: Multi-species experiments and simulations

Plutonium (Pu) in the subsurface environment can transport in different oxidation states as an aqueous solute or as colloidal particles. The transport behavior of Pu is affected by the relative abundances of these species and can be difficult to predict when they simultaneously exist. This study investigates the concurrent transport of Pu intrinsic colloids, Pu(IV) (aq) and Pu(V-VI) (aq) through a combination of controlled experiments and semi-analytical dual-porosity transport modeling. Pu transport experiments were conducted in a fractured granite at high and low flow rates to elucidate sorption processes and their scaling behavior. In the experiments, Pu(IV) (aq) was the least mobile of the Pu species, Pu(V-VI) (aq) had intermediate mobility, and the colloidal Pu, which consisted mainly of precipitated and/or hydrolyzed Pu(IV), was the most mobile. The semi-analytical modeling revealed that the sorption of each Pu species was rate-limited, as the sorption could not be described by assuming local equilibrium in the experiments. The model was able to describe the sorption of the different Pu species that occurring either on fracture surfaces, in the pores of the rock matrix, or simultaneously in both locations. While equally good fits to the data could be achieved using any of these assumptions, a fracture-dominated process was considered to be the most plausible because it provided the most reasonable estimates of sorption rate constants. Importantly, a key result of this work is that the sorption rate constant of all Pu species tends to decrease with increasing time scales, which implies that Pu will tend to be more mobile at longer time scales than observations at shorter time scales suggest. This result has important implications for predicting the environmental impacts of Pu in the safety assessments of geologic repositories for radioactive waste disposal, and we explore potential mechanistic bases for upscaling the sorption rate constants to time and distance scales that cannot be practically evaluated in experiments.

58 GEOSCIENCES↗
DOE OSTI2021

Resonant amplification of enzymatic chemical oscillations by oscillating flow

Using theory and simulation, we analyzed the resonant amplification of chemical oscillations that occur due to externally imposed oscillatory fluid flows. The chemical reactions are promoted by two enzyme-coated patches located sequentially on the inner surface of a pipe that transports the enclosed chemical solution. In the case of diffusion-limited systems, the period of oscillations in chemical reaction networks is determined by the rate of the chemical transport, which is diffusive in nature and, therefore, can be effectively accelerated by the imposed fluid flows. In this work, we first identify the natural frequencies of the chemical oscillations in the unperturbed reaction–diffusion system and, then, use the frequencies as a forcing input to drive the system to resonance. We demonstrate that flow-induced resonance can be used to amplify the amplitude of the chemical oscillations and to synchronize their frequency to the external forcing. In particular, we show that even 10% perturbations in the flow velocities can double the amplitude of the resulting chemical oscillations. Particularly, effective control can be achieved for the two-step chemical reactions where during the first half-period, the fluid flow accelerates the chemical flux toward the second catalytic patch, while during the second half-period, the flow amplifies the flux to the first patch. The results can provide design rules for regulating the dynamics of coupled reaction–diffusion processes and can facilitate the development of chemical reaction networks that act as chemical clocks.

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

Verification of the LOTUS code with C5G7 benchmark

In this study, the verification of the neutron transport code LOTUS through the well-known C5G7 benchmark is performed. LOTUS uses the current coupled collision probability method with the expansion of the flux by orthogonal polynomials (CCCPO) for the solution of the neutron transport equation. The expansion of the flux by orthogonal polynomials allows one to avoid discretization of the calculation regions and significantly decreases the simulation time. The results of the LOTUS calculations are compared with the results of the reference OpenMC Monte Carlo simulations. The results of the comparison demonstrate almost perfect agreement with Monte-Carlo for the second order of the flux expansion. Even though the results obtained in the presented study agree well with the results of reference Monte Carlo calculations, further investigations are necessary for a better understanding of the stability and limitations of the code.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

A radioisotope - enabled reactive transport model for deep vadose zone carbon

In mountainous regions, which constitute the principle source of recharge to major rivers and regional aquifers, infiltration occurs through fractured, partially saturated, weathered bedrock that acts as a boundary layer between saturated aquifers and surface soil. Commonly this deep vadose zone (DVZ) is many meters thick, and yet its role in regulating the generation, retention and mobility of reactive solutes, including nutrients, contaminants and weathering products, is largely unknown. In particular, many of the key reactions that drive the formation of the weathered DVZ and the quality of water moving through it are redox processes, regulated by the availability of organic carbon and oxygen below the soil layer. The role of the DVZ is thus also poorly constrained in the context of carbon stocks and mobility, particularly in lithologies that are naturally high in organic carbon, such as shales. The overarching hypothesis of this study is that upland regions developed in geologic settings with abundant petrogenic carbon store and actively cycle carbon in the weathered DVZ below the soil and above the water table at rates that are significant and currently unconstrained. In order to quantify this cycling, the current study combines novel instrumentation techniques allowing new direct sampling of DVZ systems with advanced numerical reactive transport simulations of carbon transport and transformation. Critically, these simulations will explicitly treat the three isotopes of carbon (the abundant 12C, the stable rare 13C and the radioactive 14C) in a unified framework, thus clearly parsing between the contributions of modern surface derived carbon and lithologic carbon sources in integrated measurements of fluid and gas phase fluxes. This novel model capability will be applied to test the role of DVZ carbon cycling as a regulator of water quality and geological weathering in two complementary field sites both located in organic carbon rich shale lithologies. The first is the Eel River Critical Zone Observatory (ERCZO) in Mendocino County, California, and the second is the Lawrence Berkeley National Laboratory Watershed Function Scientific Focus Area (SFA) in the East River watershed, near Crested Butte, Colorado. At the ERCZO, a novel vadose zone monitoring system has been installed in a 20 m thick, partially saturated, weathered shale hillslope, and preliminary data already indicate substantial CO2 flux generated many meters below the soil surface. At the SFA field site, an instrumented hillslope transect indicates a more complex multi-dimensional fluid and solute transport regime, which will serve as a key test of the calibrated models. Collectively, this project will advance understanding of the cycling of carbon belowground and in relation to transport pathways across the poorly constrained DVZ characteristic of primary water recharge areas. The key product of this work will be enhanced isotope simulation capabilities that are robust and publicly available for application across a broad diversity of systems.

58 GEOSCIENCES↗