Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solids 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.

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

At least 1,117 records · Page 62

The Role of Local Inhomogeneities on Dendrite Growth in LLZO-Based Solid Electrolytes

The majority of the ceramic solid electrolytes (LLZO, LATP) demonstrate polycrystalline grain/grain-boundary (G/GB) microstructure. Higher lithium (Li) concentration and lower mechanical stiffness result in current focusing at the GBs. Growth of Li dendrites through local inhomogeneities and subsequent short circuit of the cell is a major concern. Recent studies have revealed that bulk Li metal is a viscoplastic material that has low (~0.3 MPa) and high (~1.0 MPa) yield strength during deformation at smaller and larger rates of strain, respectively. It has been argued that during deposition at smaller current densities, due to its lower yield strength, Li metal should demonstrate plastic flow against stiff ceramic electrolytes, and Li dendrites will be prevented from penetrating through solid electrolytes. In this manuscript, a multiscale modeling framework has been developed for predicting properties of GBs and the bulk of ceramic electrolytes using atomistic calculations for input to mesoscale models. Using the parameters obtained from the atomistic simulations, the mesoscale model reveals that, given enough time, even at low charge rates, lithium dendrites can grow through the GBs of LLZO. The present multiscale model results also provide information regarding the dendrite growth velocity through LLZO.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Flow Testing of Corrugated Metal Flexhoses to Evaluate Flow-Induced Vibration and Stiffness

Corrugated metal flexhoses are used to supply fluid routing where straight rigid pipes cannot meet the design requirements due to vibrations, thermal expansion, or motion. These types of hoses are found in a wide range of applications at KSC (Exploration Ground Systems Program in particular) with various fluid commodities such as fuel, oxidizers, coolant, and cryogenics. However, flexibility of design comes at a price: due to increased levels of turbulence generated by the hose geometry, the necessary supply pressure must increase to achieve the same flow rate. Furthermore, the convolutes of the flexhoses interact with the flow field to generate areas of flow separation which leads to a phenomenon known as vortex shedding. When the frequency of vortex shedding couples with the natural frequency of the hose, this can be detrimental and cause premature failure. There are a limited amount of software that can correctly model the coupled fluid structure interactions(FSI) between the solid and fluid physics. This project aims to develop a computer model that can predict Flow-Induced Vibration (FIV) on flexhoses in one of the atypical configurations found in the State-of-the-Art (SOTA) standard for hoses in an angulated state. Our research seeks to extend the literature and expand NASA’s FIV standard. The computational rigor and resources available in the Apollo Era did not allow for flow coupling between fluids and structure interfaces to evaluate FIV of the flexhose. Our team believes that modern techniques and methods of evaluation should be used to reevaluate and extend the database of FIV and stiffness properties of flexhoses to assess the risk of failure. The core performance criteria is a computer model that is able to predict the FIV frequency within 10%.

Jared F. Congiardo↗

Opposed-Flow Flame Spreading in Reduced Gravity

Experimental results obtained in drop towers and in Space Shuttle based experiments coupled with modelling efforts are beginning to provide information that is allowing an understanding to be developed of the physics of opposed-flow flame spread at reduced gravity where the spread rate and flow velocity are comparable and of the role played by radiative and diffusive processes in flame spreading in microgravity. Here we describe one Space Shuttle based experiment on flame spreading in a quiescent environment, the Solid Surface Combustion Experiment, SSCE, one planned microgravity experiment on flame spreading in a radiatively-controlled, forced opposing flow environment, the Diffusive and Radiative Transport in Fires Experiment, DARTFire, modelling efforts to support these experiments, and some results obtained to date.

Altenkirch, Robert A.↗

A numerical study of the effects of wind tunnel wall proximity on an airfoil model

A procedure was developed for modeling wind tunnel flows using computational fluid dynamics. Using this method, a numerical study was undertaken to explore the effects of solid wind tunnel wall proximity and Reynolds number on a two-dimensional airfoil model at low speed. Wind tunnel walls are located at varying wind tunnel height to airfoil chord ratios and the results are compared with freestream flow in the absence of wind tunnel walls. Discrepancies between the constrained and unconstrained flows can be attributed to the presence of the walls. Results are for a Mach Number of 0.25 at angles of attack through stall. A typical wind tunnel Reynolds number of 1,200,000 and full-scale flight Reynolds number of 6,000,000 were investigated. At this low Mach number, wind tunnel wall corrections to Mach number and angle of attack are supported. Reynolds number effects are seen to be a consideration in wind tunnel testing and wall interference correction methods. An unstructured grid Navier-Stokes code is used with a Baldwin-Lomax turbulence model. The numerical method is described since unstructured flow solvers present several difficulties and fundamental differences from structured grid codes, especially in the area of turbulence modeling and grid generation.

Potsdam, Mark↗

Investigation of boron powder flow rates on real-time wall

The limit of boron flow rates for real-time conditioning of the first walls has been systematically investigated in the Experimental Advanced Superconducting Tokamak (EAST) with a full metal wall. Initially, solid boron injection demonstrated effective control over carbon impurities and deuterium recycling on the basis of pre-discharge boronization. A minimum flow rate, identified between 1.0 mg/s and 2.0 mg/s, was necessary for actively improving wall conditions under specific plasma operating scenarios, with this effect progressively enhancing as boron flow rates increased. Additionally, a maximum flow rate, estimated between 3.5 mg/s and 8.0 mg/s, was identified for these plasma conditions. When boron flow rates exceeded this maximum, boron-induced fueling effects influenced the plasma line-averaged density, and at excessively high flow rates, plasma disruption was observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Aerodynamic sound of flow past an airfoil

The long term objective of this project is to develop a computational method for predicting the noise of turbulence-airfoil interactions, particularly at the trailing edge. We seek to obtain the energy-containing features of the turbulent boundary layers and the near-wake using Navier-Stokes Simulation (LES or DNS), and then to calculate the far-field acoustic characteristics by means of acoustic analogy theories, using the simulation data as acoustic source functions. Two distinct types of noise can be emitted from airfoil trailing edges. The first, a tonal or narrowband sound caused by vortex shedding, is normally associated with blunt trailing edges, high angles of attack, or laminar flow airfoils. The second source is of broadband nature arising from the aeroacoustic scattering of turbulent eddies by the trailing edge. Due to its importance to airframe noise, rotor and propeller noise, etc., trailing edge noise has been the subject of extensive theoretical (e.g. Crighton & Leppington 1971; Howe 1978) as well as experimental investigations (e.g. Brooks & Hodgson 1981; Blake & Gershfeld 1988). A number of challenges exist concerning acoustic analogy based noise computations. These include the elimination of spurious sound caused by vortices crossing permeable computational boundaries in the wake, the treatment of noncompact source regions, and the accurate description of wave reflection by the solid surface and scattering near the edge. In addition, accurate turbulence statistics in the flow field are required for the evaluation of acoustic source functions. Major efforts to date have been focused on the first two challenges. To this end, a paradigm problem of laminar vortex shedding, generated by a two dimensional, uniform stream past a NACA0012 airfoil, is used to address the relevant numerical issues. Under the low Mach number approximation, the near-field flow quantities are obtained by solving the incompressible Navier-Stokes equations numerically at chord Reynolds number of 104. The far-field noise is computed using Curle's extension to the Lighthill analogy (Curle 1955). An effective method for separating the physical noise source from spurious boundary contributions is developed. This allows an accurate evaluation of the Reynolds stress volume quadrupoles, in addition to the more readily computable surface dipoles due to the unsteady lift and drag. The effect of noncompact source distribution on the far-field sound is assessed using an efficient integration scheme for the Curle integral, with full account of retarded-time variations. The numerical results confirm in quantitative terms that the far-field sound is dominated by the surface pressure dipoles at low Mach number. The techniques developed are applicable to a wide range of flows, including jets and mixing layers, where the Reynolds stress quadrupoles play a prominent or even dominant role in the overall sound generation.

Wang, Meng↗

Analysis of Off-design Performance of a 16-stage Axial-flow Compressor with Various Blade Modifications

The over-all performance of a 16-stage axial-flow compressor was determined with various stator-blade resettings and a reduction in solidity of the rotor blades in the last three stages. It was shown that little control over the sudden change in slope of the surge-limit line at intermediate speeds was obtained with the blade modifications attempted, except that some change in speed at which the change in slope occurred could be effected by stator-blade resettings. Interstage data indicated that the severe surge limit at intermediate speeds was caused by stall of the inlet stage, which, because of stage interaction effects, resulted in a simultaneous decrease in performance of the following five or six stages. Stage data are presented which indicate the flow and pressure-ratio range over which each stage is required to operate at compressor speeds from 50 to 100 percent of design speed.

COMPRESSORS - AXIAL-FLOW↗

Remobilizing the Interfaces of Thermocapillary Driven Bubbles Retarded by the Adsorption of a Surfactant Impurity on the Bubble Surface

Thermocapillary migration is a method for moving bubbles in space in the absence of buoyancy. A temperature gradient is applied to the continuous phase in which a bubble is situated, and the applied gradient impressed on the bubble surface causes one pole of the drop to be cooler than the opposite pole. As the surface tension is a decreasing function of temperature, the cooler pole pulls at the warmer pole, creating a flow which propels the bubble in the direction of the warmer fluid. A major impediment to the practical use of thermocapillarity to direct the movement of bubbles in space is the fact that surfactant impurities which are unavoidably present in the continuous phase can significantly reduce the migration velocity. A surfactant impurity adsorbed onto the bubble interface is swept to the trailing end of the bubble. When bulk concentrations are low (which is the case with an impurity), diffusion of surfactant to the front end is slow relative to convection, and surfactant collects at the back end of the bubble. Collection at the back lowers the surface tension relative to the front end setting up a reverse tension gradient. For buoyancy driven bubble motions in the absence of a thermocapillarity, the tension gradient opposes the surface flow, and reduces the surface and terminal velocities (the interface becomes more solid-like). When thermocapillary forces are present, the reverse tension gradient set up by the surfactant accumulation reduces the temperature tension gradient, and decreases to near zero the thermocapillary velocity. The objective of our research is to develop a method for enhancing the thermocapillary migration of bubbles which have been retarded by the adsorption onto the bubble surface of a surfactant impurity, Our remobilization theory proposes to use surfactant molecules which kinetically rapidly exchange between the bulk and the surface and are at high bulk concentrations. Because the remobilizing surfactant is present at much higher concentrations than the impurity, it adsorbs to the bubble much faster than the impurity when the bubble is formed, and thereby prevents the impurity from adsorbing onto the surface. In addition the rapid kinetic exchange and high bulk concentration maintain a saturated surface with a uniform surface concentrations. This prevents retarding surface tension gradients and keeps the velocity high. In our first report last year, we detailed experimental results which verified the theory of remobilization in ground based experiments in which the steady velocity of rising bubbles was measured in a continuous phase consisting of a glycerol/water mixture containing a polyethylene glycol surfactant C12E6 (CH3(CH2)11(OCH2CH2)6OH). In our report this year, we detail our efforts to describe theoretically the remobilization observed. We construct a model in which a bubble rises steadily by buoyancy in a continuous (Newtonian) viscous fluid containing surfactant with a uniform far field bulk concentration. We account for the effects of inertia as well as viscosity in the flow in the continuous phase caused by the bubble motion (order one Reynolds number), and we assume that the bubble shape remains spherical (viscous and inertial forces are smaller than capillary forces, i e. small Weber and capillary numbers). The surfactant distribution is calculated by solving the mass transfer equations including convection and diffusion in the bulk, and finite kinetic exchange the bulk and the surface. Convective effects dominate diffusive mass transfer in the bulk of the liquid (high Peclet numbers) except in a thin boundary layer near the surface. A finite volume method is used to numerically solve the hydrodynamic and mass transfer equations on a staggered grid which accounts specifically for the thin boundary layer. We present the results of the nondimensional drag as a function of the bulk concentration of surfactant for different rates of kinetic exchange, from which we develop criteria for the concentration necessary to develop a prescribed degree of remobilization. The criteria compare favorably with the experimental results.

Palaparthi, Ravi↗

Predictive models for moving contact line flows

Modeling flows with moving contact lines poses the formidable challenge that the usual assumptions of Newtonian fluid and no-slip condition give rise to a well-known singularity. This singularity prevents one from satisfying the contact angle condition to compute the shape of the fluid-fluid interface, a crucial calculation without which design parameters such as the pressure drop needed to move an immiscible 2-fluid system through a solid matrix cannot be evaluated. Some progress has been made for low Capillary number spreading flows. Combining experimental measurements of fluid-fluid interfaces very near the moving contact line with an analytical expression for the interface shape, we can determine a parameter that forms a boundary condition for the macroscopic interface shape when Ca much les than l. This parameter, which plays the role of an "apparent" or macroscopic dynamic contact angle, is shown by the theory to depend on the system geometry through the macroscopic length scale. This theoretically established dependence on geometry allows this parameter to be "transferable" from the geometry of the measurement to any other geometry involving the same material system. Unfortunately this prediction of the theory cannot be tested on Earth.

Rame, Enrique↗

Computational Simulation of Semispan Wings in Wind Tunnels

The computational modelling of experiments, with the end aim of providing sufficiently accurate simulations to assess and improve turbulence models is described. Solid wall tunnels are the only tunnels in which the boundary conditions can, in principle, be known exactly. The modelling of the tunnel walls for transonic flows requires the accurate modelling of the viscous displacement effects on the tunnel walls. This paper describes the modelling of semispan wing experiments in solid wall tunnels, with the tunnel walls modelled as inviscid walls, and with all 4 walls modelled viscously. The effect of the viscous effects is discussed, as well as the feasibility of modelling these effects in an inviscid, apriori manner.

Olsen, Mike↗

Space Launch System Base Heating Test: Environments and Base Flow Physics

The NASA Space Launch System (SLS) vehicle is composed of four RS-25 liquid oxygen-hydrogen rocket engines in the core-stage and two 5-segment solid rocket boosters and as a result six hot supersonic plumes interact within the aft section of the vehicle during flight. Due to the complex nature of rocket plume-induced flows within the launch vehicle base during ascent and a new vehicle configuration, sub-scale wind tunnel testing is required to reduce SLS base convective environment uncertainty and design risk levels. This hot-fire test program was conducted at the CUBRC Large Energy National Shock (LENS) II short-duration test facility to simulate flight from altitudes of 50 kft to 210 kft. The test program is a challenging and innovative effort that has not been attempted in 40+ years for a NASA vehicle. This paper discusses the various trends of base convective heat flux and pressure as a function of altitude at various locations within the core-stage and booster base regions of the two-percent SLS wind tunnel model. In-depth understanding of the base flow physics is presented using the test data, infrared high-speed imaging and theory. The normalized test design environments are compared to various NASA semi-empirical numerical models to determine exceedance and conservatism of the flight scaled test-derived base design environments. Brief discussion of thermal impact to the launch vehicle base components is also presented.

Mehta, Manish↗

Optimization of the moderators in the STS preliminary design

This report details the results for an optimization of the dimensions of the moderators in the preliminary design of the Spallation Neutron Source Second Target Station (STS). This study uses the optimization algorithms of Dakota and an unstructured mesh model for the moderators in MCNP. More details on the unstructured mesh model and the automated mesh generation can be found in [3]. Parallel to this effort, the same moderator geometries have been optimized using a constructive solid geometry (CSG) MCNP model. More details on this model and its results can be found in [4]. Three optimal designs are selected for each moderator: one that is optimized for maximum peak brightness, one for maximum time-integrated brightness, and one for a combination of peak and time-integrated brightness. The backbone of the optimization work flow is provided by Dakota. For each set of design parameters requested by Dakota, a new solid geometry is automatically built in Creo and SpaceClaim, and subsequently exported to Attila4MC to generate an unstructured mesh geometry for MCNP. After the MCNP calculation is finished, the objective function (e.g., brightness metric) is returned to Dakota. After the new design has been evaluated, a result-file is written, and Dakota proposes the next set of design parameters to be evaluated. The loop continues until a specified convergence criterion has been met. The design parameters of the cylindrical (upper) moderator include the hydrogen radius, the premoderator thickness (top, bottom, radial), the beryllium radius and the horizontal position of the moderator. The crucial design choice is the hydrogen radius. A radius of 62 mm is shown to provide the maximum time-integrated brightness. The maximum peak brightness occurs with a radius of 40 mm. A combined (middle) design, which balances peak and time-integrated brightnesses, is obtained with a hydrogen radius of 50 mm. The premoderator thicknesses and the beryllium radius are slightly larger in the design optimized for time-integrated brightness than in the design optimized for peak brightness. The sensitivity to these two parameters is relatively small close to the optimal configurations. The hydrogen vessel and vacuum vessel wall thicknesses are dependent on the radius of the liquid hydrogen due to structural integrity requirements. The increased wall thicknesses for larger vessels significantly penalize the time-integrated brightness, with the maximum obtainable value reduced by more than 10% relative to earlier studies which used fixed vessel wall thicknesses. The impact of the variable wall thicknesses is much less for the peak brightness and combined brightness designs. The design parameters of the tube (lower) moderator selected for the optimization are the tube length, the annular premoderator thickness, the beryllium radius and the horizontal position of the moderator. The tube length is the crucial parameter and is chosen large (210 mm) and small (125 mm) in the designs optimized for time-integrated and peak brightness respectively. A combined optimal design has a tube length of 170 mm. The premoderator thickness and the beryllium radius are chosen larger in the design optimized for time-integrated brightness.

42 ENGINEERING↗

High Solids Performance Testing in a Scaled TSCR System

The Tank Side Cesium Removal (TSCR) project is a technology demonstration that will pretreat Hanford tank waste supernatant in support of the Direct Feed Low-Activity Waste (DFLAW) mission. The TSCR system employs two key separation technologies: dead-end filtration (DEF) and ion exchange (IX) using crystalline silicotitanate (CST) media. DEF will be used to remove undissolved solids from tank waste to protect the functionality of the IX columns and the IX system will remove Cs-137 from tank waste. The separation technologies (DEF and IX) used in TSCR are technically mature and have also been successfully deployed at the Savannah River site in a similar facility known as the Tank Closure Cesium Removal (TCCR) system. While testing with simulants and real waste has been successfully performed under conditions expected during the initial operation of TSCR, test data is absent for assessing off normal high solids loading that may be in the TSCR waste feed. Normal TSCR treatment operations are expected to handle wastes with solids content on the order of 200 ppm, and off normal solids loading could be much larger than the nominal level. The testing program described in this report was conducted to understand the consequence of operating the TSCR system at elevated solids loadings up to the high-solids limit of 15,000 ppm [i.e., 1.5-wt%] identified in the TSCR design basis. Although the system is not required to make throughput above the nominal solids loading, the testing was intended to provide important information related to potential off normal operations. At off normal levels near the high-solids limit, there are potential implications for TSCR performance in the areas of throughput, DEF pressure drop, filter backflush frequency, and IX column pressure drop. In addition, intrusion of solids into the IX column was postulated to impact the Cs-137 loading behavior by promoting channeling or flow maldistribution in the column; since the magnitude of the postulated effect was unknown, assessing it was also of interest. The testing was performed using representative waste simulants and a prototypic, integrated TSCR system designed and assembled specifically to conduct the high solids performance assessment. Overall, the scaled TSCR testing demonstrated that full-scale unit operations can succeed in fulfilling their processing objectives in the presence of solids up to 3,000 ppm, but there are potential performance challenges to filtration operations at solids loadings as low as ~500 ppm. The severity of the challenge is likely to be dependent on the type and size distribution of solids, of which the current testing only examined a single type and size distribution. To provide some flexibility for future full-scale operations, the results of the testing suggest two possible risk reduction strategies that can be implemented without any changes in TSCR design or configuration. One option would be to enact an administrative limit on the solids loading to protect TSCR from feeds that are likely to require a high DEF swap frequency. Another option is to permit operation of the DEFs at differential pressures greater than 2 psid before swapping filters. The selection of a higher differential pressure target is not anticipated to adversely impact DEF backflushing efficacy and would reduce both swap frequency and the amount of waste sent to AP-108.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Particle-fluid interactions for flow measurements

Study has been made of the motion of single particle and of group of particles, emphasizing solid particles in gaseous fluid. Velocities of fluid and particle are compared for several conditions of physical interest. Mean velocity and velocity fluctuations are calculated for single particle, and some consideration is given to multiparticle systems.

Berman, N. S.↗

Solid metabolic waste transport and stowage investigation

The basic Waste Collection System (WCS) design under consideration utilized air flow to separate the stool from the WCS user and to transport the fecal material to a slinger device for subsequent deposition on a storage bowel. The major parameters governing stool separation and transport were found to be the area of the air inlet orifices, the configuration of the air inlet orifice and the transport air flow. Separation force and transport velocity of the stool were studied. The developed inlet orifice configuration was found to be an effective design for providing fecal separation and transport. Simulated urine tests and female user tests in zero gravity established air flow rates between 0.08 and 0.25 cu sm/min (3 and 9 scfm) as satisfactory for entrapment, containment and transport of urine using an urinal. The investigation of air drying of fecal material as a substitute for vacuum drying in a WCS breadboard system showed that using baseline conditions anticipated for the shuttle cabin ambient atmosphere, flow rates of 0.14 cu sm/min (5 cfm) were adequate for drying and maintaining biological stability of the fecal material.

Burt, R. A.↗

Striped Electrodes for Solid-Electrolyte Cells

Striped thick-film platinum electrodes help insure lower overall cell resistance by permitting free flow of gases in gaps between stripes. Thickfilm stripes are also easier to fabricate than porous thin-film electrodes that cover entire surface. Possible applications for improved cells include oxygen production from carbon dioxide, extraction of oxygen from air, small fluidic pumping, sewage treatment, and fuel cells.

Richter, R.↗

BETO 2021 Peer Review - Rational Design of Robust Reactor Feeding Systems for Heterogeneous Cellulosic and Agricultural Wastes Based on Biomass Quality Characteristics

Consistent and reliable preprocessing, conveyance, and reactor in-feed systems, particularly for low-cost waste feedstocks, remains a major technical challenge for the emerging Bioeconomy. By identifying critical biomass attributes and connecting them to flow and conversion behavior, science-driven system designs can address these often-overlooked solids handling challenges. The Wonderful Company (TWC) is the world's largest almond and pistachio grower, generating 250,000 dry tons/year of waste material including hulls, shells, and wood (>5 million tons/year industry wide in the U.S.). This project seeks to turn this environmental and economic liability into a sustainable and profitable resource, targeting conversion via gasification to syngas for electricity and bio-char, by addressing related material handling and feeding challenges. Based on optimized preprocessing strategies, bulk material flow, and thermal conversion properties, an overall system design will be developed. The methodology will be tested with FCIC's benchmark loblolly pine residues, to demonstrate the robustness of the overall approach and provide insight and guidance for future systems. The project will culminate in extended field trials to demonstrate an improved continuous feeding system with a commercial biomass-to-electricity gasifier vendor, and an economic analysis demonstrating a reduction in electricity production costs by maximizing on-stream time while minimizing preprocessing and CapEx costs.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Pore-Scale Modeling of Mineral Growth and Nucleation in Reactive Flow

A fundamental understanding of mineral precipitation kinetics relies largely on microscopic observations of the dynamics of mineral surfaces exposed to supersaturated solutions. Deconvolution of tightly bound transport, surface reaction, and crystal nucleation phenomena still remains one of the main challenges. Particularly, the influence of these processes on texture and morphology of mineral precipitate remains unclear. This study presents a coupling of pore-scale reactive transport modeling with the Arbitrary Lagrangian-Eulerian approach for tracking evolution of explicit solid interface during mineral precipitation. It incorporates a heterogeneous nucleation mechanism according to Classical Nucleation Theory which can be turned “on” or “off.” This approach allows us to demonstrate the role of nucleation on precipitate texture with a focus at micrometer scale. In this work precipitate formation is modeled on a 10 micrometer radius particle in reactive flow. The evolution of explicit interface accounts for the surface curvature which is crucial at this scale in the regime of emerging instabilities. The results illustrate how the surface reaction and reactive fluid flow affect the shape of precipitate on a solid particle. It is shown that nucleation promotes the formation of irregularly shaped precipitate and diminishes the effect of the flow on the asymmetry of precipitation around the particle. The observed differences in precipitate structure are expected to be an important benchmark for reaction-driven precipitation in natural environments.

58 GEOSCIENCES↗