Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model validation”

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 1,639 records · Page 91

Recovering New Type of Sources by Off-Site Source Recovery Project for WIPP Disposal - 16604

OSRP uses two types of eligible packaging configurations for Cm-244 and Cf-252 disposal to WIPP: an S-100 drum and a 12'' Pipe Overpack (POC) drum. Both types of packaging are WIPP accepted configurations for Contact Handled (CH) waste, which offer different advantages and disadvantages when used in the case of the two isotopes. The S-100 is a poly-shielded type A, 55 gallon drum, containing a 6'' diameter POC. The 12'' POC drum is a type A, 55 gallon drum without the internal poly-shield and with a pipe component of 12'' diameter. The S-100 drum presents the advantage of poly-shielding which is very important for neutron sources, such as the Cf-252, or Cm-244Be, at the expense of reduced internal volume. The 12'' POC drum presents a larger pipe component accommodating many more sources, like Cm-244 alpha emitters or shielded Cf-252, for which the CH dose rate is not an issue. All the sources prepared for disposal to WIPP are required to be special form, as defined in 49 CFR 173.403. When OSRP recovers normal form sources, these can be packaged in a field closable DOT certified Special Form Capsule (SFC) model II or model III, thus becoming special form. These capsules are then packaged in the type-A containers. The difference between the two types of packaging is that the 12'' POC can hold more sources, or a mix of sources and capsules at the expense of less shielding. This becomes important when the TRU content criteria of 100 nCi/g of TRU waste has to be met for the drums that are packaged for WIPP disposal. Measured neutron dose rates for Cm-244 packaged in S-100 agree reasonably well with calculated dose rates based on the models developed with MCNP. Future work will include validating the model for Cm-244 packaged in 12'' POC. A 1:1 ratio of neutron to gamma is expected for Cm-244Be dose rates packaged in S-100 drums. For Cf-252 packaged in S-100 and 12'' POC, the contribution of Cf-250 to the total emission neutron rate had to be taken into account due to the increasing role of Cf-250 as sources age. The uncertainty in the measurement of Cf-252 neutron dose rate for S-100 drum was quite high due to a high neutron background. A larger Cf- 252 source needs to be used for measurements in the future, in order to fully validate the model developed. A better agreement for measured versus calculated rate was obtained for a Cf-252 packaged in 12'' POC. For the WIPP disposal of these sources, OSRP will use the best packaging configuration for these isotopes in order to minimize the dose rate and maximize the TRU content. The S-100 packaging will be used for Cf-252 and Cm-244Be sources, when neutron dose reduction by the poly-shielding is important. Neutron dose rates will be estimated using the new MCNP correlations.

07 ISOTOPE AND RADIATION SOURCES↗

Linking process and structure in the friction stir scribe joining of dissimilar materials: A computational approach with experimental support

We present that friction stir welding (FSW) is a popular technique to join dissimilar materials in numerous applications. The solid state nature of the process enables joining materials with strikingly different physical properties. For welds in lap configuration, an enhancement to this technology is made by introducing a short, hard insert, referred to as a cutting-scribe, at the bottom of the tool pin. The cutting-scribe induces deformation in the bottom plate which leads to the formation of mechanical interlocks or hook like structures at the interface of two materials. A thermo-mechanical computational model employing a coupled Eulerian-Lagrangian approach is developed to quantitatively capture the morphology of these interlocks during the FSW process. Simulations using this model are validated by experimental observations. In conclusion, the identified interface morphology coupled with the predicted temperature field from this process–structure model can be used to estimate the post-weld microstructure and joint strength.

36 MATERIALS SCIENCE↗

Implementation and validation of a fully implicit accumulator model in RELAP-7

This paper presents the implementation and validation of an accumulator model in RELAP-7 under the framework of preconditioned Jacobian free Newton Krylov (JFNK) method, based on the similar model used in RELAP5. RELAP-7 is a new nuclear reactor system safety analysis code being developed at the Idaho National Laboratory (INL). RELAP-7 is a fully implicit system code. The JFNK and preconditioning methods used in RELAP-7 is briefly discussed. The slightly modified accumulator model is summarized for completeness. The implemented model was validated with LOFT L3-1 test and benchmarked with RELAP5 results. RELAP-7 and RELAP5 had almost identical results for the accumulator gas pressure and water level, although there were some minor difference in other parameters such as accumulator gas temperature and tank wall temperature. One advantage of the JFNK method is its easiness to maintain and modify models due to fully separation of numerical methods from physical models. It would be straightforward to extend the current RELAP-7 accumulator model to simulate the advanced accumulator design. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

GALVANICALLY GRADED INTERFACE: A COMPUTATIONAL MODEL FOR MITIGATING GALVANIC CORROSION BETWEEN MAGNESIUM AND MILD STEEL

The impact of a graded metallic spacer on the galvanic corrosion between magnesium and mild steel is investigated in this work using a COMSOL model based on a validated numerical model. A graded spacer of 4 mm thickness decreases the peak galvanic corrosion by 50% over a system without a spacer, and 37% over a system with an aluminum spacer. The impact of the electrochemical properties of spacer materials is also investigated.

Spies, Kurt A.↗

Optimal management of electric hotel loads in mild hybrid heavy duty truck

Here, the problem of engine idling for heavy-duty trucks has been under study for decades with Auxiliary Power Units (APUs) and Truck Stop Electrification (TSE) as the most compelling solutions. With the electrification of trucks approaching feasibility in terms of cost-effective technology, hybridization offers another “degree of freedom” to tackle the problem. Here, this work aims at exploiting a battery pack of a 48 V mild-hybrid heavy-duty truck to store sufficient onboard energy for powering the auxiliary loads during the hoteling. This problem is not trivial, as the battery packs typically cannot recover the entire energy required through regeneration alone; hence an optimal energy management strategy needs to be employed to charge the battery through the engine during drive operation. This strategy optimizes powertrain performances among the four modes: (i) Engine of Coasting (EOC), (ii) Regeneration by braking, (iii) Regeneration by engine, and (iv) engine idling. This paper presents the development of a Dynamic Programming (DP) framework that employs a multi-objective cost function to minimize the fuel consumption and maximize the regeneration using the above-mentioned four modes. A typical heavy-duty truck drive cycle is used to represent the drive phase, with mandatory hoteling stops as per regulations. A comprehensive powertrain model is developed using validated components’ model. The DP employs two state variables: battery State-of-Charge (SOC) and engine mode, and three control inputs: (i) the engine ON–OFF state, (ii) clutch engagement state, and (iii) power request at the Electric Machine (EM) for calculating optimal SOC trajectory. The framework also tackles rapid engine ON–OFF scenarios to avoid the challenges associated with DP and the compromises in fuel cost with those approaches. Finally, the effectiveness of the proposed framework is tested for potential fuel savings on two different battery packs by performing the full cycle simulations. The results show 6.47% of fuel consumption reduction as compared to traditional APU-based heavy-duty truck.

33 ADVANCED PROPULSION SYSTEMS↗

Systematic network assessment of the carcinogenic activities of cadmium

Cadmium has been defined as type I carcinogen for humans, but the underlying mechanisms of its carcinogenic activity and its influence on protein-protein interactions in cells are not fully elucidated. The aim of the current study was to evaluate, systematically, the carcinogenic activity of cadmium with systems biology approaches. From a literature search of 209 studies that performed with cellular models, 208 proteins influenced by cadmium exposure were identified. All of these were assessed by Western blotting and were recognized as key nodes in network analyses. The protein-protein functional interaction networks were constructed with NetBox software and visualized with Cytoscape software. These cadmium-rewired genes were used to construct a scale-free, highly connected biological protein interaction network with 850 nodes and 8770 edges. Of the network, nine key modules were identified and 60 key signaling pathways, including the estrogen, RAS, PI3K-Akt, NF-κB, HIF-1α, Jak-STAT, and TGF-β signaling pathways, were significantly enriched. With breast cancer, colorectal and prostate cancer cellular models, we validated the key node genes in the network that had been previously reported or inferred form the network by Western blotting methods, including STAT3, JNK, p38, SMAD2/3, P65, AKT1, and HIF-1α. These results suggested the established network was robust and provided a systematic view of the carcinogenic activities of cadmium in human. - Highlights: • A cadmium-influenced network with 850 nodes and 8770 edges was established. • The cadmium-rewired gene network was scale-free and highly connected. • Nine modules were identified, and 60 key signaling pathways related to cadmium-induced carcinogenesis were found. • Key mediators in the network were validated in multiple cellular models.

60 APPLIED LIFE SCIENCES↗

Modelling piloted ignition of wood and plastics

Highlights: Black-Right-Pointing-Pointer We model piloted ignition times of wood and plastics. Black-Right-Pointing-Pointer The model is applied on a packed bed. Black-Right-Pointing-Pointer When the air flow is above a critical level, no ignition can take place. - Abstract: To gain insight in the startup of an incinerator, this article deals with piloted ignition. A newly developed model is described to predict the piloted ignition times of wood, PMMA and PVC. The model is based on the lower flammability limit and the adiabatic flame temperature at this limit. The incoming radiative heat flux, sample thickness and moisture content are some of the used variables. Not only the ignition time can be calculated with the model, but also the mass flux and surface temperature at ignition. The ignition times for softwoods and PMMA are mainly under-predicted. For hardwoods and PVC the predicted ignition times agree well with experimental results. Due to a significant scatter in the experimental data the mass flux and surface temperature calculated with the model are hard to validate. The model is applied on the startup of a municipal waste incineration plant. For this process a maximum allowable primary air flow is derived. When the primary air flow is above this maximum air flow, no ignition can be obtained.

42 ENGINEERING↗

Electrical flicker-noise generated by filling and emptying of impurity states in injectors of quantum-cascade lasers

Free running line-widths (>100 kHz), much broader than intrinsic line-widths ∼100 Hz, of existing quantum-cascade lasers are governed by strong flicker frequency-noise originating from electrical flicker noise. Understanding of microscopic origins of the electrical flicker noises in quantum-cascade lasers is crucially important for the reduction of strength of flicker frequency-noise without assistances of any type of feedback schemes. In this article, an ad hoc model that is based on fluctuating charge-dipoles induced by electron trappings and de-trappings at indispensable impurity states in injector super-lattices of a quantum-cascade laser is proposed, developing theoretical framework based on the model. The validity of the present model is evaluated by comparing theoretical voltage-noise power spectral densities based on the model with experimental ones obtained by using mid-infrared quantum-cascade lasers with designed impurity-positioning. The obtained experimental results on flicker noises, in comparison with the theoretical ones, shed light on physical mechanisms, such as the inherent one due to impurity states in their injectors and extrinsic ones due to surface states on the ridge-walls and due to residual deep traps, for electrical flicker-noise generation in existing mid-infrared quantum-cascade lasers. It is shown theoretically that quasi-delta doping of impurities in their injectors leads to strong suppression of electrical flicker noise by minimization of the dipole length at a certain temperature, for instance ∼300 K and, in turn, is expected to result in substantial narrowing of the free running line-width down below 10 kHz.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Development and Experimental Validation of a Heat Transfer Model for Spilled Molten Salt Pools

A spill of radionuclide-bearing molten salt is one of the major postulated events that needs to be analyzed for liquid fluorine salt-cooled high-temperature reactor (FHR) or molten salt reactor licensing purposes. In this postulated event, radioactive source term materials (RSTMs) in the molten salt are discharged from the reactor vessel to the reactor building. The release of RSTMs from the spilled salt pool to the gas space in the reactor building is expected to be controlled by the cooling behavior of the spilled salt, including the growth and shrinkage of the solid crust on the surface of the spilled salt pool. This paper presents a simulation model for spilled salt pool heat transfer and validation efforts. The validation data come from two molten salt spill tests that were performed recently: the PELE2 test by the Rapid Experimental Laboratory of Kairos Power LLC (KP) and the Argonne salt cooling test conducted by Argonne National Laboratory. The former was a large-scale test involving kilograms of molten spilled FLiNaK salt, and the latter was a relatively smaller-scale test targeting various processes associated with a salt spill event. Both tests generated valuable data sets that can be used to assess salt cooling and validate evaluation models. This paper provides a new one-dimensional model that can simulate the cooling process of a spilled salt pool as well as the thermal responses of heat structures, such as the stainless steel liner and the concrete below the salt. The model has been implemented as part of KP-SAM code, which is a branch of the systems code SAM specific to KP FHR. In conclusion, the simulation results of the model are compared with the data of the PELE2 and Argonne tests, and reasonable agreements are observed between the model and test data.

heat transfer model↗

Resource Assessment for Distributed Wind Energy: An Evaluation of Best-Practice Methods in the Continental US

Current wind resources within the United States (US) indicate a potential to profitably install nearly 1,400 gigawatts of distributed wind (DW) capacity. This amount is equivalent to over half of the United States’ current energy demand from electricity, making it enough to power millions of homes and businesses and replace countless fossil fuel-based generating plants. Despite the potential growth of DW in the US, deployments are presently hindered by a lack of confidence in resource estimation methods. One potential challenge is that smaller-scale turbines, with hub heights of 40 meters or less, are disproportionately impacted by obstacles such as buildings and vegetation. These obstacles may produce complex wake effects, best modeled with high-fidelity complex fluid dynamics (CFD) models that are too computationally expensive to use for routine siting and resource assessment. Thus, installers today make use of heuristics and simple equations to approximate the impact of obstacles while also leveraging long-term resource data from commercial or publicly available atmospheric models. This study evaluates these historical and commonly used methods alongside new lower-order obstacle models produced from CFD simulations and measurement-based bias correction. The preliminary results from this study show the importance of taking care in the choice and application of mesoscale atmospheric models and the significant value of bias correction using measurements from nearby meteorological towers. Detailed obstacle modeling provides only modest additional gains in performance and, in some cases, can add error, especially at sites where turbines have already been located to avoid obvious impact from upwind obstacles. These findings reinforce the importance of collecting in situ measurements and suggest that obstacle models may be better applied in practice to automated or computer-aided siting, rather than in economic wind resource assessments.

17 WIND ENERGY↗

Hyper Suprime-Cam Year 3 results: Cosmology from galaxy clustering and weak lensing with HSC and SDSS using the minimal bias model

Here, we present cosmological parameter constraints from a blind joint analysis of three two-point correlation functions measured from the Year 3 Hyper Suprime-Cam (HSC-Y3) imaging data, covering about 416 deg 2 , and the SDSS DR11 spectroscopic galaxies spanning the redshift range [0.15, 0.70]. We subdivide the SDSS galaxies into three luminosity-cut, and therefore nearly volume-limited samples separated in redshift, each of which acts as a large-scale structure tracer characterized by the measurement of the projected correlation function, w p ⁡(R). We also use the measurements of the galaxy-galaxy weak-lensing signal Δ⁢Σ⁡(R) for each of these SDSS samples which act as lenses for a secure sample of source galaxies selected from the HSC-Y3 shape catalog based on their photometric redshifts. We combine these measurements with the cosmic shear correlation functions, ξ ± ⁡($\vartheta$) measured for our HSC source sample. We model these observables with the minimal bias model of the galaxy clustering observables in the context of a flat Λ⁢CDM cosmology. We use conservative scale cuts, R >12 and 8⁢ℎ -1 Mpc for Δ⁢Σ and w p , respectively, where the minimal bias model is valid, in addition to conservative prior on the residual bias in the mean redshift of the HSC photometric source galaxies. We present various validation tests of our model as well as analysis methods. Our baseline analysis yields S 8 = $0.77⁢5^{+0.043}_{-0.038}$ (68% C.I.) for the Λ⁢CDM model, after marginalizing over uncertainties in other parameters. Our value of S 8 is consistent with that from the Planck 2018 data, but the credible interval of our result is still relatively large. We show that various internal consistency tests based on different splits of the data are passed. Our results are statistically consistent with those of a companion paper, which extends this analysis to smaller scales with an emulator-based halo model, using Δ⁢Σ⁡(R) and w p ⁡(R) down to R >3 and 2⁢ℎ -1 Mpc, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

Development and Optimization of Voltammetric Methods for Real Time Analysis of Electrorefiner Salt with High Concentrations of Actinides and Fission Products

This project addresses the problem of achieving accurate material control and accountability (MC&A) around pyroprocessing electrorefiner systems. Spent nuclear fuel pyroprocessing poses a unique challenge with respect to reprocessing technology in that the fuel is never fully dissolved in the process fluid. In this case, the process fluid is molten, anhydrous LiCl-KCl salt. Therefore, there is no traditional input accountability tank. However, electrorefiners (ER) accumulate very large quantities of fissile nuclear material (including plutonium) and should be well safeguarded in a commercial facility. Idaho National Laboratory (INL) currently operates a pyroprocessing facility for treatment of spent fuel from Experimental Breeder Reactor-II with two such ER systems. INL implements MC&A via a mass tracking model in combination with periodic sampling of the salt and other materials followed by destructive analysis. This approach is projected to be insufficient to meet international safeguards timeliness requirements. A real time or near real time monitoring method is, thus, direly needed to support commercialization of pyroprocessing. A variety of approaches to achieving real time monitoring for ER salt have been proposed and studied to date—including a potentiometric actinide sensor for concentration measurements, a double bubbler for salt depth and density measurements, and laser induced breakdown spectroscopy (LIBS) for concentration measurements. While each of these methods shows some promise, each also involves substantial technical complexity that may ultimately limit their implementation. Yet another alternative is voltammetry —a very simple method in theory that has previously been tested for this application to a limited extent. The equipment for a voltammetry system consists of off-the-shelf components (three electrodes and a potentiostat), which results in substantial benefits relative to cost and robustness. Based on prior knowledge of electrochemical reduction potentials for each of the species of interest, voltammetry can be used to quantify concentrations of a variety of elemental species—including uranium, plutonium, minor actinides, and rare earths. Various methods have been tested by other researchers to date—including cyclic voltammetry, square wave voltammetry, normal pulse voltammetry, etc. In most cases, it has been observed that there is a very limited concentration range for which the output can be readily correlated with concentration in the salt. Furthermore, testing to date has been limited to simple ternary salts with only a single element being quantified. While incomplete for application to MC&A for pyroprocessing, these results lead us to believe that voltammetry can be optimized based on salt properties and fundamental electrochemical rate processes to yield a highly accurate and robust method. This project is divided into four tasks jointly executed by three university research groups. This includes experimental measurement of key physical data on the systems of interest, development of a predictive voltammetry model, experimental validation of the voltammetry model, and design/verification of an optimized measurement method. This project supports the goals of the US-ROK Joint Fuel Cycle Study in addition to the NA-24 Office of the National Nuclear Security Agency and the International Atomic Energy Agency (IAEA).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

University of Tennessee, Knoxville (UTK) contribution to: Deciphering the role of mixed-material deposition and temperature on lithium-coated PFCs in NSTX-U high-performance plasmas: Collaborative UIUC & UTK Proposal (Final Report)

This report describes the final product of model development to decipher the role of mixed material deposition and temperature on lithium-coated PFCs in NSTX-U. The computational modeling approach used to address these specific questions involved a pragmatic multiscale materials modeling framework, with the use of reaction-diffusion based continuum simulations, through modifications to the SDTrimSP code. This pragmatic modeling was supplemented with a few, select parametric ab-initio and molecular dynamics simulations, to investigate the effect of lithium and mixed lithium-carbon-refractory metal chemistry and surface morphology and chemistry on erosion yields and deuterium uptake, including the role of temperature and boronization. The main progress of the project was to develop the major framework and structure of the model, although validated application of the model was limited as a result of the unplanned shutdown and outage of NSTX-U.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Edge and Scrape-off-Layer Physics for NSTX-U (Final Report)

U.S. Department of Energy (DOE) supported work at Lodestar Research Corporation under grant DE-FG02-02ER54678 is summarized in this final report. The research conducted under this grant addressed key scientific issues associated with the plasma confinement and fusion energy experimental devices NSTX and its upgrade NSTX-U located at Princeton Plasma Physics Laboratory. Harnessing the energy that is released in fusion reactions would provide a safe and abundant source of power to meet the growing energy needs of the world population. Lodestar investigated the important cm-scale region near the edge of these devices where hot plasma, close to one million degrees, meets material surfaces in what is known as a scrape-off layer (SOL). This critical region is subject to plasma turbulence and plasma flows which must be understood to optimize fusion performance on future machines. The work employed theoretical and theory-based computational models for validation and applied the models to understand experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Building Initial Dynamic System Models for Digital Twins of the Cryogenic Moderator System at the ORNL Spallation Neutron Source

This work describes the initial development of dynamic system models of the cryogenic moderator system (CMS) of the Spallation Neutron Source (SNS) at ORNL as a part of the ORNL LDRD funded project Building TRANSFORM to Accelerate Digital Twin Applications for Nuclear Systems, LOIS 10563. The goal of the work is to start the dynamic system modeling effort with the end goal of using them for real-time applications as digital twins. The CMS is a cryogenic liquid hydrogen flow loop that provides moderation of the neutrons that are generated by the SNS. For optimal neutron production, the CMS needs to maintain a steady and controlled density of cryogenic hydrogen in the moderator section thus requiring precise temperature and pressure control. Due to the varied time scales and system characteristics, control of the system is complex, and diagnostics are also difficult. Difficulty in accessing the flow loop during operations, limited instrumentation and unknown design details of the equipment combine to make the case for having sophisticated digital twin models of the system. Operationally the CMS also provides a strong use case for digital twins due to the constant need of optimization and for troubleshooting/diagnostics. The large amount of data collected which are freely available for using in building the model and verifying and validating the model also makes it a great candidate for a proof-of-concept for digital twins. The project extends ORNL's capacity of development and implementation of the open-source dynamic system modeling tool TRANSFORM for engineering design and digital twin/real-time applications. Specific system configuration data for the CMS have been gathered and an initial dynamic model was created in the TRANSFORM library using Dymola as the solution platform. Models of increasing complexity are created to demonstrate the need for a multi-layered approach in digital twin modeling depending on the scale and phenomena being focused on. The dynamic modeling is shown to bring the dynamic operational aspects to the design process for systems as well as serve as a digital twin to the hardware and allow for models to be tuned and compared against real time operational data. These aims should help to push forward strategic goals of application of digital twins and increase the impact of ORNL systems modeling capabilities with TRANSFORM/Modelica for various advanced energy systems.

42 ENGINEERING↗

First Principles Model of Electric Cable Braid Penetration with Dielectrics

In this study, we report the formulation to account for dielectrics in a first principles multipole-based cable braid electromagnetic penetration model. To validate our first principles model, we consider a one-dimensional array of wires, which can be modeled analytically with a multipole-conformal mapping expansion for the wire charges; however, the first principles model can be readily applied to realistic cable geometries. We compare the elastance (i.e. the inverse of the capacitance) results from the first principles cable braid electromagnetic penetration model to those obtained using the analytical model. The results are found in good agreement up to a radius to half spacing ratio of 0.5-0.6, depending on the permittivity of the dielectric used, within the characteristics of many commercial cables. We observe that for typical relative permittivities encountered in braided cables, the transfer elastance values are essentially the same as those of free space; the self-elastance values are also approximated by the free space solution as long as the dielectric discontinuity is taken into account for the planar mode.

42 ENGINEERING↗

Towards Hyper-Dimensional Variography Using the Product-Sum Covariance Model

Modeling hyper-dimensional spatial variability is a complex task from both practical and theoretical standpoints. In this paper we develop a method for modeling hyper-dimensional covariance (variogram) structures using the product-sum covariance model initially developed to model spatio-temporal variability. We show that the product-sum model can be used recursively up to an arbitrarily large number of dimensions while preserving relative modeling simplicity and yielding valid covariance models. The method can be used to model variability in anisotropic conditions with multiple axes of anisotropy or when temporal evolution is involved, and thus is applicable to “full anisotropic 3D+time” situations often encountered in environmental sciences. It requires fewer assumptions than the traditional product-sum modeling approach. The new method also presents an alternative to classical approaches to modeling zonal anisotropy and requires fewer parameters to be estimated from data. We present an example by applying the method in conjunction with ordinary kriging to map photosynthetically-active radiation (PAR) for 2006, in Oklahoma, CA, USA and to explore effects of spatio-temporal variability in PAR on gross primary productivity (GPP).

54 ENVIRONMENTAL SCIENCES↗