Search NASA⌕ Search

SEARCH · Search NASA

Results for “Two Phase 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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Newton trust-region methods with primary variable switching for simulating high temperature multiphase porous media flow

Coupling multiphase flow with energy transport due to high temperature heat sources introduces significant new challenges since boiling and condensation processes can lead to dry-out conditions with subsequent re-wetting. The transition between two-phase and single-phase behavior can require changes to the primary dependent variables adding discontinuities as well as extending constitutive nonlinear relations to extreme physical conditions. Practical simulations of large-scale engineered domains lead to Jacobian systems with a very large number of unknowns that must be solved efficiently using iterative methods in parallel on high-performance computers. Performance assessment of potential nuclear repositories, carbon sequestration sites and geothermal reservoirs can require numerous Monte-Carlo simulations to explore uncertainty in material properties, boundary conditions, and failure scenarios. Due to the numerical challenges, standard NR iteration may not converge over the range of required simulations and require more sophisticated optimization method like trust-region. In this study, we use the open-source simulator PFLOTRAN for the important practical problem of the safety assessment of future nuclear waste repositories in the U.S. DOE geologic disposal safety assessment Framework. The simulator applies the PETSc parallel framework and a backward Euler, finite volume discretization. We demonstrate failure of the conventional NR method and the success of trust-region modifications to Newton’s method for a series of test problems of increasing complexity. Trust-region methods essentially modify the Newton step size and direction under some circumstances where the standard NR iteration can cause the solution to diverge or oscillate. Furthermore, we show how the Newton Trust-Region method can be adapted for Primary Variable Switching (PVS) when the multiphase state changes due to boiling or condensation. The simulations with high-temperature heat sources which led to extreme nonlinear processes with many state changes in the domain did not converge with NR, but they do complete successfully with the trust-region methods modified for PVS. This implementation effectively decreased weeks of simulation time needing manual adjustments to complete a simulation down to a day. Finally, we show the strong scalability of the methods on a single node and multiple nodes in an HPC cluster.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

MERIDIONAL FLOW IN THE SOLAR CONVECTION ZONE. II. HELIOSEISMIC INVERSIONS OF GONG DATA

Meridional flow is thought to play a very important role in the dynamics of the solar convection zone; however, because of its relatively small amplitude, precisely measuring it poses a significant challenge. Here we present a complete time–distance helioseismic analysis of about 2 years of ground-based Global Oscillation Network Group (GONG) Doppler data to retrieve the meridional circulation profile for modest latitudes in an attempt to corroborate results from other studies. We use an empirical correction to the travel times due to an unknown center-to-limb systematic effect. The helioseismic inversion procedure is first tested and reasonably validated on artificial data from a large-scale numerical simulation followed by a test to broadly recover the solar differential rotation found from global seismology. From GONG data, we measure poleward photospheric flows at all latitudes with properties that are comparable with earlier studies and a shallow equatorward flow about 65 Mm beneath the surface, in agreement with recent findings from Helioseismic and Magnetic Imager (HMI) data. No strong evidence of multiple circulation cells in depth or latitude is found, yet the whole phase space has not yet been explored. Tests of mass flux conservation are then carried out on the inferred GONG and HMI flows and compared to a fiducial numerical baseline from models, and we find that the continuity equation is poorly satisfied. While the two disparate data sets do give similar results for about the outer 15% of the interior radius, the total inverted circulation pattern appears to be unphysical in terms of mass conservation when interpreted over modest time scales. We can likely attribute this to both the influence of realization noise and subtle effects in the data and measurement procedure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Effect of Aluminum Alloying on the Hot Deformation Behavior of Nano-bainite Bearing Steel

Interest in using aluminum in nano-bainite steel, especially for high-carbon bearing steel, is gradually growing. In this study, GCr15SiMo and GCr15SiMoAl steels are introduced to investigate the effect of Al alloying on the hot deformation behavior of bearing steel. Results show that the addition of Al not only notably increases the flow stress of steel due to the strong strengthening effect of Al on austenite phase, but also accelerates the strain-softening rates for its increasing effect on stacking fault energy. Al alloying also increases the activation energy of deformation. Two constitutive equations with an accuracy of higher than 0.99 are proposed. The constructed processing maps show the expanded instability regions for GCr15SiMoAl steel as compared with GCr15SiMo steel. This finding is consistent with the occurrence of cracking on the GCr15SiMoAl specimens, revealing that Al alloying reduces the high-temperature plasticity of the bearing steel. On the contrary, GCr15SiMoAl steel possesses smaller grain size than GCr15SiMo steel, manifesting the positive effect of Al on bearing steel. Attention should be focused on the hot working process of bearing steel with Al.

36 MATERIALS SCIENCE↗
DOE OSTI2015

A moist Boussinesq shallow water equations set for testing atmospheric models

The shallow water equations have long been used as an initial test for numerical methods applied to atmospheric models with the test suite of Williamson et al. being used extensively for validating new schemes and assessing their accuracy. However the lack of physics forcing within this simplified framework often requires numerical techniques to be reworked when applied to fully three dimensional models. In this paper a novel two-dimensional shallow water equations system that retains moist processes is derived. This system is derived from three-dimensional Boussinesq approximation of the hydrostatic Euler equations where, unlike the classical shallow water set, we allow the density to vary slightly with temperature. This results in extra (or buoyancy) terms for the momentum equations, through which a two-way moist-physics dynamics feedback is achieved. The temperature and moisture variables are advected as separate tracers with sources that interact with the mean-flow through a simplified yet realistic bulk moist-thermodynamic phase-change model. This moist shallow water system provides a unique tool to assess the usually complex and highly non-linear dynamics–physics interactions in atmospheric models in a simple yet realistic way. The full non-linear shallow water equations are solved numerically on several case studies and the results suggest quite realistic interaction between the dynamics and physics and in particular the generation of cloud and rain. - Highlights: • Novel shallow water equations which retains moist processes are derived from the three-dimensional hydrostatic Boussinesq equations. • The new shallow water set can be seen as a more general one, where the classical equations are a special case of these equations. • This moist shallow water system naturally allows a feedback mechanism from the moist physics increments to the momentum via buoyancy. • Like full models, temperature and moistures are advected as tracers that interact through a simplified yet realistic phase-change model. • This model is a unique tool to test numerical methods for atmospheric models, and physics–dynamics coupling, in a very realistic and simple way.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A physics-based model for frost buildup under turbulent flow using direct numerical simulations

We present a new model for frost buildup under turbulent (and laminar) flow using direct numerical simulations. The physical model consists of two layers, the air and the frost. The air layer is fully resolved and consists of solving for the velocity, temperature, and vapor mass fraction fields. The frost layer thickness is resolved using conservation of mass and energy. Both phases are dynamically coupled using the immersed boundary method. Three-dimensional simulations are conducted in an open-channel configuration. A number of challenges need to be overcome to make these simulations feasible. First, to enforce far-field conditions of zero gradient and prescribed mean temperature and humidity, a source term is added to the energy and transport equations in the flow solver. Second, the mean frost thickness is subtracted after each time step to ensure a constant mean flow thickness and level of turbulence in the numerical domain. Third, a slow-time acceleration approach, which accelerates the frost buildup by a predetermined factor, is employed to bridge the gap between the fast turbulent and slow frost buildup time scales. Finally, a frost densification scheme is used to overcome the difficulties of vertically varying frost properties. The model is validated by comparing the frost thickness and frost thickness buildup rate over a period of one hour from a cooled flat plate experiment. As a result, both quantities compare favorably with experiments.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2019

Dominant two-fluid magnetohydrodynamic instabilities in CFETR upgrade phase-I scenario in presence of perfect conducting wall

Here, the dominant two-fluid magnetohydrodynamic (MHD) stabilities of China Fusion Engineering Test Reactor upgrade phase-I baseline equilibria in presence of perfect conducting wall have been evaluated using the initial value code NIMROD. Our MHD simulations show that for equilibrium with q axis < 1, high-n edge localized modes are slightly stabilized by toroidal shear flow, and n = 1 mode is found to be unstable in the core region with the characteristic of sawtooth instability. For q axis > 2 equilibrium, all modes are found to be unstable but edge localized both in single-fluid and two-fluid MHD models. Specifically for n ≤ 13 modes, the two-fluid MHD model gives a slightly higher growth rate than the single-fluid MHD model, while for n > 13 modes, this trend becomes the opposite, which means the two-fluid MHD model is needed for high-n mode analysis. In addition, n = 1–10 modes are found to be more unstable with increasing wall position and eventually their growth rates approach values in the no-wall limit.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Multiscale characterization of phase change materials for building thermal energy storage applications

Phase change materials (PCMs) store and release large amounts of thermal energy because of their high latent energy storage capacity. However, long-term cyclic stability, supercooling and performance-scalability are some of the major challenges for their use in building thermal energy storage (TES) applications. Here, in this study, we present a comprehensive multiscale characterization of two commercially available organic PCMs, Puretemp 18 and Puretemp 23. At the microscale, differential scanning calorimetry (DSC) was used to characterize phase change temperature, specific heat, and latent heat. At the mesoscale, a heat flow meter apparatus (HFMA), following the ASTM C1784 standard, was employed to measure the phase change temperature, specific heat, and latent heat properties. A comparative analysis of latent heat as a function of temperature was conducted by integrating the DSC and HFMA results. At the macroscale, the thermal performance and cyclic stability of the TES system was evaluated using Puretemp 23. The TES system consisted of a finned tube heat exchanger with a storage volume of 0.0189 m 3 (5 gal), which represents a compact, real-world TES solution suitable for building energy storage. The results showed consistent thermal stability of the PCM over 200 cycles, and the supercooling temperature remained within 0.2 °C, which was not detected in smaller-scale characterization methods. Additionally, the macroscale testing methodology of the PCM revealed that the TES is able to charge and discharge stored latent energy within 2 h under a temperature differential of 16.67 °C measured between the inlet water temperature and the phase transition temperature of the PCM. The proposed multiscale PCM characterization method provides a systematic basis for comparing important thermal storage properties while also investigating the scalability, reliability and integration challenges in large scale TES applications.

Latent heat↗
DOE OSTI2020

Development and validation of the cavitation-induced erosion risk assessment tool

This work presents the development of a cavitation-induced erosion risk assessment (CIERA) tool that links multiphase flow simulation predictions with the progress towards material erosion. To develop a robust erosion modeling tool, the cavitation and erosion predictions for pressurized diesel fuel flow within a channel geometry were validated over a range of Reynolds and cavitation number conditions in two different aluminum channel geometries, one featuring a rounded inlet corner and the other with a sharp inlet corner. The multiphase flow development within the channel was modeled using a compressible mixture model, where phase change was represented with the homogeneous relaxation model and the turbulent flow evolution was modeled using a dynamic structure approach for Large Eddy Simulations. To improve representation of the incubation period before material rupture over existing approaches, a physics-based metric was derived based on the cumulative energy absorbed by the solid material from repeated hydrodynamic impacts. When the average peak pressure was related to the incubation period, the incubation period and its sensitivity to changes in flow conditions were found to be overpredicted. In contrast, predictions from CIERA provided a more accurate means to qualitatively and quantitatively predict the influence of flow conditions on the incubation period before material erosion. When the predicted stored energy was related to the solid material properties to estimate the incubation period, multiphase flow simulations demonstrated accurate representation of the sensitivity of erosion severity to changes in flow conditions. The use of CIERA led to quantitative agreement of the predicted incubation period within 5% of the experimentally measured incubation period.

42 ENGINEERING↗
DOE OSTI2015

Richtmyer-Meshkov evolution under steady shock conditions in the high-energy-density regime

This work presents direct experimental evidence of long-predicted nonlinear aspects of the Richtmyer-Meshkov (RM) process, in which new modes first arise from the coupling of initially-present modes, and in which shorter-wavelength modes are eventually overtaken by longer-wavelength modes. This is accomplished using a technique we developed employing a long driving laser pulse to create a strong (Mach ~ 8) shock across a well-characterized material interface seeded by a two-mode sinusoidal perturbation. Furthermore, this technique further permits the shock to be sustained, without decay of the high-energy-density flow conditions, long enough for the system to evolve into the nonlinear phase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Cross-Plane Acquisitions in Electrical Capacitance Volume Tomography

Electrical capacitance tomography (ECT) is a widely used imaging modality to image two-dimensional cross sections of multiphase flows. Recent developments in electrical capacitance volume tomography (ECVT) have made it possible to directly obtain volumetric images from measured data. An ECVT is instrumental for obtaining accurate phase hold up information and velocity information that are needed for the optimization of certain flow processes. However, compared to ECT, the high correlation between the measurements in ECVT exacerbates the ill-conditioning of the associated image reconstruction problem. Previous studies have suggested that neglecting mutual capacitance data between the ECVT electrodes located at cross-planes that are well separated along the sensor axis can be done without significantly affecting the reconstructed image. In addition, this may help constrain the ill-conditioning of the reconstruction problem. Here, we examine in detail and quantify the effect of reduced cross-plane acquisition strategies for optimizing the image reconstruction and constraining the ill-conditioning of typical ECVT settings.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2019

Asymmetries in adaptive optics point spread functions

An explanation for the origin of asymmetry along the preferential axis of the point spread function (PSF) of an AO system is developed. When phase errors from high-altitude turbulence scintillate due to Fresnel propagation, wavefront amplitude errors may be spatially offset from residual phase errors. These correlated errors appear as asymmetry in the image plane under the Fraunhofer condition. In an analytic model with an open-loop AO system, the strength of the asymmetry is calculated for a single mode of phase aberration, which generalizes to two dimensions under a Fourier decomposition of the complex illumination. Other parameters included are the spatial offset of the AO correction, which is the wind velocity in the frozen flow regime multiplied by the effective AO time delay and propagation distance or altitude of the turbulent layer. In this model, the asymmetry is strongest when the wind is slow and nearest to the coronagraphic mask when the turbulent layer is far away, such as when the telescope is pointing low toward the horizon. A great emphasis is made about the fact that the brighter asymmetric lobe of the PSF points in the opposite direction as the wind, which is consistent analytically with the clarification that the image plane electric field distribution is actually the inverse Fourier transform of the aperture plane. Validation of this understanding is made with observations taken from the Gemini Planet Imager, as well as being reproducible in end-to-end AO simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Two dimensional topology optimization of heat exchangers with the volume fraction method

We perform a comparison study of two topology optimizations methods applied to the design of a two fluids heat exchanger modeled with a coupled thermal-flow problem. The flow follows an isothermal and incompressible Stokes-Brinkman equation and the heat transfer is governed by a convection-diffusion equation without internal generation and high Peclet number. To keep the two fluid phases separated, we solve two Stokes-Brinkman equations, where the Brinkman term models the other phase as a solid. These two velocity fields are then fed to the heat transfer equation. Our goal is to maximize the enthalpy at the cold outlet while constraining the pressure drop. We first solve the design modeling the solid and fluid phases with a volume fraction variable. A SIMP-like penalization in the Brinkman term drives the optimization to a discrete design. The cost and constraint function derivatives are calculated with the library pyadjoint and the optimization is performed by IPOPT. We present optimized designs in two dimensions and discuss the influence of the parameters.

Beck, VictorA.↗
DOE OSTI2017

Generation Mechanisms of Quasi-parallel and Quasi-circular Flare Ribbons in a Confined Flare

We analyze a confined multiple-ribbon M2.1 flare (SOL2015-01-29T11:42) that originated from a fan-spine coronal magnetic field configuration, within active region NOAA 12268. The observed ribbons form in two steps. First, two primary ribbons form at the main flare site, followed by the formation of secondary ribbons at remote locations. We observe a number of plasma flows at extreme-ultraviolet temperatures during the early phase of the flare (as early as 15 minutes before the onset) propagating toward the formation site of the secondary ribbons. The secondary ribbon formation is co-temporal with the arrival of the pre-flare generated plasma flows. The primary ribbons are co-spatial with Ramaty High Energy Spectroscopic Imager ( RHESSI ) hard X-ray sources, whereas no enhanced X-ray emission is detected at the secondary ribbon sites. The (E)UV emission, associated with the secondary ribbons, peaks ∼1 minute after the last RHESSI hard X-ray enhancement. A nonlinear force-free model of the coronal magnetic field reveals that the secondary flare ribbons are not directly connected to the primary ribbons, but to regions nearby. Detailed analysis suggests that the secondary brightenings are produced due to dissipation of kinetic energy of the plasma flows (heating due to compression), and not due to non-thermal particles accelerated by magnetic reconnection, as is the case for the primary ribbons.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Modeling and Optimization of a Rotating Packed Bed Contactor with a Tetraamine-Appended Metal−Organic Framework for CO 2 Capture

A potential contactor technology for sorbent-based CO 2 capture is the rotating packed bed that contains separate sections for continuous adsorption and desorption. A heat exchanger can be embedded to remove heat in the adsorption section and add heat in the desorption section. In this work, we develop a two-dimensional (2D) model of a rotating packed bed for use in CO 2 capture applications. Mass and energy balances for the model are developed based on a Ljungström-type air preheater, which accounts for the counter-current axial flow of gas phases in separate sections of the bed and the rotation of a solid sorbent, which cycles between adsorption and desorption sections. The sorbent used for this analysis is the tetraamine-appended metal−organic framework Mg 2 (dobpdc)(3−4− 3), chosen for its stability and affinity for CO 2 at low partial pressures, such as those from a natural gas power plant source. An optimization problem is solved that considers the trade-off between maximizing the productivity of the bed and minimizing energy consumption. Maximum productivity and minimum energy are found to be 8.53 kg/h/m 3 and 3.84 MJ/kg, respectively, when these objectives are optimized independently. It is observed that the flue gas pressure and bed rotational speed are the desired operating variables to vary for model-based design of experiments to reduce uncertainty in parameter estimation, as these two variables yielded the most information content based on the Fisher information matrix.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2015

Performance of a Cross-Flow Humidifier with a High Flux Water Vapor Transport Membrane

Water vapor transport (WVT) flux across a composite membrane that consists of a very thin perfluorosulfonic acid (PFSA) ionomer layer sandwiched between two expanded polytetrafluoroethylene (PTFE) microporous layers is investigated. Static and dynamic tests are conducted to measure WVT flux for different composite structures; a transport model shows that the underlying individual resistances for water diffusion in the gas phase and microporous and ionomer layers and for interfacial kinetics of water uptake at the ionomer surface are equally important under different conditions. A finite-difference model is formulated to determine water transport in a full-scale (2-m2 active membrane area) planar cross-flow humidifier module assembled using pleats of the optimized composite membrane. In agreement with the experimental data, the modeled WVT flux in the module increases at higher inlet relative humidity (RH) of the wet stream and at lower pressures, but the mass transfer effectiveness is higher at higher pressures. The model indicates that the WVT flux is highest under conditions that maintain the wet stream at close to 100% RH while preventing the dry stream from becoming saturated. The overall water transport is determined by the gradient in RH of the wet and dry streams but is also affected by vapor diffusion in the gas layer and the microporous layer.

combined thermal and mass transfer↗
DOE OSTI2020

Development of PFLOTRAN Transport Capability for Use in the Waste Isolation Pilot Plant Performance Assessment - 20545

Waste Isolation Pilot Plant (WIPP) performance assessment (PA) calculations estimate the probability of radionuclide release from the repository to the land surface and across the land withdrawal boundary for a regulatory period of 10,000 years after facility closure. Simulations of flow and transport in the repository and the surrounding Salado Formation are foundational to the PA. Because proposed additional waste emplacement panels would result in an asymmetric repository layout, the US Department of Energy (DOE) is preparing to transition to use of a three-dimensional (3-D) model domain for simulation of flow and transport instead of the two-dimensional (2-D) flared grid domain currently used. DOE has charged Sandia National Laboratories with developing the capability necessary to simulate processes affecting flow and transport in the WIPP in PFLOTRAN, an open-source massively parallel multi-phase flow and reactive transport code. The new flow and transport capabilities developed in PFLOTRAN incorporate WIPP-specific process models and will replace the 2-D simulators (BRAGFLO and NUTS) that are currently utilized for Salado flow and transport calculations in WIPP PA. The focus of this paper is on the development of a new Nuclear Waste Transport (NWT) mode in PFLOTRAN that has all of the capabilities necessary for Salado transport simulations, including the ability to handle complete dry-out (100% gas saturation) of arbitrary cells in the model domain, radionuclide mass conservation at step changes in porosity associated with borehole intrusion, and the ability to calculate fluxes on a flared grid. The new PFLOTRAN transport capability and a suite of verification tests were designed around a list of functional requirements for WIPP PA calculations. (authors)

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

Self-diffusion in the non-Newtonian regime of shearing liquid crystal model systems based on the Gay-Berne potential

The self-diffusion coefficients of nematic phases of various model systems consisting of regular convex calamitic and discotic ellipsoids and non-convex bodies such as bent-core molecules and soft ellipsoid strings have been obtained as functions of the shear rate in a shear flow. Then the self-diffusion coefficient is a second rank tensor with three different diagonal components and two off-diagonal components. These coefficients were found to be determined by a combination of two mechanisms, which previously have been found to govern the self-diffusion of shearing isotropic liquids, namely, (i) shear alignment enhancing the diffusion in the direction parallel to the streamlines and hindering the diffusion in the perpendicular directions and (ii) the distortion of the shell structure in the liquid whereby a molecule more readily can escape from a surrounding shell of nearest neighbors, so that the mobility increases in every direction. Thus, the diffusion parallel to the streamlines always increases with the shear rate since these mechanisms cooperate in this direction. In the perpendicular directions, these mechanisms counteract each other so that the behaviour becomes less regular. In the case of the nematic phases of the calamitic and discotic ellipsoids and of the bent core molecules, mechanism (ii) prevails so that the diffusion coefficients increase. However, the diffusion coefficients of the soft ellipsoid strings decrease in the direction of the velocity gradient because the broadsides of these molecules are oriented perpendicularly to this direction due the shear alignment (i). The cross coupling coefficient relating a gradient of tracer particles in the direction of the velocity gradient and their flow in the direction of the streamlines is negative and rather large, whereas the other coupling coefficient relating a gradient in the direction of the streamlines and a flow in the direction of the velocity gradient is very small.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Simulation of mineral dissolution at the pore scale with evolving fluid-solid interfaces: review of approaches and benchmark problem set

This manuscript presents a benchmark problem for the simulation of single-phase flow, reactive transport, and solid geometry evolution at the pore scale. The problem is organized in three parts that focus on specific aspects: flow and reactive transport (part I), dissolution-driven geometry evolution in two dimensions (part II), and an experimental validation of three-dimensional dissolution-driven geometry evolution (part III). Five codes are used to obtain the solution to this benchmark problem, including Chombo-Crunch, OpenFOAM-DBS, a lattice Boltzman code, Vortex, and dissolFoam. These codes cover a good portion of the wide range of approaches typically employed for solving pore-scale problems in the literature, including discretization methods, characterization of the fluid-solid interfaces, and methods to move these interfaces as a result of fluid-solid reactions. A short review of these approaches is given in relation to selected published studies. Results from the simulations performed by the five codes show remarkable agreement both quantitatively--based on upscaled parameters such as surface area, solid volume, and effective reaction rate--and qualitatively--based on comparisons of shape evolution. This outcome is especially notable given the disparity of approaches used by the codes. Therefore, these results establish a strong benchmark for the validation and testing of pore-scale codes developed for the simulation of flow and reactive transport with evolving geometries. They also underscore the significant advances seen in the last decade in tools and approaches for simulating this type of problem.

58 GEOSCIENCES↗