Search NASA⌕ Search

SEARCH · Search NASA

Results for “Process simulation”

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

BCARS Simulated Phantom Dataset for Evaluation of Processing Pipelines

Broadband coherent anti-Stokes Raman scattering (BCARS) microscopy is a powerful label-free biological imaging technique, but the raw signal requires careful processing. The vibrationally resonant (Raman) fingerprint signal is usually small compared with instrumental noise sources and the nonresonant background (NRB) inherent in the BCARS signal. Fortunately, the NRB exhibits a systematic phase relationship with the coherent Raman response, acting as a heterodyne amplifier for the weak fingerprint signal. Due to this heterodyne effect, the Raman response can be recovered quantitatively and invariantly across different instruments, provided the NRB shape is known. Even with heterodyne amplification, the amplitudes of fingerprint signal components are often comparable to system noise. Singular value decomposition (SVD), which utilizes spatial information, is often employed for additional noise filtering. Consequently, finding optimal processing parameters to properly distinguish the NRB and Raman responses and suppress noise in the complex BCARS signal requires a reference system that realistically represents the spectral and spatial properties of BCARS signals obtained from biological samples. We present a digital tissue phantom that meets these criteria as a tool for testing candidate signal processing pipelines. The digital phantom is generated with simulated hyperspectral Raman images having system-specific noise and background characteristics. Here, we analyze phantom datasets with differing background and signal-to-noise conditions to evaluate their impact on the performance of multiple signal processing pipelines. Specifically, we investigate the application of a Butterworth filter-based routine to directly estimate the NRB from the BCARS signal. Additionally, we evaluate a Lorentzian wavelet transform as an alternative to the Hilbert transform for extracting the Raman spectrum from the BCARS signal. While we demonstrate this phantom for BCARS, it can be used for any spectroscopic Raman imaging approach.

Dixon, Jessica Z. [Georgia Institute of Technology↗
DOE OSTI2024

Noble Metals Testing on the DM1200 Pilot Melter System with AZ-101 HLW Simulants, VSL-03R3740-2, Rev. 0

This report documents the results obtained with respect to the behavior of noble metals in tests on the DM1200 HLW Pilot Melter during processing of simulated HLW AZ-101 feed. The principal objective of these tests was to collect data to determine the extent of retention of noble metals in the melter since that is an important uncertainty in determining the expected life of the WTP HLW melter. The present results constitute one part of an integrated testing strategy that has been developed to address this risk.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

On the relationship between precipitation extreme and local temperature over eastern China based on convection permitting simulations: roles of different moisture processes and precipitation types

The Clausius–Clapeyron (CC) scaling, which indicates a roughly 7% increase in saturated water vapor per 1 °C increase in temperature, can serve as a strong constraint linking the intensity of precipitation extremes and local temperature. However, the relationship between precipitation extreme and local temperature (referred to as the PE-T relationship) does not always follow the CC scaling and is highly dependent on climate regimes. In this study, we investigated the impacts of different moisture processes and precipitation types on the PE-T relationship over eastern China during the summertime based on convection-permitting model simulations. Consistent with observations, the simulated intensity of precipitation extremes increases with temperature at a rate close to CC (double-CC) scaling below (above) 20 °C. When the temperature exceeds 25 °C, precipitation intensity starts to drop. Precipitation extremes are mainly contributed by the stratiform, MCS (i.e., mesoscale convective system) convective, and non-MCS convective precipitation at low (< 20 °C), medium (20–25 °C), and high (> 25 °C) temperatures, respectively, suggesting that the double-CC scaling occurs when convective types become dominant, while the negative scaling at high temperatures is attributed to the reduced horizontal scale of convection. Corresponding to the reduced intensity of precipitation at high temperatures, there are stronger divergence and subsidence in the low-level atmosphere, which is probably caused by the net cooling associated with the enhanced melting and evaporation of falling hydrometeors due to the lower relative humidity in the low-level atmosphere. Overall, our findings contribute to a deeper understanding of the temperature dependence of precipitation extremes in eastern China.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Membrane-based solvent extraction for the recovery of rare earths from phosphate mining process streams

This study reports on the capture of rare earth elements (REEs) from phosphate industry process streams, including phosphoric acid (PA) sludge and phosphogypsum (PG), using a membrane solvent extraction (MSX) process. While MSX has been proven effective for a relatively concentrated feed, its effectiveness for dilute REEs solutions remains unexplored. Investigated PA-sludge and PG particles contain total REEs concentrations of ∼1100 and ∼320 ppm, respectively. Acid leaching, implemented to dissolve the REEs, significantly dilutes the REEs concentration to ∼210 ppm for PA-sludge leachate and ∼60 ppm for PG leachate. These low concentrations, compounded by the higher levels of non-REE ions and radioactive species, uranium (U) and thorium (Th), poses challenges to the MSX process. Here, we demonstrated that N,N,N′,N′-tetraoctyl-diglycolamide (TODGA) selectively binds REEs from a >3 M nitric-acid leachate while effectively rejecting U and Th. Concentrations of light REEs in strip solution were doubled compared to the feed, while heavy REEs were preferentially extracted. Furthermore, >99% purity gypsum, free of U and Th, was precipitated during the acid leaching process, aiding separation by removing significant amounts of non-REEs species (e.g., calcium) prior to the MSX process. Molecular simulations support the experimental data, suggesting preferential separation of heavy over light REEs. Based on these results, a cost-effective integrated process including pretreatment, acid leaching, MSX, and wastewater treatment is proposed for the co-recovery of REEs, phosphoric acid, gypsum, and U. This study shows MSX as a technically and economically feasible process for the recovery of REEs from low-concentration process streams, offering advantages over conventional solvent extraction.

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

Using MARCUS, MICRE, and COMBLE data to improve understanding and modeling of cloud, aerosol, and boundary layer processes at high-latitudes

Because it is believed general circulation (GCM) and numerical weather prediction (NWP) models underestimate shortwave radiation over the Southern Ocean (SO) due to inadequate representations of boundary layer (BL) and cloud processes, it is critical to improve our understanding of key aerosol, cloud, precipitation, and BL processes. Over the north Atlantic Ocean (NA), cold-air outbreaks are common, yet few studies of the aerosol and environmental controls of the associated convective BL clouds exist as needed to develop and evaluate GCM representations. Although processes cannot be observed, cloud and aerosol properties can be measured in-situ or remotely retrieved, which when combined with numerical simulations enable process level understanding required to improve model representations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Viscosity Measurements of the HLW APPS Glasses

The U.S. Department of Energy is responsible for building the Hanford Waste Treatment and Immobilization Plant (WTP) at the Hanford site in Washington to remediate 56 million gallons of radioactive waste historically stored in 177 underground tanks. The Office of River Protection has requested that the Savannah River National Laboratory (SRNL) contribute in areas of recognized capabilities and expertise for glass waste form development to support successful startup of the WTP. SRNL support of this work is defined in the Task Technical and Quality Assurance Plan.1 This report provides results from viscosity measurements on a series of simulated nuclear waste glasses designed and fabricated at Pacific Northwest National Laboratory (PNNL). The glasses were designated the high-level waste – Aspen Process Performance Simulation (HLW-APPS) study glasses. The data provided in this report is to be used in the development, validation, and implementation of enhanced property/composition models for nuclear waste glasses.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Electronic effects in radiation-induced collision cascades in nickel

The accurate treatment of electronic effects in multi-million-atom simulations of radiation-induced collision cascades is crucial for reliable predictions of primary radiation damage. In this work, we explore the fidelity of a recently developed two-temperature molecular dynamics model implementing an electron density-dependent coupling of electronic and atomic subsystems for cascade simulations in nickel. We show that the parameter-free model realistically captures the instantaneous energy losses during all stages of the highly nonequilibrium cascade process. Our simulations predict two distinct coupling regimes, corresponding to the rapid energy losses through electronic stopping in the early stages of the cascade and to the slow equilibration through the electron-phonon coupling mechanism in the later stages, without the use of separate models or coupling terms. The intermediate stage of the cascade dynamics displays a complex energy transfer between the subsystems, which cannot be interpreted by comparison to either electronic stopping or electron-phonon coupling theories. We therefore compare the predicted atomic mixing, which is sensitive to the energy losses during the intermediate cascade stage, with experimental ion beam mixing measurements. We find good agreement with the experiments, validating the coupling model for the intermediate stage of the cascade. Predictions of final defect numbers and cluster sizes are found in line with the results from conventional electronic stopping-based methods, while significantly reducing the theoretical uncertainty in the outcomes of conventional models stemming from arbitrary choices of thresholds for different coupling terms. Our results represent a notable improvement in cascade damage predictions in nickel, providing validation of the electron density-dependent coupling model for radiation damage simulations in general. The results lead us to propose an interpretation of the electronic energy losses in the intermediate regime of velocities, where we find an effectively nonlinear dissipation.

Crystal defects↗
DOE OSTI2025

Numerical simulations of three-dimensional ion crystal dynamics in a Penning trap using the fast multipole method

We simulate the dynamics, including laser cooling, of three-dimensional (3-D) ion crystals confined in a Penning trap using a newly developed molecular dynamics-like code. The numerical integration of the ions’ equations of motion is accelerated using the fast multipole method to calculate the Coulomb interaction between ions, which allows us to efficiently study large ion crystals with thousands of ions. In particular, we show that the simulation time scales linearly with ion number, rather than with the square of the ion number. By treating the ions’ absorption of photons as a Poisson process, we simulate individual photon scattering events to study laser cooling of 3-D ellipsoidal ion crystals. Initial simulations suggest that these crystals can be efficiently cooled to ultracold temperatures, aided by the mixing of the easily cooled axial motional modes with the low frequency planar modes. In our simulations of a spherical crystal of 1000 ions, the planar kinetic energy is cooled to several millikelvin in a few milliseconds while the axial kinetic energy and total potential energy are cooled even further. This suggests that 3-D ion crystals could be well suited as platforms for future quantum science experiments.

Zaris, John (ORCID:0009000196476323)↗
DOE OSTI2025

Analysis of a Disturbance Event with Inverter-Based Resources Using EMT Simulations

Increasing penetration of inverter-based resources (IBRs) necessitates newer methods of planning and analysis of disturbances. The existing phasor-domain transient stability (TS) analysis may not capture the dynamics of IBRs during fault events. Here, in this paper, electromagnetic transient (EMT) simulations using high-fidelity detailed model of power grid and one of the affected photovoltaic (PV) plants during the Angeles Forest disturbance in 2018 are performed. In these simulations, the processes to develop EMT models of power grid from traditional phasor-domain TS data and PV plant from collected data are described. Thereafter, using these simulations, the response of the PV plant during the fault event in 2018 is replicated and a sensitivity analysis is performed. The sensitivity analysis consists of making changes to the components within the PV plant and in the power grid to evaluate the impact they have on the response observed by the PV plant during the fault event. This analysis provides an understanding of the components that impact the operation of a PV plant during fault events and provide guidance to system planners on the studies that need to be performed to maintain a reliable power grid as new IBR plants are integrated.

42 ENGINEERING↗
DOE OSTI2025

2025 Workshop on Envisioning Frontiers in AI and Computing for Biological Research: Position Papers

This workshop aims to identify key research directions for transforming biology using artificial intelligence (AI), machine learning (ML) and computational methods to facilitate the discovery of new behaviors, mechanisms, and designs of biological processes relevant to DOE missions, underpinning a broader U.S. bioeconomy. By developing novel AI/ML technologies to analyze and interpret complex biological data, researchers can organize and simulate biological processes at various scales as well as advance predictive understanding and manipulation of biological systems. This integration of computation, experimentation, and next-generation experimental technologies can lead to discoveries in new biological behaviors and mechanisms relevant to DOE missions. The focus is on how advanced computational and mathematical methods can impact this mission by exploring digital twins, foundation models, automated laboratory experiments, modeling of complex living systems, and data-driven approaches for the biodesign of plants and microbial systems. While data management is important, it is not the primary focus of this workshop, which will assess the current state, trends, and AI/ML challenges at the interface between biology and computational science to identify opportunities for high-impact research at their intersection. The goal is to define research needs and opportunities that align with biological sciences, computational sciences, and applied mathematics research.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Direct Feed High-Level Waste APPS Model Glass Testing (DFHLW APPS) Matrix, Phase 2

This report summarizes the data collected during the batching and melting of a second matrix of Direct Feed High-Level Waste (DFHLW) glasses generated using the preliminary enhanced waste glass models (EWG2.5) and the Britton and Anderson (2024) preliminary DFHLW feed vector. The purpose of these glasses is two-fold: 1. Validate EWG2.5 glass calculations being used in the Aspen Process Performance Simulation (APPS) model. 2. Evaluate and ultimately improve the glass property models and formulation methods used for design of DFHLW glasses as part of an iterative process of data collection and model refinement. Some of the 16 APPS2 glasses tested did not satisfy all target property constraints due to the limited data on DFHLW glass supporting the EWG2.5 models. • One glass, APPS2-10, formed nepheline on canister centerline cooling (CCC) heat-treatment and failed the product consistency test (PCT) response limits. This glass also had high B and Cr release rates for the toxicity characteristic leaching procedure (TCLP). All other glasses were found to satisfy the PCT and TCLP constraints for both quenched and CCC samples. • One glass, APPS2-08, had higher than acceptable viscosity due to magnetite crystallization. • One glass, APPS2-09, formed greater than 2 vol% crystals at 950 °C. As the glass design criterion was that the temperature at 2 vol% crystal (T 2% ) be less than 950 °C, only one glass failed the criteria. However, this criterion is being reevaluated. Four additional glasses formed crystal fractions between 1 and 2 vol% at 950 °C (APPS2-03, -08, -12, and -14). • Four glasses – APPS2-01, -02, -04, and -16 – failed the Monofrax K-3 refractory neck corrosion (k neck ) design limit of 0.04 in. at 1208 °C for 6 d. This is another criterion being reevaluated. Four additional glasses (APPS2-05, -06, -11, and -13) exhibited 0.025 = k neck = 0.04 in. • All 16 glasses passed the sulfur solubility and TCLP constraints. The measured property values were compared to predicted values using EWG2.5 and a selection of other existing models. A few models (e.g., electrical conductivity, TCLP) were found to be adequate for designing DFHLW glasses in the near future, while others require refits or offsets. It is recommended that new property models be developed for EWG3.0, as a large amount of DFHLW glass property data (> 14 × existing data) is expected to be collected in the compositional spaces where no data was previously available. To enable near-term calculations and formulations for designing DFHLW glasses and processing rate estimations, a formulation algorithm with minor modifications will be developed, EWG2.6.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Comparative study of machine learning techniques for post-combustion carbon capture systems

Computational analysis of countercurrent flows in packed absorption columns, often used in solvent-based post-combustion carbon capture systems (CCSs), is challenging. Typically, computational fluid dynamics (CFD) approaches are used to simulate the interactions between a solvent, gas, and column's packing geometry while accounting for the thermodynamics, kinetics, heat, and mass transfer effects of the absorption process. These simulations can then be used explain a column's hydrodynamic characteristics and evaluate its CO 2 -capture efficiency. However, these approaches are computationally expensive, making it difficult to evaluate numerous designs and operating conditions to improve efficiency at industrial scales. In this work, we comprehensively explore the application of statistical ML methods, convolutional neural networks (CNNs), and graph neural networks (GNNs) to aid and accelerate the scale-up and design optimization of solvent-based post-combustion CCSs. We apply these methods to CFD datasets of countercurrent flows in absorption columns with structured packings characterized by several geometric parameters. We train models to use these parameters, inlet velocity conditions, and other model-specific representations of the column to estimate key determinants of CO 2 -capture efficiency without having to simulate additional CFD datasets. We also evaluate the impact of different input types on the accuracy and generalizability of each model. We discuss the strengths and limitations of each approach to further elucidate the role of CNNs, GNNs, and other machine learning approaches for CO 2 -capture property prediction and design optimization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Numerical Mixing Suppresses Submesoscale Baroclinic Instabilities Over Sloping Bathymetry

Abstract The impacts of spurious numerical salinity mixing on the larger‐scale flow and tracer fields are characterized using idealized simulations. The idealized model is motivated by realistic simulations of the Texas‐Louisiana shelf and features oscillatory near‐inertial wind forcing. can exceed the physical mixing from the turbulence closure in frontal zones and within the mixed layer. This suggests that simulated mixing processes in frontal zones are driven largely by . Near‐inertial alongshore wind stress amplitude is varied to identify a base case that maximizes the ratio of to in simulations with no prescribed horizontal mixing. We then test the sensitivity of the base case with three tracer advection schemes (MPDATA, U3HC4, and HSIMT) and conduct ensemble runs with perturbed bathymetry. Instability growth is evaluated using the volume‐integrated eddy kinetic energy and available potential energy . While all schemes have similar total mixing, the HSIMT simulations have over double the volume‐integrated and 20% less relative to other schemes, which suppresses the release of and reduces the by roughly 25%. This results in reduced isohaline variability and steeper isopycnals, evidence that enhanced suppresses instability growth. Differences in and between the MPDATA and U3HC4 simulations are marginal. However, the U3HC4 simulations have 25% more . Experiments with variable horizontal viscosity and diffusivity coefficients show that small amounts of prescribed horizontal mixing improve the representation of the ocean state for all advection schemes by reducing the and increasing the .

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A Versatile Simulated Data Transport Layer for in Situ Workflows Performance Evaluation

In situ processing does not only allow scientific applications to face the explosion in data volume and velocity but also to address the time constraints of many simulation-analysis workflows by providing scientists with early insights about their applications at runtime. Multiple frameworks implement the concept of a data transport layer (DTL) to enable such in situ workflows. These tools are very versatile, directly or indirectly access the data generated on the same node, another node of the same compute cluster, or a completely distinct node, and allow data publishers and subscribers to run on the same computing resources or not. This versatility puts on researchers the onus of taking key decisions related to resource allocation and how to transport data to ensure the most efficient execution of their in situ workflows. However, domain scientists and workflow practitioners lack the appropriate tools to assess the respective performance of particular design and deployment options. In this paper we introduce a versatile simulated DTL designed to provide researchers with insights on the respective performance of different execution scenarios of in situ workflows. This open-source, standalone library builds on the SimGrid toolkit and can be linked to any SimGrid-based simulator. It facilitates the evaluation of the performance behavior, at scale, of different data transport configurations and the study of the effects of resource allocation strategies. We demonstrate the scalability, versatility, and accuracy of this simulated DTL by reproducing the execution of two synthetic benchmarks and of a real-world in situ workflow composed of an MPI application and a parallel data analysis. Results of simulations run on a single core show that the proposed library can simulate the interactions of tens of thousands of simulated processes deployed on two interconnected commodity clusters in a few seconds, and the execution by a thousand simulated processes of an in situ workflow in less than three minutes.

Suter, Fred [ORNL] (ORCID:0000000319021955)↗
DOE OSTI2025

Effect of toolpath in large-format additive manufacturing with bio-derived composites

There has been growing interest in integrating bio-derived composites into Large-format additive manufacturing (LFAM) feedstocks to reduce the use of petroleum-derived materials and reduce the overall carbon footprint of LFAM. However, these materials present unique challenges during manufacturing due to their variability, which can lead to unintended deformations and failures attributable to suboptimal process conditions. While numerical modelling has been extensively employed to simulate numerous manufacturing processes, its application in LFAM with bio-derived composites remains limited. This study addresses this gap by systematically developing a numerical model to simulate the LFAM process using bio-based materials, specifically wood fiber-reinforced polylactic acid (PLA/WF). Experimental investigations were conducted to characterise the thermal and mechanical properties of additively manufactured PLA/WF specimens. Numerical simulations were performed to predict temperature profiles and deformations during LFAM. The effect of varying infill patterns, internal structures, and tool paths on the temperature distribution and deformation of printed parts was explored using the developed model. This article aims to advance the utilisation of bio-derived composites in LFAM systems and provide a comprehensive understanding of the LFAM process. The findings offer valuable insights for optimising process parameters and enhancing the performance of LFAM with bio-based composites.

Large-format additive manufacturing↗
DOE OSTI2026

Accelerating laser ray tracing in high fidelity physics simulations of laser melting using squeeze U-net

Laser melting is a core component of the ongoing industrial revolution, dubbed Industry 4.0, as lasers facilitate fast and precise melting and fusion in advanced manufacturing. There is a strong need to optimize the laser process using simulations. However, this has proven challenging as high fidelity simulations are needed for predictive modeling and this is currently prohibitively expensive even when run on hundreds of processors on high performance computers. The challenge is capturing complex physics of laser material interaction, fluid dynamics, thermal physics and material phase transformations at various length and time scales. To close this technological gap, we modified a squeeze U-net to accelerate the laser ray tracing component of such high fidelity models by ~4x–40x while preserving the core physics principle of conservation of energy with 97% accuracy. This approach enables the accurate modeling of global laser energy absorption as a function of local surface temperatures and complex surface topologies, which govern the reflection directions and energy losses of laser rays upon interacting with the material surface.

Computer science↗
DOE OSTI2026

Dynamic evolution of ultracold electron beams from nanostructured cathodes

Dynamic Coulomb expansion of dense, nonuniform particle beamlets can lead to causticlike formations in the charge density profile of the resultant beam. Such features only occur for cold beams where thermal diffusion of the particles is minimized. These phenomena have been observed experimentally for select cases with ion beams but not yet for electron beams as the electron dynamics occur over far faster timescales. We simulate the interaction of electron bunches while varying the initial transverse temperature and density profiles to determine the thresholds that characterize this pattern formation. Furthermore, we extend the simulations to geometric arrangements that are more pertinent to the nanotip cathode application and consider the effects of asymmetries and the impact of a low-density halo on the overall process. The simulations are conducted using a previously documented high-fidelity algorithm for collisional particle dynamics. These studies identify the temperature and density thresholds required for caustic formation and provide guidance for optimizing the design of nanoscale structured cathodes. We also demonstrate that the correlation between the final density patterns and the underlying beam parameters enables a virtual diagnostic tool for measuring ultralow emittance beams solely based on a transverse charge density profile.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Modeling Mycorrhizal Carbon Costs in Temperate Forests: The Impacts of Functional Diversity and Global Change Factors

Mycorrhizal fungi form symbiotic relationships with most plant species, facilitating nutrient acquisition while consuming a significant fraction of the plant's photosynthetic carbon (C), which we define as the mycorrhizal C cost. Drivers of the mycorrhizal C cost, which is crucial for predicting environmental impacts on plant productivity, remain under-explored and difficult to quantify. Ecosystem models that incorporate mycorrhizae can offer insights into mycorrhizal C cost dynamics, but their predictions have rarely been validated against empirical data. Here, in this study, we used the Myco-CORPSE model, which explicitly simulates mycorrhizal processes alongside soil carbon and nitrogen cycling, to investigate the drivers of mycorrhizal C cost in temperate forests. Applying this model to over 1,800 forest inventory plots across the eastern United States, we found that the simulations matched published data, showing higher C allocation to ectomycorrhizal (ECM) fungi (16.0% of net primary production (NPP)) compared to arbuscular mycorrhizal (AM) fungi (5.8% of NPP). Further analysis showed that mixed forests, co-dominated by both AM and ECM trees, allocated less C to mycorrhizal fungi compared to forests dominated by either AM or ECM fungi alone, due to complementary nutrient acquisition strategies. Elevated Nitrogen (N) deposition and higher temperatures reduce mycorrhizal C costs, favoring AM strategies. Conversely, elevated CO 2 (eCO 2 ) increased plant N demand and mycorrhizal C costs, favoring ECM strategies that access organic N sources. These findings underscore the critical role of mycorrhizal functional diversity in plant nutrient acquisition and C dynamics, providing new insights into how mycorrhizal symbioses respond to global change.

Shao, Siya [Dartmouth College, Hanover, NH (United↗