Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mass diffusivity”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

An extension of the localized artificial diffusivity method for immiscible and high density ratio flows

The localized artificial diffusivity (LAD) method is widely regarded as the preferred multi-material regularization scheme for the compact finite difference method, because it is conservative, easy to implement, and generally robust for a wide range of multi-material problems. However, traditional LAD methods face significant challenges when applied to flows with large density ratios and when maintaining thermodynamic equilibrium across material interfaces. These limitations arise from the formulation of the artificial diffusivity flux and the reliance on enthalpy diffusion for interface regularization. Additionally, traditional LAD methods struggle to ensure stability under large density ratio conditions, fail to maintain a finite interface thickness, and are therefore unsuitable for modeling immiscible interfaces. Here, in this work, we discuss the origins of these issues in traditional LAD methods and propose modifications which enable the simulation of large density ratio and immiscible flows. The proposed method targets the artificial diffusion fluxes at gradients and ringing in the volume fraction, rather than the mass fraction in traditional methods, to consistently regularize large density ratio interfaces. Furthermore, the proposed method introduces an artificial bulk density diffusion term to enforce equilibrium conditions across interfaces. To address the challenge of modeling immiscible flows, a conservative diffuse interface term is incorporated into the formulation to ensure a finite interface thickness. Specific consideration is taken in the design of the method to ensure that these crucial properties are maintained for N -material flows. The effectiveness of the proposed method is demonstrated through a series of canonical test cases, and its accuracy is validated by comparison with experimental data on micro-bubble collapse in water. These results highlight the method’s robustness and its ability to overcome the limitations of traditional LAD approaches.

Artificial diffusivity↗
DOE OSTI2020

Fission Product Diffusion in Silicon Carbide: As Revealed by Computer Simulations and Experimental Measurements

The diffusion and release of fission products through silicon carbide in the tri-structural isotropic (TRISO) fuel particles remain unsolved for decades. The underlying mechanism is quite challenging to be determined. To help unveil the mysterious story, the current work applies molecular dynamics method to show the stability of silver, palladium, ruthenium and iodine as an interstitial and their atomic diffusion along coincident site lattice (CSL) boundary, especially ?3 grain boundary (GB). The major finding presents a much faster diffusion along GB than in bulk for all elements considered. The reasonably close estimate to experiments and simulations where available has confirmed the important role of grain boundary diffusion of Ag and Pd in SiC. However, the discrepancy addressed in Ag with measurements from fuel studies suggest a more complicated mechanism, which might be in correlation with high energy grain boundaries or the presence of crack. The subsequent characterization of Ru and I distribution in SiC-PyC-SiC diffusion couples, which have been ion irradiated at 900?C to 10 dpa and 20 dpa, has performed by secondary ion mass spectrometry (SIMS) analysis. The experiment measurements correlate well with the grain boundary diffusion by simulation, which provide further evidences that the grain boundary diffusion cannot be neglected once the fission products are accessible at the grain boundary.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

A phase-field diffraction model for thermo-hydro-mechanical propagating fractures

This paper introduces a novel diffraction based thermo-hydraulic–mechanical (THM) model for fracture propagation using a phase-field fracture (PFF) approach. The key innovation of the THM-PFF model lies in its integrated treatment of four solution variables—displacements, phase-field, pressure, and temperature—each governed by a combination of conservation of momentum (mechanics problem), a variational inequality (constrained minimization problem), mass conservation (pressure problem), and energy conservation (temperature problem). This leads to a new formulation of a coupled variational inequality system. A major advancement is the development of an extended fixed-stress algorithm, where displacements, phase-field, pressures, and temperatures are solved in a staggered sequence. An important aspect of this work is the global coupling of pressures and temperatures across the domain using diffraction systems, with diffraction coefficients defined by material parameters weighted by the diffusive phase-field variable. To ensure robust local mass conservation, we employ enriched Galerkin finite elements (EG) for both pressure and temperature diffraction equations. By enriching the continuous Galerkin basis functions with discontinuous piecewise constants, EG accurately represents solution and parameter discontinuities while preserving local mass and energy conservation—crucial aspects for THM problems and realistic behavior. Moreover, the use of a predictor–corrector local mesh adaptivity scheme is employed, allowing the model to handle small phase-field length-scale parameters while maintaining high numerical accuracy and reasonable computational cost. Furthermore, these new model and algorithmic developments represent significant advances in the field and have been substantiated through rigorous numerical tests.

Diffraction systems↗
DOE OSTI2021

Large-eddy simulation of a multi-injection flame in a diesel engine environment using an unsteady flamelet/progress variable approach

In this work, large-eddy simulations (LESs) are conducted for a multiple-injection flame in a diesel engine environment using an unsteady flamelet/progress variable (UFPV) approach in which differential diffusion is considered. The suitability of the UFPV tabulation approach is first evaluated through a priori analyses using the state-of-the-art direct numerical simulation (DNS) dataset. Both the instantaneous data and the conditional values for the major and minor species' mass fractions are compared between the UFPV and the DNS. The comparisons show that the proposed UFPV tabulation approach can give good predictions for the multiple-injection flame at different injection phases. While the gas temperature and major species mass fractions can be accurately predicted with or without differential diffusion being considered in the UFPV flamelet library, the prediction accuracy for the highly diffusive species (e.g., hydrogen) in the main injection phase can be noticeably improved when differential diffusion is taken into account. The fully coupled LES/UFPV simulations show that the overall structure of the multiple-injection flame can be predicted, and the conditional thermo-chemical values are close to the filtered DNS dataset. The reasons for the remaining discrepancies found in the a priori analyses and the a posteriori simulations using the UFPV approach are analyzed.

Wen, Xu (ORCID:0000000331334292)↗
DOE OSTI2024

Modeling the Progenitor Stars of Observed Type IIP Supernovae

Abstract Type IIP supernovae (SNe IIP) are thought to originate from the explosion of massive stars >10M ⊙ . Their luminosity is primarily powered by the explosion energy and the radioactive decay energy of 56 Co, with the photosphere location regulated by hydrogen recombination. However, the physical connections between SNe IIP and their progenitor stars remain unclear. This paper presents a comprehensive study of SNe IIP and their progenitor stars by using the one-dimensional stellar evolution code,MESA. Our model grids consider the effects of stellar metallicity, mass, and rotation in the evolution of massive stars, as well as the explosion energy and 56 Ni production in modeling supernovae. To elucidate the observed SNe IIP and their origins, we compare their light curves (LCs) with our models. Furthermore, we investigate the impact of stellar parameters on LCs by considering stellar mass, metallicity, rotation, explosion energy, and 56 Ni production. We find that more massive stars exhibit longer plateaus due to increased photon diffusion time caused by massive ejecta. Higher metallicity leads to increased opacity and mass loss of progenitor stars. Rapid rotation affects internal stellar structures, enhancing convective mixing and mass loss, potentially affecting the plateau’s brightness and duration. Higher explosion energy results in brighter but shorter plateaus due to faster-moving ejecta. 56 Ni mass affects late-time luminosity and plateau duration, with larger masses leading to slower declines.

Astronomy & Astrophysics↗
DOE OSTI2021

Bulk and Short-Circuit Anion Diffusion in Epitaxial Fe 2 O 3 Films Quantified Using Buried Isotopic Tracer Layers

Self-diffusion is a fundamental physical process that, in solid materials, is intimately correlated with both microstructure and functional properties. In this work, a universal approach is presented to precisely characterize self-diffusion in ionic solids by isotopically enriching anions and/or cations at specific locations within an epitaxial film stack, and characterize their redistribution at high spatial resolution with atom probe tomography. Nanoscale anion diffusivity is quantified in epitaxial α-Fe 2 O 3 thin films deposited by molecular beam epitaxy with a thin (10 nm) buried layer highly enriched in 18 O. The isotopic sensitivity of the atom probe allows precise measurement of 18 O distribution across the sharp interfaces between this layer and the surrounding Fe 2 O 3 after annealing. Short-circuit anion diffusion through 1D and 2D structural defects in Fe 2 O 3 are also directly visualized in 3D. This versatile approach to study precisely tailored thin film samples at high spatial and mass fidelity will facilitate a deeper understanding of atomic-scale diffusion phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Phosphorus diffusion and deactivation during SiGe oxidation

Dopant profiles near the semiconductor–oxide interface are critical for microelectronic device performance. As the incorporation of Si 1−x Ge x into transistors continues to increase, it is necessary to understand the behavior of dopants in Si 1−x Ge x . In this paper, the diffusion and electrical activation of phosphorus within a strained, single-crystal Si 0.7 Ge 0.3 layer on Si during oxidation are reported. Both layers were uniformly doped, in situ, with an average phosphorus concentration of 4 × 10 19 atoms/cm 3 . After high-temperature oxidation, secondary ion mass spectrometry measurements revealed that the bulk of the phosphorus diffuses out of only the SiGe layer and segregates at the oxidizing SiGe–SiO 2 interface. Hall effect measurements corroborate the observed phosphorus loss and show that the phosphorus diffusing to the oxidizing interface is electrically inactive. Through density functional theory (DFT) calculations, it is shown that phosphorus interstitials prefer sites near the SiGe–SiO 2 interface. Finally, based on a combination of experimental data and DFT calculations, we propose that the phosphorus atoms are displaced from their lattice sites by Ge interstitials that are generated during SiGe oxidation. The phosphorus atoms then migrate toward the SiGe–SiO 2 interface through a novel mechanism of hopping between Ge sites as P–Ge split interstitials. Once they reach the interface, they are electrically inactive, potentially in the form of interstitial clusters or as part of the reconstructed interface or oxide.

Thornton, Chappel S.↗
DOE OSTI2023

Kinetic model describing self-limiting CO 2 diffusion in supported amine adsorbents

A reaction–diffusion shrinking core model describing the decay in diffusivity of supported amine sorbents upon CO 2 sorption under both simulated direct air capture and point source capture conditions is described. The decay in CO 2 diffusivity is associated with crosslinking in the aminopolymer samples and general pore blockage in the amino-silane derived samples, which occurs as CO 2 is adsorbed. The model is used to extract four kinetic parameters that govern the CO 2 uptake kinetics and working capacity: an apparent reaction rate constant, an initial effective diffusivity, and two dimensionless decay parameters. Ideally, an initially reaction limited system would allow for direct determination of the intrinsic reaction rate constant; however, sorption experiments suggest mass transfer resistances related to gas mixing, external boundary layers and intraparticle diffusion are present. Reaction rate constants are determined and agree well with theoretical values predicted with the Eyring equation parameterized using density functional theory energies from literature sources. The kinetic performance is expressed as the average effective diffusivity as a function of average conversion, which can be correlated to the dispersion of sorption sites on the support and the morphology of the active sorbent phase. Four supports are impregnated or grafted with amines, SBA-15, single-walled zeolite nanotubes (ZNT), Syloid SiO 2 , and γ-Al 2 O 3 . Due to its pore structure, γ-Al 2 O 3 supported amines can be modeled at the μm scale or at the nm scale, where the shell balance is on the μm-sized macroporous particle aggregate or on the nm-sized amine film on the surface of the Al 2 O 3 nanoparticles, which comprise the spherical particle aggregates. Faster diffusion rates are maintained under 400 ppm rather than 10% CO 2 due to a slower reaction rate giving a slower decay in diffusivity. In conclusion, this work provides a first principles kinetic analysis of CO 2 sorption where previous models are semi-empirical and use arbitrary kinetic parameters.

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

Designing a Zn–Ag Catalyst Matrix and Electrolyzer System for CO 2 Conversion to CO and Beyond

CO 2 emissions can be transformed into high-added-value commodities through CO 2 electrocatalysis; however, efficient low-cost electrocatalysts are needed for global scale-up. Inspired by other emerging technologies, the authors report the development of a gas diffusion electrode containing highly dispersed Ag sites in a low-cost Zn matrix. Here, this catalyst shows unprecedented Ag mass activity for CO production: –614 mA cm –2 at 0.17 mg of Ag. Subsequent electrolyte engineering demonstrates that halide anions can further improve stability and activity of the Zn–Ag catalyst, outperforming pure Ag and Au. Membrane electrode assemblies are constructed and coupled to a microbial process that converts the CO to acetate and ethanol. Combined, these concepts present pathways to design catalysts and systems for CO 2 conversion toward sought-after products.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Temperature-dependent diffuse reflectance measurements of ceramic powders in the near- and mid-infrared spectra

Radiative properties are critical to quantify radiative energy fluxes between surfaces and in participating media. However, there is limited experimental data on temperature-dependent radiative properties of materials. Herein, this work focuses on experimentally measuring temperature-dependent diffuse reflectance in the near- and mid-infrared spectra (1–20 μm) for ceramic particles with applications as heat-transfer and thermal-storage media in concentrated solar power (CSP) plants. Specifically, a commercially available sintered bauxite proppant, ACCUCAST ID80, and its primary chemical constituents—alumina (Al 2 O 3 ) and silica (SiO 2 )—are measured in powder form using a Fourier transform infrared spectrometer (FTIR) coupled with a specialized diffuse reflectance accessory and a heated stage. Room-temperature diffuse reflectance measurements show increased absorption in tests with greater mass fractions of the ceramic samples. There is a strong correlation in the measured reflectance spectra of ACCUCAST with alumina and silica in the spectral range 2000–500 cm –1 (5–20 μm). Whereas, for shorter wavelengths (< 5 μm), the absorptance for ACCUCAST is greater than the absorptance for alumina and silica, indicating contributions from other chemical species present in the composite material. For the first time, temperature-dependent diffuse reflectance measurements are reported for ACCUCAST up to 1000 °C. These results are compared against those of alumina and silica through the calculation of a thermal emittance. All three materials exhibit a calculated emittance of ~0.9 at room temperature. However, this value decreases to ~0.6 for ACCUCAST and drops to less than 0.4 for alumina and silica at 1000 °C. Thermal cycling in air at 1000 °C resulted in a visible color change from dark grey to light orange for ACCUCAST and a subsequent larger increase in reflectance for wavelengths less than 5 μm as compared to ACCUCAST thermally cycled at 1000 °C in vacuum. Alumina and silica spectra proved to be largely unaffected by thermal cycling under atmospheric conditions. Overall, this study establishes a powerful technique for the characterization of radiative properties of particulate materials as a function of temperature and presents a detailed case study of ACCUCAST, a candidate for next-generation particle-based CSP.

14 SOLAR ENERGY↗
DOE OSTI2024

Quantifying local and global mass balance errors in physics-informed neural networks

Physics-informed neural networks (PINN) have recently become attractive for solving partial differential equations (PDEs) that describe physics laws. By including PDE-based loss functions, physics laws such as mass balance are enforced softly in PINN. This paper investigates how mass balance constraints are satisfied when PINN is used to solve the resulting PDEs. We investigate PINN’s ability to solve the 1D saturated groundwater flow equations (diffusion equations) for homogeneous and heterogeneous media and evaluate the local and global mass balance errors. We compare the obtained PINN’s solution and associated mass balance errors against a two-point finite volume numerical method and the corresponding analytical solution. We also evaluate the accuracy of PINN in solving the 1D saturated groundwater flow equation with and without incorporating hydraulic heads as training data. We demonstrate that PINN’s local and global mass balance errors are significant compared to the finite volume approach. Tuning the PINN’s hyperparameters, such as the number of collocation points, training data, hidden layers, nodes, epochs, and learning rate, did not improve the solution accuracy or the mass balance errors compared to the finite volume solution. Mass balance errors could considerably challenge the utility of PINN in applications where ensuring compliance with physical and mathematical properties is crucial.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Modeling the Enceladus dust plume based on in situ measurements performed with the Cassini Cosmic Dust Analyzer

We analyzed data recorded by the Cosmic Dust Analyzer on board the Cassini spacecraft during Enceladus dust plume traversals. Our focus was on profiles of relative abundances of grains of different compositional types derived from mass spectra recorded with the Dust Analyzer subsystem during the Cassini flybys E5 and E17. The E5 profile, corresponding to a steep and fast traversal of the plume, has already been analyzed. In this paper, we included a second profile from the E17 flyby involving a nearly horizontal traversal of the south polar terrain at a significantly lower velocity. Additionally, we incorporated dust detection rates from the High Rate Detector subsystem during flybys E7 and E21. We derived grain size ranges in the different observational data sets and used these data to constrain parameters for a new dust plume model. This model was constructed using a mathematical description of dust ejection implemented in the software package DUDI. Further constraints included published velocities of gas ejection, positions of gas and dust jets, and the mass production rate of the plume. Our model employs two different types of sources: diffuse sources of dust ejected with a lower velocity and jets with a faster and more colimated emission. From our model, we derived dust mass production rates for different compositional grain types, amounting to at least 28 kg s –1 . Previously, salt-rich dust was believed to dominate the plume mass based on E5 data alone. The E17 profile shows a dominance of organic-enriched grains over the south polar terrain, a region not well constrained by E5 data. By including both E5 and E17 profiles, we find the salt-rich dust contribution to be at most 1% by mass. This revision also results from an improved understanding of grain masses of various compositional types that implies smaller sizes for salt-rich grains. Our new model can predict grain numbers and masses for future mission detectors during plume traversals.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

High-efficiency single-junction p-i-n GaAs solar cell on roll-to-roll epi-ready flexible metal foils for low-cost photovoltaics

We demonstrate the fabrication and characterization of high-efficiency, single-junction p-i-n GaAs solar cells, on flexible metal foil with epi-ready buffer via roll-to-roll fabrication. Single-junction p-i-n GaAs solar cells were fabricated using metal–organic chemical vapor deposition (MOCVD). Furthermore, an efficiency greater than 13% was obtained at 1 sun, which is the highest reported efficiency on GaAs photovoltaics directly deposited on metal tapes. This exceeds our previously reported study showcasing 11.5% efficiency on single-junction p-n solar cell structure. Improved morphology of p-i-n structure compared with p-n is explained by atomic force microscopy (AFM), scanning electron microscopy (SEM), and helium ion microscopy (HIM) measurements. Time-of-flight secondary ion mass spectrometry (TOF-SIMS) analysis showed reduced Zn diffusion in p-i-n cell compared with the p-n cell. We attribute the improvement in efficiency of p-i-n cells on flexible metal tapes to the quality of the junction, surface morphology, controlled diffusion of species within the active layers, and increase in absorption due to the optimized intrinsic layer thickness.

14 SOLAR ENERGY↗
DOE OSTI2024

Geochemical-hydromechanical couplings during water-serpentinized harzburgite interactions at 20° C and 30 bars

Interaction of groundwater with the serpentinized harzburgites of the upper mantle plays an important role in the geologic carbon cycle and serpentinization processes. Here, in this study, an intact serpentinized harzburgite with more than 60% lizardite from the Semail ophiolite in Oman was reacted with water at near neutral pH and low PCO 2 to observe its dissolution behavior and reaction path within the CaO-MgO-SiO 2 -CO 2 -H 2 O system at 20 °C and 30 bars. In the experiment, mass (i.e. H 2 O, ions) transfer is by diffusion and the selected temperature and pressure conditions mimic water-rock interactions at the shallow portions of most ultramafic rock hosted aquifers. Equilibrium activity-activity diagrams and thermodynamic data were used to compare the experimentally determined reaction path with the theoretical stability boundaries for minerals involved to investigate the evolution of bulk chemical composition, mineralogy, and water composition. The experiments demonstrated that (i) the water-rock interactions increased the pH of the aqueous phase from 5.9 to 7.5 over a period of 18 weeks; (ii) the dissolution of lizardite, orthopyroxene (enstatite), and clinopyroxene (diopside and augite) increased the concentration of Mg, Ca, and Si in the aqueous phase; (iii) the dissolution was incongruent with respect to Mg and Si, favoring Si release at pH > 6 due to (a) breaking of more reactive cation‑oxygen bonds that is consistent with the stoichiometry of magnesium for proton exchange reaction which favors a surface that is enriched in Mg at basic conditions, and (b) preferential dissolution of clinopyroxene Mg 1.30 Ca 0.52 Fe 0.06 Si 2 O 6 ; (iv) the aqueous phase was undersaturated with respect to carbonate (e.g. calcite, magnesite, and hydromagnesite) and hydrous (e.g. lizardite, chrysotile, brucite, and talc) minerals; (v) the bulk chemical composition and mineralogy of the intact serpentinized harzburgite matrix did not change; (vi) the thermodynamic data can successfully predict water-intact serpentinized harzburgite behavior if the water chemical composition can be constrained; and (vii) no detectable macroscopic form of damage (e.g. microcracks forming as a result of differential stress fields due to changes in volume during chemical reaction) was observed.

Carbonation↗
DOE OSTI2021

Mechanistic Understanding of Methane Combustion over Ni/CeO 2 : A Combined Experimental and Theoretical Approach

Catalytic oxidation of methane (CH 4 ) over nonprecious Ni/CeO 2 catalysts has received a lot of attention due to the large natural gas reserves found in North America and the prohibitive cost of palladium-based catalysts, commonly used for CH 4 oxidation. However, the catalytic mechanism of CH 4 oxidation over Ni/CeO 2 still remains unclear. Moreover, the parameters affecting the reaction rates, the interaction between nickel and CeO 2 , and the reaction intermediates are still not well understood. In this study, kinetic model fitting, CH 4 temperature-programmed reduction-mass spectroscopy (CH 4 TPR-MS), in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS), and density functional theory (DFT) calculations were combined to elucidate the mechanism of complete oxidation of CH 4 over Ni/CeO 2 . CH 4 TPR-MS showed that the complete oxidation of CH 4 over Ni/CeO 2 requires 55–120 °C lower compared to bare CeO 2 or Ni/quartz sand; complete oxidation of CH 4 took place when the surface oxygen species were abundant, while partial oxidation products (CO, H 2 ) were formed when the oxygen species were depleted. In situ DRIFTS showed that CH 3 , CH 2 , CO, and CO 2 were formed after CH 4 activation over Ni/CeO 2 , while CH 3 O species were not observed. Combining those findings with kinetic model fitting, a redox Mars–van Krevelen (MvK) mechanism showed the best description of the experimental observations. The MvK mechanism involves the reaction of dissociated oxygen species with gas-phase CH 4 while water inhibits the reaction rate by adsorbing on the oxidized sites. Moreover, CH 4 activation leads to the reduction of the active sites and oxygen vacancy formation followed by reoxidation of the active sites by gas-phase O 2 . A CH 4 oxidation reaction pathway over Ni/CeO 2 is proposed by DFT calculations. In summary, the findings shown here suggest that CH 4 oxidation over Ni/CeO 2 follows a redox MvK mechanism and provides guidance for the rational design of non-precious-metal catalysts for CH 4 oxidation reactions.

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

Superionic iron oxide–hydroxide in Earth’s deep mantle

Water ice becomes a superionic phase under the high pressure and temperature conditions of deep planetary interiors of ice planets such as Neptune and Uranus, which affects interior structures and generates magnetic fields. The solid Earth, however, contains only hydrous minerals with a negligible amount of ice. Here we combine high pressure and temperature electrical conductivity experiments, Raman spectroscopy and first-principles simulations to investigate the state of hydrogen in the pyrite-type FeO 2 H x (x ≤ 1), which is a potential H-bearing phase near the core–mantle boundary. We find that when the pressure increases beyond 73 GPa at room temperature, symmetric hydroxyl bonds are softened and the H + (or proton) becomes diffusive within the vicinity of its crystallographic site. Increasing temperature under pressure, the diffusivity of hydrogen is extended beyond the individual unit cell to cover the entire solid, and the electrical conductivity soars, indicating a transition to the superionic state, which is characterized by freely moving protons and a solid FeO 2 lattice. Furthermore, the highly diffusive hydrogen provides fresh transport mechanisms for charge and mass, which dictate the geophysical behaviours of electrical conductivity and magnetism, as well as geochemical processes of redox, hydrogen circulation and hydrogen isotopic mixing in Earth’s deep mantle.

58 GEOSCIENCES↗
DOE OSTI2024

Biopolymer-supramolecular polymer hybrids for photocatalytic hydrogen production

Solar generation of H 2 is a promising strategy for dense energy storage. Supramolecular polymers composed of chromophore amphiphile monomers containing perylene monoimide (PMI) have been reported as crystalline light-harvesting assemblies for aqueous H 2 -evolving catalysts. Gelation of these supramolecular polymers with multivalent ions creates hydrogels with high diffusivity but insufficient mechanical stability and catalyst retention for reusability. We report here on using sodium alginate (SA) biopolymer to both induce supramolecular polymerization of PMI and co-immobilize them with catalysts in a robust hydrogel with high diffusivity that can also be 3D-printed. Faster mass transfer was achieved by controlling the material macrostructure by reducing gel diameter and microstructure by reducing biopolymer loading. Optimized gels produce H 2 at rates rivaling solution-based PMI and generate H 2 for up to 6 days. The PMI assemblies in the SA matrix create a percolation network capable of bulk-electron transfer under illumination. These PMI–SA materials were then 3D-printed on conductive substrates to create 3D hydrogel photoelectrodes with optimized porosity. The design of these versatile hybrid materials was bioinspired by the soft matter environment of natural photosynthetic systems and opens the opportunity to carry out light-to-fuel conversion within soft matter with arbitrary shapes and particular local environments.

36 MATERIALS SCIENCE↗