Search NASA⌕ Search

SEARCH · Search NASA

Results for “Process Flow”

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 253 records · Page 14

Solid State Solar Thermochemical Fuel (SoFuel) for Long Duration Storage

Efficient thermal storage systems, when coupled with renewable energy, enable the decarbonization of numerous industrial processes requiring high temperature steam or air, and provide a path for seasonal building heating, especially for colder climates. Existing thermal storage systems face a significant challenge due to losses inherent to all high temperature systems. A viable route to long-term storage is to use thermochemical reactions to convert concentrated solar energy to a fuel that is shelf-stable and can be stored at room temperature, thus eliminating losses associated with high temperature storage. The Solid-State Solar Thermochemical Fuel (SoFuel) technology developed by Michigan State University, Oregon State University, and Mississippi State University provides reactors and processes with minimal sensible heat losses and allows storing solar energy as a solid-state fuel at room temperature for long duration. The production of SoFuel occurs within a cylindrical cavity reduction chemical reactor that captures concentrated solar radiation from a solar field. Reactive magnesium manganese oxide (Mg-Mn-O) resides within the cylindrical cavity chemical reactor and undergoes thermal reduction as the temperature exceeds 1350°C. The thermally reduced Mg-Mn-O pellets (the SoFuel) are cooled down through a recuperative process and stored within a bin until used. The SoFuel can directly supply up to 1100C heat to an adjacent power plant for electricity generation or industrial heating. Oxidation of SoFuel pellets occurs in a counter flow reactor and supplies heat to the user for electricity generation or industrial processing, after which the fuel is returned to the concentrating solar field where it is regenerated for re-use. Both reactors can be controlled well using a variety of strategies. With the low cost of the material, its cyclability, and the possibility of using the pelletized with on-sun reactors, or with electricity that would be curtailed, this project offers a viable option of medium- and long-term thermal energy storage.

25 ENERGY STORAGE↗

Solid State Solar Thermochemical Fuel (SoFuel) for Long Duration Storage

Efficient thermal storage systems, when coupled with renewable energy, enable the decarbonization of numerous industrial processes requiring high temperature steam or air, and provide a path for seasonal building heating, especially for colder climates. Existing thermal storage systems face a significant challenge due to losses inherent to all high temperature systems. A viable route to long-term storage is to use thermochemical reactions to convert concentrated solar energy to a fuel that is shelf-stable and can be stored at room temperature, thus eliminating losses associated with high temperature storage. The Solid-State Solar Thermochemical Fuel (SoFuel) technology developed by Michigan State University, Oregon State University, and Mississippi State University provides reactors and processes with minimal sensible heat losses and allows storing solar energy as a solid-state fuel at room temperature for long duration. The production of SoFuel occurs within a cylindrical cavity reduction chemical reactor that captures concentrated solar radiation from a solar field. Reactive magnesium manganese oxide (Mg-Mn-O) resides within the cylindrical cavity chemical reactor and undergoes thermal reduction as the temperature exceeds 1350°C. The thermally reduced Mg-Mn-O pellets (the SoFuel) are cooled down through a recuperative process and stored within a bin until used. The SoFuel can directly supply up to 1100°C heat to an adjacent power plant for electricity generation or industrial heating. Oxidation of SoFuel pellets occurs in a counter flow reactor and supplies heat to the user for electricity generation or industrial processing, after which the fuel is returned to the concentrating solar field where it is regenerated for re-use. Both reactors can be controlled well using a variety of strategies. With the low cost of the material, its cyclability, and the possibility of using the pelletized with on-sun reactors, or with electricity that would be curtailed, this project offers a viable option of medium- and long-term thermal energy storage.

14 SOLAR ENERGY↗

De-Risking High-Recovery Brackish Water Desalination via Flow Reversal and Feed Flushing Using Techno-Economic Assessment

Novel desalination technologies have demonstrated enhanced performance and improved financial metrics over existing processes adopted by industry. Establishing quantitative performance targets is essential for achieving financial benefits over the current state of the art. Herein, we demonstrate how WaterTAP, a techno-economic assessment (TEA) tool, can be used to identify minimum performance metrics necessary to achieve financial benefit over using existing processes. This study evaluates the feasibility of increasing water recovery at the Chino Desalter I above 90 % through the addition of a third variable configuration reverse osmosis (VCRO) stage. Sensitivity analyses revealed flow reversal frequency, feed flushing volume (used as a cleaning step), and membrane lifespan are key factors influencing the financial viability of the VCRO process. The TEA analysis demonstrated that the system must achieve a recovery of 84 % and a 1-year membrane lifespan to have a breakeven LCOW, while achieving 90 % recovery can ultimately reduce the LCOW by 16 %. Notably, a trade-off between decreasing frictional losses and increased osmotic pressure across the recovery range, resulted in a stable specific energy consumption across the recovery range, enabling meaningful LCOW reductions without an energy penalty, a key finding that contrasts with conventional RO. This work demonstrates how TEA can guide system design by identifying key performance targets and exploring trade-offs, enabling data-driven decisions to de-risk innovative processes. These findings underscore the importance of leveraging TEA to evaluate scaling mitigation strategies and optimize inland desalination systems for sustainable and cost-effective operation.

14 SOLAR ENERGY↗

Scalable Risk Assessment of Rare Events in Power Systems With Uncertain Wind Generation and Loads

Risk assessment of rare events has become increasingly important in power system planning and operation with the increasing integration of renewable energy and the presence of system uncertainties. However, quantifying the risk posed by rare events via the traditional method, i.e., Monte Carlo sampling (MCS), incurs substantial computational expense stemming from the vast ensemble of power flow simulations. To accelerate the assessment, this paper proposes a Deep Neural Network (DNN)-kernelized vector-valued Gaussian Process (VVGP) approach with excellent computational efficiency while maintaining high accuracy. Consequently, serving as a surrogate model for the power flow solver, the DNN-kernelized VVGP enables significantly faster but accurate risk assessment compared to the power flow solver. The developed surrogate model evaluates low-order N - k events that contain more than 90% instances by adeptly capturing the topological features while the high-order N - k events are assessed via a power flow solver, thereby striking a balance between computational efficiency and uncertainty quantification accuracy. Moreover, the model incorporates a Support Vector Machine (SVM) classifier to resample concerning low-probability tail events to counteract the biases potentially introduced during the DNN-kernelized VVGP evaluations. Simulations conducted on the modified IEEE 24-bus, 118-bus, and European 1354-bus systems demonstrate that the proposed method maintains the accuracy benchmark set by MCS while significantly reducing computational demands in large-scale power systems as compared to other state-of-the-art methods.

17 WIND ENERGY↗

Polymer coatings on zirconia microspheres via rotating flow fluid dynamics

Conventional tristructural isotropic (TRISO) coatings for nuclear fuels require multi-step and expensive formation processes; any breakage of the coatings may lead to fission species release. Polymer derived ceramic (PDC) coatings can be a suitable alternative to address these issues. In this study, allylhydridopolycarbosilane (SMP-10) coatings were created on yttria stabilized zirconia (YSZ) microspheres using a Rotating Flow Fluid Dynamics (RFFD) coating method. The effects of curing temperature, rotation speed, and coating cycle/time on the coating were analyzed. Surface functionalization of YSZ microspheres with NaOH resulted in good adhesion between the polymer precursor and the YSZ kernel particles. Spectroscopy analysis revealed complete curing of SMP-10 coated YSZ at 160 °C. Rotation speed and coating time significantly affect the coating characteristics. After 3 cycles at lower rotation speed (50 rpm) for 10 min each, the coating obtained was uniform and homogeneous. In comparison, the coating showed almost 10 times less eccentricity at a high rotation speed of 300 rpm. Compared with the coatings prepared at 300 rpm condition, the coatings have higher sphericity and lower eccentricity at 50 rpm. Overall, this study provides a novel and effective route for fabricating and curing SMP-10 precursor coatings on YSZ microspheres.

Ravi, Nivetha [University of Alabama, Birmingham]↗

Stochastic fracture generation and thermo-hydro-mechanical modeling in an equivalent continuum framework for enhanced geothermal systems

Enhanced geothermal systems (EGS) involve fracturing low permeability material to establish well connectivity and then injecting and circulating fluid into the fractured subsurface for geothermal power production. Changes in fracture aperture from contraction of the cooling matrix rock may alter network connectivity and risk thermal short-circuiting. Thermo-hydro-mechanical (THM) models are a useful tool to study these processes. However, as fracture networks are complex, and data may be limited, fracture networks in THM models are often stochastically generated. Given reliance on stochastic fracture networks and THM modeling to represent the subsurface and assess productivity of EGS, increased understanding of the influence of such statistically derived fracture networks on flow and heat transport in THM models is needed. Here, a new fracture process model is developed in the reactive transport code PFLOTRAN to stochastically generate fracture families and simulate changes in fracture aperture over time due to temperature changes of the rock matrix. Sixty-four different fracture networks ranging from well to poorly-connected, are modeled in PFLOTRAN with and without mechanical processes (THM vs TH). Results indicate that for well-connected fracture networks, thermal short-circuiting is less of a concern due to the abundance of available alternative flowpaths. For poorly-connected fracture networks, inclusion of mechanical processes showed steep thermal drawdown coincident with increase in fracture aperture along developing colder flowpaths, demonstrating the risk of thermal short-circuiting. Simulations with additional, larger fractures engineered to establish connectivity in a poorly-fractured subsurface, indicate that while stochastic variation of fracture orientation of the background network had limited influence, such variation in the engineered fractures significantly affected flow and heat transport.

Discrete fracture networks (DFN)↗

Mechanically Accelerated Depolymerization of Entangled Linear Polymer Melts

Mechanical forces can enhance the chemical depolymerization of synthetic polymers when shear flow accelerates chain scission. To quantify the extent of mechanically-accelerated scission, the effect of simple shear flow (duration and strength) with low Weissenberg and Deborah numbers was investigated by considering the impact of applied work in both simple shear and shear dominated mixed flows. Hydrogenated polyisoprene was chosen as a model linear, entangled system. The conditions (strain amplitude, frequency, and shearing time) necessary to increase chain scission were assessed in the rubbery melt. Shear flow accelerated chain scission at higher temperatures, suggesting an activated process. Isothermal scission versus work curves were superposed by applying shift factors a T,S , whose Arrhenius-like temperature dependence gave an apparent activation energy for chain scission of ~ 110 kJ/mol, which is likely a combination of the activation energy of viscosity and bond energy. This work provides a base for quantifying the impact of shear on depolymerization of polymer melts and highlight the connection between viscous dissipation and scission chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanistic Insights for Plasma-Catalytic CO 2 Reduction over TiO 2 in a Dielectric Barrier Discharge Reactor

Reaction kinetics experiments coupled with phenomenological kinetic modeling and parameter estimation are used to elicit insights into the mechanism and active sites for the plasma-catalytic dissociation of CO 2 on TiO 2 . Experimental and model insights showed that gas-phase reactions contribute at least two-thirds of the overall product formation at explored conditions; weak temperature dependence, strong sensitivity to specific energy input (SEI), apparent first order in CO 2 , and positive influence of cofed argon (Ar) and oxygen (O 2 ) for the gas-phase contributions all suggest that expected plasma reaction steps such as electron-impact and high-energy collisions are the dominant modes for CO 2 dissociation. The Arrhenius-like expression for gas contributions resulted in a preexponential of 4.40 × 10 –3 s –1 , an E SEI,g of 7.90 × 10 –4 mol/kJ, and an E a,g of 1.00 × 10 –3 J/mol. For surface contributions, the small apparent barrier of 16.3 kJ/mol, relatively weaker dependence on SEI, first-order dependence on CO 2 , and insensitivity to cofed Ar and O 2 all point to CO 2 dissociation on TiO 2 surface facets without vacancies and aided by plasma (leading to vibrationally excited CO 2 and/or a reactive surface with significant surface charge accumulation). The Arrhenius-like expression resulted in a preexponential of 7.81 × 10 –2 s –1 , an E SEI,s of 1.90 × 10 –3 mol/kJ, and an E a,s of 1.63 × 10 4 J/mol. The derived kinetic model further enabled a systematic evaluation of the effect of inputs (plasma power, flow rate, CO 2 inlet concentration, and temperature) to identify process trends and optimal operating conditions.

catalyst↗

Experimental analysis of conductive drying of paper and board

The primary heat transfer source in the conventional multi-cylinder drying of paper and board is conduction. Conduction is facilitated by high-tension contact with steam-heated cylinders, while convective drying, the main mass transfer source, operates as heated air flows over the paper web in the pockets. This conduction process occurs at elevated temperatures and contact pressures to ensure effective contact between the wet paper web and the heated cylinders. The contact pressures and temperatures of the hot surface significantly influence the conductive drying characteristics. This research paper presents an experimental study that involves designing a simple lab-scale setup with in-situ and continuous sensing capabilities for various commercially available grades of paper and board. Embedded thermocouples measure the temperatures of the heated platen and the sheet, allowing the collection of flux data to determine heat transfer characteristics. Here, the goal of this study is to ascertain the instantaneous contact coefficients for different basis weights as a function of moisture content. The acquired data will provide valuable insights and information towards process development, design and simulation of the paper drying process.

42 ENGINEERING↗

Continued Testing and Modeling of the AirSquared All-Metal Tritium Scroll Pump

The implementation of fusion energy requires the processing of the deuterium-tritium mixture used to fuel the reaction. Many gas processing operations require either a gas-motive force or a pressure gradient across a system; both are achieved using vacuum pumps. However, because of the material challenges of working with tritium, specialized pumps are required that use only metal on wetted surfaces. This prevents tritium retention and degradation of materials such as oils and polymers and reduces the generation of tritiated waste. Typical gas handling systems for tritium operations involve metal bellows pumps, all-metal scroll pumps, and metal turbomolecular pumps, depending on the pressure requirements. With the industry standard Normetex Type 15 m 3 /h (typically referred to as the Normetex 15) pump no longer available, the fusion market requires viable alternatives. In this study, we test an American-made pump, the AirSquared V16H034A-C01, and compare its performance to the Normetex. A performance model of the AirSquared pump is also produced, which can be applied to general process models to generate an accurate description of pressures and flows through various unit operations.

08 HYDROGEN↗

Distributed Sapphire Fiber Bragg Gratings Based Thermal Profiling of Submerged Entry Nozzles

This research focuses on the application of sapphire fiber Bragg gratings (FBGs) for instrumentation in submerged entry nozzles (SEN) within the steelmaking industry. The SEN is pivotal for transferring molten steel from a tundish to a mold, while preventing the infiltration of oxygen and nitrogen from the surrounding environment. Maintaining optimal flow conditions in the mold is crucial for ensuring casting process stability and maintaining high quality steel. Sapphire FBG sensors have been instrumented in SENs to enable distributed thermal mapping for monitoring the health of the SEN. The optical sensor comprises three cascaded sapphire FBGs inscribed using femtosecond laser technology into a one-meter-long sapphire crystalline fiber. The sensor underwent characterization in a laboratory setting up to 1600°C and was tested for long-term stability over 40 hours under extreme environmental conditions. The coupling between silica and sapphire fibers was investigated and implemented during sensor packaging. Here, the sensor successfully captured the pre-heat sequence of the SEN in real-world steelmaking operations. Compared to conventional thermocouples, sapphire FBG sensors demonstrated exceptional efficiency and precision. They offer potential benefits such as increased productivity, reduced energy consumption, and minimized carbon footprint in the steel industry.

Sapphire Fiber Bragg Grating↗

Particle Handling in Particle-Based Pumped Thermal Energy Storage

A particle-based pumped thermal electricity storage system stores high-temperature heat (~1000 degrees C) in low-cost silica sand and generates power through an efficient power cycle. Central to this system is a counterflow direct-contact gas/particle fluidized-bed heat exchanger, which can significantly improve the heat exchange process due to the large heat transfer surface area of particles. To showcase a lab-scale 10- to 20-kWe heat exchange process with particles heating up to 300 degrees C, a comprehensive hydrodynamic analysis of the fluidization condition inside the heat exchanger was conducted. The heat exchanger was designed to target for heat exchange process of 10-20 kW and accommodate an air mass flow rate of 0.3-0.7 kg/s. The fluidization condition was set to maintain stable bubbling fluidization, thereby maximizing the particle-to-air heat transfer. An air distributor was designed to equally distribute air over the bed, which can avoid defluidized zones in the heat exchanger. Additionally, particle handling systems including L-valves, screw conveyors, and pneumatic conveyors were developed for the prototype heat exchanger, efficiently conveying high-temperature particles at 0.2-0.5 kg/s. This work lays the foundation for scaling up the system and integrating it into larger energy storage applications, demonstrating its potential for efficient, high-temperature thermal energy storage and power generation.

25 ENERGY STORAGE↗

Monitoring Sulfuric Acid and Temperature Using Raman Spectroscopy and Multivariate Chemometrics

Multivariate regression models were optimized for the quantification of sulfuric acid (H 2 SO 4 ) [0–8 M] and temperature (20 °C–80 °C) in the presence of ammonium sulfate ((NH 4 ) 2 SO 4 [0–0.6 M]) using Raman spectroscopy. Optical vibrational spectroscopy is a useful nondestructive technique for the in situ analysis of complex chemical systems notoriously difficult to monitor in situ and in real-time. Multivariate analysis, a chemometrics method, can be paired with these nondestructive optical methods for determining analyte concentration and speciation in complex solutions, such as dissociated species in polyprotic acids, e.g., H 2 SO 4 . The effect of temperature is often overlooked although it can have a major influence on speciation and the corresponding Raman spectra. Here, in this study, partial least squares regression models were optimized for the quantification of H 2 SO 4 and its two deprotonated forms as a function of temperature. Measuring bisulfate as a function of temperature is particularly challenging owing to changes in the second dissociation constant. A designed training set effectively minimized the sample set size and trained a robust predictive model with percent root mean square error of <3% for H 2 SO 4 . The practical strategy employed here was demonstrated to be effective for building chemometric models that directly account for dynamic temperatures with static samples and is shown to be amenable to flow cell analysis applications with a simple calibration transfer for process monitoring applications.

D-optimal design↗

Granular Temperature Modeling with Large-Scale CFD-DEM Data

The granular temperature is a measure of solids-phase fluctuating kinetic energy in particle-laden gas-solid flows. This report focuses on reduced-order models of granular temperature in triply-periodic cluster induced turbulent flow. The granular temperature data was generated by post-processing large-scale CFD-DEM simulations that were carried out as part of an ALCC awarded study. The algebraic model of Tang et al, originally fit to DNS data, was found to be amenable to the present data. The results of several regression analyses are presented building on the original expression. This work was performed in part through the National Nuclear Security Administration (NNSA) Minority Serving Institutions Internship Program (MSIIP).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Code for the manuscript "Mori-Zwanzig Modal Decomposition"

We would like to create an open source repository in LANL's github on code written in Julia, in which we implement and extend the data-driven Mori-Zwanzig method for extracting large-scale spatio-temporal structures from data, which we call MZMD. This method is an extension of Dynamic Mode Decomposition (DMD) in which Mori-Zwanzig memory kernels are included into the associated companion matrix. In the code we would like to release, we apply MZMD to a flow over a cylinder with Reynolds number 100 rather than the much larger data set used in the associated manuscript. DMD is used extensively in the fluid dynamics community mainly for extracting large scale spatio-temporal structures (patters) from flow data. This is useful for understanding the key mechanisms that generate certain complex dynamical process relevant in engineering design. In MZMD, we improve upon DMD by adding the Mori-Zwanzig memory kernels, and show this improvement is especially important in strongly nonlinear regions of the flow.

Woodward, Michael↗

Quantifying concentration distributions in redox flow batteries with neutron radiography

Abstract The continued advancement of electrochemical technologies requires an increasingly detailed understanding of the microscopic processes that control their performance, inspiring the development of new multi-modal diagnostic techniques. Here, we introduce a neutron imaging approach to enable the quantification of spatial and temporal variations in species concentrations within an operating redox flow cell. Specifically, we leverage the high attenuation of redox-active organic materials (high hydrogen content) and supporting electrolytes (boron-containing) in solution and perform subtractive neutron imaging of active species and supporting electrolyte. To resolve the concentration profiles across the electrodes, we employ an in-plane imaging configuration and correlate the concentration profiles to cell performance with polarization experiments under different operating conditions. Finally, we use time-of-flight neutron imaging to deconvolute concentrations of active species and supporting electrolyte during operation. Using this approach, we evaluate the influence of cell polarity, voltage bias and flow rate on the concentration distribution within the flow cell and correlate these with the macroscopic performance, thus obtaining an unprecedented level of insight into reactive mass transport. Ultimately, this diagnostic technique can be applied to a range of (electro)chemical technologies and may accelerate the development of new materials and reactor designs.

Science & Technology - Other Topics↗

Anisotropic Turbulent Flows Observed in Above-the-loop-top Regions during Solar Flares

Abstract Solar flare above-the-loop-top (ALT) regions are vital for understanding solar eruptions and fundamental processes in plasma physics. Recent advances in three-dimensional (3D) magnetohydrodynamic (MHD) simulations have revealed unprecedented details on turbulent flows and MHD instabilities in flare ALT regions. Here, for the first time, we examine the observable anisotropic properties of turbulent flows in ALT by applying a flow-tracking algorithm on narrow-band extreme-ultraviolet images that are observed from the face-on viewing perspective. First, the results quantitatively confirm the previous observation that vertical motions dominate and that the anisotropic flows are widely distributed in the entire ALT region with the contribution from both upflows and downflows. Second, the anisotropy shows height-dependent features, with the most substantial anisotropy appearing at a certain middle height in ALT, which agrees well with the MHD modeling results where turbulent flows are caused by Rayleigh–Taylor-type instabilities in the ALT region. Finally, our finding suggests that supra-arcade downflows (SADs), the most prominently visible dynamical structures in ALT regions, are only one aspect of turbulent flows. Among these turbulent flows, we also report the antisunward-moving underdense flows that might develop due to MHD instabilities, as suggested by previous 3D flare models. Our results indicate that the entire flare fan displays group behavior of turbulent flows where the observational bright spikes and relatively dark SADs exhibit similar anisotropic characteristics.

Xie, Xiaoyan (ORCID:0009000705827807)↗

Hydrologic Model Data for the East Fork Poplar Creek Watershed Simulated with the Advanced Terrestrial Simulator (ATS): Streamflow and Network Expansion–Contraction Dynamics

This dataset supports hydrologic modeling and stream network expansion–contraction analysis for the East Fork Poplar Creek (EFPC) Watershed in Tennessee. It includes a Jupyter notebook for model setup, model configuration files, simulation outputs, and derived products used to evaluate model performance and investigate stream dynamics under varying hydrologic conditions. The dataset was generated using the Watershed Workflow Python package and the Advanced Terrestrial Simulator (ATS), enabling integrated surface–subsurface hydrologic simulations using a stream-aligned mesh. Outputs include high-resolution time series of streamflow, active network length, water table depth, and related hydrologic variables. Also included are spatially explicit stream persistency indices and classifications of reaches as perennial or non-perennial. These data facilitate reproducibility and support further research on stream intermittency and variability in network extent.The model data archive is organized in following directories:1) model_setup_inputsContains the Watershed Workflow Jupyter notebooks (accessed through any open source code editor), selected input datasets, and resulting ATS input files, including XML files (access through any open source code editor), computational mesh (.exo files can be viewed using Paraview), and meteorological forcing files (.h5 files can be accessed through h5py python package and HDFView open source software). 2) model_outputsIncludes ATS simulation outputs relevant to this study. Time series of spatially integrated or averaged variables (e.g., streamflow, water table depth) are provided as CSV files. Select spatial fields (e.g., ponded depth and water table depth) are saved as pickled Python objects to reduce file size, and can be accessed through pickle package in Python. Key geometry objects from Watershed Workflow—such as the surface mesh and river tree—are also included to support analysis of streamflow persistency and expansion–contraction dynamics. These files can also be accessed through Watershed Workflow Python package.3) model_evaluationProvides observed streamflow time series and field survey-based flow regime classifications used to evaluate model performance. Jupyter notebooks for processing ATS outputs and comparing model predictions with observations to build confidence in the model prior to scientific analysis are also included.4) Q_L_relationshipsContains workflows for generating time series of discharge, active network length, and related hydrologic variables used in the stream network expansion–contraction analysis. Includes routines for delineating baseflow-dominated periods. For each catchment, notebooks and processed data (as pickled DataFrames accessed through Pandas Python package) are provided. 5) figure_scriptsProvides the Jupyter notebooks used to generate the figures presented in the paper.

54 ENVIRONMENTAL SCIENCES↗