Search NASA⌕ Search

SEARCH · Search NASA

Results for “process simulation”

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

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

At least 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

In-situ tube burst testing and high-temperature deformation behavior of candidate materials for accident tolerant fuel cladding

The high resistance of cladding to plastic deformation and burst failure is one of the most essential properties of accident tolerant fuel (ATF) for maintaining structural integrity during a loss-of-coolant accident (LOCA) since the deformation and burst behavior governs the cooling efficiency of flow channels and process of fission product release. To simulate and evaluate such deformation and burst process of thin-walled cladding, an in-situ testing and evaluation method has been developed on the basis of visual imaging and image analysis techniques. The method uses a specialized optics system consisted of a high-resolution video camera, light filtering unit, and monochromatic light sources, and the in-situ testing is performed using a 50 mm long pressurized thin-walled tubular specimen set in a programmable furnace. In this study eleven (11) candidate cladding materials for ATF, i.e., 6 FeCrAl alloys and 5 nanostructured steels, were tested using the newly developed method, and the time-dependent images were analyzed to produce detailed deformation and burst data such as true hoop stress, strain (creep) rate, and failure stress. Relatively soft FeCrAl alloys deformed and burst below 800°C while negligible strain rates were measured for higher strength alloys and/or for relatively thick wall specimens.

Gussev, Maxim N.↗
DOE OSTI2022

Avoiding excess computation in asynchronous evolutionary algorithms

Abstract Asynchronous evolutionary algorithms are becoming increasingly popular as a means of making full use of many processors while solving computationally expensive search and optimization problems. These algorithms excel at keeping large clusters fully utilized, but may sometimes inefficiently sample an excess of fast‐evaluating solutions at the expense of higher‐quality, slow‐evaluating ones. We have previously introduced a steady‐state parent selection strategy, SWEET (“Selection whilE EvaluaTing”), that sometimes selects individuals that are still being evaluated and allows them to reproduce early. We perform a takeover‐time analysis that confirms that this strategy gives slow‐evaluating individuals that have higher fitnesses an increased ability to multiply in the population. We also find that SWEET appears effective at improving optimization performance on problems in which solution quality is positively correlated with evaluation time. We evaluate our approach on six simulated real‐valued optimization problems and three real‐world applications: an autonomous vehicle controller problem that involves tuning a spiking neural network and two adversarial EA problems. We further evaluate SWEET versus a basic asynchronous process in a simulated setting. We present evidence that SWEET outperforms basic asynchronous processes in a use‐case in which performance is positively correlated with evaluation time, and performs comparably (and often better) than basic asynchronous processes in several use‐cases where performance is negatively correlated with evaluation time. That said, in the cases where performance and evaluation time are negatively correlated the variance of outcomes for SWEET is notably high.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Long-term measurements in a mixed-grass prairie reveal a change in soil organic carbon recalcitrance and its environmental sensitivity under warming

Soil respiration, the major pathway for ecosystem carbon (C) loss, has the potential to enter a positive feedback loop with the atmospheric CO2 due to climate warming. For reliable projections of climate-carbon feedbacks, accurate quantification of soil respiration and identification of mechanisms that control its variability are essential. Process-based models simulate soil respiration as functions of belowground C input, organic matter quality, and sensitivity to environmental conditions. However, evaluation and calibration of process-based models against the long-term in situ measurements are rare. Here, we evaluate the performance of the Terrestrial ECOsystem (TECO) model in simulating total and heterotrophic soil respiration measured during a 16-year warming experiment in a mixed-grass prairie; calibrate model parameters against these and other measurements collected during the experiment; and explore whether the mechanisms of C dynamics have changed over the years. Calibrating model parameters against observations of individual years substantially improved model performance in comparison to pre-calibration simulations, explaining 79–86% of variability in observed soil respiration. Interannual variation of the calibrated model parameters indicated increasing recalcitrance of soil C and changing environmental sensitivity of microbes. Overall, we found that (1) soil organic C became more recalcitrant in intact soil compared to root-free soil; (2) warming offset the effects of increasing C recalcitrance in intact soil and changed microbial sensitivity to moisture conditions. Furthermore, these findings indicate that soil respiration may decrease in the future due to C quality, but this decrease may be offset by warming-induced changes in C cycling mechanisms and their responses to moisture conditions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Benchmark Problems of the Geothermal Technologies Office Code Comparison Study

A diverse suite of numerical simulators is currently being applied to predict or understand the performance of enhanced geothermal systems (EGS). To build confidence and identify critical development needs for these analytical tools, the United States Department of Energy, Geothermal Technologies Office has sponsored a Code Comparison Study (GTO-CCS), with participants from universities, industry, and national laboratories. A principal objective for the study was to create a community forum for improvement and verification of numerical simulators for EGS modeling. Teams participating in the study were those representing U.S. national laboratories, universities, and industries, and each team brought unique numerical simulation capabilities to bear on the problems. Two classes of problems were developed during the study, benchmark problems and challenge problems. The benchmark problems were structured to test the ability of the collection of numerical simulators to solve various combinations of coupled thermal, hydrologic, geomechanical, and geochemical processes. This class of problems was strictly defined in terms of properties, driving forces, initial conditions, and boundary conditions. Study participants submitted solutions to problems for which their simulation tools were deemed capable or nearly capable. Some participating codes were originally developed for EGS applications whereas some others were designed for different applications but can simulate processes similar to those in EGS. Solution submissions from both were encouraged. In some cases, participants made small incremental changes to their numerical simulation codes to address specific elements of the problem, and in other cases participants submitted solutions with existing simulation tools, acknowledging the limitations of the code. The challenge problems were based on the enhanced geothermal systems research conducted at Fenton Hill, near Los Alamos, New Mexico, between 1974 and 1995. The problems involved two phases of research, stimulation, development, and circulation in two separate reservoirs. The challenge problems had specific questions to be answered via numerical simulation in three topical areas: 1) reservoir creation/stimulation, 2) reactive and passive transport, and 3) thermal recovery. Whereas the benchmark class of problems were designed to test capabilities for modeling coupled processes under strictly specified conditions, the stated objective for the challenge class of problems was to demonstrate what new understanding of the Fenton Hill experiments could be realized via the application of modern numerical simulation tools by recognized expert practitioners.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Long-time integration of parametric evolution equations with physics-informed DeepONets

Ordinary and partial differential equations (ODEs/PDEs) play a paramount role in analyzing and simulating complex dynamic processes across all corners of science and engineering. In recent years machine learning tools are aspiring to introduce new effective ways of simulating such equations, however existing approaches are not able to reliably return stable and accurate predictions across long temporal horizons. We aim to address this challenge by introducing an effective framework for learning evolution operators that map random initial conditions to associated ODE/PDE solutions within a short time interval. Such operators can be parametrized by deep neural networks that are trained in an entirely self-supervised manner without requiring one to generate any paired input-output observations. Global long-time predictions across a range of initial conditions can be then obtained by iteratively evaluating the trained model using each prediction as the initial condition for the next evaluation step. Here, this introduces a new approach to temporal domain decomposition that is shown to be effective in performing accurate long-time simulations for a wide range of parametric ODE and PDE systems, from wave propagation, to reaction-diffusion dynamics and stiff chemical kinetics, introducing a new way of rapidly emulating non-equilibrium processes in science and engineering.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Predicting the structural basis of targeted protein degradation by integrating molecular dynamics simulations with structural mass spectrometry

Targeted protein degradation (TPD) is a promising approach in drug discovery for degrading proteins implicated in diseases. A key step in this process is the formation of a ternary complex where a heterobifunctional molecule induces proximity of an E3 ligase to a protein of interest (POI), thus facilitating ubiquitin transfer to the POI. In this work, we characterize 3 steps in the TPD process. (1) We simulate the ternary complex formation of SMARCA2 bromodomain and VHL E3 ligase by combining hydrogen-deuterium exchange mass spectrometry with weighted ensemble molecular dynamics (MD). (2) We characterize the conformational heterogeneity of the ternary complex using Hamiltonian replica exchange simulations and small-angle X-ray scattering. (3) We assess the ubiquitination of the POI in the context of the full Cullin-RING Ligase, confirming experimental ubiquitinomics results. Differences in degradation efficiency can be explained by the proximity of lysine residues on the POI relative to ubiquitin.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Modeling the Capacitive Deionization Process in Dual-Porosity Electrodes

In many areas of the world, there is a need to increase water availability. Capacitive deionization (CDI) is an electrochemical water treatment process that can be a viable alternative for treating water and for saving energy. A model is presented to simulate the CDI process in heterogeneous porous media comprising two different pore sizes. It is based on a theory for capacitive charging by ideally polarizable porous electrodes without Faradaic reactions or specific adsorption of ions. A two steps volume averaging technique is used to derive the averaged transport equations in the limit of thin electrical double layers. A one-equation model based on the principle of local equilibrium is derived. The constraints determining the range of application of the one-equation model are presented. The effective transport parameters for isotropic porous media are calculated solving the corresponding closure problems. The source terms that appear in the average equations are calculated using theoretical derivations. The global diffusivity is calculated by solving the closure problem.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2018

When managers forage for pests: Implementing the functional response in pest management

In this work, we explore how the functional response framework can be implemented in pest management. Here, managers take the role of predators foraging on pests and facing monetary costs for survey and control in a spatial domain where the pest distribution and control strategy do not have to be random. To investigate this framework quantitatively, we simulated various management processes on different pest spatial distributions using a spatially-explicit individual-based model and Monte-Carlo simulations, and also confirmed some of the results analytically. By graphing the number of pests controlled versus pest density, we obtained management functional response curves. Whether the management functional response was shaped like a type I, type II or type III functional response depended on the management costs and the search area. Yet, the management spatial strategy and the pest spatial distribution had little effect on the functional response. We applied our model to the management of mountain pine beetle epidemic in Cypress Hills, Saskatchewan, Canada, with simulations matching the real number of attacked trees controlled by managers. We showed how to make an analogy between functional responses in predator–prey interactions and in human–pest interactions and thereby, apply insights from the functional response framework to pest management.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Interinstitutional Study of the New EURO-GANEX Process Resistance by Gamma Irradiation Test Loops

As part of the homogeneous actinides recycling strategy, the EURO-GANEX process is one of the most promising options to achieve the goal of minor actinides recovery. However, EURO-GANEX also has various drawbacks that are being addressed. Improvements made to EURO-GANEX system have resulted in the emergence of the so-called New EURO-GANEX process, where the composition of the solvent has been modified by replacing TODGA and DMDOHEMA with cis-mTDDGA in the organic phase and SO3-Ph-BTP with PyTri-Diol in the aqueous phase in order to resolve important issues. The objective of this work is twofold: evaluate the gamma radiolysis resistance of the new EURO-GANEX process by dynamic irradiation conditions simulating the three main steps of the process and validate the design of CIEMAT Náyade, CEA Marcel and INL irradiation loop devices since each of them mimics different aspects of the real process. Náyade and the INL loops could irradiate together the organic and aqueous phases, whereas in the CEA loop, the irradiated solvent is recycled continuously inside a platform with several stages of mixer-settlers containing aqueous flows simulating the three main steps of the process. The extraction performances and changes in the composition of the solvent have been analysed during the irradiation experiment by different techniques: gamma spectrometry and ICP-MS/OES for cations or tracers extraction, and HPLC-MS to quantify the degradation and evaluate the degradation compounds. Additionally, in spite of some differences between the three-irradiation facilities, this inter-institutional study shows that they are three comparative tools, which provide results on the stability towards radiolysis of a liquid-liquid extraction system. Favourable extraction results for the different steps are obtained according to the static irradiation studies found in literature. However, the degradation of cis-mTDDGA is appreciable leading to degradation compounds, some of which form precipitates and produce important changes in viscosity, important aspects that must be addressed prior to the successful industrial application of the new EURO-GANEX process.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

DOE-GENIORS Dynamic Irradiation Testing

As part of the homogeneous actinides recycling strategy, the EURO-GANEX process is one of the most promising options to achieve the goal of minor actinides recovery. Improvements made to EURO-GANEX system have resulted in the emergence of the so-called New EURO-GANEX process, where the composition of the solvent has been modified by replacing TODGA and DMDOHEMA with cis-mTDDGA in the organic phase and SO3-Ph-BTP with PyTri-Diol in the aqueous phase in order to resolve important issues. The objective of this work is twofold: evaluate the gamma radiolytic resistance of the new EURO-GANEX process by dynamic irradiation conditions simulating the three main steps of the process, and validate the design of CIEMAT Náyade, CEA Marcel and INL irradiation loop devices since each of them mimics different aspects of the real process. Náyade and the INL loops could irradiate the organic and aqueous phases together, whereas in the CEA loop, the irradiated solvent is recycled continuously inside a platform with several stages of mixer-settlers containing aqueous flows simulating the three main steps of the process. The extraction performances and changes in the composition of the solvent have been analyzed during the irradiation experiment by different techniques: gamma spectrometry and ICP-MS/OES for cations or radioactive tracer extraction, and HPLC-MS to identify and quantify the degradation compounds. Despite some differences between the three-irradiation facilities, this inter-institutional study shows that these three comparative tools provide similar trends in the radiolytic stability of a liquid-liquid extraction system. Favourable extraction results for the different steps are obtained according to the static irradiation studies found in literature. However, the degradation of cis-mTDDGA is appreciable leading to degradation compounds, some of which form precipitates and produce important changes in viscosity, important aspects that must be addressed prior to the successful industrial application of the new EURO-GANEX process.

38 - RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCL↗
DOE OSTI2023

Modeling of Chemical Slurry Rheology in DWPF Sludge Batch (SB) 10 Simulants

The Defense Waste Processing Facility (DWPF) treats high-activity radionuclides from sludge through a process called vitrification. This process converts radioactive liquid waste currently stored in tank farms into a solid glass form that is suitable for long-term storage and disposal. Due to the complexities involved in vitrifying this waste within each operation of the Chemical Processing Cell (CPC), waste rheology is studied to characterize the fluid-mechanical properties as it passes through the CPC and into the Melter. To better understand the waste and validate flow behavior, slurry rheology of simulants that represents the waste was studied at various acid stoichiometry percentages and solids concentrations to determine the simulant’s yield stress and viscosity. This research work has been supported by the DOE-FIU Science & Technology Workforce Development Initiative, an innovative program developed by the U.S. Department of Energy’s Office of Environmental Management (DOE-EM) and Florida International University’s Applied Research Center (FIU-ARC). During the spring of 2022, a DOE Fellow intern, Brendon Cintas, spent 10 weeks doing a summer internship at Savannah River National Laboratory (SRS) under the supervision and guidance of Dan Lambert, Chemical Flowsheet Development. The intern’s project was initiated on June 6, 2022, and continued through August 11, 2022 with the objective of assisting scientists at SRNL’s Rheology and Grout Laboratory at Aiken Country Technology Lab (ACTL) better understand the sludge composition on the rheology of a simulant slurry using a HAAKE RheoStress 6000 rheometer and extrapolate the results to the real-waste data.

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

Sensitivities of subgrid-scale physics schemes, meteorological forcing, and topographic radiation in atmosphere-through-bedrock integrated process models: a case study in the Upper Colorado River basin

Abstract. Mountain hydrology is controlled by interacting processes extending from the atmosphere through the bedrock. Integrated process models (IPMs), one of the main tools needed to interpret observations and refine conceptual models of the mountainous water cycle, require meteorological forcing that simulates the atmospheric process to predict hydroclimate then subsequently impacts surface–subsurface hydrology. Complex terrain and extreme spatial heterogeneity in mountainous environments drive uncertainty in several key considerations in IPM configurations and require further quantification and sensitivity analyses. Here, we present an IPM using the Weather Research and Forecasting (WRF) model which forces an integrated hydrologic model, ParFlow-CLM, implemented over a domain centered over the East River watershed (ERW), located in the Upper Colorado River basin (UCRB). The ERW is a heavily instrumented 300 km2 region in the headwaters of the UCRB near Crested Butte, CO, with a growing atmosphere-through-bedrock observation network. Through a series of experiments in the water year 2019 (WY19), we use four meteorological forcings derived from commonly used reanalysis datasets, three subgrid-scale physics scheme configurations in WRF, and two terrain shading options within WRF to test the relative importance of these experimental design choices for key hydrometeorological metrics including precipitation and snowpack, as well as evapotranspiration, groundwater storage, and discharge simulated by the ParFlow-CLM. Our hypothesis is that uncertainty from synoptic-scale forcings produces a much larger spread in surface–subsurface hydrologic fields than subgrid-scale physics scheme choice. Results reveal that the WRF subgrid-scale physics configuration leads to larger spatiotemporal variance in simulated hydrometeorological conditions, whereas variance across meteorological forcing with common subgrid-scale physics configurations is more spatiotemporally constrained. Despite reasonably simulating precipitation, a delay in simulated discharge peak is due to a systematic cold bias across WRF simulations, suggesting the need for bias correction. Discharge shows greater variance in response to the WRF simulations across subgrid-scale physics schemes (26 %) rather than meteorological forcing (6 %). The topographic radiation option has minor effects on the watershed-average hydrometeorological processes but adds profound spatial heterogeneity to local energy budgets (±30 W m−2 in shortwave radiation and 1 K air temperature differences in late summer). This is the first presentation of sensitivity analyses that provide support to help guide the scientific community to develop observational constraints on atmosphere-through-bedrock processes and their interactions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

The Catalogue for Astrophysical Turbulence Simulations (CATS)

Turbulence is a key process in many fields of astrophysics. Advances in numerical simulations of fluids over the last several decades have revolutionized our understanding of turbulence and related processes such as star formation and cosmic ray propagation. However, data from numerical simulations of astrophysical turbulence are often not made public. We introduce a new simulation-oriented database for the astronomical community: the Catalogue for Astrophysical Turbulence Simulations (CATS), located at www.mhdturbulence.com. CATS includes magnetohydrodynamic (MHD) turbulent box simulation data products generated by the public codes athena++, arepo, enzo, and flash. CATS also includes several synthetic observational data sets, such as turbulent HI data cubes. We also include measured power spectra and three-point correlation functions from some of these data. We discuss the importance of open-source statistical and visualization tools for the analysis of turbulence simulations such as those found in CATS.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Adaptive simulations enable computational design of electron beam processing of nanomaterials with supersonic micro-jet precursor

Focused Electron Beam Induced Processing (FEBIP) is a powerful tool for the “direct-write” of nanomaterials with the possibility of atomistic control on suspended 2D material substrates. FEBIP capabilities have been significantly expanded by using a localized jet-based delivery of precursors, and especially when a thermally energized supersonic micro-jet enhances the delivery of mass flux along with controlling the far-from-equilibrium thermodynamic state of adsorbed adatoms. The possibilities of growing nanomaterials with “dialed-in” composition with ultra-high growth rates and aspect ratios up to 100:1 have been demonstrated using the supersonic micro-jet-FEBIP. Bringing this scientific discovery to the level of maturity required for practical applications in additive nanomanufacturing requires simulation tools that are capable of capturing the complex flow physics of micro-jet-substrate interactions that bridge a wide range of flow regimes from the high-density gas micro-jet expanding into the vacuum environment of FEBIP. To address this significant computational challenge, a new approach has been developed and described in this work for a multiscale adaptive DSMC (Direct Simulation Monte Carlo) algorithm to predict the micro-jet gas dynamics in an FEBIP environment, spanning the full range of flow regimes from low Knudsen (Kn) number O(0.01) continuum flow to high Kn of O(10) for the molecular flow within a unified computational framework. Here, the fundamental principles of the adaptive DSMC algorithm are described, its viability as a DSMC technique is demonstrated, and computational improvements for the benchmark cases are discussed in the context of advancing 3D nanofabrication with FEBIP. Ultimately, combining the first principle simulations via adaptive DSMC with the complementary experimental data will enable the creation of the powerful CAD tools for in silico design and optimal operation of micro-jet-FEBIP.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A detailed experimental and modeling comparison of molecular radiative heat loss in a spark-ignition engine

Radiative heat transfer has been chiefly considered negligible in internal combustion engines, except for Diesel engines where soot radiation was recognized as a significant radiative transfer source. Only more recently, detailed simulations have shown that molecular radiation can be substantial as well. In extension to this, molecular radiative heat transfer can reach detectable levels of about 5–10% of the total heat transfer in spark-ignited engines. The broadband radiative nature of the significant emitting molecules, carbon dioxide and water, makes it necessary to address whether radiative trapping plays a substantial role in either total heat loss or energy redistribution within the cylinder. An experimental setup that allows measurements of the infrared emissions at a 2-crank-angle-degree resolution was developed to determine the importance of radiative trapping by carbon dioxide and water. Measurements were conducted in the well-characterized and documented TCC-III engine at the University of Michigan. The engine operated on a stoichiometric propane/air mixture at a speed of 1,300 RPM. Large Eddy Simulations with added line-by-line photon Monte-Carlo simulations of the molecular radiation were an integral part of this study to plan and devise the measurements. Post-processing of the simulation data included an accurate representation of the experimental volume from where infrared signals are collected. Additionally, the experimental data were used for validation of the photon Monte-Carlo simulations. Joint analysis of experimental and simulated spectra allowed quantifying the significance of radiative trapping. Results show the significant role of radiative trapping when predicting radiative heat transfer in the TCC-III engine.

42 ENGINEERING↗
DOE OSTI2021

Modeling fission product diffusion in TRISO fuel particles with BISON

Diffusion of fission products in intact TRISO particles depends on particle geometry, fission product source rates, time, temperature, and temperature-dependent diffusion coefficients. Simulating this diffusion process requires models for source rates and diffusion coefficients, plus computation of the temperature field if not prescribed. In addition, simulation quality depends on discretization of the geometry, appropriate time stepping, and the accuracy of the solution method. In this paper, we explore the simulation of fission product diffusion in TRISO fuel particles using the finite element method via the fuel performance code Bison. Recent material model development has occurred in Bison for each material present in tri-structural isotropic (TRISO) fuel particles: the buffer, inner pyrolytic carbon, silicon carbide, and outer pyrolytic carbon layers, as well as the fuel kernel. Also, new mesh generation and fission product release fraction capabilities have been added. Diffusion capabilities are shown to converge to the correct solution via formal verification tests. A large number of code benchmarking problems are also given, with good results, showing that Bison’s computed release fractions closely match those of other software tools. Finally, a significant validation effort is detailed in which fission product release, measured as part of the AGR-1 capsule experiments, is compared to Bison outputs. Bison outputs compare very well to the experimental data and to PARFUME results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Observation and Estimation of Evapotranspiration from an Irrigated Green Roof in a Rain-Scarce Environment

While the rain-driven evapotranspiration (ET) process has been well-studied in the humid climate, the mixed irrigation and rain-driven ET process is less understood for green roof implementations in dry regions, where empirical observations and model parameterizations are lacking. This paper presents an effort of monitoring and simulating the ET process for an irrigated green roof in a rain-scarce environment. Annual ET rates for three weighing lysimeter test units with non-vegetated, sedums, and grass covers were 2.01, 2.52, and 2.69 mm d -1 , respectively. Simulations based on the three Penman–Monteith equation-derived models achieved accuracy within the reported range of previous studies. Compared to the humid climate, the overestimation of high ET rates by existing models is expected to cause a larger error in dry environments, where the enhanced ET process caused by repeated irrigations overlapped with hot, dry conditions often occurs during summer. The studied sedum species did not show significantly lower ET rates than native species, and could not effectively take advantage of the deep moisture storage. Therefore, native species, instead of the shallow-rooted species commonly recommended in humid climates, might be a better choice for green roofs in rain-scarce environments.

54 ENVIRONMENTAL SCIENCES↗