Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal system modeling”

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 505 records · Page 28

Modelling Nuclear Thermal Propulsion Reactor Expander Cycle Startup Transients

As the interest in nuclear thermal propulsion grows, the need for high fidelity and comprehensive models of the system increases. An expander cycle model was developed using RELAP-7 for the thermal-hydraulic performance and Griffin to provide the power level using a point kinetics approximation, for a reference NTP system. The tank and pump are approximated using boundary conditions. Flow is split between the moderator cooling channels flow path and the regenerative cooling/reflector cooling channels flow path. The flow recombines before entering the turbine and the fuel cooling channels. The turbine pressure ratio is maintained at 1.44 using a PID controller. A startup transient was simulated and shows the capability of the model to run from low temperature and pressures to nominal conditions.

33 ADVANCED PROPULSION SYSTEMS↗

Modelling Nuclear Thermal Propulsion Reactor Expander Cycle Startup Transients

As the interest in nuclear thermal propulsion grows, the need for high fidelity and comprehensive models of the system increases. An expander cycle model was developed using RELAP-7 for the thermal-hydraulic performance and Griffin to provide the power level using a point kinetics approximation, for a reference NTP system. The tank and pump are approximated using boundary conditions. Flow is split between the moderator cooling channels flow path and the regenerative cooling/reflector cooling channels flow path. The flow recombines before entering the turbine and the fuel cooling channels. The turbine pressure ratio is maintained at 1.44 using a PID controller. A startup transient was simulated and shows the capability of the model to run from low temperature and pressures to nominal conditions.

33 ADVANCED PROPULSION SYSTEMS↗

An Overview of NASA Efforts on Zero Boiloff Storage of Cryogenic Propellants

Future mission planning within NASA has increasingly motivated consideration of cryogenic propellant storage durations on the order of years as opposed to a few weeks or months. Furthermore, the advancement of cryocooler and passive insulation technologies in recent years has substantially improved the prospects for zero boiloff storage of cryogenics. Accordingly, a cooperative effort by NASA's Ames Research Center (ARC), Glenn Research Center (GRC), and Marshall Space Flight Center (MSFC) has been implemented to develop and demonstrate "zero boiloff" concepts for in-space storage of cryogenic propellants, particularly liquid hydrogen and oxygen. ARC is leading the development of flight-type cryocoolers, GRC the subsystem development and small scale testing, and MSFC the large scale and integrated system level testing. Thermal and fluid modeling involves a combined effort by the three Centers. Recent accomplishments include: 1) development of "zero boiloff" analytical modeling techniques for sizing the storage tankage, passive insulation, cryocooler, power source mass, and radiators; 2) an early subscale demonstration with liquid hydrogen 3) procurement of a flight-type 10 watt, 95 K pulse tube cryocooler for liquid oxygen storage and 4) assembly of a large-scale test article for an early demonstration of the integrated operation of passive insulation, destratification/pressure control, and cryocooler (commercial unit) subsystems to achieve zero boiloff storage of liquid hydrogen. Near term plans include the large-scale integrated system demonstration testing this summer, subsystem testing of the flight-type pulse-tube cryocooler with liquid nitrogen (oxygen simulant), and continued development of a flight-type liquid hydrogen pulse tube cryocooler.

Hastings, Leon J.↗

Asymptotic Giant Branch stars as a source of short-lived radioactive nuclei in the solar nebula

We carried out a theoretical evaluation of the contribution of Asymptotic Giant Branch (AGB) stars to some short-lived (10(exp 6) less than or equal to Tau-bar less than or equal to 2 x 10(exp 7) yr) isotopes in the Interstellar Medium (ISM) and in the early solar system using stellar model calculations for thermally pulsing evolutionary phases of low-mass stars. The yields of s-process nuclei in the convective He-shell for different neutron exposures tau(sub 0) were obtained, and AGB stars were shown to produce several radioactive nuclei (especially Pd-107, Pb-205, Fe-60, Zr-93, Tc-99, Cs-135, and Hf-182) in diferent amounts. Assuming either contamination of the solar nebula from a single AGB star or models for continuous injection and mixing from many stars into the ISM, we calculate the ratios of radioactive to stable nuclei at the epoch of the Sun's formation. The dilution factor between the AGB ejecta and the early solar system matter is obtained by matching the observed Pd-107/Pd-108 and depends on the value of tau(sub 0). It is found that small masses M(sub He) of He-shell material (10(exp -4)-10(exp -7) solar mass) enriched in s-process nuclei are sufficient to contaminate 1 solar mass of the ISM to produce the Pd-107 found in the early solar system. Predictions are made for all of the other radioactive isotopes. The optimal model to explain several observed radioactive species at different states of the proto-solar nebula involves a single AGB star with a low neutron exposure (tau(sub 0) = 0.03 mbarn(sup -1)) which contaminated the cloud with a dilution factor of M(sub He)/solar mass approximately 1.5 x 10(exp -4). This will also contribute newly synthesized stable s-process nuclei in the amount of approximately 10(exp -4) of their abundances already present in the proto-solar cloud. Variations in the degree of homogenization (approximately 30%) of the injected material may account for some of the small general isotopic anomalies found in meteorites. It is also found that Fe-60 is produced in small but significant quantities that may be sufficient to explain the observations if the time elapsed delta from the contamination of the ISM to the formation of protoplanetary bodies is not higher than delta = 5 x 10(exp 6) yr. If delta is longer, up to 10 x 10(exp 6) yr, this would require the single AGB star to experience enhanced neutron densities (n(sub n) approximately 3 x 10(exp 9)n/cu cm) in the s-processing zone in order to compensate for the branching at Fe-59. The alternative model of long-term continuous ejection of matter from many AGB stars does not appear to match the observations. We also estimate the Al-26 production from the H-shell and find that the Al-26 abundance in the early solar system may be readily explained in a self-consistent manner. Moreover, Al-26 from AGB stars may contribute substantially to the galactic Al-26 gamma-source, while no significant gamma-flux from Co-60 (deriving from Fe-60 decay) is to be expected.

Wasserburg, G. J.↗

Thermomechanics coupling to Monte Carlo particle transport on unstructured mesh geometries using Cardinal

Geometry deformation due to thermal expansion influences neutron transport in many systems. Studying this phenomenon involves coupling models for neutronics, thermal hydraulics, and solid mechanics. To enable high fidelity modeling of these coupled physics, new capabilities were introduced in Cardinal, coupling OpenMC Monte Carlo particle transport models with MOOSE thermomechanical physics on unstructured moving-mesh geometries. In this work, we present a fully open-source capability leveraging on-the-fly mesh skinning to automatically regenerate OpenMC geometry, which allows multiphysics feedback from temperature, density, and geometry changes. The new capability is verified using an analytic benchmark slab problem, which couples S 2 neutron transport with thermal conduction, convective boundary conditions, Doppler-broadened cross sections, and nonlinear thermal expansion effects along the heated slab. Cardinal reproduces the analytic solutions for the neutron flux, heating, k eff , and temperature with demonstrated convergence in various error terms including mesh resolution and cross section temperature library spacing. For the nominal benchmark conditions and with a fine mesh, maximum relative errors for neutron flux, temperature, and heating are lower than 1%, while errors in integral quantities such as k eff and slab length are within 1 pcm and 48 µm, respectively. This work (i) presents a new numerical approach to thermomechanics coupling with OpenMC models, (ii) is the first (to our knowledge) to utilize a mechanical partial differential equation (PDE) solution to solve the (Griesheimer and Kooreman, 2022) analytic benchmark, and (iii) develops this verified capability within an open-source package.

97 - MATHEMATICS AND COMPUTING↗

Correlation of predicted and measured thermal stresses on an advanced aircraft structure with similar materials

A laboratory heating test simulating hypersonic heating was conducted on a heat-sink type structure to provide basic thermal stress measurements. Six NASTRAN models utilizing various combinations of bar, shear panel, membrane, and plate elements were used to develop calculated thermal stresses. Thermal stresses were also calculated using a beam model. For a given temperature distribution there was very little variation in NASTRAN calculated thermal stresses when element types were interchanged for a given grid system. Thermal stresses calculated for the beam model compared similarly to the values obtained for the NASTRAN models. Calculated thermal stresses compared generally well to laboratory measured thermal stresses. A discrepancy of signifiance occurred between the measured and predicted thermal stresses in the skin areas. A minor anomaly in the laboratory skin heating uniformity resulted in inadequate temperature input data for the structural models.

Jenkins, J. M.↗

Validation Database Based Thermal Analysis of an Advanced RPS Concept

Advanced RPS concepts can be conceived, designed and assessed using high-end computational analysis tools. These predictions may provide an initial insight into the potential performance of these models, but verification and validation are necessary and required steps to gain confidence in the numerical analysis results. This paper discusses the findings from a numerical validation exercise for a small advanced RPS concept, based on a thermal analysis methodology developed at JPL and on a validation database obtained from experiments performed at Oregon State University. Both the numerical and experimental configurations utilized a single GPHS module enabled design, resembling a Mod-RTG concept. The analysis focused on operating and environmental conditions during the storage phase only. This validation exercise helped to refine key thermal analysis and modeling parameters, such as heat transfer coefficients, and conductivity and radiation heat transfer values. Improved understanding of the Mod-RTG concept through validation of the thermal model allows for future improvements to this power system concept.

numerical model↗

Component Selection, Accelerated Testing, and Improved Modeling of AMTEC Systems for Space Power

Alkali metal thermal to electric converter (AMTEC) designs for space power are numerous, but the selection of materials for construction of long-lived AMTEC devices has been limited to electrodes, current collectors, and the solid electrolyte. AMTEC devices with lifetimes greater than 5 years require careful selection and life testing of all hot-side components. The likely selection of a remote condensed design for initial flight test and probable use with a radioisotope heat source in AMTEC powered planet probes requires the device to be constructed to tolerate operating T greater than 1150K, as well as exposure to Na (g) , and Na (liq) on the high pressure side. The temperatures involved make the characterization of high strength and chemical resistance to Na containing Na 2 O critical. Selection among materials which can be worked should not be driven by ease of fabricablity, as high temperature stability is the critical issue. These concepts drive the selecton of Mo alloys for Na (liq) containment in AMTEC cells for T to 1150K operation, as they are significantly stronger than comparable Nb or Ta alloys, are less soluble in Na (liq) containing dissolved Na 2 O, are workable compared with W alloys (which might be used for certain components), and are ductile at the T greater than 500K of proposed AMTEC modules in space applications.

AMTEC↗

Engineering-Scale Integrated Energy System Data Projection Demonstration via the Dynamic Energy Transport and Integration Laboratory

The objective of this study is to demonstrate and validate the Dynamic Energy Transport and Integration Laboratory (DETAIL) preliminary scaling analysis using Modelica language system-code Dymola. The DETAIL preliminary scaling analysis includes a multisystem integral scaling package between thermal-storage and hydrogen-electrolysis systems. To construct the system of scaled equations, dynamical system scaling (DSS) was applied to all governing laws and closure relations associated with the selected integral system. The existing Dymola thermal-energy distribution system (TEDS) facility and high-temperature steam electrolysis (HTSE) facility models in the Idaho National Laboratory HYBRID repository were used to simulate a test case and a corresponding scaled case for integrated system HYBRID demonstration and validation. The DSS projected data based on the test-case simulations and determined scaling ratios were generated and compared with scaled case simulations. The preliminary scaling analysis performance was evaluated, and scaling distortions were investigated based on data magnitude, sequence, and similarity. The results indicated a necessity to change the normalization method for thermal storage generating optimal operating conditions of 261 kW power and mass flow rate of 6.42 kg/s and the possibility of reselecting governing laws for hydrogen electrolysis to improve scaling predictive properties. To enhance system-scaling similarity for TEDS and HTSE, the requirement for scaling validation via physical-facility demonstration was identified.

08 HYDROGEN↗

Technology Development and Infusion by NASA's Entry Systems Modeling Project

This paper describes recent development of modeling and simulation technologies for entry systems and their infusion into NASA's exploration missions. Technology development is organized and prioritized using a system-level perspective, resulting in four broad technical areas of investment: (1) Thermal protection material modeling, (2) Shock layer kinetics and radiation, (3) Computational and experimental aerosciences, and (4) Guidance, navigation, and control. The paper will highlight key contributions from each of these areas, their impacts from a spacecraft and mission design perspective, and discuss planned future investment. Aspects of each technical area are only briefly summarized here. Thermal protection material modeling is geared toward high-fidelity, predictive models capable of optimizing design performance, post-flight reconstruction, and quantifying thermal protection system reliability. New computational tools and experimental techniques have been applied to Orion, MSL/Mars 2020, Mars InSight, and Mars Sample Return missions. Research and development in the area of shock layer kinetics has focused on air and CO2-based atmospheres. In both cases, substantial improvements in model uncertainty have directly impacted the development of mission margin policies, flight instrumentation design and analysis (Orion and Mars 2020), and have even revealed the importance of neglected phenomena like mid-wave infrared radiation of CO2. Aerosciences is a very broad area of interest in entry systems, yet a number of important challenges are being addressed: Coupled fluid-structure simulations of parachute inflation and dynamics affecting Orion, Commercial Crew, and Mars programs; Experimental and computational studies of vehicle dynamics; Multi-phase flow with dust particles to simulate augmentation of aerothermal environments at Mars during dust storms; and studies of roughness-induced heating augmentation relevant to tiled (Orion, Mars 2020) and woven (Mars Sample Return) thermal protection systems. Guidance and control in the context of entry systems has focused on development of methods for multi-axis control (i.e. pitch and yaw, rather than bank angle alone) of spacecraft during entry and descent, with precision landing requirements driven by Mars human exploration goals.

Barnhardt, Michael D.↗

Renewable Thermal Energy Systems: Characterization of the Most Important Thermal Energy Applications in Buildings and Industry (Report 1)

This report is the first in a three-report series that evaluates the provision of renewable heat for industry and buildings via current and prospective renewable thermal energy system (RTES) technologies. The RTES project has undertaken initial research focused on technologies that could be suited for industrial process heat applications at different temperature levels, and, where possible, gathered performance and cost data for these technologies. This project does not directly evaluate RTES for distributed residential or commercial applications, nor does it yet include documented cases or modeling of RTES using geothermal, biomass, waste heat, renewable fuels like renewable natural gas, or hydrogen production. The three technical reports are summarized as follows: Renewable Thermal Energy Systems: Characterization of the Most Important Thermal Energy Applications in Buildings and Industry (Report 1), this report: summary of thermal demands of U.S. industry and buildings, and relevant hybrid RTES configurations; Renewable Thermal Energy Systems: Systemic Challenges and Transformational Policies (Report 2): discussion of socio-technical characteristics of RTES, innovation challenges, and supporting policies. Available at: https://www.nrel.gov/docs/fy23osti/83020.pdf; Renewable Thermal Energy Systems: Modeling Developments and Future Directions (Report 3): Energy yield and performance modeling of RTES, techno-economic analysis via case studies, and proposed development of a user decision support tool. Available at: https://www.nrel.gov/docs/fy23osti/83021.pdf.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Techno-economic design of a linear Fresnel reflector for industrial process heat

A techno-economic model of a Concentrating Solar Thermal (CST) system using a Linear Fresnel Reflector (LFR) has been developed. LFRs can deliver process heat suitable for a range of industries, including food and beverage production. This model uses an adaptive algorithm to calculate the optimal secondary reflector shape given the geometry and optical properties of the rest of the system. A ray-tracing program is used to calculate optical efficiency over a wide range of longitudinal and transversal incidence angles and subsequently evaluate the annual efficiency at a given geographical location. A specific LFR design developed by Hyperlight Energy was modelled, and this industrial partner provided a detailed cost breakdown which was used as the basis of an economic model. Combining the technical and economic data facilitates the calculation of the Levelized Cost of Heat (LCOH). The influence of a number of parameters on the annual efficiency and LCOH is explored; notable parameters include the absorber height, the number, width, and spacing of the primary mirrors, the aim point, and the secondary reflector shape and width. By identifying an optimal combination of these parameters, we reduce the LCOH of the industry partner’s system design by 9.2%, from 14.4 $\$/MWh_{th}$ to 13.0 $\$/MWh_{th}$. In comparison, the LCOH of a natural gas boiler delivering the same annual quantity of heat is 29 $\$/MWh_{th}$, which indicates that LFRs can be a competitive heat source for industrial processes.

14 SOLAR ENERGY↗

Evolution of a Geological Model for Co-Producing Electricity at the Blackburn Oil Field, Nevada: Preprint

Oil fields characterized by significant proportions of hot water in the reservoir offer a unique opportunity for geothermal energy co-production and utilization. Development of an accurate geological model is essential in understanding the subsurface. However, this is a challenge at the preliminary stages of resource assessment when data on the subsurface parameters and subterranean fault network is scarce. Our focus in this paper is to explore a methodology for developing an accurate conceptual geological model of the complexly faulted reservoir at the Blackburn Field, Nevada, to create a numerical simulation that will eventually be used to assess the potential of co-producing electricity. There are four steps to the evolution of this model: (1) collection of public data and research on the geological setting of the system; (2) digitization of seismic time-interpretation and cross-sections to create the most representative model of the subsurface; and (3) creation of a natural-state numerical model to simulate the thermal gradient through the subsurface. Utilizing publicly available data of the geological system and Leapfrog conceptual modeling software, we refined a preliminary model to be more representative of the subsurface geology of an active geothermal exploration stage project. The conceptual geological model developed has been imported into VOLSUNG reservoir modeling software to create a natural-state model, further characterizing the Blackburn Field reservoir's thermal gradient. In the future, new seismic data will be utilized to further evolve the conceptual geological model for numerical simulations.

Basin and Range province↗

Techno-Economic Analysis and Market Potential of Geological Thermal Energy Storage (GeoTES) Charged With Solar Thermal and Heat Pumps

In this project, we developed a techno-economic analysis (TEA) model that can be used to evaluate the viability of a proposed Geological Thermal Energy Storage (GeoTES) design. This MATLAB-based model integrates distinct subsystem models for the reservoir, wells, power cycle, and solar field to capture their distinct characteristics. It applies this approach in simulating GeoTES storage and dispatch operations for durations ranging from hourly to seasonal. Using cases studies based on GeoTES designs provided by industry partners - Premier Resource Management (PRM) and EarthBridge Energy - we validated the TEA model estimations of system performance and costs (such as thermal and electrical power/energy inflow and outflow, capital costs, and levelized costs of energy and storage) for both concentrating solar thermal (CST) and Carnot Battery (CB) pairings with GeoTES (CST-GeoTES and CB-GeoTES). For the CST-GeoTES case, the model was validated against the proposed system designed by PRM. It showed good agreement with PRM's estimations when well and pump costs derived from PRM's estimations were used. When GETEM-based costs were used, there was a slight overprediction due to GETEM's project/site agnostic assumption of these costs. From a sensitivity analysis perspective, the levelized cost of electricity (LCOE) of the CST-GeoTES case was most sensitive to well flow rate and the charging temperature. An optimal design scenario resulted in an LCOE of 0.11 $\$$/kWhe. CST-GeoTES can also provide a source of heat to meet seasonal demands. With 12-hour and 24-hour levelized cost of heat (LCOH) of 0.018 $\$$/kWhth and 0.022 $\$$/kWhth, respectively, CST-GeoTES could be competitive in the California market with an average industrial price of natural gas in California between 0.041-0.047 $\$$/kWhth. The levelized cost of storage (LCOS) for CST-GeoTES depends on the energy storage duration. Although the LCOS is relatively higher for shorter durations (e.g., ~0.50 $\$$/kWhe for 1 hour of storage), it is an order of magnitude lower (0.06 $\$$/kWhe) for longer storage durations and competitive with lithium-ion batteries (beyond 12 hours of storage) and molten-salt thermal energy storage (beyond 32 hours). Energy. Three options were explored and applied to the EarthBridge case study: (1) A Carnot Battery design using R125 working fluid with both hot and cold storage; (2) A Carnot Battery design using R125 working fluid with only hot storage; (3) A Carnot Battery using a commercially available heat pump with carbon dioxide (CO2) working fluid and hot storage only. The CB-GeoTES with cold storage only had a slight (round-trip) efficiency advantage over the system without (43.4% vs. 42.8%). This is because the cold storage is limited by the freezing point of water, so the cold storage is not much colder than the environment. The system using commercially available technologies was the least efficient - partly because different cycles were used in the heat pump (CO2) and heat engine (binary cycle) which leads to some inefficiencies. Using the commercially available design, the levelized cost of energy (LCOS) from the model (0.10 $\$$/kWhe) was higher than that estimated by EarthBridge (0.068 $\$$/kWhe). This is because of the low round-trip (38.7%) efficiency of the commercially available design. Sensitivity analysis reveals that the model is most sensitive to electricity price. Including electricity price in the TEA for CB-GeoTES leads to an increase in LCOS from the base value to 0.25 $\$$/kWhe. To determine storage sites suitable for GeoTES, we gathered and analyzed geological, petrophysical, and geophysical data of oil and gas reservoir and aquifers in California and Texas. We down-selected possible sites based on cut-off values for site characteristics (e.g., reservoir temperature, formation thickness, permeability, porosity, depth, and brine salinity) and preliminary costs. Using this approach, the Carrizo-Wilcox, Yegua-Jackson, and Dockum brackish aquifers in Texas were identified as having the highest suitability. Similarly, in the central California region, the White Wolf, Belridge South Tulare, and Belridge South Reef Ridge were the most suitable. Going further, we assessed the storage potential in the selected sites. To do this we developed distributions of reservoir characteristic data and applied a Monte Carlo-based analysis to account for intrinsic uncertainty in the acquired data. The analysis revealed that the Carrizo-Wilcox aquifer had the highest storage potential with a mean capacity of 554 TWhth (i.e., 63 TWhe). The estimated capacity serves as an upper limit of storage potential given that not all fields in the basin will be developed. We participated in multiple outreach activities including conference presentations, panel session discussions, and the facilitation of a GeoTES workshop at the NREL Golden campus.

15 GEOTHERMAL ENERGY↗

Development of a Three-Dimensional, Unstructured Material Response Design Tool

A preliminary verification and validation of a new material response model is presented. This model, Icarus, is intended to serve as a design tool for the thermal protection systems of re-entry vehicles. Currently, the capability of the model is limited to simulating the pyrolysis of a material as a result of the radiative and convective surface heating imposed on the material from the surrounding high enthalpy gas. Since the major focus behind the development of Icarus has been model extensibility, the hope is that additional physics can be quickly added. This extensibility is critical since thermal protection systems are becoming increasing complex, e.g. woven carbon polymers. Additionally, as a three-dimensional, unstructured, finite-volume model, Icarus is capable of modeling complex geometries. In this paper, the mathematical and numerical formulation is presented followed by a discussion of the software architecture and some preliminary verification and validation studies.

Thermal Protection System↗

Use of MSC/NASTRAN for the thermal analysis of the Space Shuttle Orbiter braking system

A description is given of the thermal modeling and analysis effort being conducted to investigate the transient temperature and thermal stress characteristics of the Space Shuttle Orbiter brake components and subsystems. Models are constructed of the brake stator as well as of the entire brake assembly to analyze the temperature distribution and thermal stress during the landing and braking process. These investigations are carried out on a UNIVAC computer system with MSC/NASTRAN Version 63. Analytical results and solution methods are presented and comparisons are made with SINDA results.

Shu, James↗

Analysis of GaAs and Si solar energy hybrid systems

Various silicon hybrid systems are modeled and compared with a gallium arsenide hybrid system. The hybrid systems modeled produce electric power and also thermal power which can be used for heating or air conditioning. Various performance indices are defined and used to compare the system performance: capital cost per electric power out; capital cost per total power out; capital cost per electric power plus mechanical power; annual cost per annual electric energy; and annual cost per annual electric energy plus annual mechanical work. These performance indices indicate that concentrator hybrid systems can be cost effective when compared with present day energy costs.

Heinbockel, J. H.↗

Field testing of thermal canopy models in a spruce-fir forest

Recent advances in remote sensing technology allow the use of the thermal infrared region to gain information about vegetative surfaces. Extending existing models to account for thermal radiance transfers within rough forest canopies is of paramount importance. This is so since all processes of interest in the physical climate system and biogeochemical cycles are thermally mediated. Model validation experiments were conducted at a well established boreal forest; northern hardwood forest ecotone research site located in central Maine. Data was collected to allow spatial and temporal validation of thermal models. Emphasis was placed primarily upon enhancing submodels of stomatal behavior, and secondarily upon enhancing boundary layer resistance submodels and accounting for thermal storage in soil and vegetation.

Source record↗