Predicting PEMFC performance from a volumetric image of catalyst layer structure using pore network modeling
Not Available
SEARCH · Search NASA
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.
Not Available
A new concept for both transport and storage of heat at high temperatures and heat fluxes is introduced and the first steps in analysis of its feasibility is taken. The concept utilizes the high energy storage capability of materials undergoing change of phase. The phase change material, for example a salt, is encapsulated in corrosion resistant sealed pellets and transported in a carrier fluid to heat source and storage. Calculations for heat transport from a typical solar collector indicate that the pellet mass flow rates are relatively small and that the required pumping power is only a small fraction of the energy transport capability of the system. Salts and eutectic salt mixtures as candidate phase change materials are examined and discussed. Finally, the time periods for melting or solidification of sodium chloride pellets is investigated and reported.
There is considerable interest in developing high-performance electrolytes for rechargeable lithium batteries. For practical applications, the electrolyte must support large dc currents. However, the parameters most often reported in the literature, conductivity, κ, and current fraction, ρ + , reflect ion transport in the limit of infinitesimal currents. In this limit, the efficacy of an electrolyte is given by the product κρ + . The limiting current density, i lim , is the maximum current density that can be applied across an electrolyte; the cell voltage diverges if the applied current density exceeds ilim. This parameter reflects ion transport in the limit of large dc currents and is therefore of practical interest. It would therefore be convenient if i lim could be predicted from measurements of κρ + . In order to explore this possibility, we studied six malonate-based polymers and PEO at a fixed salt concentration (r = 0.08) and temperature (90°C) using symmetric cells with planar electrodes. Unfortunately, there is no correlation between ilim and κρ + . When the applied current density, i, is less than ilim, the cell voltage approaches a stable plateau, ϕ plateau . Here, we found a linear dependence between i and thickness-normalized plateau potential, ϕ plateau L –1 , irrespective of the magnitude of the applied current. In all seven polymer electrolytes, we found a linear correlation between ilim and the slopes of these lines, σ. In other words, measurements of σ can be used to predict the limiting current.
Plants play a critical role in the cycling of potassium (K) and the fractionation of its isotopes. However, little is known about K stable isotopic compositions in natural soil–plant systems and possible fractionation during intra-plant transport and root-soil uptake of K. Herein this study focuses on K isotopic fractionation within a humid and an arid soil–plant system sampled on the windward and leeward sides of Kohala Mountain, Hawai‘i. We determined the K isotopic compositions of < 2-mm bulk soil, soil saturation extraction, and selected plant tissues by multi-collector inductively coupled plasma mass spectrometry and X-ray absorption spectroscopy. We studied soils and individual tissue samples such as roots, stems, barks, shoots (a sum of stems and fresh leaves), leaves (fresh and dead), seeds, and flowers of trees and grasses. The results demonstrated that: (i) tissue δ 41 K values ranged from –1.06 ± 0.06 to 1.15 ± 0.09‰; (ii) within the same plant, stems (barks), dead leaves, and reproductive tissues (flowers and seeds) were isotopically lighter compared to fresh leaves, and to a lesser extent, roots; (iii) δ 41 K values of the humid soil (-0.54 ± 0.07 to –0.49 ± 0.06‰) were lower than those of the arid soil (-0.24 ± 0.07 to –0.14 ± 0.06‰); and (iv) soil bioavailable pool δ 41 K (saturation extracts) ranged from –0.63 ± 0.08 to 0.34 ± 0.08‰ and 0.48 ± 0.08 to 0.54 ± 0.10‰ in the humid and arid soils, respectively. From synchrotron-based analysis of K atoms, we identified two major K-bearing phases co-existing as ionic K + and K-pectate association of different fractions. Based on K isotopic and synchrotron data, we conclude that two dominant processes are responsible for plant-mediated K isotopic fractionation, including (1) K redistribution during intra-plant circulation and (2) uptake at the root-soil interface. For intra-plant circulation of K, there is a high affinity of isotopically lighter K to organic complexes as K-pectate, and K-pectate is particularly enriched in roots and fresh leaves. For K uptake at root-soil interface, isotopically lighter K is preferentially taken by roots from soil bioavailable pools following a low-affinity (passive) transport path. Soil K budget in two sites reflects strong source mixing effects with limited plant imprints. This work provides exploratory data on the biogeochemical fractionation of K isotopes in the soil–plant system.
The discrepancy between the observed and computed flux of solar neutrinos can be eliminated if the temperature gradient in the core of the model, and hence its central temperature, is reduced. This could in principle be accomplished by reducing the core opacity; alternatively, it has been proposed that a fraction of the energy transport in the solar core results from the motion of hypothetical particles (the so-called WIMPs). The resulting changes in solar structure have measurable effects on the frequencies of the solar 5-minute oscillations. Here, models are considered where the opacity in the solar core has been decreased in a manner which roughly simulates the effects of the WIMPs. The analysis of the models and their frequencies provides insight into the consequences of modifications to the physics of solar models. In particular, it is found that models with the observed neutrino flux are inconsistent with observations of low-degree solar oscillations.
In this work, we present a site-specific three-dimensional analysis of porosity evolution in neutron-irradiated U–10Zr metallic fuel using high-resolution synchrotron X-ray tomography. Focused ion beam cubic lift-outs from four radial positions—fuel center, middle, edge, and fuel-cladding interaction (FCCI) zones—were imaged, reconstructed, and segmented to quantify pore volume, density, morphology, and pore connectivity. Porosity increased modestly from 5.5% at the center to 10.5% at the edge, yet pore number density increased by over two orders of magnitude in the fuel portion near the FCCI interface (from 4.8 x 10 4 to 6.0 x 10 6 pores/mm 3 ). Morphological classification revealed a progression from small spherical pores at the center to equiaxed and tortuous networks at the periphery, with enhanced orientation along the radial direction. Connectivity and permeability analysis reveal that the FCCI region maintains dense, highly interconnected pores despite a moderate volume fraction, enabling rapid fission gas transport and lanthanide migration. Effective thermal conductivity models incorporating sodium logging confirm that these site-specific pore features critically influenced heat transport during reactor irradiation. These findings demonstrate that pore topology—not just porosity fraction—controls thermal performance and thermal pathways in U–10Zr fuels, impacting fuel thermal performance in a reactor.
Glassy solids evolve towards lower-energy structural states by physical aging. This can be characterized by structural relaxation times, the assessment of which is essential for understanding the glass’ time-dependent property changes. Conducted over short times, a continuous increase of relaxation times with time is seen, suggesting a time-dependent dissipative transport mechanism. By focusing on micro-structural rearrangements at the atomic-scale, we demonstrate the emergence of sub-diffusive anomalous transport and therefore temporal fractional diffusion in a metallic glass, which we track via coherent x-ray scattering conducted over more than 300,000 s. At the longest probed decorrelation times, a transition from classical stretched exponential to a power-law behavior occurs, which in concert with atomistic simulations reveals collective and intermittent atomic motion. Our observations give a physical basis for classical stretched exponential relaxation behavior, uncover a new power-law governed collective transport regime for metallic glasses at long and practically relevant time-scales, and demonstrate a rich and highly non-monotonous aging response in a glassy solid, thereby challenging the common framework of homogeneous aging and atomic scale diffusion.
Abstract This work presents a new reactive transport framework that combines a powerful geochemistry engine with advanced numerical methods for flow and transport in subsurface fractured porous media. Specifically, the PhreeqcRM interface (developed by the USGS) is used to take advantage of a large library of equilibrium and kinetic aqueous and fluid-rock reactions, which has been validated by numerous experiments and benchmark studies. Fluid flow is modeled by the Mixed Hybrid Finite Element (FE) method, which provides smooth velocity fields even in highly heterogenous formations with discrete fractures. A multilinear Discontinuous Galerkin FE method is used to solve the multicomponent transport problem. This method is locally mass conserving and its second order convergence significantly reduces numerical dispersion. In terms of thermodynamics, the aqueous phase is considered as a compressible fluid and its properties are derived from a Cubic Plus Association (CPA) equation of state. The new simulator is validated against several benchmark problems (involving, e.g., Fickian and Nernst-Planck diffusion, isotope fractionation, advection-dispersion transport, and rock-fluid reactions) before demonstrating the expanded capabilities offered by the underlying FE foundation, such as high computational efficiency, parallelizability, low numerical dispersion, unstructured 3D gridding, and discrete fraction modeling.
Vertical transport by moist sub-grid scale processes such as deep convection is a well-known source of uncertainty in CO2 source/sink inversion. However, a dynamical link between vertical transport, satellite based retrievals of column mole fractions of CO2, and source/sink inversion has not yet been established. By using the same offline transport model with meteorological fields from slightly different data assimilation systems, we examine sensitivity of frontal CO2 transport and retrieved fluxes to different parameterizations of sub-grid vertical transport. We find that frontal transport feeds off background vertical CO2 gradients, which are modulated by sub-grid vertical transport. The implication for source/sink estimation is two-fold. First, CO2 variations contained in moist poleward moving air masses are systematically different from variations in dry equatorward moving air. Moist poleward transport is hidden from orbital sensors on satellites, causing a sampling bias, which leads directly to small but systematic flux retrieval errors in northern mid-latitudes. Second, differences in the representation of moist sub-grid vertical transport in GEOS-4 and GEOS-5 meteorological fields cause differences in vertical gradients of CO2, which leads to systematic differences in moist poleward and dry equatorward CO2 transport and therefore the fraction of CO2 variations hidden in moist air from satellites. As a result, sampling biases are amplified and regional scale flux errors enhanced, most notably in Europe (0.43+/-0.35 PgC /yr). These results, cast from the perspective of moist frontal transport processes, support previous arguments that the vertical gradient of CO2 is a major source of uncertainty in source/sink inversion.
Turbulent thermal convection is of considerable importance in fluid dynamical transport phenomena occurring, for example, in the planetary boundary layer of the Earth, the interiors of stars, and accretion disks. In particular, during a significant portion of the evolutionary phase of many stars having convectively unstable cores or outer envelopes, a substantial fraction of energy is transported from the central layers to the outer layers by thermal convection. Moreover, as much of the interior of a star is in highly turbulent motion, a complete theory of stellar structure and evolution requires the explicit consideration of turbulence in order to have expressions for the turbulent quantities arising in the stellar structure equations, and particularly, the turbulent fluxes that appear in the total flux conservation equation, such as the convective flux, kinetic energy flux, etc. A reliable quantification of these fluxes continues to present a challenge in astrophysical fluid dynamics, primarily because astrophysical turbulence is almost always fully-developed and nearly inviscid, and therefore governed by strong nonlinear interactions that distribute the energy among a very wide spectrum of eddies with scales ranging from the characteristic dimension of the flow to those sufficiently small to be affected by viscous dissipation. Furthermore, astrophysical flows are invariably compressible, anisotropic, and inhomogeneous, which requires the consideration of the dynamics of longitudinal modes and their interaction with the transverse modes, as well as complicated boundary conditions. In order to reach a compromise between analytical and numerical tractability and the basic physics of turbulent convection, we have constructed a model of stationary turbulent convection that yields various turbulence statistics, including the convective flux, that are required in stellar evolution models.
Membranes are essential components of photoelectrochemical CO 2 reduction cells (PEC-CRC), as they regulate solute transport and impact device efficiency. These cells are a promising approach for converting carbon dioxide into valuable chemicals. However, electrochemical reactions involve not just a single electrolyte; membrane transport typically includes multiple electrolytes. When a mixture of solutes is present, the presence of one solute in the membrane can impact the diffusion and sorption of other solutes. In complex systems where multiple solutes are diffusing simultaneously through a hydrated, dense polymer membrane, the resulting transport behavior is poorly understood. This study highlights the individual and two solute transport behavior of potassium formate, acetate, and bicarbonate ions in neutral membranes and quaternary ammonium-functionalized anion exchange membranes (AEMs) based on phenyl acrylate polymer backbones. By isolating the effect of water content, we emphasize how fixed charge density and solute-membrane interactions influence transport behavior. In neutral membranes, the diffusivity follows the order of potassium formate > potassium acetate > potassium bicarbonate. However, in AEMs, the order shifts due to changes in ion-membrane interactions, with potassium bicarbonate diffusivity surpassing that of potassium acetate. These results underscore the importance of ion exchange membranes and ion-membrane interaction in modulating transport properties, providing valuable insights for optimizing membranes in electrochemical applications.
The OMFIT STEP (Meneghini et al. , Nucl. Fusion , vol. 10, 2020, p. 1088) workflow has been used to develop inductive and steady-state H-mode core plasma scenario use cases for a $B_0 = 8 \, {\rm T}$ , $R_0 = 4 \, {\rm m}$ machine to help guide and inform future higher-fidelity studies of core transport and confinement in compact tokamak reactors. Both use cases are designed to produce 200 MW or more of net electric power in an up-down symmetric plasma with minor radius $a = 1.4 \, {\rm m}$ , elongation $\kappa = 2.0$ , triangularity $\delta = 0.5$ and effective charge $Z_{{\rm eff}} \simeq 2$ . Additional considerations based on the need for compatibility of the core with reactor-relevant power exhaust solutions and external actuators were used to guide and constrain the use case development. An extensive characterization of core transport in both scenarios is presented, the most important feature of which is the extreme sensitivity of the results to the quantitative stiffness level of the transport model used as well as the predicted critical gradients. This sensitivity is shown to arise from different levels of transport stiffness exhibited by the models, combined with the gyroBohm-normalized fluxes of the predictions being an order of magnitude larger than other H-mode plasmas. Additionally, it is shown that although heating in both plasmas is predominantly to the electrons and collisionality is low, the plasmas remain sufficiently well coupled for the ions to carry a significant fraction of the thermal transport. As neoclassical transport is negligible in these conditions, this situation inherently requires long-wavelength ion gyroradius-scale turbulence to be the dominant transport mechanism in both plasmas. These results are combined with other basic considerations to propose a simple heuristic model of transport in reactor-relevant plasmas, along with simple metrics to quantify coupling and core transport properties across burning and non-burning plasmas.
Prompt loss of beam injected fast ions approaches 100% in lithium tokamak experiment-beta (LTX-β) discharges, though significantly improved confinement is expected for the higher current plasmas made available by a recent upgrade to the Ohmic heating power supply. Modeling of fast ions using TRANSP/NUBEAM finds a maximum coupled beam fraction of 76% at the near-term limits of the LTX-β operating space. The full ion orbit code POET is employed to validate NUBEAM results against possible non-adiabatic effects on fast ion orbits, but corrections to the prompt loss fraction due to collisionless transport are found to be small. The graphical method code CONBEAM is used to investigate the topology of fast ion phase space as it relates to neutral beam deposition, and counter-injected NBI is considered as a way to access a region of high field side beam deposition. A metric is developed within the CONBEAM using a beam filament model to estimate the prompt loss fraction and shown to agree well with both POET and NUBEAM, enabling near real-time analysis and potential feedback to operators between plasma discharges.
The data set consists of lead (Pb) isotopic ratios (206Pb/207Pb and 208Pb/207Pb, as well as 206Pb/204Pb, 207Pb/204Pb, 208Pb/204Pb) for particulate samples collected from snow and dry deposition. This was part of a study (O’Day et al. 2020) that was looking at the speciation of phosphorous in splits of these samples to evaluate the flux and bio-availability of phosphorous to the mountain sites in the study. The Pb isotopic data allow sourcing of long-range transported Pb, specifically the fraction of Pb sourced in Asia, and so a comparison of the relative distances of PM sources between the high and low elevation sites. Samples were collected at the East River Watershed, CO (dry deposition collected over Sept.-Oct. 2016, show sampling winter 2017) and at the Southern Sierra Critical Zone Observatory (dry deposition collected over Sept-Nov 2016, July-Sept 2017, June-July 2017). The data includes Pb concentration data, as well as calculated fraction of the Pb sourced in Asia. This data set was published in part in O’Day et al. (2020) Phosphorus Speciation in Atmospherically Deposited Particulate Matter and Implications for Terrestrial Ecosystem Productivity. ES&T, https://dx.doi.org/10.1021/acs.est.9b06150
Abstract This study aimed to better understand how tillage row orientation with respect to dominant flow‐pathway along hillslope impacts runoff and the transport of different sediment size fractions. Experimental plots were constructed in contour ridge till (CRT) and parallel ridge till (PRT) sites to monitor runoff and sediment fluxes. Particle size fractions of the sediment, along with organic C and N contents were measured to quantify physical and chemical enrichment ratios for the two tillage orientations. In the CRT plot, tillage produced large oriented roughness elements along the contours, which acted as little check dams, while in the PRT site, the roughness elements helped confine and concentrate the runoff. In the CRT, the “check dams” resulted in a runoff coefficient of .03, while in the PRT, the flow confinement between rows produced a runoff coefficient of .82. Moreover, the erosion rates at the CRT site were 97% less than those in the PRT plot. Large contours produced finer sediment fractions due to the selective sorting in ponded furrows. More aggregate sediment fractions were present in the PRT site, which was dominated by rill erosion. Physical enrichment revealed selective entrainment of finer sediment particles during erosion. Finer‐grain particles with higher specific surface areas, attached more organic C resulting in chemical enrichment. Physical and chemical enrichment methods were found to be in good agreement. These findings suggest transport models that can simulate size fraction updates to the soil active layer can be used to estimate SOC redistribution and hence more accurate C budgets.
The large gradual solar-energetic-particle (SEP) events, where abundances are commonly measured, are produced when coronal mass ejections (CMEs) drive shock waves through the corona and the interplanetary medium. The shock accelerates particles from the highly-ionized, approximately 1.5 MK, plasma in a manner that depends only weakly upon the Q/A of the ion, except at very high energies. Averaging the approximately 1 MeV/amu abundances over many events compensates for the acceleration effects to produce abundances that appear to correspond directly to those in the coronal source for all observed elements, including H. The resulting abundances reflect the 4 x enhancement of ions with low values of first ionization potential (FIP) arising from ion-neutral fractionation that occurs as the atoms are transported up from the photosphere. A different pattern of fractionation is found for ions that are shock-accelerated from the high speed solar wind emerging from coronal holes.
A regional-scale estimate of the fraction of biomass burning emissions that are transported to the free troposphere by deep convection is presented. The focus is on CO and the study region is a part of Brazil that underwent intensive deforestation in the 1980s. The method of calculation is stepwise, scaling up from a prototype convective event, the dynamics of which are well-characterized, to the vertical mass flux of carbon monoxide over the region. Given uncertainties in CO emissions from biomass burning and the representativeness of the prototype event, it is estimated that 10-40 percent of CO emissions from the burning region may be rapidly transported to the free troposphere over the burning region. These relatively fresh emissions will produce O3 efficiently in the free troposphere where O3 has a longer lifetime than in the boundary layer.
Alkenes formed, during the refining of crude oil, by cracking the heavier fractions are present in transportation fuels in significant amounts, up to as much as 15–20% in gasoline. Moreover, alkenes are also the major intermediate products of the oxidation of alkanes, which play a significant role in autoignition chemistry. This review has assessed the recent progress in gas-phase detailed kinetic model development for species with C=C double bond, mostly C2–C10 alkenes and 1,3-butadiene. The compiled knowledge on alkene combustion chemistry enabled a better understanding of the influence of the number and the position of the C=C double bond on the chemical kinetics and hence combustion behavior of alkenes in engines. At first, the article gives an extensive overview of fundamental combustion experiments by considering studies of C2–C10 alkenes and 1,3-butadiene in shock tubes, rapid compression machines, laminar flames, and jet-stirred and flow reactors. The value of the data from such experiments is critically discussed. Secondly, this article highlights the important reaction classes involved in alkene oxidation over low-, intermediate- and high-temperature ranges. Combustion chemistry covering C2 to C10 alkenes, with a special emphasis on C2 to C7 isomers is discussed by presenting a large body of experimental and modeling investigations. Detailed chemistry differences between alkene isomers and also between alkenes and alkanes are also addressed. Thirdly, the article presents important reaction pathways for PAH precursor formation in different alkenes. Finally, a summary of the distinguishing features of alkene combustion chemistry and an outlook towards future research in this area are presented. In conclusion, this review is focused on linear and branched chain alkenes, and the chemistry of cyclo-alkenes is not included.