Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “drops”

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 163 records · Page 9

Solid surface wetting and the deployment of drops in microgravity

The complete or partial deployment of liquid samples in low gravity is primarily influenced by the interfacial properties of the specific liquid and solid materials used because the overwhelming bias of the Earth gravitational acceleration is removed. This study addresses the engineering aspects of injecting and deploying drops of prescribed volume into an acoustic positioning chamber in microgravity. The specific problems of interest are the design, testing, and implementation of injector tips to be used in a simultaneously retracting dual-injector system in the Drop Physics Module microgravity experiment facility. Prior to release, the liquid to be deployed must be retained within a restricted area at the very end of the injectors under dynamic stimuli from the continuous injection flow as well as from the stepped motion of the injectors. The final released drop must have a well determined volume and negligible residual linear or angular momentum. The outcome of Earth-based short-duration low gravity experiments had been the selection of two types of injector tips which were flown as back-up parts. They were successfully utilized during the USML-1 Spacelab mission as the primary tips. The combination of a larger contact surface, liquid pinning with a sharp edge, and selective coating of strategic tip surfaces with a non-wetting compound has allowed a significant increase in the success rate of deployment of simple and compound drops of aqueous solutions of glycerol and silicone oil. The diameter of the samples studied in the Drop Physics Module range between 0.3 and 2.7 cm. The tests conducted on-orbit with a manually operated small device have allowed the calibration of the volume deployed for a few drop sizes. The design for improved tips to be used during the next USML flight is based on these results.

Trinh, E. H.↗

Computations of Drop Collision and Coalescence

Computations of drops collisions, coalescence, and other problems involving drops are presented. The computations are made possible by a finite difference/front tracking technique that allows direct solutions of the Navier-Stokes equations for a multi-fluid system with complex, unsteady internal boundaries. This method has been used to examine the various collision modes for binary collisions of drops of equal size, mixing of two drops of unequal size, behavior of a suspension of drops in linear and parabolic shear flows, and the thermal migration of several drops. The key results from these simulations are reviewed. Extensions of the method to phase change problems and preliminary results for boiling are also shown.

Tryggvason, Gretar↗

System for Manipulating Drops and Bubbles Using Acoustic Radiation Pressure

The manipulation and control of drops of liquid and gas bubbles is achieved using high intensity acoustics in the form of and/or acoustic radiation pressure and acoustic streaming. generated by a controlled wave emission from a transducer. Acoustic radiation pressure is used to deploy or dispense drops into a liquid or a gas or bubbles into a liquid at zero or near zero velocity from the discharge end of a needle such as a syringe needle. Acoustic streaming is useful in manipulating the drop or bubble during or after deployment. Deployment and discharge is achieved by focusing the acoustic radiation pressure on the discharge end of the needle, and passing the acoustic waves through the fluid in the needle. through the needle will itself, or coaxially through the fluid medium surrounding the needle. Alternatively, the acoustic waves can be counter-deployed by focusing on the discharge end of the needle from a transducer axially aligned with the needle, but at a position opposite the needle, to prevent premature deployment of the drop or bubble. The acoustic radiation pressure can also be used for detecting the presence or absence of a drop or a bubble at the tip of a needle or for sensing various physical characteristics of the drop or bubble such as size or density.

Oeftering, Richard C.↗

Flow Visualization in Evaporating Liquid Drops and Measurement of Dynamic Contact Angles and Spreading Rate

A new hybrid optical system, consisting of reflection-refracted shadowgraphy and top-view photography, is used to visualize flow phenomena and simultaneously measure the spreading and instant dynamic contact angle in a volatile-liquid drop on a nontransparent substrate. Thermocapillary convection in the drop, induced by evaporation, and the drop real-time profile data are synchronously recorded by video recording systems. Experimental results obtained from this unique technique clearly reveal that thermocapillary convection strongly affects the spreading process and the characteristics of dynamic contact angle of the drop. Comprehensive information of a sessile drop, including the local contact angle along the periphery, the instability of the three-phase contact line, and the deformation of the drop shape is obtained and analyzed.

Zhang, Neng-Li↗

Spin-Up Instability of a Levitated Molten Drop in MHD-Flow Transition to Turbulence

When an alternating magnetic field interacts with induced eddy currents in a conducting body, there will be a repulsive force between the body and the driving coil system generating the field. This repulsive force is the basis of electromagnetic levitation, which allows containerless processing of different materials. The eddy currents in the conducting body also generate Joule heating. Axial rotation of electromagnetically levitated objects is a common observation in levitation systems and often an undesirable side effect of such experiments on 1-g and -g. There have been recent efforts to use magnetic damping and suppress this tendency of body rotation. The first report of rotation in EML drops was attributed to a slight asymmetry of the shape and location of the levitation coils could change the axis and speed of rotation. Other theories of sample rotation include a frequency difference in the traveling electromagnetic waves and a phase difference in two different applied fields of the same frequency. All of these different mechanisms share the following characteristics: the torque is small, constant for constant field strength, and very weakly dependent on the sample's temperature and phase (solid or liquid). During experiments on the MSL-1 (First Microgravity Science Laboratory) mission of the Space Shuttle (STS-83 and STS-94, April and July 1997), a droplet of palladium-silicon alloy was electromagnetically levitated for viscosity measurements. For the non-deforming droplet, the resultant MHD flow inside the drop is inferred from motion of impurities on the surface. These observations indicate formation of a pair of co-rotating toroidal flow structures inside the spheroidal levitated drop that undergo secondary flow instabilities. As rise in the fluid temperature rises, the viscosity falls and the internal flow accelerates and becomes oscillatory; and beyond a point in the experiments, the surface impurities exhibit non-coherent chaotic motion signifying emergence of turbulence inside the drop. In this work, a background of these set of observations will be given followed by a presentation of our results on the digital particle tracking analysis that has been performed on a number of available videos. The analysis indicates that the levitated drop attains a constant rotational speed during the melting phase and formation of the co-rotating axi-symmetric laminar toroidal structures. However, the rate of axial rotation increases dramatically during the deformation of the toroidal structures anti their breakup into chaotic entities. This new data suggests an interaction between the flow inside the levitated molten drop and the driving coils in the experiments. Possible mechanisms for this interaction will be reviewed. The data will also be used to make an assessment of existing theories on droplet rotation.

Abedian, B.↗

Solid Surface Wetting and the Deployment of Drops in Microgravity

The complete or partial deployment of liquid samples in low gravity is primarily influenced by the interfacial properties of the specific liquid and solid materials used because the overwhelming bias of the Earth gravitational acceleration is removed. This study addresses the engineering aspects of injecting and deploying drops of prescribed volume into an acoustic positioning chamber in microgravity. The specific problems of interest are the design, testing, and implementation of injector tips to be used in a simuttaneously retracting dual-injector system used in the Drop Physics Module microgravity experiment facility. Prior to release, the liquid to be deployed must be retained within a restricted area at the very end of the injectors even under dynamic stimuli due to continuous injection flow as well as to the stepped motion of the injectors, and the final released drop must have a well determined volume as well as negligible residual linear or angular momentum from the deployment process. The outcome of Earthbased short-duration low gravity experiments had been the selection of two types of injector tips which were flown as back-up parts and were successfully utilized during the USML-1 Spacelab mission. The combination of a larger contact surface, liquid pinning with a sharp edge, and selective coating of strategic tip surfaces with a non-wetting compound has allowed a significant increase in the success rate of deployment of simple and compound drops of aqueous solutions of glycerol and silicone oil. The diameter of the samples studied in the Drop Physics Module ranged between 0.3 and 2.7 cm. The tests conducted onsrbit with a manually operated small device have allowed the calibration of the volume deployed for a few drop sizes. The design for improved tips to be used during the next USML flight is based on these results.

Trinh, E. H.↗

Experiments on the Motion of Drops on a Horizontal Solid Surface due to a Wettability Gradient

Results from experiments performed on the motion of drops of tetraethylene glycol in a wettability gradient present on a silicon surface are reported and compared with predictions from a recently developed theoretical model. The gradient in wettability was formed by exposing strips cut from a silicon wafer to decyltrichlorosiland vapors. Video images of the drops captured during the experiments were subsequently analyzed for drop size and velocity as functions of position along the gradient. In separate experiments on the same strips, the static contact angle formed by small drops was measured and used to obtain the local wettability gradient to which a drop is subjected. The velocity of the drops was found to be a strong function of position along the gradient. A quasi-steady theoretical model that balances the local hydrodynamic resistance with the local driving force generally describes the observations; possible reasons for the remaining discrepancies are discussed. It is shown that a model in which the driving force is reduced to accomodate the hysteresis effect inferred from the data is able to remove most of the discrepancy between the observed and predicted velocities.

Moumen, Nadjoua↗

Hanging drop crystal growth apparatus and method

An apparatus (10) is constructed having a cylindrical enclosure (16) within which a disc-shaped wicking element (18) is positioned. A well or recess (22) is cut into an upper side (24) of this wicking element, and a glass cover plate or slip (28) having a protein drop disposed thereon is sealably positioned on the wicking element (18), with drop (12) being positioned over well or recess (22). A flow of control fluid is generated by a programmable gradient former (16), with this control fluid having a vapor pressure that is selectively variable. This flow of control fluid is coupled to the wicking element (18) where control fluid vapor diffusing from walls (26) of the recess (22) is exposed to the drop (12), forming a vapor pressure gradient between the drop (12) and the control fluid vapor. Initially, this gradient is adjusted to draw solvent from the drop (12) at a relatively high rate, and as the critical supersaturation point is approached (the point at which crystal nucleation occurs), the gradient is reduced to more slowly draw solvent from the drop (12). This allows discrete protein molecules more time to orient themselves into an ordered crystalline lattice, producing protein crystals which, when processed by X-ray crystallography, possess a high degree of resolution.

Carter, Daniel C.↗

Multi-terrain Vertical Drop Tests of a Composite Fuselage Section

A 5-ft-diameter composite fuselage section was retrofitted with four identical blocks of deployable honeycomb energy absorber and crash tested on two different surfaces: soft soil, and water. The drop tests were conducted at the 70-ft. drop tower at the Landing and Impact Research (LandIR) Facility of NASA Langley. Water drop tests were performed into a 15-ft-diameter pool of water that was approximately 42-in. deep. For the soft soil impact, a 15-ft-square container filled with fine-sifted, unpacked sand was located beneath the drop tower. All drop tests were vertical with a nominally flat attitude with respect to the impact surface. The measured impact velocities were 37.4, and 24.7-fps for soft soil and water, respectively. A fuselage section without energy absorbers was also drop tested onto water to provide a datum for comparison with the test, which included energy absorbers. In order to facilitate this type of comparison and to ensure fuselage survivability for the no-energy-absorber case, the velocity of the water impact tests was restricted to 25-fps nominal. While all tests described in this paper were limited to vertical impact velocities, the implications and design challenges of utilizing external energy absorbers during combined forward and vertical impact velocities are discussed. The design, testing and selection of a honeycomb cover, which was required in soft surface and water impacts to transmit the load into the honeycomb cell walls, is also presented.

Sotirios, Kellas↗

Multi-Terrain Vertical Drop Tests of a Composite Fuselage Section

A 5-ft-diameter composite fuselage section was retrofitted with four identical blocks of deployable honeycomb energy absorber and crash tested on two different surfaces: soft soil, and water. The drop tests were conducted at the 70-ft. drop tower at the Landing and Impact Research (LandIR) Facility of NASA Langley. Water drop tests were performed into a 15-ft-diameter pool of water that was approximately 42-in. deep. For the soft soil impact, a 15-ft-square container filled with fine-sifted, unpacked sand was located beneath the drop tower. All drop tests were vertical with a nominally flat attitude with respect to the impact surface. The measured impact velocities were 37.4, and 24.7-fps for soft soil and water, respectively. A fuselage section without energy absorbers was also drop tested onto water to provide a datum for comparison with the test, which included energy absorbers. In order to facilitate this type of comparison and to ensure fuselage survivability for the no-energy-absorber case, the velocity of the water impact tests was restricted to 25-fps nominal. While all tests described in this paper were limited to vertical impact velocities, the implications and design challenges of utilizing external energy absorbers during combined forward and vertical impact velocities are discussed. The design, testing and selection of a honeycomb cover, which was required in soft surface and water impacts to transmit the load into the honeycomb cell walls, is also presented.

Kellas, Sotiris↗

Part 2 of a Computational Study of a Drop-Laden Mixing Layer

This second of three reports on a computational study of a mixing layer laden with evaporating liquid drops presents the evaluation of Large Eddy Simulation (LES) models. The LES models were evaluated on an existing database that had been generated using Direct Numerical Simulation (DNS). The DNS method and the database are described in the first report of this series, Part 1 of a Computational Study of a Drop-Laden Mixing Layer (NPO-30719), NASA Tech Briefs, Vol. 28, No.7 (July 2004), page 59. The LES equations, which are derived by applying a spatial filter to the DNS set, govern the evolution of the larger scales of the flow and can therefore be solved on a coarser grid. Consistent with the reduction in grid points, the DNS drops would be represented by fewer drops, called computational drops in the LES context. The LES equations contain terms that cannot be directly computed on the coarser grid and that must instead be modeled. Two types of models are necessary: (1) those for the filtered source terms representing the effects of drops on the filtered flow field and (2) those for the sub-grid scale (SGS) fluxes arising from filtering the convective terms in the DNS equations. All of the filtered-sourceterm models that were developed were found to overestimate the filtered source terms. For modeling the SGS fluxes, constant-coefficient Smagorinsky, gradient, and scale-similarity models were assessed and calibrated on the DNS database. The Smagorinsky model correlated poorly with the SGS fluxes, whereas the gradient and scale-similarity models were well correlated with the SGS quantities that they represented.

Okongo, Nora↗

Construction of 3D MHD pressure drop correlation and flow characterization in the contraction region of a fusion blanket manifold

Inlet and outlet manifolds are typical components of liquid metal (LM) blanket designs of a fusion power reactor to be used to distribute the LM flow into breeding channels and collect it at the exit of the blanket. High pressure loss in the magnetohydrodynamic (MHD) flows featuring abrupt geometrical changes is one of the main feasibility issues of such designs. Recently, optimization studies were conducted to construct 3D MHD pressure drop correlations for a LM flow in an electrically insulating manifold with gradual expansion. Here, the 3D computational approach developed in that study is applied to the outlet manifold featuring gradual contraction. A systematic analysis was performed with a total number of 135 flow cases computed with COMSOL Multiphysics for Hartmann numbers 1000 < Ha < 10,000, Reynolds numbers 100 < Re < 12,000, and contraction angles 45° < θ < 75° for a fixed contraction ratio of 4. The effects of Ha, Re and θ on the flow recirculation, development length and the total pressure drop were carefully examined. A linear regression analysis was used to determine the power rule of pressure drop coefficient k related to Ha and Re, demonstrating a good match with the Ludford layer theory. Eventually, a correlation for the 3D MHD pressure drop coefficient was constructed as a function of Ha, Re and θ. Further, the results were compared against the inlet manifold. It was found that the flow in the inlet manifold exhibits larger recirculation zones. In the investigated range of Ha, Re and θ, the pressure drop coefficient k of the LM MHD flow in the gradual contraction is only slightly lower (< 8 %) than that in the gradual expansion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Computing the viscous effect in early-time drop impact dynamics

The impact of a liquid drop on a solid surface involves many intertwined physical effects, and is influenced by drop velocity, surface tension, ambient pressure and liquid viscosity, among others. Experiments by Kolinski et al. ( Phys. Rev. Lett. , vol. 112, no. 13, 2014 b , p. 134501) show that the liquid–air interface begins to deviate away from the solid surface even before contact. They found that the lift-off of the interface starts at a critical time that scales with the square root of the kinematic viscosity of the liquid. To understand this, we study the approach of a liquid drop towards a solid surface in the presence of an intervening gas layer. We take a numerical approach to solve the Navier–Stokes equations for the liquid, coupled to the compressible lubrication equations for the gas, in two dimensions. With this approach, we recover the experimentally captured early time effect of liquid viscosity on the drop impact, but our results show that lift-off time and liquid kinematic viscosity have a more complex dependence than the square-root scaling relationship. We also predict the effect of interfacial tension at the liquid–gas interface on the drop impact, showing that it mediates the lift-off behaviour.

42 ENGINEERING↗

Pressure drop in a prototypical 3D magnetohydrodynamic flow across contraction of a fusion blanket manifold

Predicting 3D magnetohydrodynamic (MHD) pressure losses associated with liquid metal flows in complex geometry ducts is required to design breeding blankets for fusion reactors. Of such components, manifolds exhibit major pressure losses. Recent correlations [T. Rhodeset al. Magnetohydrodynamic pressure drop and flow balancing of liquid metal flow in a prototypic fusion blanket manifold. Phys. Fluids. 2018; 30: 057101.] demonstrate promise for predicting pressure drops in electrically insulated inlet manifolds. In the present work, the extend to which these correlations can be applied to outlet manifolds, which feature sudden contractions, is investigated in both viscous-electromagnetic (VE) and inertial-electromagnetic (IE) regimes for Reynolds numbers 50 < Re < 1500, Hartmann numbers 2500 < Ha < 5475, and expansion/contraction ratio 10. Numerical computations have shown that the MHD pressure drops in the eutectic lead-lithium (PbLi) flows in contractions are almost identical to those in the flows featuring expansion with slightly higher magnitudes in the contraction cases. The small discripancy in the MHD pressure drop between contractions and expansions (<8%) suggests that the earlier obtained correlations for the 3D MHD pressure drop in a duct flow with a sudden expansion can also be applied to flows with a sudden contraction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evaluation and optimization of flow boiling frictional pressure drop correlations using the data from traditional and next-generation refrigerants in a micro-fin tube

This study presents an experimental evaluation and optimization of flow boiling frictional pressure drop correlations for conventional and next-generation refrigerants in a horizontal micro-fin tube, with particular emphasis on the newly emerging refrigerant blends R-454C and R-455A, for which pressure-drop data in enhanced tubes remain limited. Experiments were conducted with R-410A, R-454C, R-455A, R-134a, R-1234yf, and R-1234ze(E) in a copper micro-fin tube with an inner diameter of 8.468 mm, over mass fluxes ranging from 100 to 300 kg/(m²·s) depending on the refrigerant, and evaporation temperatures of 7, 12, and 14 °C. Frictional pressure gradients were determined from measured total pressure drops after subtracting acceleration pressure drop, and the resulting database was used to assess four existing models: Kuo and Wang (1996), Cavallini et al. (1997), Goto et al. (2001), and Diani et al. (2014). The measured frictional pressure gradient increased with vapor quality and mass flux for all refrigerants and increased further at lower evaporation temperatures, with the overall trend strongly related to liquid viscosity. Among the four correlations, the Goto et al. (2001) model provided the best overall agreement with the measured data before optimization. To further improve prediction accuracy, the Kuo and Wang (1996) and Goto et al. (2001) models were optimized using the complete experimental database. After optimization, both models reduced the overall mean absolute deviation to below 15%, while the optimized Goto et al. (2001) model maintained the best and most consistent overall performance. The results provide new pressure-drop data for next-generation refrigerants and demonstrate that parameter optimization can significantly enhance the applicability of existing micro-fin-tube correlations.

Hu, Yifeng [ORNL] (ORCID:0000000242875185)↗

Density drop at the divertor target in the prototype material plasma exposure eXperiment (Proto-MPEX)

The steady-state linear device “Material Plasma Exposure eXperiment” (MPEX) is currently under construction at Oak Ridge National Laboratory with the goal of enabling Plasma-Material Interaction studies at future fusion reactor relevant plasma conditions. In this work, a newly in-house developed hybrid Particle-In-Cell code-PICOS++ is applied to understand the experimental results obtained from the prototype of MPEX referred to as the “Proto-MPEX” during its helicon-only and helicon with ion cyclotron resonance heating (ICRH) experiments. This study explains the physics of the experimentally observed plasma density-drop at the divertor target in Proto-MPEX device during ICRH. In contrast to previous work on ICRH in MPEX, this study demonstrates that the mirror force plays a central role in the Proto-MPEX plasma transport during ICRH, which has new features not previously explored. Force balance analyses reveal that the temperature anisotropy produced by ICRH leads to a significant increase in the mirror force downstream of the resonance where the magnetic field is diverging. This force accelerates ions toward the target and leads to a drop in plasma density to ensure conservation of particle flux. Simulations with ICRH where the magnetic field divergence downstream of the resonance has been removed, do not produce plasma acceleration nor density drop at the target despite efficient ion heating at the resonance. Moreover, simulation results demonstrate that for a given ICRH power, lowering the source rate produces ions with increased perpendicular energy which interact with the mirror force to produce higher plasma acceleration which increases the strength of the density-drop at the target. The strength of the density drop appears to reach an asymptotic limit at a certain threshold ICRH power. Simulations show that this threshold power increases with increasing particle source rate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

AN ESTIMATE OF SPENT NUCLEAR FUEL MECHANICAL LOADS IN THE GENERAL 30 CM PACKAGE DROP SCENARIO

The US Department of Energy Spent Fuel and Waste Science and Technology (SFWST) program is performing research to determine the mechanical loading conditions applied to spent nuclear fuel (SNF) during normal conditions of transport to inform mechanical tests of SNF and close an important knowledge gap related to the practical disposition of SNF in the US. Researchers at Pacific Northwest National Laboratory (PNNL) have completed an extensive finite element study to characterize and estimate the potential mechanical loads on SNF during a hypothetical 30 cm drop of an SNF transportation package. This modeling study is validated with test data collected by the SFWST program during a physical test campaign that included one-third scale package drop tests and full-scale single fuel assembly drop tests. The test campaign was led by Sandia National Laboratories (SNL) and included international collaboration with Equipos Nucleares S.A, S.M.E (ENSA) and Bundesanstalt für Materialforschung und -prüfung (BAM). The key contribution of the modeling study is to go beyond the limitations of the limited number of physical tests to estimate the impact response to variations in impact angle, initial gap conditions, fuel assembly design, burnup and other parameters that affect the mechanical loads. The methodology of this study included a classic parametric study to calculate the impact response of highly detailed fuel assemblies over many combinations of parameters. Models of a 17x17 pressurized water reactor fuel assembly and a generic 10x10 boiling water reactor fuel assembly were both used in this study to cover the major fuel assembly types in the US inventory. Over 2,000 impact responses were calculated. The results of the parametric study were evaluated using traditional methods and basic statistics. The results were also used to construct a damage model using multiple nonlinear regression techniques to predict the mechanical loads over the full range of all input parameters. The damage model was found to work very well for all impact angle cases where the cask came to rest on its side. It was concluded that end drop cases where the cask remained vertical (instead of tipping over onto its side) were not sufficiently characterized by the current set of parametric study cases to include in the damage model, but it was not a priority to fully investigate that range because the highest mechanical loads were observed in the broader range of side impact cases. This modeling work provides sufficient insight into the mechanical loads on SNF during a hypothetical 30 cm package drop that, when considered along with the physical test data collected by the SNL-led team, the SFSWT program can consider the knowledge gap closed.

Klymyshyn, Nicholas A.↗

Factors That Influence Variability in Stress-Drop Measurements Using Spectral Decomposition and Spectral-Ratio Methods for the 2019 Ridgecrest Earthquake Sequence

Stress drop is a fundamental parameter related to earthquake source physics, but is hard to measure accurately. To better understand how different factors influence stress-drop measurements, we compare two different methods using the Ridgecrest stress-drop validation data set: spectral decomposition (SD) and spectral ratio (SR), each with different processing options. Here, we also examine the influence of spectral complexity on source parameter measurement. Applying the SD method, we find that frequency bandwidth and time-window length could influence spectral magnitude calibration, while depth-dependent attenuation is important to correctly map stress-drop variations. For the SR method, we find that the selected source model has limited influence on the measurements; however, the Boatwright model tends to produce smaller standard deviation and larger magnitude dependence than the Brune model. Variance reduction threshold, frequency bandwidth, and time-window length, if chosen within an appropriate parameter range, have limited influence on source parameter measurement. For both methods, wave type, attenuation correction, and spectral complexity strongly influence the result. The scale factor that quantifies the magnitude dependence of stress drop show large variations with different processing options, and earthquakes with complex source spectra deviating from the Brune-type source models tend to have larger scale factor than earthquakes without complexity. Based on these detailed comparisons, we make a few specific suggestions for data processing workflows that could help future studies of source parameters and interpretations.

58 GEOSCIENCES↗