Search NASA⌕ Search

SEARCH · Search NASA

Results for “chemical reaction networks”

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.

611 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Growth of deuterium supersaturated surface layer with increasing ion flux and fluence in plasma-exposed tungsten

Deuterium supersaturated surface layer (DSSL) in tungsten, a few nm thick layer exhibiting extremely high deuterium content (> 5%), has been studied as a function of deuterium plasma ion flux and fluence. Tungsten samples were exposed at 400 K and deuterium ion energy of ∼ 60 eV. The deuterium ion flux spanned over an order of magnitude, being 7.3 × 10 20 , 4.2 × 10 21 , and 3.8 × 10 22 D/m 2 . The fluence ranged from 3.4 × 1024 to 3.4 × 1025 D/m2. The samples and their deuterium content were analyzed by nuclear reaction analysis (NRA), scanning transmission electron microscopy (STEM), and thermal desorption spectroscopy (TDS). The largest thickness of DSSL was found to be around 10.5 nm and was observed in the case of the highest exposure flux and fluence, whereas the layer was only about 3.8 nm thick in the case of the lowest value of the two parameters. The thickness of the DSSL was found to monotonically increase with both deuterium flux and fluence. Finally, similar to the thickness, the estimated D concentration seems to follow the same trend, increasing from the lowest value of around 3 at.% to the highest value of around 5 at.%.

Deuterium↗
DOE OSTI2024

Simulating hindered grain boundary diffusion using the smoothed boundary method

Abstract Grain boundaries can greatly affect the transport properties of polycrystalline materials, particularly when the grain size approaches the nanoscale. While grain boundaries often enhance diffusion by providing a fast pathway for chemical transport, some material systems, such as those of solid oxide fuel cells and battery cathode particles, exhibit the opposite behavior, where grain boundaries act to hinder diffusion. To facilitate the study of systems with hindered grain boundary diffusion, we propose a model that utilizes the smoothed boundary method to simulate the dynamic concentration evolution in polycrystalline systems. The model employs domain parameters with diffuse interfaces to describe the grains, thereby enabling solutions with explicit consideration of their complex geometries. The intrinsic error arising from the diffuse interface approach employed in our proposed model is explored by comparing the results against a sharp interface model for a variety of parameter sets. Finally, two case studies are considered to demonstrate potential applications of the model. First, a nanocrystalline yttria-stabilized zirconia solid oxide fuel cell system is investigated, and the effective diffusivities are extracted from the simulation results and are compared to the values obtained through mean-field approximations. Second, the concentration evolution during lithiation of a polycrystalline battery cathode particle is simulated to demonstrate the method’s capability.

Materials Science↗
DOE OSTI2025

Blockage effects in the chemotaxis of diffusiophoretic particles

Transport mechanisms at the micro- and nano-scale play an essential role in regulating intracellular organization. Recent work indicates that directed motion of constituents inside cells can emerge through diffusiophoretic transport, in which colloidal particles move under the influence of chemical gradients. Here, we examine how blockers—passive or actively consuming—reshape those gradients and thereby influence the motion of diffusiophoretic particles. By combining analytical solutions with finite element simulations, we first show that a single blocker can distort a background gradient enough to create or eliminate stagnation points, significantly modifying particle transport. We then introduce a second, explicitly sized blocker at one of these stagnation points and measure how its finite radius alters the diffusiophoretic velocity field for a test particle. Even moderate changes in the second blockers size can cause noticeable shifts in the substrate distribution, highlighting the importance of accounting for explicit particle radii under crowded or consumption-driven conditions. Our findings underscore that subtle geometric variations—such as the radii and positions of two or more blockers—can profoundly affect diffusiophoretic motion, providing a more complete picture of how blocking and crowding phenomena shape intracellular transport.

Song, Zehao [Northwestern Univ., Evanston, IL (Uni↗
DOE OSTI2026

Cave-based chronology of valley incision in the glaciated Northern Carpathians

While glaciers leave a mark on mountain landscapes, the importance of glacial erosion in mountains at a global scale remains insufficiently constrained despite over a century of research aimed at quantifying the rates and magnitudes of such erosion. Measuring rates of landscape change across the Quaternary in repeatedly glaciated mountain landscapes requires the ability to accurately partition the relative role of glacial versus fluvial processes in forming relief; however, the modification of landforms by surface erosion hinders quantification of process rates over the last few million years. This study presents 32 paired 26 Al and 10 Be cosmogenic nuclide burial ages from allochthonous, quartz-bearing sediment deposited within a paleophreatic karst system in the northern Carpathians. These data quantify Pliocene-to-present timing of valley incision and karst system reactivation, as well as erosion rates over time. Post-burial, muon production-correction yields 24 fully resolved ages, six minimum ages, and two unresolvable ages. Burial ages are interpreted to document the tempo of valley incision via phreatic passage formation and reveal that ∼80% of current valley relief (∼700 m) had developed by approximately 1.6 Ma. Paleoerosion rates derived from this dataset suggest little change in the processes generating sediment washed into the karst networks since the Pliocene. These results suggest that river incision dominated relief creation in Carpathian valleys that bear the imprint of glaciation. Our findings suggest the importance of investigating other mid-latitude mountainous regions worldwide to determine if they have similar erosion and topographic histories.

Carpathians↗
DOE OSTI2025

CASM Monte Carlo: Calculations of the thermodynamic and kinetic properties of complex multicomponent crystals

Monte Carlo techniques play a central role in statistical mechanics approaches that connect macroscopic thermodynamic and kinetic properties to the electronic structure of a material. This paper describes the implementation of Monte Carlo techniques for the study of multicomponent crystalline materials within the Clusters Approach to Statistical Mechanics (CASM) software suite, and demonstrates their use in model systems to calculate free energies and kinetic coefficients, study phase transitions, and construct phase diagrams from first principles. Many crystal structures are complex, with multiple sublattices occupied by differing sets of chemical species, along with the presence of vacancies or interstitial species. This imposes constraints on concentration variables, the form of thermodynamic potentials, and the values of kinetic transport coefficients. The framework used by CASM to formulate thermodynamic potentials and kinetic transport coefficients accounting for arbitrarily complex crystal structures is presented and demonstrated with examples of increasing complexity. Additionally, an overview of the capabilities of the CASM software specific to Monte Carlo methods is given, and a new CASM software package is introduced, casm-flow, which helps automate the setup, submission, management, and analysis of Monte Carlo simulations.

Cluster expansion↗
DOE OSTI2026

Contrasting structural reversibility and magnetic correlations in isostructural honeycomb magnets CrCl3 and 𝛼−RuCl3

We report a comparative neutron single crystal diffraction study of the structural and magnetic properties of layered halides CrCl3 and 𝛼−RuCl3. They host a honeycomb arrangement of transition metal ions with distinct electronic configurations and undergo a first-order structural transition between high-temperature 𝐶⁢2/𝑚 and low-temperature 𝑅⁢‾‾‾3. Both compounds show a step-like change in the 𝑐-lattice parameter across the structure transition. In contrast, the in-plane lattice response is quite different: 𝛼−RuCl3 exhibits an abrupt hysteretic change across the transition accompanied by progressive crystalline degradation upon thermal cycling, whereas CrCl3 shows a smooth in-plane lattice evolution and remains structurally robust. Magnetically, CrCl3 orders into an A-type antiferromagnetic structure at 𝑇𝑁=14 K and exhibits pronounced diffuse magnetic scattering extending up to about 40 K. 𝛼−RuCl3 shows no observable magnetic diffuse scattering above its zigzag antiferromagnetic ordering temperature 𝑇𝑁=7.6 K. These results suggest that the contrasting structural responses arise from an interplay between interlayer sliding energetics, stacking-strain coupling, and elastic accommodation of the stacking transition. The distinct chemical bonding and electronic configurations of the two compounds provide a microscopic basis for their different lattice responses to the structure transition and magnetic correlations.

Morgan, Zachary [ORNL] (ORCID:0000000243625911)↗
NASA NTRSDate not supplied

Discrete Rotational Energy for Polyatomic Molecules in Direct Simulation Monte Carlo

The environments experienced by vehicles upon entry into planetary atmospheres generally result in severe aerothermodynamic loading. These flows are characterized by the formation of strong shock waves, behind which high temperatures and non-equilibrium environments are generated, and the accurate prediction of internal energy relaxation and chemical kinetics becomes important. This manuscript focuses on the internal relaxation of the rotational energy of polyatomic molecules. Previous work has described the implementation of a discrete rotational energy model for diatomic molecules (linear rotators). Later, Gimelshein et. al extended this model to include polyatomic molecular internal energies. However, they did not include details on how to generate the particles at equilibrium conditions and centrifugal distortion was not included in their model. The purpose of the present study is to extend the afore-mentioned models to the treatment of discrete rotational energy in polyatomic molecules with the addition of centrifugal distortion. An overview of the background theory needed for the implementation of the model will be given followed by the model itself and verification of the model. For example, sampled versus Boltzmann distributions are compared in Figure 1 for CH3 (oblate rotator – two quantum numbers) and examples of thermal relaxation behavior for linear and non-linear rotators are presented in Figure 2. In addition, comparisons of thermophysical properties will be presented, as well as examples of adiabatic relaxation and application to hypersonic flow.

Rarefied Gas Dynamics↗
DOE OSTI2026

Metal–bicarbonate ion pairing in alkaline aqueous solutions from multilevel embedded correlated wavefunction theory and molecular dynamics

In this work, we examine ion-pairing mechanisms of bicarbonates in alkaline aqueous solutions with the divalent metal ions most abundantly present in seawater, namely, Ca 2+ and Mg 2+ . We employ a rare-event enhanced sampling approach within first-principles molecular dynamics to explore regions of phase space spanning solvent-shared to contact ion pairs. Second-order Møller–Plesset perturbation theory (MP2) corrections are subsequently applied in an embedding framework (EMB) to refine the electronic structure of stationary states and associated reaction barriers along the free-energy profiles while retaining the extended solvent effects at the density functional theory (DFT) level. Ca 2+ –HCO 3 − was previously hypothesized to exist in a solvent-shared ion pair (SSHIP) by DFT studies with an endergonic contact ion pair (CIP) formation; however, our EMB-MP2 refinement of the DFT ion-pairing pathways reveals that Ca 2+ and HCO 3 − form a virtually barrier-free CIP in alkaline solutions, with even more energetic ease than the widely studied Ca–CO 3 ion pair. We find qualitative agreement between DFT and EMB-MP2 for Mg 2+ —unlike Ca 2+ , Mg 2+ refuses to shed its strong hydration shell, thereby preferring a SSHIP state with a significant activation barrier to crossover to the CIP forms—a trait reminiscent of ion pairing in Mg–CO 3 and closely related to the kinetic limitations underlying the famous subject of the dolomite problem. Our study highlights the importance of improved electronic structure descriptions of liquids, modeled as a condensed phase of matter lacking in long-range crystalline order. It also strongly suggests that Ca 2+ –HCO 3 − CIPs are likely precursors involved in prenucleation of CaCO 3 mineral formation in seawater.

Sharma, Vidushi [Princeton Plasma Physics Laborato↗
DOE OSTI2024

Mechanochemically responsive polymer enables shockwave visualization

Abstract Understanding the physical and chemical response of materials to impulsive deformation is crucial for applications ranging from soft robotic locomotion to space exploration to seismology. However, investigating material properties at extreme strain rates remains challenging due to temporal and spatial resolution limitations. Combining high-strain-rate testing with mechanochemistry encodes the molecular-level deformation within the material itself, thus enabling the direct quantification of the material response. Here, we demonstrate a mechanophore-functionalized block copolymer that self-reports energy dissipation mechanisms, such as bond rupture and acoustic wave dissipation, in response to high-strain-rate impacts. A microprojectile accelerated towards the polymer permanently deforms the material at a shallow depth. At intersonic velocities, the polymer reports significant subsurface energy absorption due to shockwave attenuation, a mechanism traditionally considered negligible compared to plasticity and not well explored in polymers. The acoustic wave velocity of the material is directly recovered from the mechanochemically-activated subsurface volume recorded in the material, which is validated by simulations, theory, and acoustic measurements. This integration of mechanochemistry with microballistic testing enables characterization of high-strain-rate mechanical properties and elucidates important insights applicable to nanomaterials, particle-reinforced composites, and biocompatible polymers.

Science & Technology - Other Topics↗
DOE OSTI2026

Observation of N-rich solid-electrolyte interphase by ToF-SIMS.

Formation of a stable solid electrolyte interphase (SEI) between lithium electrodes and electrolyte upon multiple charge/discharge cycles is crucial to a long-term lithium-ion battery performance. Addition of LiNO3 to lithium bis (fluorosulfonyl) imide/poly(ethylene oxide) (LiFSI/PEO) electrolyte leads to a durable SEI that is electrically insulating yet highly conductive to Li ions, chemically and electrochemically stable, physically uniform, and mechanically robust. ToF-SIMS was used here in combination with sputtering by a gaseous cluster ion beam (GCIB) to examine how the addition of a small proportion of LiNO3 to the LiFSI/PEO electrolyte affects the SEI composition. Negative ion ToF-SIMS spectra of the cycled samples display an intense m/z 26 peak associated with the SEI. Exact mass assignments and isotopic ratios indicate that this peak should be assigned as (CN-)-C-12, with little to no negative secondary ion signal arising from (LiF-)-Li-7. This CN- signal appears to arise from an N-rich portion of the SEI adjacent to the Li electrode that is depleted in LiF relative to the bulk electrolyte. The dearth of LiF- (and LiF+ from the positive ion spectra) is unexpected because LiF has been identified in the SEI in similar samples. Finally, GCIB sputtering indicates that the SEI adheres more strongly to the Li electrode than to the LiFSI/PEO electrolyte.

Shavandi, Seyedeh Reyhaneh↗
DOE OSTI2026

Treatment of Lagoon Dairy Manure Wastewater via Iron Electrocoagulation, Microfiltration, and Adsorption

Dairy manure wastewater generated by flushing barn cow waste contains nutrients, pathogens, and organic and inorganic contaminants. This study utilized a process consisting of iron electrocoagulation (Fe-EC), microfiltration (MF), and activated carbon (AC) adsorption to treat farm wastewater and explore the reclamation of clean water for irrigation and livestock consumption. Significant removal (>99.9%) of chemical oxygen demand (COD), total organic carbon (TOC), phosphorus (P), turbidity, and microorganisms, as well as ions including magnesium, calcium, sulfur, and silica was achieved by the combined EC-MF-AC process. Specifically, a charge loading of ∼37,500 C/L in a continuous-flow EC configuration, followed by MF, achieved more than 95% removal of TOC and COD. Characterization of produced flocs and foam via scanning-electron microscopy with energy-dispersal spectroscopy and Fourier transform infrared spectroscopy confirmed the removal of ions, including calcium, sulfur, and silica. A key finding was the electrocatalytic conversion of nitrogen species to ammonia gas through the intermediate reduction of nitrate/nitrite, which led to ∼60% total nitrogen (TN) removal. AC treatment further improved TN removal to ∼70%. The Fe-EC process also eradicated >99.9% of bacteria. Preliminary process cost assessment, based on recycled materials for EC electrodes, showed significant cost savings (∼2 times) compared to commercial electrodes.

Dutta, Swapnamoy [ORNL]↗
DOE OSTI2023

Fabrication of Large-Area Metal-on-Carbon Catalytic Condensers for Programmable Catalysis

Catalytic condensers stabilize charge on either side of a high-k dielectric film to modulate the electronic states of a catalytic layer for the electronic control of surface reactions. Here, carbon sputtering provided for fast, large-scale fabrication of metal–carbon catalytic condensers required for industrial application. Carbon films were sputtered on HfO 2 dielectric/p-type Si with different thicknesses (1, 3, 6, and 10 nm), and the enhancement of conductance and capacitance of carbon films was observed upon increasing the carbon thickness following thermal treatment at 400 °C. After Pt deposition on the carbon films, the Pt catalytic condenser exhibited a high capacitance of ∼210 nF/cm 2 that was maintained at a frequency ∼1000 Hz, satisfying the requirement for a dynamic catalyst to implement catalytic resonance. Temperature-programmed desorption of carbon monoxide yielded CO desorption peaks that shifted in temperature with the varying potential applied to the condenser (−6 or +6 V), indicating a shift in the binding energy of carbon monoxide on the Pt condenser surface. A substantial increase in capacitance (∼2000 nF/cm 2 ) of the Pt-on-carbon devices was observed at elevated temperatures of 400 °C that can modulate ∼10% of charge per metal atom when 10 V potential was applied. A large catalytic condenser of 42 cm 2 area Pt/C/HfO 2 /Si exhibited a high capacitance of 9393 nF with a low leakage current/capacitive current ratio (<0.1), demonstrating the practicality and versatility of the facile, large-scale fabrication method for metal–carbon catalytic condensers.

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

Novel PtNi single-atom–nanocluster (SA–NC) ensembles promote Tafel kinetics and ampere-class AEM hydrogen evolution

The development of efficient, durable, and low-PGM electrocatalysts for the hydrogen evolution reaction (HER) in alkaline media is critical for next-generation electrolysis technologies. We report a facile two-step synthesis of highly dispersed PtNi and PtNi-nitride nanoclusters (NCs) (~2.3 nm) with ultralow Pt content (0.5 at.%) anchored on N-doped Vulcan carbon. Structural and compositional characterization via XAS, XPS, HAADF-STEM, HRTEM, and EDS mapping established key structure–activity relationships across varying Pt/Ni ratios and pyrolysis temperatures. The Pt0.5Ni0.5/C-750 catalyst, an ensemble of PtNi M-N-C type single-atom (SA) moieties with neighboring PtNi nanoclusters (NC), exhibited superior HER performance in alkaline media, achieving overpotentials of 30, 115, and 210 mV at 10, 100, and 500 mA cm−2, respectively. Despite at a lower Pt content, this novel SA–NC ensemble outperformed commercial Pt/C by ~36%. A standardized literature comparison with contemporary Pt- and Ru-doped analogues reveals the as-prepared Pt0.5Ni0.5/C-750 to sit at the apex of Tafel-limited kinetics and low overpotential at 100 mA cm−2. Tafel-limited Tafel slopes in both alkaline and acidic regimes confirm favorable proton recombination kinetics. Mass activities at 200 mV reached 13.8 and 18.84 A mgPt−1 in alkaline and acidic media, respectively. However, excessive nitridation (e.g., at 650 °C) adversely altered Pt electronic structure and HER kinetics. While Ni enhanced alkaline HER, acidic HER favored Ni-free analogues. Pt0.5Ni0.5/C-750 also demonstrated robust temperature responsiveness and 300-h operational stability at high current densities (0.5–1.0 A cm−2) in MEA tests. This work presents a scalable strategy for designing thermally responsive, durable, and compositionally tunable NC catalysts with neighboring SA moieties for alkaline electrolysis.

AEMWE↗
DOE OSTI2026

Accelerating actinium-225 purification by high-pressure ion chromatography

Actinium-225 (t1/2 = 9.92 days) is an important radioisotope for targeted alpha therapy applications. The limited supply obtained through the decay of thorium-229 has motivated accelerator-based production routes, including irradiation of thorium targets. Irradiated targets can produce useful quantities of actinium-225, but the product requires final purification from chemically similar lanthanide contaminants. This work describes an automated high-pressure ion chromatography method for this final polishing step. The method uses a reusable strong-acid cation-exchange column bearing sulfonic acid functional groups. α-Hydroxyisobutyric acid (α-HIBA), adjusted to pH 4.3 with lithium hydroxide, complexes and elutes lanthanides, a dilute hydrochloric acid matrix-exchange step removes residual α-HIBA, and concentrated hydrochloric acid then elutes retained actinium(III). The protocol purified actinium-225 to >99% radiopurity across tracer-level samples and samples containing >150 µCi (5.6 MBq) of activity. A 10 min, 0.1 M hydrochloric acid matrix exchange substantially reduced organic eluent carryover, and in-line sodium iodide detection enabled real-time monitoring of actinium and lanthanide elution. The developed method can be completed in <1 h and provides a basis for automated purification workflows for accelerator-produced actinium-225.

Gaddis, Kevin [ORNL] (ORCID:0000000183398314)↗
NASA NTRSDate not supplied

NASA's RPS Program Status Update and Future Outlook

Radioisotope Power Systems (RPS) have been safely providing the United States with the power to explore space for almost 65 years. NASA established the RPS Program in 2010 to support specialized science missions and potential future opportunities more effectively and efficiently. The RPS Program ensures the availability of RPS for the exploration of the solar system in environments where conventional solar or chemical power generation is impractical or impossible. NASA missions have utilized space nuclear power to explore planets, moons, and interstellar space, enabling missions to some of the deepest, darkest, and dustiest regions in the solar system and beyond. This scientific exploration has deepened our understanding of the solar system and our role within. The RPS Program, in partnership with the Department of Energy (DOE) Office of Nuclear Energy continues to operate as an interagency partnership to provide robust radioisotope power system solutions to spacecraft that conduct missions to extreme environments for science and exploration. This Program invests in systems and technologies to ensure that NASA maintains this capability well into the future and it manages processes that ensure the safe use and launch of these systems. This paper provides a synopsis of current activities of the RPS Program, and it introduces a focus on the commercialization of RPS to enable a variety of future mission applications.

radioisotope power systems↗
DOE OSTI2025

Interleaved bond frustration in a triangular lattice antiferromagnet

Frustration of long-range order via lattice geometry amplifies fluctuations and generates ground states that are highly sensitive to perturbations. Traditionally, geometric frustration is used to engineer unconventional magnetic states; however, the charge degree of freedom and bond order can be similarly frustrated. Finding materials that host both frustrated magnetic and bond networks holds promise for engineering structural and magnetic states with the potential of coupling to one another via either magnetic or strain fields. Here we identify an unusual instance of this coexistence in the triangular lattice antiferromagnets LnCd 3 P 3 (Ln = lanthanides). These compounds feature two-dimensional planes of unique trigonal planar CdP 3 units with an underlying bond instability that is frustrated via emergent kagome ice correlations. This bond instability is interleaved in between layers of frustrated magnetic moments. Furthermore, our results establish LnCd 3 P 3 as a rare materials class in which frustrated magnetism is embedded within a dopable semiconductor with a frustrated bond order instability.

Gomez Alvarado, S. J. [University of California, S↗
DOE OSTI2025

Large language model-driven database for thermoelectric materials

Thermoelectric materials have the ability to convert waste heat into electricity, offering a valuable solution for energy harvesting. However, their widespread use is hindered by low conversion efficiency, the reliance on expensive rare earth elements, and the environmental and regulatory concerns associated with lead-based materials. A fast and cost-effective way to identify highly efficient thermoelectric materials is through data-driven methods. These approaches rely on robust and comprehensive datasets to train models. Although there are several databases on thermoelectric materials, there is still a need to collect and integrate experimental data from peer-reviewed research articles to capture diverse compositions and properties of materials. Here, in this work, we developed a comprehensive database of 7,123 thermoelectric compounds, containing key information such as chemical composition, structural detail, seebeck coefficient, electrical and thermal conductivity, power factor, and figure of merit (ZT). We used the GPTArticleExtractor workflow, powered by large language models (LLM), to extract and curate data automatically from the scientific literature published in Elsevier journals. This process enabled the creation of a structured database that addresses the challenges of manual data collection. The open access database could stimulate data-driven research and advance thermoelectric material analysis and discovery.

Database↗
DOE OSTI2025

High-voltage water-scarce hydrogel electrolytes enable mechanically safe stretchable Li-ion batteries

Soft Li-ion batteries, based on conventional organic electrolytes, face performance degradation challenges due to moisture penetration and safety concerns due to possible leakage of toxic fluorine compounds and flammable solvents under mechanical damage. We design a water-scarce hydrogel electrolyte with fluorine-free lithium salt to achieve wide electrochemical stability window (up to 3.11 volts) in ambient air without hermetic packaging while balancing high stretchability (1348%), ion conductivity (41 millisiemens per centimeter), and self-healing capabilities for mechanically and chemically safe stretchable Li-ion batteries. Molecular synergy between hydrophilicity and lithiophilicity of zwitterionic polymer backbone is revealed by molecular dynamics simulations. The battery exhibits capacity retention under harsh mechanical stresses—enduring stretching, twisting, folding, and multiple through-punctures by a needle—while self-healing from repeated through cuts by a razor blade. Stable ambient operation for 1 month over 500 charge-discharge cycles (average coulomb efficiency, 95%) is achieved. A prototype self-healing electronic system with embedded soft batteries demonstrates practical application as a durable embodied energy source.

Science & Technology - Other Topics↗