SEARCH · Search NASA

Results for “MIXED FLOW”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Challenge Problem 1: Preliminary Model Development and Assessment of Flexible Heat Transfer Modeling Approaches

This report presents the modeling progress of a university consortium studying flexible modeling for turbulent heat transfer. In the context of advanced nuclear reactor design, where challenges include non-unity Prandtl fluids, free and mixed convection, and transitional flow, computational fluid dynamics methods are underutilized due to low confidence in modeling approaches and high computational cost. This report evaluates the performance of engineering CFD models in diabatic, buoyant turbulent flow. It finds that all turbulence models including several variants of the k-ε and k-ω models struggle to predict accurate turbulent momentum and heat transfer in such flows. The Nusselt numbers have been compared between the models and the DNS data, where calculations have been performed for each case and trends have shown a good agreement between DNS estimated Nusselt numbers and available correlations as well as experimental data. A novel DNS correlation for high Pr cases as they are transitioning from mixed convection to natural convection has been proposed. While Nusselt number errors relative to DNS range from 20% to 50%, the models capture similar trends to DNS with respect to Nusselt suppression and amplification under varying levels of buoyancy effect. This report also highlights model form error as a significant contributor to CFD predictions and proposes a framework for quantifying model error and improving confidence in CFD calculations. To improve the predictive capability of engineering CFD models, data-driven approaches for turbulence models are investigated. Theoretical frameworks based on the invariant tensor / vector basis neural networks for prediction of Reynolds stresses and turbulent heat fluxes are employed. The models are developed using direct numerical simulations data for forced convection flows of different fluids in vertical planar channel domain. The framework is implemented in spectral element solvers Nek5000 / nekRS and has shown a potential for future development and consideration of mixed convection flows.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Information loss, mixing and emergent type III 1 factors

A manifestation of the black hole information loss problem is that the two-point function of probe operators in a large Anti-de Sitter black hole decays in time, whereas, on the boundary CFT, it is expected to be an almost periodic function of time. We point out that the decay of the two-point function (clustering in time) holds important clues to the nature of observable algebras, states, and dynamics in quantum gravity. We call operators that cluster in time “mixing” and explore the necessary and sufficient conditions for mixing. The information loss problem is a special case of the statement that in type I algebras, there exists no mixing operators. We prove that, in a thermofield double state (KMS state), if mixing operators form an algebra (close under multiplication), the resulting algebra must be a von Neumann type III 1 factor. In other words, the physically intuitive requirement that all nonconserved operators should exponentially mix is so strong that it fixes the observable algebra to be an exotic algebra called a type III1 factor. More generally, for an arbitrary out-of-equilibrium state of a general quantum system (von Neumann algebra), we show that if the set of operators that mix under modular flow forms an algebra, it is a type III1 von Neumann factor.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Flexible dynamic boundary microgrid operation considering network and load unbalances

Flexible microgrids with dynamic boundaries have recently been introduced in the literature. With the ability to reconfigure the topology of the microgrids dynamically through remotely controlled switches, flexible microgrids with dynamic boundaries can further improve the resiliency and energy efficiency of microgrids with distributed energy resources (DERs). This paper focuses on the optimal operation considering one of the predominant characteristics of microgrids and distribution systems – unbalanced networks and loads. In existing literature, balanced modeling of microgrids is more common due to its attractive simplicity. The three-phase power unbalance has not been considered as a constraint on the generation units in a microgrid. Further, negative sequence constraints have also been neglected. In this article, we propose a set of constraints that is specifically related to the capabilities of inverter interfaced resources to supply unbalanced current/power when the microgrid is islanded from the main distribution grid. We incorporate the new set of constraints into two optimization formulations leveraging two convex relaxations of the three-phase power flow equations: mixed-integer linear programming (MILP) and mixed-integer semidefinite programming (MISDP) that optimize the dispatch of controllable switches and DERs in the microgrid. The algorithms are then extended to networked microgrids with grid-forming sources. We test the algorithms on a realistic community microgrid model in Puerto Rico as well as standardized IEEE distribution test feeders. The testing results demonstrate the performance of the proposed algorithms. The MILP is fast and scalable, and the MISDP enforces the negative sequence voltage constraints.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

Improvements on the Diagnostic Residual Gas Analyzer at Wendelstein 7-X

Exhaust gas analysis provides key information on fusion processes, divertor operation, and wall state in fusion experiments and future reactors. The diagnostic residual gas analyzer (DRGA) concept has been developed for ITER with a focus on fast helium and hydrogen isotope detection. The first operation of the prototype DRGA (P-DRGA) at the stellarator Wendelstein 7-X showed potential for improvement in terms of magnetic sensor shielding, data acquisition automation, and potential new additions to the cluster of sensors on the P-DRGA. More recently, a Monte Carlo simulation of the flow of the mixed gas species effluent from the pressure-reducing orifice, down to about 8-m sampling tube and into the analysis region of the sensors, has been found to generally agree with previous calculations and measurements but revealed potential back-streaming effects for light gases, with impact on detection limits both for the prototype and for the ITER DRGA currently in design. For the upcoming campaign of the prototype, an enhanced soft iron shield will safeguard the gauges against magnetic stray field influence. The newly introduced shielding has been tested for its effect on magnetic stray fields and found to reduce the inside residual field by about two orders of magnitude.

Schlisio, G↗

DOE OSTI2026

Watershed response of the Feather River Basin, California, United States of America to future climate changes

Increasing mean annual temperatures under climate change are expected to reduce seasonal snowpack, increase evapotranspiration (ET), and alter summer baseflow in headwater watersheds worldwide. Strong regional variability in hydrologic responses highlights the need for catchment-scale, physically based models to assess future flood risk and water availability. This study examines climate-driven changes in the hydrologic response of the Upper Feather River watershed in the Sierra Nevada Mountains, California. Four priority climate models and two representative concentration pathways (RCP4.5 and RCP8.5) are used to evaluate future hydrologic responses of the watershed. Hydrologic processes are simulated using an objectively calibrated Soil and Water Assessment Tool Plus model for a historical baseline (1986–2005) and a future period 2070–2099), driven by observed and projected precipitation and temperature. Results indicate a declining contribution of snowfall to annual precipitation, with peak snowfall and water yield shifting 1–3 months earlier. Long-term annual maximum flows are projected to increase considerably, whereas low-flow responses are mixed, with both increases and decreases projected by the end of the century. These findings highlight the need for adaptive watershed management to enhance flood protection, water storage, and drought resilience. Future water resource planning should also account for one-to-three-month shifts in peak water yield and surface runoff due to changes in snowmelt timing and a lower snowfall-to-rainfall ratio under climate change. The modeling framework and insights are transferable to other snow-dominated headwater watersheds experiencing climate-driven change.

Tigabu, T [UC Davis]↗

DOE OSTI2020

RELAP5-3D Modeling of High Temperature Test Facility (HTTF) Test PG-26

The High Temperature Test Facility (HTTF) at Oregon State University (OSU) is a scaled integral effects experiment designed to investigate transient behavior in high-temperature gas-cooled nuclear reactors with prismatic fuel and reflector blocks. Several tests have been completed, and more are still planned to at the HTTF, including depressurized conduction cooldown (DCC) and pressurized conduction cooldown (PCC) transients. This report analyses test PG-26, a progression of the Double Ended Inlet-Outlet Crossover Duct Break transient that is referred to as a DCC. PG-26 has been performed at the HTTF between May 30 and June 30, 2019. Core initial conditions (i.e., before the DCC started) have been met using low power (<100 kW) and two of ten available electric heaters. The DCC transient was initiated during the 50 th hour of the test. The break valves were opened, and hot helium from the core and cold helium from the reactor cavity simulation tank (RCST) started mixing. The gases flowed in a countercurrent fashion, where the top half of the hot duct contained hot helium that flowed in one direction and cold helium that flowed in the other direction in the bottom half of the duct. After the pressure and density reached equilibrium, the event entered a diffusion mode. The onset of a reverse natural circulation was not observed during the DCC period of the test. Version 4.4.2ie of the RELAP5-3D computer code has been used to model the HTTF PG-26 test, and results have been compared to available high-quality measured data. The model used in this study is the quality-controlled HTTF RELAP5-3D model (HTTF base 2018-04-19 QA), originally developed by P. Bayless. The report includes RELAP5-3D results of the “base calculations” as well as some sensitivities to important uncertain model inputs, such as primary helium mass flow rate, core ceramic thermal properties, as well as heat evacuation and loop friction models. Using the base RELAP5-3D model predicts a countercurrent helium flow in the hot duct observed at the beginning of the DCC, but instead of going into a molecular diffusion mode, the model predicts the onset of natural convection. Increasing friction in the core and hot duct prevents the natural convection from happening in some of the simulations. Although some temperatures are well predicted (and even overpredicted), the general tendency is to underpredict the ceramic and helium temperatures and heat removal rates during the DCC, resulting in many of the assessment findings being in minimal or insufficient agreement with the data. It is worth noting that the described discrepancies between measured data and RELAP5-3D predictions are not RELAP5-3D code limitations. More so, they reflect limitations in boundary condition and thermal property knowledge. While the RELAP5-3D calculations of the test provide some insights into what happens during the transient, and point to missing or potentially uncertain data to which the experimenters can direct their attention, the principal conclusion is that the PG-26 test data are insufficient for a system code assessment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Compact Counter Flow Fluidized Bed Heat Exchanger (Final Report)

Generation 3 concentrating solar power systems which utilize solid granular media as the heat transfer and storage media hold promise in reducing the levelized cost of electricity relative to generation 2 systems which utilize molten salt. The particles can be heated to temperatures exceeding 800 °C which enables the use of high efficiency supercritical CO 2 Brayton power cycles. A key challenge associate with the use of solid media is its low bulk thermal conductivity which limits the ability to transfer heat from the particles to a working fluid, such as supercritical CO 2 . Tradition shell and tube and shell and plate heat exchangers feature gravity driven particle flow which limits the maximum particle side heat transfer coefficient in these systems to approximately 400 W/m 2 K. A fluidized bed heat exchanger does not rely solely on gravity for particle flow. Instead, enhanced mixing and horizontal conveyance is achieved through fluidization induced through air injection at the bottom of a bed of particles. A tube bundle within the bed conveys the supercritical CO 2 which is heated via the turbulent particle motion at the tube wall.

14 SOLAR ENERGY↗

DOE OSTI2021

Modeling of HTTF test PG-26 using RELAP5-3D and SAM

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Experimental measurement of the rapid mixing of fuel and air in a multi-element diffusion (Hencken) burner

Flat flame burners are valuable tools for studying complex combustion phenomena like particle synthesis due to the well-controlled and well-behaved initial and boundary conditions. In this study, the methane/air flame temperature and tracer density profiles of a multi-element diffusion burner equipped with a central tube were characterized using krypton (Kr) x-ray fluorescence at the Advanced Photon Source at Argonne National Laboratory. The experiments were performed using a fuel-to-oxygen equivalence ratio of 0.8 and 3.3% Kr by volume was added to the reactant gases. One- and two-dimensional scans of non-combusting (i.e., cold flow) and combusting conditions, both with and without flow through the central tube were performed. Spatially-resolved Kr number density profiles were obtained at different planar locations and heights above the surface of the burner. The results show rapid mixing occurred near the exit-plane of the burner and increased spatial uniformity was achieved from 0.1 mm to 0.5 mm above the surface of the burner. Further, the data also showed the impact of the central tube flow on the main methane/air flame, including evidence of the central tube flow mixing with the main flow. Flame temperatures were calculated from the fluorescence measurements at various heights above the burner surface and were in good agreement with model predictions for a burner stabilized flame. These high-fidelity measurements provide accurate baseline data that are particularly valuable for understanding the effects of fuel and air mixing on downstream phenomena, such as the formation of condensed-phase products.

42 ENGINEERING↗

DOE OSTI2020

Hydrodynamic behavior of wrapped coil water heater tank

In this work, a numerical investigation using computational fluid dynamics (CFD) has been carried out to analyze the thermal behavior of a water tank integrated to a heat pump water heating system (HPWH). The heat pump condenser is a long coil wrapped around the water tank and insulated from outside from which heat is transferred from the heat pump cycle into the water stored inside the tank. The simulated water tank has a capacity of about 170 L with the charging inlet located at the bottom segment of the tank and the discharging outlet located at the top segment of the tank like commercially available water heaters. The tank was simulated during charging/discharging process with isothermal boundary condition at the tank wall to represent the majority of the process for refrigerant condensation inside the heat pump condenser (excluding superheated and sub-cooled section). Three flow rates were selected to represent low, medium and high flow pattern. Both mean and turbulent flow properties were calculated and analyzed. Heat transfer mechanisms inside the tank are quantified by calculating Richardson and Péclet numbers to show the relative contribution of driving forces. The calculations show that advection is dominant over diffusion. Mixing factor is calculated and presented. It is found that the rate of increase of mixing is nonlinear with the increase in flow rate. Mixing increased by a factor of four when the flow rate increased by a factor of two.

stratification↗

DOE OSTI2021

Development of five-moment two-fluid modeling for Z-pinch physics

The Z-pinch m = 0 instability as well as its stabilization by radially sheared axial flow is studied using the nonlinear ideal five-moment twofluid (5M2F) model with an extension of that model to include Braginskii heat and momentum transport. Using the ideal 5M2F model, linear growth rate results are compared with prior work using MHD and Hall MHD. At small normalized wavenumber, 1 < k a < 4 , where a is the effective pinch radius, 5M2F results agree with Hall MHD within 20% in scenarios without radially sheared axial flow. With the sheared flow and focusing on ka = 10/3, agreement with Hall MHD is excellent. In the limit of small ion inertial length, results also match with MHD. A comparison with PIC modeling of shear-free m = 0 stability focuses on a plasma scenario based on recent experimental results. In a scan of mode wavenumber, ideal 5M2F results are qualitatively similar to PIC: the growth rate rises to a peak at a moderate wavenumber and declines at a large wavenumber in contrast to MHD results, which show the saturation of the growth rate with the increasing wavenumber rather than a decline. The peak normalized 5M2F growth rate is γ τ A = 1.5 , where sA is the Alfven transit time across the pinch. The peak occurs at normalized wavenumber ka = 10. For comparison, PIC results have a peak growth of γ τ A = 0.77 at ka = 5. Including Braginskiibased closure of the 5M2F model does not qualitatively change the ideal results in this particular case. Nonlinear saturation is studied using the 5M2F model with the dissipative Braginskii-based closure in cases with pinch-edge sheared-flow speed equal to half the Alfven speed. Nonlinear mixing due to the sheared flow yields a quasi-steady state after modest losses of pinch ion inventory and pinch thermal energy, approximately 30% and 10%, respectively. 5M2F modeling captures the essential physics of m = 0 instability and offers a computationally tractable route to high-fidelity modeling of 3D Z-pinch behavior, including m = 1 instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Dual Modification of MOFs Improves Dispersion and Ionic Conductivity of Mixed Matrix Membranes

Nonaqueous redox flow batteries (NARFBs) are a promising class of energy storage devices, but the lack of a chemically stable, conductive membrane that exhibits size-selectivity over redox-active species prevents their broader implementation. Recently, metal–organic frameworks (MOFs) have been implemented into mixed-matrix membranes (MMMs) for NARFBs, but the effects of the MOF linker functionality on membrane properties are not well-understood. In this work, we develop a series of MOF-based MMMs by blending postsynthetically modified variants of UiO-66-NH2 with poly(ethylene-co-vinyl acetate). The modification of UiO-66-NH2 with sulfate groups initially resulted in poor dispersion throughout the MMMs, but when dual-modified with poly(N-isopropylacrylamide), MOF dispersion throughout the MMM was improved, and ionic conductivity was significantly higher than the UiO-66-NH2 MMMs. Furthermore, the dual-modified MMMs demonstrated excellent size-selectivity by blocking redox active species transport. This work demonstrates a synergy between the MOF functional groups to improve MMM properties critical for the development of practical NARFBs.

Carey, Cassidy A.↗

DOE OSTI2024

Fluid inertia controls mineral precipitation and clogging in pore to network-scale flows

Mineral precipitation caused by fluid mixing presents complex control and predictability challenges in a variety of natural and engineering processes, including carbon mineralization, geothermal energy, and microfluidics. Precipitation dynamics, particularly under the influence of fluid flow, remain poorly understood. Combining microfluidic experiments and three-dimensional reactive transport simulations, we demonstrate that fluid inertia controls mineral precipitation and clogging at flow intersections, even in laminar flows. We observe distinct precipitation regimes as a function of Reynolds number (Re). At low Reynolds numbers (Re < 10), precipitates form a thin, dense layer along the mixing interface, which shuts precipitation off, while at high Reynolds numbers (Re > 50), strong three-dimensional flows significantly enhance precipitation over the entire intersection, resulting in rapid clogging. When injection rates from two inlets are uneven, flow symmetry-breaking leads to unexpected flow bifurcation phenomena, which result in enhanced concurrent precipitation in both downstream channels. Finally, we extend our findings to rough channel networks and demonstrate that the identified inertial effects on precipitation at the intersection scale are also present and even more dramatic at the network scale. This study sheds light on the fundamental mechanisms underlying mixing-induced mineral precipitation and provides a framework for designing and optimizing processes involving mineral precipitation.

Science & Technology - Other Topics↗

DOE OSTI2024

A method for generating quantitative vapor-phase infrared spectra of solids: results for phenol, camphor, menthol, syringol, dicyclopentadiene and naphthalene

Here, a method is presented to generate quantitative vapor-phase infrared spectra from substances that naturally occur as solids with moderate volatility. The solid is gravimetrically dissolved into a solvent that has few infrared spectral features, typically CS 2 and CCl 4 separately. The solution is flowed at a constant rate from a linearly pumped syringe into a metered stream of nitrogen carrier gas regulated by a mass flow controller. The analyte/solvent mix is flash vaporized by volatilizing the solution across a heated stainless-steel surface as it emanates from the syringe tip. The N 2 gas-solution mixture is flowed into a long-path White cell thermostatted at a desired temperature, the long optical path compensating for the modest analyte mixing ratio. A composite spectrum is generated from typically ten or more 760-Torr pressure-broadened spectra over the 600 to 6500 cm -1 spectral range at 0.1 cm -1 spectral resolution. The solid analytes reported here using this novel technique include dicyclopentadiene, menthol, syringol, phenol, camphor, and naphthalene.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2025

Tightest Mixed-Integer Programming Formulations for Quadratic SCUC Optimization

In this project, we developed new, tighter Mixed-Integer Programming (MIP) formulations for the combined Alternating Current (AC) Security-Constrained Unit Commitment (SCUC) and Security-Constrained Optimal Power Flow (SCOPF). The work addresses a critical challenge in power system operations: efficiently determining which generation units to commit and how to optimally dispatch them while maintaining network reliability constraints for both normal and contingency scenarios. Our efforts: 1. Advance the Understanding of SCUC/SCOPF Modeling: By introducing tighter MIP formulations and leveraging cutting-edge optimization tools (Julia/JuMP, PowerModels.jl), this project has pushed forward the state of the art in efficient power systems scheduling. 2. Enhance Technical and Economic Feasibility: The methods developed provide more accurate and potentially faster solutions to large-scale, realistic scheduling and dispatch problems in electric power systems, which can translate into improved reliability and potentially lower costs for grid operations. 3. Benefit to the Public: Greater efficiency in power system operations leads to cost savings for utilities and end-users. Improved reliability and integration of advanced modeling approaches can facilitate the adoption of clean energy resources and better accommodate uncertainties in renewable generation. Because this technology could impact bulk power markets and reliability, these innovations have far-reaching public benefits in terms of cost savings, reliability, and sustainability.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2024

Magneto-Stokes flow in a shallow free-surface annulus

In this study, we analyse ‘magneto-Stokes’ flow, a fundamental magnetohydrodynamic (MHD) flow that shares the cylindrical-annular geometry of the Taylor–Couette cell but uses applied electromagnetic forces to circulate a free-surface layer of electrolyte at low Reynolds numbers. The first complete, analytical solution for time-dependent magneto-Stokes flow is presented and validated with coupled laboratory and numerical experiments. Three regimes are distinguished (shallow-layer, transitional and deep-layer flow regimes), and their influence on the efficiency of microscale mixing is clarified. The solution in the shallow-layer limit belongs to a newly identified class of MHD potential flows, and thus induces mixing without the aid of axial vorticity. We show that these shallow-layer magneto-Stokes flows can still augment mixing in distinct Taylor dispersion and advection-dominated mixing regimes. The existence of enhanced mixing across all three distinguished flow regimes is predicted by asymptotic scaling laws and supported by three-dimensional numerical simulations. Mixing enhancement is initiated with the least electromagnetic forcing in channels with order-unity depth-to-gap-width ratios. If the strength of the electromagnetic forcing is not a constraint, then shallow-layer flows can still yield the shortest mixing times in the advection-dominated limit. Our robust description of momentum evolution and mixing of passive tracers makes the annular magneto-Stokes system fit for use as an MHD reference flow.

58 GEOSCIENCES↗

DOE OSTI2022

A 1D–0D–3D coupled model for simulating blood flow and transport processes in breast tissue

In this work, we present mixed dimensional models for simulating blood flow and transport processes in breast tissue and the vascular tree supplying it. These processes are considered, to start from the aortic inlet to the capillaries and tissue of the breast. Large variations in biophysical properties and flow conditions exist in this system necessitating the use of different flow models for different geometries and flow regimes. In total, we consider four different model types. First, a system of 1D nonlinear hyperbolic partial differential equations (PDEs) is considered to simulate blood flow in larger arteries with highly elastic vessel walls. Second, we assign 1D linearized hyperbolic PDEs to model the smaller arteries with stiffer vessel walls. The third model type consists of ODE systems (0D models). It is used to model the arterioles and peripheral circulation. Finally, homogenized 3D porous media models are considered to simulate flow and transport in capillaries and tissue within the breast volume. Sink terms are used to account for the influence of the venous and lymphatic systems. Combining the four model types, we obtain two different 1D–0D–3D coupled models for simulating blood flow and transport processes: The first model results in a fully coupled 1D–0D–3D model covering the complete path from the aorta to the breast combining a generic arterial network with a patient specific breast network and geometry. The second model is a reduced one based on the separation of the generic and patient specific parts. The information from a calibrated fully coupled model is used as inflow condition for the patient specific sub-model allowing a significant computational cost reduction. Several numerical experiments are conducted to calibrate the generic model parameters and to demonstrate realistic flow simulations compared to existing data on blood flow in the human breast and vascular system. Moreover, we use two different breast vasculature and tissue data sets to illustrate the robustness of our reduced sub-model approach.

multiscale blood flow models↗

DOE OSTI2025

Model Scripts for "Old-Aged Groundwater Contributes to Mountain Hillslope Hydrologic Dynamics"

The partitioning of water inputs between deep and shallow groundwater flow paths is a fundamental processes, yet is challenging to observe. Numerical models provide a valuable tool to further develop insights on these groundwater mixing processes. This repository contains scripts to run the ParFlow-CLM and EcoSLIM integrated hydrologic models along the Pumphouse Hillslope in the East River Watershed and associated python scripts to process model outputs. Files includes input decks for the associated models, python scripts, and text and css files to support the model runs and processing. The 2-D hillslope model simulates from 2000 to 2021 using transient forcing conditions from the publicly available NLDAS-2 dataset. We develop an ensemble of models with variable hydrogeologic and soil parameters. Model outputs from the ensemble of runs are compared to to water level observations at the PLM1 and PLM6 wells published at (https://data.ess-dive.lbl.gov/view/doi:10.15485/1866836) and groundwater mean ages published at (https://data.ess-dive.lbl.gov/datasets/doi:10.15485/1960042. The model results are used to evaluate hillslope water mass-balance transience and mixing dynamics between groundwater flow paths with young (<10 year) and old (>10 year) ages.

54 ENVIRONMENTAL SCIENCES↗