Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Extreme scale”

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

Is the Local Ion Density Sufficient to Drive NaCl Nucleation from the Melt and Aqueous Solution?

Even though nucleation is ubiquitous in different science and engineering problems, investigating nucleation is extremely difficult due to the complicated ranges of time and length scales involved. In this work, we simulate NaCl nucleation in both molten and aqueous environments using enhanced sampling of all-atom molecular dynamics with deep-learning-based estimation of reaction coordinates. By incorporating various structural order parameters and learning the reaction coordinate as a function thereof, we achieve significantly improved sampling relative to traditional ad hoc descriptions of what drives nucleation, particularly in an aqueous medium. Our results reveal a one-step nucleation mechanism in both environments, with reaction coordinate analysis highlighting the importance of local ion density in distinguishing solid and liquid states. However, although fluctuations in the local ion density are necessary to drive nucleation, they are not sufficient. Our analysis shows that near the transition states, descriptors such as enthalpy and local structure become crucial. Furthermore, our protocol proposed here enables robust nucleation analysis and phase sampling and could offer insights into nucleation mechanisms for generic small molecules in different environments.

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

Extreme ultraviolet microscope characterization using photomask surface roughness

We demonstrate a method for characterizing the field-dependent aberrations of a full-field synchrotron-based extreme ultraviolet microscope. The statistical uniformity of the inherent, atomic-scale roughness of readily-available photomask blanks enables a self-calibrating computational procedure using images acquired under standard operation. We characterize the aberrations across a 30-um field-of-view, demonstrating a minimum aberration magnitude of smaller than λ/21 rms averaged over the center 5-um area, with a measurement accuracy better than λ/180 rms. The measured field variation of aberrations is consistent with system geometry and agrees with prior characterizations of the same system. In certain cases, it may be possible to additionally recover the illumination wavefront from the same images. Our method is general and is easily applied to coherent imaging systems with steerable illumination without requiring invasive hardware or custom test objects; hence, it provides substantial benefits when characterizing microscopes and high-resolution imaging systems in situ.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Efficient numerical method for evaluating normal and anomalous time-domain equilibrium Green's functions in inhomogeneous systems

In this work we develop EPOCH (equilibrium propagator by orthogonal polynomial chain), a computationally efficient method to calculate the time-dependent equilibrium Green's functions, including the anomalous Green's functions of superconductors, to capture the time evolution in large inhomogeneous systems. The EPOCH method generalizes the Chebyshev wave-packet propagation method from quantum chemistry and efficiently incorporates the Fermi-Dirac statistics that is needed for equilibrium quantum condensed matter systems. The computational cost of EPOCH scales only linearly in the system degrees of freedom, generating an extremely efficient algorithm also for very large systems. We demonstrate the power of the EPOCH method by calculating the time evolution of an excitation near a superconductor–normal metal interface in two and three dimensions, capturing transmission as well as normal and Andreev reflections.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Uniqueness criteria for the Virasoro-Shapiro amplitude

The scattering amplitudes of string theory exhibit many extraordinary properties. But are they the unique mathematical objects to do so? Recently, it has been shown how the spectrum and amplitudes of open string theory follow directly from the assumptions of faster than power-law falloff at high energies and a property dubbed level truncation. At present there is no analogous principle for closed string scattering, which is famously rigid and naively impervious to modification. In this paper we analytically bootstrap the spectrum and four-point amplitudes of the closed string—together with a parametrized space of deformations—from conditions on high-energy falloff and level truncation. While these deformations exhibit the same Regge scaling as pure gravity, in the tensionless limit they reproduce remarkable extremal amplitudes that have appeared in bottom-up studies of positivity.

Field & string theory models & techniques↗
DOE OSTI2025

The Origins of the ParaView and VisIt Scientific Visualization Tools

ParaView and VisIt play a key role in the visual understanding of scientific simulation data. These tools are open source, designed to handle extremely large datasets, and can run on supercomputers and large-scale display walls. They are in daily use by scientists, practitioners, and students at supercomputing centers, in industry, and at universities. In conclusion, this article gives personal accounts of the origins of these visualization tools.

Ahrens, James [Los Alamos National Laboratory (LAN↗
DOE OSTI2023

Resilience Analytics for Interdependent Power and Water Distribution Systems

The occurrence of large-scale outages in power distribution systems (PDSs) caused, for instance, by extreme weather events has raised concerns on the cascading impact of the power outages on the resilience of other critical infrastructures, including water distribution systems (WDSs), whose healthy operation highly relies on the availability of power from PDSs. Here, this paper proposes an analytics model for quantifying the interdependence between the resilience of power and water distribution systems. The proposed model first performs spatio-temporal outage analysis on the PDS and WDS to determine the PDS power serving capability and the resulting cascading impact on the WDS water serving capability during the occurrence of disruptions. The resilience interdependence is then determined by a proposed set of six metrics calculated using the power and water serving capability profiles, which trace and quantify the cascading impacts of power outages in the operation of WDSs during the degradation and recovery stages. The proposed analytics model is implemented on the IEEE 33-bus power distribution system supplying power for pumping stations of a 15-node test WDS, considering multiple cases of water storage capacity, single and multiple outage scenarios, and availability of distributed generation. The numerical results demonstrate that the proposed metrics track the temporal and spatial interdependence between the resilience of power and water distribution systems, providing a mechanism to identify weak spots, and assess options to enhance the resilience of the critical water infrastructure.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Soil microbiome interventions for carbon sequestration and climate mitigation

Mitigating climate change in soil ecosystems involves complex plant and microbial processes regulating carbon pools and flows. Here, we advocate for the use of soil microbiome interventions to help increase soil carbon stocks and curb greenhouse gas emissions from managed soils. Direct interventions include the introduction of microbial strains, consortia, phage, and soil transplants, whereas indirect interventions include managing soil conditions or additives to modulate community composition or its activities. Approaches to increase soil carbon stocks using microbially catalyzed processes include increasing carbon inputs from plants, promoting soil organic matter (SOM) formation, and reducing SOM turnover and production of diverse greenhouse gases. Marginal or degraded soils may provide the greatest opportunities for enhancing global soil carbon stocks. Among the many knowledge gaps in this field, crucial gaps include the processes influencing the transformation of plant-derived soil carbon inputs into SOM and the identity of the microbes and microbial activities impacting this transformation. As a critical step forward, we encourage broadening the current widespread screening of potentially beneficial soil microorganisms to encompass functions relevant to stimulating soil carbon stocks. Moreover, in developing these interventions, we must consider the potential ecological ramifications and uncertainties, such as incurred by the widespread introduction of homogenous inoculants and consortia, and the need for site-specificity given the extreme variation among soil habitats. Incentivization and implementation at large spatial scales could effectively harness increases in soil carbon stocks, helping to mitigate the impacts of climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Observation of time-dependent internal charge amplification in a planar germanium detector at cryogenic temperature

For the first time, time-dependent internal charge amplification through impact ionization has been observed in a planar germanium (Ge) detector operated at cryogenic temperature. In a time period of 30 and 45 min after applying a bias voltage, the charge energy corresponding to a baseline of the 59.54 keVγ rays from a 241 Am source is amplified for a short period of time and then decreases back to the baseline. The amplification of charge energy depends strongly on the applied positive bias voltage with drifting holes across the detector. No such phenomenon is visible with drifting electrons across the detector. We find that the observed charge amplification is dictated by the impact ionization of charged states, which has a strong correlation with impurity level and applied electric field. We analyze the dominant physics mechanisms that are responsible for the creation and the impact ionization of charged states. Our analysis suggests that the appropriate level of impurity in a Ge detector can enhance charge yield through the impact ionization of charged states to achieve extremely low-energy detection threshold (< 10 meV) for MeV-scale dark matter searches if the charge amplification can be stabilized.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Machine Learning for Adaptive Model Refinement to Bridge Scales

This whitepaper is responsive to focal area (2) Predictive modeling through the use of AI techniques and AI-derived model components: the use of AI and other tools to design a prediction system comprising of a hierarchy of models (e.g., AI driven model/component/parameterization selection). Here we describe scale-aware ML models for adaptive model refinement that allow us to bridge the spatial and/or temporal scales in simulation models and observation data for capturing and predicting extreme water cycles.

58 GEOSCIENCES↗
DOE OSTI2024

Low-Temperature Annealing of Nanoscale Defects in Polycrystalline Graphite

Polycrystalline graphite contains multi-scale defects, which are difficult to anneal thermally because of the extremely high temperatures involved in the manufacturing process. In this study, we demonstrate annealing of nuclear graphite NBG-18 at temperatures below 28 °C, exploiting the electron wind force, a non-thermal stimulus. High current density pulses were passed through the specimens with a very low-duty cycle so that the electron momentum could mobilize the defects without heating the specimen. The effectiveness of this technique is presented with a significant decrease in electrical resistivity, defect counts from X-ray computed tomography, Raman spectroscopy, and nanoindentation-based mechanical characterization. Such multi-modal evidence highlights the feasibility of nanoscale defect control at temperatures about two orders of magnitude below the graphitization temperature.

Liu, Gongyuan↗
DOE OSTI2025

Extreme Heating of Minor Ions in Imbalanced Solar-wind Turbulence

Minor ions in the solar corona are heated to extreme temperatures, far in excess of those of the electrons and protons that comprise the bulk of the plasma. These highly nonthermal distributions make minor ions sensitive probes of the collisionless processes that heat the corona and power the solar wind. The recent discovery of the "helicity barrier" offers a mechanism in which imbalanced Alfvénic turbulence in low-β plasmas preferentially heats protons over electrons, generating high-frequency, proton-cyclotron-resonant fluctuations. We use the hybrid-kinetic particle-in-cell code Pegasus++ to drive imbalanced Alfvénic turbulence in a 3D low-β plasma with additional passive ion species, He 2+ and O 5+ . A helicity barrier naturally develops, followed by clear phase-space signatures of oblique proton-cyclotron-wave heating and Landau-resonant heating from the imbalanced Alfvénic fluctuations. The former results in characteristically arced ion velocity distribution functions, whose non-bi-Maxwellian features are shown by linear ALPS calculations to be critical to the heating process. Additional features include a steep transition-range electromagnetic spectrum, proton-cyclotron waves propagating in the direction of the imbalance, significantly enhanced proton-to-electron heating ratios, ion temperatures that are considerably more perpendicular with respect to magnetic field, and extreme heating of heavier species in a manner consistent with mass scalings inferred from spacecraft measurements. None of these features are realized in an otherwise equivalent simulation of balanced turbulence. If seen simultaneously in the fast solar wind, these signatures of the helicity barrier would testify to the necessity of incorporating turbulence imbalance in a complete theory for the evolution of the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Joint X-Ray, Kinetic Sunyaev–Zeldovich, and Weak Lensing Measurements: Toward a Consensus Picture of Efficient Gas Expulsion from Groups and Clusters

There is no consensus on how baryon feedback shapes the underlying matter distribution from either simulations or observations. We confront the uncertain landscape by jointly analyzing new measurements of the gas distribution around groups and clusters—DESI+ACT kinetic Sunyaev–Zel’dovich (kSZ) effect profiles and eROSITA X-ray gas masses—with mean halo masses characterized by galaxy–galaxy lensing. Across a wide range of halo masses ( M 500 = 10 13−14 M ⊙ ) and redshifts (0 < z < 1), we find evidence of more efficient gas expulsion beyond several R 500 than predicted by most state-of-the-art simulations. A like-with-like comparison reveals all kSZ and X-ray observations are inconsistent with the fiducial 1 Gpc 3 hydrodynamical FLAMINGO simulation, which was calibrated to reproduce pre-eROSITA X-ray gas fractions: eROSITA X-ray gas fractions are 2 × lower than the simulation, and the kSZ measurements are combined >8σ discrepant. The FLAMINGO simulation variant with the most gas expulsion, and therefore the most suppression of the matter power spectrum relative to a dark-matter-only simulation, provides a good description of how much gas is expelled and how far it extends; the enhanced gas depletion is achieved by more powerful but less frequent AGN outbursts. Joint kSZ, X-ray, and lensing measurements form a consistent picture of gas expulsion beyond several R 500 , implying a more suppressed matter power spectrum than predicted by most recent simulations. Complementary observables (e.g., thermal Sunyaev–Zel’dovich effect and fast radio bursts) and next-generation simulations are critical to understanding the physical mechanism behind this extreme gas expulsion and mapping its impact on the large-scale matter distribution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

IM3/HyperFACETS Thermodynamic Global Warming (TGW) Simulation Datasets

Publication For a thorough description of the methods, see the peer-reviewed paper: Jones, A.D., Rastogi, D., Vahmani, P. et al. Continental United States climate projections based on thermodynamic modification of historical weather. Sci Data 10, 664 (2023). https://doi.org/10.1038/s41597-023-02485-5 Overview The IM3 / HyperFACETS climate simulations provide 40-year historical (1980-2019) as well as four 80-year future simulations (2020-2099) over the U.S. The future simulations are split into near (2020-2059) and far future (2060-2099) segments. The future scenarios span a range of plausible changes in future climate (both Global Circulation Model (GCM) and Representative Concentration Pathways/Shared Socioeconomic Pathway (RCP/SSP) dimensions). The simulations provide climate variables with high spatiotemporal resolution (25 hourly variables and 207 3-hourly variables at 12 km2). The datasets are generated using dynamical downscaling with the WRF (Weather Research and Forecasting) model (version 4.2.1) and therefore preserve physical consistency across variables. WRF is a state-of-the-art, fully compressible, non-hydrostatic, mesoscale numerical weather prediction model. WRF is coupled with an urban canopy model (UCM), which resolves urban surfaces. The future scenarios were developed using a thermodynamic global warming approach where past events are replayed under a range of future warming conditions. These scenarios therefore provide a perspective on potential increases in extreme event intensity, geographic scope, and duration, with previously non-extreme conditions potentially crossing new thresholds to be considered extreme by today's standards. This approach is not intended to estimate future changes in extreme event frequency that might result from changes in large-scale atmospheric dynamics. This dataset has NOT been bias corrected. A bias corrected version of selected variables is under development and will be released here when available. Scenarios Files Data for each scenario is provided in weekly NetCDF files. 25 variables are available at hourly resolution, and 207 variables are available at three-hourly resolution. Spatial resolution is 12km and spans the conterminous United States (CONUS), including some areas of Canada and Mexico, resulting in a grid of 424 by 299 cells. The spatial projection is a Lambert Conformal Conic with the following proj-string: "+proj=lcc +lat_0=40.0000076293945 +lon_0=-97 +lat_1=30 +lat_2=45 +x_0=0 +y_0=0 +R=6370000 +units=m +no_defs". The available scenarios and simulation periods are listed below: historical | 1980 - 2019 rcp45cooler | 2020 - 2059 rcp45cooler | 2060 - 2099 rcp45hotter | 2020 - 2059 rcp45hotter | 2060 - 2099 rcp85cooler | 2020 - 2059 rcp85cooler | 2060 - 2099 rcp85hotter | 2020 - 2059 rcp85hotter | 2060 - 2099 * The first year (1979, 2019, and 2059) of data within each scenario represents a model warmup period and should not be used. These are located in the `spinup_files` directory. Historical year 2020 is considered an extra year of data beyond the simulation period and can be found in the `additional_files` directory. For information on specific variables and a more in-depth discussion of methodology, please refer to the data landing page at https://tgw-data.msdlive.org. Delta Warming Files The global and CONUS warming deltas for each scenario are provided in degrees Celsius annually and monthly. Restart Files Yearly restart files are provided for each scenario which can be used to restart the WRF model at a particular point in time. Spinup Files The first year of data within each simulation period represents a model warmup period and should not be used. The files are provided here for the sake of reproducibility. Additional Files Additional years of data are provided as an extension of the historic simulation.

Jones, Andrew D.↗
DOE OSTI2024

Cold Climate Degradation: An Analysis of Double-Axis Tracked, E-W Vertical, and Fixed-Tilt Photovoltaic Deployments in Alaska

As countries around the world transition towards renewable energy, there is increasing interest in using photovoltaic (PV) technologies to help decarbonize remote northern communities due to their scalability and affordability. However, a major barrier towards large-scale adoption of PV in cold climates is performance uncertainty under extreme environmental conditions including snowfall, freeze-thaw cycles, and high wind loads. Existing literature on PV degradation rates in the North is relatively limited, with published degradation rates varying between -0.2%/year (Sweden) to -1.3%/year (Scotland). At this workshop, we will present preliminary results on the long-term performance of two diverse photovoltaic sites located in Fairbanks, Alaska at 64.8 degrees N: a monofacial Al-BSF double-axis tracking site maintained by the Cold Climate Housing Research Center (CCHRC), and a bifacial PERC/SHJ E-W vertical and south-facing fixed-tilt site maintained by the Alaska Center Energy and Power (ACEP). CCHRC data has been collected over a period of 15 years, while ACEP site data has been collected over 4 years. Using the degradation analysis tool, RdTools, we will present annual system degradation rates, seasonal performance ratio, and identify potential cold climate failure mechanisms for commercially available PV technologies. This analysis will add to existing literature by directly comparing the performance of multiple PV configurations in Alaska.

bifacial↗
DOE OSTI2026

Robotics for Energy Infrastructure Construction on Earth and Beyond: Pilot Study Exploring Level 3-4 Semi-Autonomous Installation of a Reduced-Scale Photovoltaic Power Plant Prototype: Preprint

There is a growing need to build energy infrastructure at massive scale on both Earth and in space. On Earth, increasing electricity demand from applications such as AI data centers is driving rapid expansion of power systems, while in space, long-duration missions such as lunar operations require reliable and resilient energy infrastructure that can be deployed quickly. These needs outpace the capabilities of conventional construction methods. Robotic construction offers a path to faster, lower-cost, and more consistent deployment, yet most existing approaches focus on automating individual tasks rather than end-to-end assembly. The Industrialized Construction Innovation (ICI) team at the National Laboratory of the Rockies (NLR), in collaboration with the Colorado School of Mines, has been investigating robotics and automation strategies for rapid end-to-end energy infrastructure construction. This paper presents a pilot study demonstrating level 3-4 semi-autonomous, end-to-end robotic assembly of a reduced-scale photovoltaic power plant prototype under terrestrial conditions. The approach combines hardware redesign of photovoltaic components for robotic compatibility, discrete-event and physics-based robotic simulations in synthetic environments, and sim-to-real deployment to execute end-to-end robotic assembly of[NM1.1][NM1.2] tracking components, structural members, photovoltaic modules, and embedded electrical connections. The pilot identifies key limitations like tolerance management, perception and localization, end effector capability, fastening methods, and multi-physics simulation fidelity, that must be addressed to enable end-to-end robotic assembly of full-scale energy infrastructure. The study highlights how these insights inform scaling strategies for terrestrial energy infrastructure and adaptation to unstructured and extreme environments, including off-world scenarios such as robotic assembly of vertical tall lunar towers.

14 SOLAR ENERGY↗
DOE OSTI2026

Integrating agrivoltaics into smallholder farming systems to enhance food production and irrigation efficiency under climate stress

Smallholder agricultural growers and entrepreneurs (SAGE) have a crucial role in sustaining food security; their cultivation accounts for roughly 30% of global crop production and contributes substantially to the diversity and nutritional value of food systems. Despite their awareness of and flexibility to adapt to productivity challenges, SAGE communities remain highly vulnerable to heat and drought stress, which threatens both livelihoods and local food availability. While intensification and adaptive farm practices offer partial resilience, their capacity to buffer sudden climatic extremes remains limited. Agrivoltaics, the dual use of land for solar electricity generation and crop cultivation, has emerged as a promising strategy to enhance agricultural productivity and its resilience to climate extremes by altering microclimatic conditions. Although proven effective in hot and arid regions, the benefits of agrivoltaics for temperate climates and small-scale farms remain minimally explored. To address this gap, we evaluated the performance of a small-scale agrivoltaics system in a temperate climate for high-value crops, including leafy vegetables and garlic. We investigate whether the system can (i) protect crops during extreme heat events and (ii) enhance productivity and reduce irrigation requirements during hot and dry periods unsuitable for conventional production. Our findings provide evidence that agrivoltaics is a climate-resilient farming strategy under current and projected climate scenarios, capable of improving yields (by +43% to +127% for the leafy vegetables grown) and reducing water consumption, while creating complementary economic opportunities through decentralized energy generation systems. This work supports the integration of agrivoltaics into small-scale agricultural systems as an innovative pathway to strengthen food security, bolster farmer livelihoods, and enable multiple co-benefits from broader solar energy adoption.

14 SOLAR ENERGY↗
DOE OSTI2020

Coupling of Spark Plasma Sintering with Advanced Modeling to Enable Process Scale-Up: Presentation to DOE-NE [Slides]

The research goal of this project is to develop at Idaho National Laboratory (INL) a first-of-its-kind Multiphysics Object-Oriented Simulation Environment (MOOSE)-based, multiscale, multiphysics spark plasma sintering (SPS) modeling and simulation code application, termed “Freya.” Freya will simulate the thermo-mechanical-electrical aspects of the SPS fabrication process and will be paired with lower length scale sub-models, such as phase-field, to predict the resulting microstructure. SPS is an advanced manufacturing process that can be used to solve a variety of material manufacturing challenges; however, this process is an extremely challenging problem for modeling and simulation. The SPS process is inherently multiphysics and multi-scale, with the macroscale electro-thermo-mechanical behavior linked intricately to the microstructure evolution of the part being sintered. Accurate modeling and simulation tools, specifically geared towards the SPS process, are needed to predict the influence of the multiple variables involved in the manufacturing process. Modeling and simulation accuracy is achieved and demonstrated through comparison to multiple validation experiments. The validation efforts for Freya include both separate effects and complete multiphysics SPS process experiments. One of the key benefits this Laboratory Directed Research & Development (LDRD) project offers stems from the emphasis placed on experimental validation of the Freya models, both on the individual length scales and of the final coupled multiscale multiphysics simulations. Experimental validation of Freya’s multiscale coupling capability provides the technical credibility necessary for potential future industry and research partners to accept the simulation predictions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Enzyme selection, optimization, and production toward biodegradation of post–consumer poly(ethylene terephthalate) at scale

Poly(ethylene terephthalate) (PET) is one of the world's most widely used polyester plastics. Due to its chemical stability, PET is extremely difficult to hydrolyze in a natural environment. Recent discoveries in new polyester hydrolases and breakthroughs in enzyme engineering strategies have inspired enormous research on biorecycling of PET. This study summarizes our research efforts toward large-scale, efficient, and economical biodegradation of post-consumer waste PET, including PET hydrolase selection and optimization, high-yield enzyme production, and high-capacity enzymatic degradation of post-consumer waste PET. First, genes encoding PETase and MHETase from Ideonella sakaiensis and the ICCG variant of leaf-branch compost cutinase (LCC ICCG ) were codon-optimized and expressed in Escherichia coli BL21(DE3) for high-yield production. To further lower the enzyme production cost, a pelB leader sequence was fused to LCC ICCG so that the enzyme can be secreted into the medium to facilitate recovery. To help bind the enzyme on the hydrophobic surface of PET, a substrate-binding module in a polyhydroxyalkanoate depolymerase from Alcaligenes faecalis (PBM) was fused to the C-terminus of LCC ICCG . The resulting four different LCC ICCG variants (LCC, PelB-LCC, LCC-PBM, and PelB-LCC-PBM), together with PETase and MHETase, were compared for PET degradation efficiency. A fed-batch fermentation process was developed to produce the target enzymes up to 1.2 g L –1 . Finally, the best enzyme, PelB-LCC, was selected and used for the efficient degradation of 200 g L –1 recycled PET in a well-controlled, stirred-tank reactor. Furthermore, the results will help develop an economical and scalable biorecycling process toward a circular PET economy.

59 BASIC BIOLOGICAL SCIENCES↗