Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical transport”

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

CONNECT Neutronics Initial Report

This report represents the first status update from the Creation of Next-generation Nuclear Energy Computational Technology (CONNECT) effort within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program and is intended to satisfy the reporting requirements of the following milestones: • M3MS-24OR0101321: Assess needs, requirements and opportunities for high-fidelity neutronics and transport as used in DOE-NE program. • M3MS-24AN0101301: Generate a community report on the needs, requirements, and opportunities for high-fidelity neutronics and transport as used in DOE-NE programs and industry. • M3MS-24OR0202433: Implement and assess numerical strategies leveraging Monte Carlo neutron transport on GPUs for production analysis. In particular, Sections 1, 2, and 3 detail the assessment of the current state of practice and the landscape of the broader advanced computing world to meet the objectives of M3MS-24OR0101321 and M3MS-24AN0101301, while Sections 4 and 5 describe new work during FY24 to extend current GPU Monte Carlo capabilities in areas relevant to active NEAMS development to satisfy M3MS-24OR0202433. Section 6 offers some thoughts about the potential future impact of Monte Carlo neutronics on NEAMS and the nuclear industry.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Numerical investigations to identify environmental factors for field-scale reactive transport of pathogens at riverbank filtration sites

While induced bank filtration is a proven method for facilitating sustainable drinking water production, it is at risk from surface water contaminations (e.g., pathogens). Induced bank filtration and pathogen transport in groundwater have been studied extensively. However, long-term studies that consider real-world conditions are missing. These conditions include seasonal changes to environmental conditions and waterworks operations. Therefore, to analyze the effect of seasonal changes on the transport of human pathogenic viruses and their indicators in induced bank filtration, concentrations of adenoviruses and pathogen indicators were monitored over 16 months at an active bank filtration plant at the Rhine River, in Düsseldorf (Germany). Based on this data, a 2D groundwater model was created in PFLOTRAN that simulated flow, heat transport, conservative transport of chloride and the resulting electrical conductivity, reactive transport of oxygen and nitrate, and colloid-based transport of coliforms, somatic coliphages, and adenoviruses. The results show that reduced travel time was the key factor determining periods with a low removal of coliforms and somatic coliphages in the aquifer. Travel time was controlled by river level variations during rainy seasons, and the waterworks extraction rates during dry seasons. Further, for adenovirus transport, travel distance in the subsurface appeared to be the key factor, while travel time had no significant impact. Coliform removal increased when the colmation layer permeability decreased, while coliphage and adenovirus removal was unaffected by the colmation layer permeability. Seasonal changes in temperature and oxygen content did not significantly impact the removal of coliphages and adenoviruses in groundwater. Denitrifying conditions correlated with a lowered coliform removal, but the modelling could not establish a connection between denitrifying conditions and coliform removal. Our study showed that removal of pathogens and pathogen indicators at induced bank filtration plants varies greatly in time and space (e.g., for coliforms from 1 to 4 log-levels at 20m travel distance), and that adenovirus transport differs considerably from transport of coliforms and somatic coliphages.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Reactive Transport Modeling of Mineral Precipitation and Carbon Trapping in Discrete Fracture Networks

In this study we use numerical experiments to analyze reactive flow and transport behavior in discrete intersecting fracture networks, focusing on (a) how reaction-induced changes in physical and chemical properties affect flow connectivity and (b) how fracture networks developed in the Earth's critical zone contribute to carbon sequestration via mineral weathering reactions. In the first part of the study, we used two-dimensional reactive flow and transport simulations to analyze the impacts of mixing in a natural discrete fracture network. We concluded that reaction-induced changes can substantially alter the flow connectivity, especially at fracture intersections. The second set of simulations considered the problem of natural weathering of fractured mafic and ultramafic rocks in the partially saturated Earth's critical zone as a function of infiltration rates, fracture permeability, and partially saturated flow parameters. As a model system, we considered an incongruent reaction network with dissolution of forsterite and precipitation of magnesite. The behavior is complex in terms of the rate-controlling processes because of the multicomponent nature of the system as shown by the grid Peclet number: the CO 2 behavior is gas diffusion-controlled in the partially saturated zone, while the rate of water flow via the Damkӧhler number controls Mg 2+ transport through the fracture network. The amounts of carbon that can be trapped are modest, but the naturally fractured domain considered here provides a useful “base case” against which various engineered solutions can be compared.

58 GEOSCIENCES↗
DOE OSTI2024

A Novel Workflow for Coupled Simulation of Hydraulic Stimulation with Simultaneous Injection of Proppant

Enhanced Geothermal Systems (EGS) are located in impermeable host rock. Therefore, to achieve promising circulation, hydraulic stimulation is a key factor. Stimulation involves injecting a combination of fluids and proppant particles into existing fractures under high pressure. The effectiveness of stimulation depends on a variety of different factors, including the pressurization scheme and the optimal choice of carrying fluid and proppant particles to be used. For example, high-density particles can settle out of the injected slurry quickly, whereas lighter particles as well as a dense carrying fluid allow for a long transport into the reservoir. A numerical tool is therefore essential to take all these factors into account and to make robust predictions about the performance of the hydraulic stimulation. We present a novel computational application for coupled numerical simulation of hydraulic stimulation with concurrent injection of proppant in a fractured geothermal reservoir. This application is based on the open-source MOOSE framework that has been enriched with a set of equations to account for the circulation of proppant-laden slurries. It accounts for both the fluid and proppant components of the mixture, as well as particle-driven processes that include gravity settling, particle-particle interactions, and strong density and viscosity contrasts. A particular challenge is the mutual coupling with fracture mechanical processes in a propagating fracture using the FE code Defmod. Aperture and fluid pressure on one side and slurry density and viscosity changes on the other side interact and affect each other. Their basic understanding is therefore essential for the success of the stimulation. We compare our application to common numerical benchmark problems for both slurry flow and hydraulic stimulation.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Enhancement of ECCD by the current condensation effect for stabilizing large magnetic islands caused by neoclassical tearing modes in tokamak plasmas

The radio frequency current condensation effect reported in Reiman and Fisch (2018 Phys. Rev. Lett. 121 225001) is modeled in the nonlinear resistive magnetohydrodynamic code. A series of numerical investigations have been performed to investigate the enhancement of electron cyclotron current drive (ECCD) by the current condensation effect during the control of neoclassical tearing mode (NTM) in tokamak plasmas. In the numerical model, both the parallel transport and the perpendicular transport of electron temperature are considered. The EC driven current and driven perturbed electron temperature can nonlinearly evolve within the given magnetic configuration and eventually reach saturation states. The input power threshold of ECCD and the fold bifurcation phenomenon are numerically verified via nonlinear simulations. The numerical results show good agreements with the analytical results. Moreover, spatial distributions of EC current for the two solutions at different condensed level are displayed. The control effectiveness of ECCD for large NTM islands has been evaluated while considering the current condensation effect. While taking into account current condensation effect, for a sufficiently large input power, a larger island can be more effectively stabilized than a smaller one, which suggests a reassessment of the previous idea that the ECCD should always be turned on as early as possible. The potential physics mechanism behind the ECCD control have all been discussed in detail. Furthermore, the condensation effect is found to have favorable effects on the radial misalignment of ECCD. In the consideration of the situation for extremely localized control needs, a highly peaked heating profile is adopted to verify that the fold bifurcation phenomenon still exists and the current condensation effect can still take effect in this extreme condition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Progress Report on Model Development for the Transport of Aerosol Through Microchannels

Accuracy of consequence assessments could be improved by accounting for the leak path deposition of aerosol in the source term. Recent experimental measurements from Sandia (2018) provided evidence of partial microchannel plugging due to particulate deposition. This filtration effect is relevant to the source term assessments and it is not unreasonable to anticipate that aerosol deposition in more realistic geometries could even lead to complete plugging of the leak path. This report summarizes current progress (as of end of Q3 FY 2020) on the development of a phenomenological model of aerosol transport, deposition, and plugging through microchannels. The purpose is to introduce a generic, reliable numerical model for the prediction of aerosol transport, deposition, and plugging in leak paths while accounting for potential plugging formation. This report includes (1) an overview of recent additions and improvements on the recently developed numerical model to analyze the various deposition processes in leak paths, to provide quantitative estimates of penetration factors, and to gain an understanding of the variables that affect them and (2) a summary of recent benchmarking results of the model’s validity with recent experimental measurements from Sandia.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Zonally dominated dynamics and Dimits threshold in curvature-driven ITG turbulence

The saturated state of turbulence driven by the ion-temperature-gradient instability is investigated using a two-dimensional long-wavelength fluid model that describes the perturbed electrostatic potential and perturbed ion temperature in a magnetic field with constant curvature (a $Z$ -pinch) and an equilibrium temperature gradient. Numerical simulations reveal a well-defined transition between a finite-amplitude saturated state dominated by strong zonal-flow and zonal temperature perturbations, and a blow-up state that fails to saturate on a box-independent scale. We argue that this transition is equivalent to the Dimits transition from a low-transport to a high-transport state seen in gyrokinetic numerical simulations (Dimits et al. , Phys. Plasmas , vol. 7, 2000, 969). A quasi-static staircase-like structure of the temperature gradient intertwined with zonal flows, which have patch-wise constant shear, emerges near the Dimits threshold. The turbulent heat flux in the low-collisionality near-marginal state is dominated by turbulent bursts, triggered by coherent long-lived structures closely resembling those found in gyrokinetic simulations with imposed equilibrium flow shear (van Wyk et al. , J. Plasma Phys. , vol. 82, 2016, 905820609). The breakup of the low-transport Dimits regime is linked to a competition between the two different sources of poloidal momentum in the system – the Reynolds stress and the advection of the diamagnetic flow by the $\boldsymbol {E}\times \boldsymbol {B}$ flow. By analysing the linear ion-temperature-gradient modes, we obtain a semi-analytic model for the Dimits threshold at large collisionality.

Physics↗
DOE OSTI2021

Numerical Modeling of CO2 Sequestration within a Five-Spot Well Pattern in the Morrow B Sandstone of the Farnsworth Hydrocarbon Field: Comparison of the TOUGHREACT, STOMP-EOR, and GEM Simulators

The objectives of this study were (1) to assess the fate and impact of CO2 injected into the Morrow B Sandstone in the Farnsworth Unit (FWU) through numerical non-isothermal reactive transport modeling, and (2) to compare the performance of three major reactive solute transport simulators, TOUGHREACT, STOMP-EOR, and GEM, under the same input conditions. The models were based on a quarter of a five-spot well pattern where CO2 was injected on a water-alternating-gas schedule for the first 25 years of the 1000 year simulation. The reservoir pore fluid consisted of water with or without petroleum. The results of the models have numerous broad similarities, such as the pattern of reservoir cooling caused by the injected fluids, a large initial pH drop followed by gradual pH neutralization, the long-term persistence of an immiscible CO2 gas phase, the continuous dissolution of calcite, very small decreases in porosity, and the increasing importance over time of carbonate mineral CO2 sequestration. The models differed in their predicted fluid pressure evolutions; amounts of mineral precipitation and dissolution; and distribution of CO2 among immiscible gas, petroleum, formation water, and carbonate minerals. The results of the study show the usefulness of numerical simulations in identifying broad patterns of behavior associated with CO2 injection, but also point to significant uncertainties in the numerical values of many model output parameters.

STOMP↗
DOE OSTI2022

A Brief Overview of Reactive Transport Codes Used in CO 2 Applications

Reactive transport modelling consists of computational and numerical models that describe the coupled physical, chemical, mechanical, and biological processes interacting with each other over a broad range of spatial and temporal scales. This report reviews some of the reactive transport codes available in the literature for carbon dioxide (CO 2 ) applications and can assist the scientific community with the applicability of each code to solve a specific problem. Approaches implemented in these codes can be categorized in three groups: (i) continuum-scale, (ii) pore-scale, and (iii) hybrid-scale approaches. Theoretical foundations, numerical implementations, and application examples using the models are described. In the continuum approach, flow and transport are formulated in terms of a representative control volume of the medium and coupled with geochemical reactions. The medium is characterized by bulk parameters such as porosity, permeability, or reactive surface area. In the pore-scale approach, each point of space in the pore network is occupied by either a fluid or solid phase. The pore-scale approach requires an exact knowledge of the spatial and temporal phase distribution. In the hybrid-scale approach, media are described by multiple characteristic length-scales, with some regions using a pore-scale approach while others used a continuum approach.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

On the implementation of dislocation reactions in continuum dislocation dynamics modeling of mesoscale plasticity

The continuum dislocation dynamics framework for mesoscale plasticity is intended to capture the dislocation density evolution and the deformation of crystals when subjected to mechanical loading. It does so by solving a set of transport equations for dislocations concurrently with crystal mechanics equations, with the latter being cast in the form of an eigenstrain problem. Incorporating dislocation reactions in the dislocation transport equations is essential for making such continuum dislocation dynamics predictive. A formulation is proposed to incorporate dislocation reactions in the transport equations of the vector density-based continuum dislocation dynamics. This formulation aims to rigorously enforce dislocation line continuity using the concept of virtual dislocations that close all dislocation loops involved in cross slip, annihilation, and glissile and sessile junction reactions. The addition of virtual dislocations enables us to accurately enforce the divergence free condition upon the numerical solution of the dislocation transport equations for all slip systems individually. A set of tests were performed to illustrate the accuracy of the formulation and the solution of the transport equations within the vector density-based continuum dislocation dynamics. Comparing the results from these tests with an earlier approach in which the divergence free constraint was enforced on the total dislocation density tensor or the sum of two densities when only cross slip is considered shows that the new approach yields highly accurate results. Bulk simulations were performed for a face centered cubic crystal based on the new formulation and the results were compared with discrete dislocation dynamics predictions of the same. The microstructural features obtained from continuum dislocation dynamics were also analyzed with reference to relevant experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Investigation of Biomass Fouling on Screw Feeder in Preconversion of Pyrolysis

Biomass is a source of renewable energy and can undergo thermochemical processes to be converted into gaseous and liquid fuels. Thermal decomposition-induced biomass fouling on the screw feeder is a major challenge in preconversion of pyrolysis, leading to plugging of the feed line. Here, this work intends to gain fundamental understanding of the biomass fouling phenomena across the heat gradient on the screw feeder via a combination of materials characterization and thermal simulation. An actual screw feeder that failed because of biomass fouling was examined to reveal the deposit’s morphology, composition, and mechanical properties. Thermal simulations were developed for both the screw feeder and the wood particles being transported using a combined analytical and numerical approach. The simulated temperature profiles correlated well with the deposit observations, based on which potential mitigations were proposed from both the operation and screw design perspectives and supported by preliminary experimental validation.

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

Progress Report on Model Development for the Transport of Aerosol through Microchannels

This report summarizes the current progress in the development of a phenomenological model of aerosol transport, deposition, and plugging through microchannels. The purpose is to introduce a generic, reliable numerical model for the prediction of aerosol transport, deposition, and plugging in leak paths while accounting for potential plugging formation, to a user community involving researchers, regulators, and industry. In that regard, a Graphical User Interface (GUI) was generated by integrating the individual MATLAB scripts that make up the model and adding additional features to aid the general understanding of the problem to a user. This report focuses on the model development and features included. Further, predictions from the GUI are compared with experimental data for validation. The strength of the GUI is the ability of a user to plug in basic parameters such as the initial pressure conditions and canister model specifications in order to obtain a first principal approximation of vital information such as the blowdown pressure differential, aerosol penetration, and deposition as a function of time. The user can obtain this without the know-how of the underlying MATLAB scripts and thus enables the model to be readily applied by regulators, industry, and shareholders to reduce the uncertainty in the off-site radiological consequences. The report also lists future tasks under this work scope. This includes ongoing efforts to include additional crack geometries (divergent and divergent-convergent slots) to the model; GUI development and improvement based on feedback from users, and integration of aerosol source term data from Sibling Pin tests into the model as we approach realistic canister stress corrosion cracking–induced aerosol release scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Development of the Table of Initial Isolation and Protective Action Distances for the 2024 Emergency Response Guidebook

The transportation of hazardous materials creates numerous opportunities for the release of toxic substances into the environment, whether caused by traffic accidents, train derailments, equipment failures, or human error. Such releases can pose acute hazards to the general public and to emergency response personnel who are the first to arrive at the scene. To help first responders determine whether a shipment is potentially hazardous and decide what actions should be taken if a toxic spill does occur, the Emergency Response Guidebook (ERG) is published by the U.S. Department of Transportation (DOT), Transport Canada, and the Secretariat of Transport and Communications of Mexico; with contributions from Centro de Informaciòn Quìmica para Emergencias of Argentina. The most recent version is the 2024 edition of the ERG (ERG 2024), titled 2024 Emergency Response Guidebook (ERG2024). The ERG provides essential information about firefighting, spill response, and potential public health effects. For chemicals that are toxic by inhalation (TIH) and chemicals that produce TIH gases upon reaction with water (TIH by water reactivity or TIHWR), the ERG provides initial isolation distances (IIDs) and protective action distances (PADs). The IID defines the radius of the zone around the spill that should be accessed solely by people who are directly involved in emergency response. The PAD is the distance downwind of the source of the release within which persons should be either evacuated or sheltered in place, depending on the severity of the incident and the nature of the population (e.g., density, age, health).

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2025

Realizability-preserving discontinuous Galerkin method for spectral two-moment radiation transport in special relativity

Here we present a realizability-preserving numerical method for solving a spectral two-moment model to simulate the transport of massless, neutral particles interacting with a steady background material moving with relativistic velocities. The model is obtained as the special relativistic limit of a four-momentum-conservative general relativistic two-moment model. Using a maximum-entropy closure, we solve for the Eulerian-frame energy and momentum. The proposed numerical method is designed to preserve moment realizability, which corresponds to moments defined by a nonnegative phase-space density. The realizability-preserving method is achieved with the following key components: (i) a discontinuous Galerkin phase-space discretization with specially constructed numerical fluxes in the spatial and energy dimensions; (ii) a strong stability-preserving implicit-explicit time-integration method; (iii) a realizability-preserving conserved to primitive moment solver; (iv) a realizability-preserving implicit collision solver; and (v) a realizability-enforcing limiter. Component (iii) is necessitated by the closure procedure, which closes higher order moments nonlinearly in terms of primitive moments. The nonlinear conserved to primitive and the implicit collision solves are formulated as fixed-point problems, which are solved with custom iterative solvers designed to preserve the realizability of each iterate. With a series of numerical tests, we demonstrate the accuracy and robustness of this discontinuous-Galerkin-implicit-explicit method.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A low-rank power iteration scheme for neutron transport criticality problems

Computing effective eigenvalues for neutron transport often requires a fine numerical resolution. Here, the main challenge of such computations is the high memory effort of classical solvers, which limits the accuracy of chosen discretizations. In this work, we derive a method for the computation of effective eigenvalues when the underlying solution has a low-rank structure. This is accomplished by utilizing dynamical low-rank approximation (DLRA), which is an efficient strategy to derive time evolution equations for low-rank solution representations. The main idea is to interpret the iterates of the classical inverse power iteration as pseudo-time steps and apply the DLRA concepts in this framework. In our numerical experiment, we demonstrate that our method significantly reduces memory requirements while achieving the desired accuracy. Analytic investigations show that the proposed iteration scheme inherits the convergence speed of the inverse power iteration, at least for a simplified setting.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Vortex guidance and transport in channeled pinning arrays

Here, we numerically examine vortices interacting with pinning arrays where a portion of the pinning sites have been removed in order to create coexisting regions of strong and weak pinning. The region without pinning sites acts as an easy-flow channel. For driving in different directions with respect to the channel, we observe distinct types of vortex flow. When the drive is parallel to the channel, the flow first occurs in the pin free region followed by a secondary depinning transition in the pinned region. At high vortex densities there is also an intermediate plastic flow phase due to the coupling between the weak and strong pinning regions. For driving applied perpendicular to the channel, we observe a jammed phase in which vortices accumulate on the boundary of the pinned region due to the vortex-vortex repulsion, while at higher drives the vortices begin to flow through the pinning array. For driving at an angle to the channel, depending on the filling we observe a drive-induced reentrant pinning effect as well as negative differential mobility which occurs when vortices move from the unpinned to the pinned portion of the sample.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Plant SWEETs: from sugar transport to plant–pathogen interaction and more unexpected physiological roles

Abstract Sugars Will Eventually be Exported Transporters (SWEETs) have important roles in numerous physiological mechanisms where sugar efflux is critical, including phloem loading, nectar secretion, seed nutrient filling, among other less expected functions. They mediate low affinity and high capacity transport, and in angiosperms this family is composed by 20 paralogs on average. As SWEETs facilitate the efflux of sugars, they are highly susceptible to hijacking by pathogens, making them central players in plant–pathogen interaction. For instance, several species from the Xanthomonas genus are able to upregulate the transcription of SWEET transporters in rice (Oryza sativa), upon the secretion of transcription-activator-like effectors. Other pathogens, such as Botrytis cinerea or Erysiphe necator, are also capable of increasing SWEET expression. However, the opposite behavior has been observed in some cases, as overexpression of the tonoplast AtSWEET2 during Pythium irregulare infection restricted sugar availability to the pathogen, rendering plants more resistant. Therefore, a clear-cut role for SWEET transporters during plant–pathogen interactions has so far been difficult to define, as the metabolic signatures and their regulatory nodes, which decide the susceptibility or resistance responses, remain poorly understood. This fuels the still ongoing scientific question: what roles can SWEETs play during plant–pathogen interaction? Likewise, the roles of SWEET transporters in response to abiotic stresses are little understood. Here, in addition to their relevance in biotic stress, we also provide a small glimpse of SWEETs importance during plant abiotic stress, and briefly debate their importance in the particular case of grapevine (Vitis vinifera) due to its socioeconomic impact.

Breia, Richard↗
DOE OSTI2023

Model data for numerical evaluation of photosensitive tracers as a strategy for separating surface and subsurface transient storage in streams

This model-data archive pertains to a study aimed at the numerical evaluation of the photosensitive tracers as a potential strategy for separating the effects of surface and hyporheic storage zones (SSZs and HSZs, respectively). Separating effects of SSZ and HSZ are important for accurately representing stream function as HSZs and SSZs expose solutes to significantly different biogeochemical conditions. We perform numerical experiments using a multiscale reactive transport model for stream corridors implemented in ATS code, which allows for representing multiple storage zones with their respective travel time distributions and biogeochemistry. For each of the numerical experiment, we provide in this model-data archive python wrapper to drive ATS (forward_model.py), synthetic observation (synthetic_btc.py, BTC_observed.csv), forward model files including ATS input (multiscale_transport.tpl), PFLOTRAN inputs (reactions_channel.tpl, reactions_hz.tpl, reactions_sz.tpl), MCMC (mcmc_run.py), predictive uncertainty (pred_uncert.py, BTCs_simulated_day.csv, BTCs_simulated_night.csv), MCMC outputs (tracer_test_GR.npy, tracer_test_logps.npy, tracer_test_parameters.npy) and Jupyter notebook for post-processing and visualization (post-processing.ipynb). For the denitrification application, ATS input file (denitrification_multisubgrid.tpl), PFLOTRAN input files (denitrification_channel.in, denitrification_hz.in), predictive uncertainty (pred_uncert_denitrification.py, BTCs_simulated_DO.csv, BTCs_simulated_DOC.csv, BTCs_simulated_Nitrate.csv).

54 ENVIRONMENTAL SCIENCES↗