Fractal textured surfaces for high temperature corrosion mitigation in molten salts
Not provided.
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 provided.
This effort employs and adapts an existing, rigorous multiscale electricity modeling platform at the National Renewable Energy Laboratory (NREL) to evaluate carbon capture and sequestration (CCS) and negative emissions technologies (NET) from the ARPA-E FLECCS program. NREL's modeling platform includes the Regional Energy Deployment System (ReEDS) electric sector capacity expansion model, which projects future electricity generation mixes at sub-state-level resolution that are downscaled to the unit-level to enable hourly, zonal or nodal electricity production cost modeling in the PLEXOS model. The resulting hourly price data from PLEXOS is provided to technology developers under the ARPA-E FLECCS program to enable technology-specific economic analysis. The ReEDS-PLEXOS modeling suite is well-established for examining electric sector futures with high renewable energy penetrations, energy storage, electrification, and distributed generation. The key advancement proposed herein utilizes collaboration with CCS experts at the University of Wyoming and the FLECCS teams to create innovative methods for representing CCS and NET in the ReEDS and PLEXOS models. Expanded technology options, new operational parameterizations, and detailed data defining CO2 capture, transportation, and storage systems are being integrated into these models to allow an unprecedented combination of scope and resolution for exploring the future of CCS and NET. The resulting capabilities will take advantage of high-performance computing resources to permit wide-ranging scenario analysis to assess CCS and NET deployment under alternative CO2 prices, fossil fuel prices, electricity demand growth, and other electric sector characteristics. Final outcomes will include publicly available hourly grid operation and price data for any U.S. region of interest along with open-access capacity expansion tools for evaluating CCS/NET systems. These products will help an emerging CCS/NET industry in the United States by providing economics-driven guidance to technology developers while informing policy and investment decisions in the public and private sectors.
The tortuosity factor of composite cathodes significantly affects the rate performance of all-solid-state batteries (ASSBs) and has significant differences from systems with liquid electrolytes. Here, in this work, we report a simulation-experiment combined approach that quantifies the effective Li-ion transport in an ASSB composite cathode, which links tortuosity factor on ∼ μm scale to terminal voltage during cycling at the cell level (on ∼ cm scale). Two independent approaches of tortuosity factor quantification are considered: fitting electrochemical cycling data and verifying at different cycling rates and calculating from segmented tomography images, with the tortuosity factor quantified from both methods reaching self-consistency. The simulated terminal voltage using the quantified tortuosity factor has a small relative error of <3% compared to the experimental measurements. We find a significantly reduced value of the Bruggeman exponent of the catholyte phase (1.75), and using shape analysis, we show that rod-shaped catholyte particles play an important role in lowering the tortuosity factor.
Post combustion carbon capture by solvent absorption in a structured packed column is a promising technology for mitigating greenhouse gas emissions. Computational fluid dynamic (CFD) modelling of such a column is a challenging multiscale problem due to the range in length and time scales. The microscale hydrodynamics play a key role in the overall column efficiency with the interfacial area significantly influencing the mass transfer between the gas and liquid phases. In this context, multiphase flow simulations using the volume of fluid (VOF) method in a representative elementary unit (REU) of the packed column can provide fundamental insights into the microscale hydrodynamics, such as, interfacial area and liquid holdup. The present paper systematically examines the impact of various factors (e.g., physical properties and contact angle) on the interfacial area. The results are compared with existing correlations and a scaling analysis is also performed. The solvent physical properties are characterized by the Kapitza number (Ka), a dimensionless number that depends only on fluid properties. So that the Ka number decreases with increasing viscosity. At a fixed liquid load, the interfacial area and liquid holdup are observed to increase with decreasing Ka number. The impact of contact angle (i.e., solid surface characteristics) is effectively investigated by modifying the wall boundary conditions. The interfacial area and liquid holdup are found to decrease with increasing contact angle. Subsequently, a phenomenological correlation for interfacial area is proposed that includes the impact of these parameters. Finally, this correlation may be used to predict the interfacial area for gas-liquid flow in a structured packing for rivulet to fully wetted flow regimes.
Purpose: The interactions of neoplastic cells with each other and the microenvironment are complex. To understand intratumoral heterogeneity, subtle differences should be quantified. Main factors contributing to heterogeneity include the gradient ischemic level within neoplasms, action of microenvironment, mechanisms of intercellular transfer of genetic information, and differential mechanisms of modifications of genetic material/proteins. This may reflect on the expression of biomarkers in the context of prognosis/stratification. Hence, a rigorous approach for assessing the spatial intratumoral heterogeneity of histological biomarker expression with accuracy and reproducibility is required, since patterns in immunohistochemical images can be challenging to identify and describe. Methods: A quantitative method that is useful for characterizing complex irregular structures is lacunarity; it is a multiscale technique that exhaustively samples the image, while the decay of its index as a function of window size follows characteristic patterns for different spatial arrangements. In histological images, lacunarity provides a useful measure for the spatial organization of a biomarker when a sampling scheme is employed and relevant features are computed. The proposed approach quantifies the segmented proliferative cells and not the textural content of the histological slide, thus providing a more realistic measure of heterogeneity within the sample space of the tumor region. The aim is to investigate in whole sections of primary pancreatic neuroendocrine neoplasms (pNENs), using whole-slide imaging and image analysis, the spatial intratumoral heterogeneity of Ki-67 immunostains. Unsupervised learning is employed to verify that the approach can partition the tissue sections according to distributional heterogeneity. Results: The architectural complexity of histological images has shown that single measurements are often insufficient. Inhomogeneity of distribution depends not only on percentage content of proliferation phase but also on how the phase fills the space. Lacunarity curves demonstrate variations in the sampled image sections. Since the spatial distribution of proliferation in each case is different, the width of the curves changes too. Image sections that have smaller numerical variations in the computed features correspond to neoplasms with spatially homogeneous proliferation, while larger variations correspond to cases where proliferation shows various degrees of clumping. Grade 1 (uniform/nonuniform: 74%/26%) and grade 3 (uniform: 100%) pNENs demonstrate a more homogeneous proliferation with grade 1 neoplasms being more variant, while grade 2 tumor regions render a more diverse landscape (50%/50%). Hence, some cases show an increased degree of spatial heterogeneity comparing to others with similar grade. Whether this is a sign of different tumor biology and an association with a more benign/malignant clinical course needs to be investigated further. The extent and range of spatial heterogeneity has the potential to be evaluated as a prognostic marker. Conclusions: The association with tumor grade as well as the rationale that the methodology reflects true tumor architecture supports the technical soundness of the method. This reflects a general approach which is relevant to other solid tumors and biomarkers. Drawing upon the merits of computational biomedicine, the approach uncovers salient features for use in future studies of clinical relevance.
Root anatomy is an important determinant of root metabolic costs, soil exploration, and soil resource capture. Root anatomy varies substantially within and among plant species. RootSlice is a multicellular functional-structural model of root anatomy developed to facilitate the analysis and understanding of root anatomical phenotypes. RootSlice can capture phenotypically accurate root anatomy in three dimensions of different root classes and developmental zones, of both monocotyledonous and dicotyledonous species. Several case studies are presented illustrating the capabilities of the model. For maize nodal roots, the model illustrated the role of vacuole expansion in cell elongation; and confirmed the individual and synergistic role of increasing root cortical aerenchyma and reducing the number of cortical cell files in reducing root metabolic costs. Integration of RootSlice for different root zones as the temporal properties of the nodal roots in the whole-plant and soil model OpenSimRoot/maize enabled the multiscale evaluation of root anatomical phenotypes, highlighting the role of aerenchyma formation in enhancing the utility of cortical cell files for improving plant performance over varying soil nitrogen supply. Such integrative in silico approaches present avenues for exploring the fitness landscape of root anatomical phenotypes.
Spectral analysis of local atomic environments has become a powerful tool for studying solute atom segregation and interactions at grain boundaries in nanocrystalline alloys. When applied to individual grain boundaries, the spectral analysis has shown that solute-solute interaction can be either attractive or repulsive, with long-range relative attraction enhancing the likelihood for solute atoms to begin clustering. In this article, we combine this analysis with a new grain-boundary structure descriptor based on grain-boundary atom coordination, to investigate the impact of grain-boundary junctions on solute atom segregation in polycrystals. Specifically, we systematically characterize the tendency of solute clusters to begin forming at various types of ordinary grain boundaries, triple junctions, and high-order junctions in Ag polycrystals containing either Ni or Cu solute atoms. Our findings demonstrate that the formation of solute clusters at grain boundaries is primarily driven by long-range relative solute attraction, rather than short-range solute-solute interactions. Furthermore, this effect is most pronounced near grain-boundary junctions. Our study highlights the multiscale nature of solute segregation at crystalline interfaces and provides new insights into the complex phenomena governing heterogeneous solute segregation in grain-boundary networks.
Many science and engineering problems exhibit scale disparity and high contrast. The small scale features cannot be omitted in the physical models because they can affect the macroscopic behavior of the problems. However, resolving all the scales in these problems can be prohibitively expensive. As a consequence, some types of model reduction techniques are required to design efficient solution algorithms. For practical purpose, we are interested in mixed finite element problems as they produce solutions with certain conservative properties. Existing multiscale methods for such problems include the mixed multiscale finite element methods. We show that for complicated problems, the mixed multiscale finite element methods may not be able to produce reliable approximations. This motivates the need of enrichment for coarse spaces. Two enrichment approaches are proposed, one is based on generalized multiscale finte element metthods (GMsFEM), while the other is based on spectral element-based algebraic multigrid (rAMGe). The former one, which is called mixed GMsFEM, is developed for both Darcy’s flow and linear elasticity. Application of the algorithm in two-phase flow simulations are demonstrated. For linear elasticity, the algorithm is subtly modified due to the symmetry requirement of the stress tensor. The latter enrichment approach is based on rAMGe. The algorithm differs from GMsFEM in that both of the velocity and pressure spaces are coarsened. Due the multigrid nature of the algorithm, recursive application is available, which results in an efficient multilevel construction of the coarse spaces. Stability, convergence analysis, and exhaustive numerical experiments are carried out to validate the proposed enrichment approaches. iii
Collisionless shocks are a grand challenge problem in plasma physics and have been the subject of study for more than six decades. Shocks are ubiquitous phenomena in the universe and are responsible for transforming high energy and momentum flows into thermal energy and energetic particles. Understanding how shocks operate is of primary importance to understand the sun-earth coupling, protecting manned missions and spacecrafts from high energy particles, inertial confinement fusion, and radiation observed from astrophysical plasmas, such as supernova remnants and astrophysical jets. A major unanswered question on this frontier is, how does a collisionless plasma transform flow energy and momentum into electron thermal energy? Many potential mechanisms have been proposed to perform the conversion between flow and thermal energy in a collisionless plasma, but the answer has been elusive. To solve the mystery, non-relativistic perpendicular and quasi-perpendicular shocks are studied via a coordinated approach employing cutting-edge numerical simulations with the continuum kinetic code Gkeyll and spacecraft observations supplied by NASA's Magnetospheric Multiscale (MMS) mission. In both sets of data, we employ novel techniques to identify how and where energy is exchanged between electromagnetic fields and plasma particles. By combining cutting-edge spacecraft data and numerical simulations with these diagnostics, we have been able to identify some of the primary mechanisms responsible electron heating and dissipation in heliospheric shocks. In addition to electron dynamics, we have complete access to proton simulation and in situ data, enabling us to characterize proton heating and acceleration in quasi-perpendicular shocks. Many astrophysical quasi-perpendicular shocks have similar parameters to those in this study, including supernova remnants and shocks in galaxy clusters. Additionally, the simulation code, Gkeyll, and associated analysis scripts are openly available for community use. This award has supported a collaborative effort between Princeton University, the University of Arizona - Tucson, and the University of Maryland - College Park. The award funded three early career scientists, including the funding and training of one female postdoctoral researcher, who was involved in all aspects of the MMS data analysis and some portions of the simulation analysis.
Climate and population change are known to influence electricity demand, but what is the impact of uncertainty in climate and population projections on electricity demand in the United States (U.S.) over the next 30 years? This question has important implications for investment decisions in the energy sector, which are typically made using a 15- to 30-year time horizon. Simply put, if future climate and population scenarios do not lead to meteorological conditions and subsequent energy demands that are distinctly different within the first 30 years then, for the purposes of investment decisions, it may not matter which potential pathway we are most likely on. The Integrated, Multisector, Multiscale Modeling (IM3) project has generated a wide, yet plausible range of 21st century high-resolution climate and population scenarios for the U.S. The IM3 projections span two climate scenarios (Representative Concentration Pathways 4.5 and 8.5) and two population scenarios (Shared Socioeconomic Pathways 3 and 5). For each of the climate scenarios, we reflect a range of climate model uncertainty by using boundary conditions from a set of climate models that are hotter and colder than the multi-model mean. In total there are eight (2 x 2 x 2) IM3 scenarios. This work explores a basic question: When do future climate and energy outcomes start to meaningfully diverge across the eight scenarios? The analysis is done on both the raw meteorology time-series from each scenario as well as a time-series of projected total electricity demand generated using IM3's Total ELectricity Loads (TELL) model. We also explore the spatial heterogeneity of the divergence signal to identify regions of the country that may be more or less path dependent.
In Part I of this study, it was shown that the Eurasian cold anomalies related to Arctic warming depend strongly on the quasi stationarity and persistence of the Ural blocking (UB). The analysis here revealed that under weak mean westerly wind (MWW) and vertical shear (VS) (quasi barotropic) conditions with weak synoptic-scale eddies and a large planetary wave anomaly, the growth of UB is slow and its amplitude is small. For this case, a quasi-stationary and persistent UB is seen. However, under strong MWW and VS (quasi baroclinic) conditions, synoptic-scale eddies are stronger and the growth of UB is rapid; the resulting UB is less persistent and has large amplitude. In this case, a marked retrogression of the UB is observed. The dynamical mechanism behind the dependence of the movement and persistence of UB upon the background conditions is further examined using a nonlinear multiscale model. The results show that when the blocking has large amplitude under quasi-baroclinic conditions, the blocking-induced westward displacement greatly exceeds the strong mean zonal-wind-induced eastward movement and hence generates a marked retrogression of the blocking. By contrast, under quasi-barotropic conditions because the UB amplitude is relatively small the blocking-induced westward movement is less distinct, giving rise to a quasi-stationary and persistent blocking. It is further shown that the strong mid–high-latitude North Atlantic mean zonal wind is the quasi-barotropic condition that suppresses UB’s retrogression and thus is conducive to the quasi stationarity and persistence of the UB. The model results show that the blocking duration is longer when the mean zonal wind in the blocking region or eddy strength is weaker.
Two magnetopause reconnection events of the Magnetospheric Multiscale mission with whistler wave activity are presented. The whistler mode around half of the electron cyclotron frequency is excited near the magnetospheric separatrix. In both events, there are positive correlations between the whistler wave and the lower hybrid drift instability (LHDI) activities, suggesting a possible role of LHDI in the whistler wave generation. A sudden change in the electron pitch angle distribution (PAD) function of energetic electrons is observed right after intense LHDI activity. This change in the PAD leads to temperature anisotropy in energetic electrons which is responsible for the whistler wave excitation. The measured dispersion relation demonstrates that the whistler wave propagates toward the X line nearly parallel to the magnetic field line. Furthermore, a linear analysis with the measured distribution function verifies that the whistler mode is excited by the temperature anisotropy in energetic electrons.
Mantis shrimp are aggressive marine crustaceans well known for their rapid and powerful hunting strategies. Less well known, yet, is the ability of some species of mantis shrimp to defend themselves from the repeated blows of conspecifics during ritualized fighting using a shield-like segment of abdominal armor called the telson. Multiscale structure–mechanical property relationships of this damage-tolerant biological composite is examined in order to reveal strategies that nature uses for resisting failure from repeated high-energy impacts. The telson structures of the smashing-type species, Odontodactylus scyllarus , and the less aggressive spearing-type species, Lysiosquillina maculata , are compared in order to better understand the ecological pressures driving the formation and use of the telson as a biological shield. A higher bulk compressive stiffness is identified within the smasher telson, which is attributed to its concave macromorphology, thicker cuticle, and higher degree of mineralization within its exocuticle. The presence of ridges at the dorsal surface indicates a role in imparting compliance for energy absorption. Fracture analysis identifies an enhanced toughening mechanism of crack twisting within the smasher telson, attributed to its well-defined pitch-graded helicoidal fibrous micro-architecture. Such findings may prove useful for the design of lightweight composite materials with potential flexibility and improved damage tolerance.
Abstract Photovoltaic (PV) modules are subjected to environmental stressors (UV exposure, temperature, and humidity) that cause degradation within the encapsulant and its interfaces with adjacent glass and cell substrates. To save experimental time and to enable long‐term assessment with intensive degradation only taking place after many years, the development of predictive models is indispensable. Previous works have modeled the delamination of the ethylene vinyl acetate (EVA) encapsulant/glass and encapsulant/cell interfaces under field aging conditions with fundamental photochemical degradation reactions that lead to molecular scission and loss of interfacial adhesion, characterized by the fracture resistance, G c . However, these models were fundamentally limited in that the following aspects were not incorporated: (i) molecular crosslinking in the field, (ii) synergistic autocatalytic interactions of degradation mechanisms, (iii) connection between degraded encapsulant structure and its mechanical properties, and (iv) rigorous treatment of the plasticity contribution to G c with finite element models. Here, we present a time‐dependent multiscale model that addresses these limitations and is applicable to a wide range of encapsulants and interfaces. For the reference EVA encapsulant and its interfaces with the glass and cell, the presented model predicts an initial rise in G c in the first 3 years of field aging from crosslinking, then a subsequent sharp decline from degradation mechanisms. We used nanoindentation to measure the changes in EVA mechanical properties over exposure time to tune the model parameters. The model predictions of G c and mechanical properties match with experimental data and show an improvement compared to previous models. The model can even predict switches in failure interfaces, such as the observed EVA/cell to EVA/glass transition. We also conducted a sensitivity analysis study by varying the degradation and crosslinking kinetic parameters to demonstrate their effects on G c . Model extensions to polyolefin elastomer‐ and silicone‐encapsulants and their interfaces are also demonstrated.
In this paper, we study the multiscale Boltzmann equation with multi-dimensional random parameters by a bi-fidelity stochastic collocation (SC) method developed in [52], [70], [71]. By choosing the compressible Euler system as the low-fidelity model, we adapt the bi-fidelity SC method to combine computational efficiency of the low-fidelity model with high accuracy of the high-fidelity (Boltzmann) model. With only a small number of high-fidelity asymptotic-preserving solver runs for the Boltzmann equation, the bi-fidelity approximation can capture well the macroscopic quantities of the solution to the Boltzmann equation in the random space. Additionally, a priori estimate on the accuracy between the high- and bi-fidelity solutions together with a convergence analysis is established. Finally, we present extensive numerical experiments to verify the efficiency and accuracy of our proposed method.
Bipolar membrane (BPM) electrochemical processes are a promising platform for carbon dioxide (CO 2 ) separations, but the molecular level thermodynamic and kinetic understanding of CO 2 -to-bicarbonate (HCO 3 − ) transformation remain poorly understood. This study employs a multiscale computational approach to systematically explore the adsorption and reactive transformation of CO 2 in five anion exchange ionomer systems. Classical molecular dynamics (MD) simulation results demonstrate that polymers with imidazolium groups significantly reduce CO 2 diffusion and enhance (OH − )–CO 2 interactions due to stronger electrostatic and π-interactions. Compared to the commonly used quaternary ammonium ionomers, imidazolium-functionalized ionomers show improved CO 2 proximity and interaction strength. Ab initio MD and density functional theory (DFT) calculations reveal that the benzyl-substituted imidazolium (IM-Ben) substantially reduces the energy barrier for HCO 3 − formation (∼72 meV lower) compared to the alkyl-substituted IM-nBu, while also mitigating imidazolium deprotonation under moderate hydration conditions. Transition state analysis shows IM-Ben forms more extensive hydrogen-bonding networks, which stabilize the transition state structure and contribute to a lower energy barrier for bicarbonate formation. These findings highlight the advantage of the adjacent benzyl moiety in enabling efficient CO 2 -to-bicarbonate transformation via hydrated hydroxide ion counterions, offering mechanistic insights and clear molecular design principles for optimizing anion exchange ionomers at bipolar membrane interfaces for electrochemical CO 2 separation applications.
Explore the source record for details and available documents.
The mapper graph is a popular tool from topological data analysis that provides a graphical summary of point cloud data. It has been used to study data from cancer research, sports analytics, neurosciences, and machine learning. In particular, mapper graphs have been used recently to visualize the topology of high-dimensional artificial neural activations from convolutional neural networks and large language models. However, a key question that arises from using mapper graphs across applications is how to compare mapper graphs to study their structural differences. In this paper, we introduce a distance between mapper graphs using tools from optimal transport. We demonstrate the utility of such a distance by studying the topological changes of neural activations across convolutional layers in deep learning, as well as by capturing the loss of structural information for multiscale mapper.