Search NASA⌕ Search

SEARCH · Search NASA

Results for “deformation function”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Calibration of RAFM Micromechanical Model for Creep Using Bayesian Optimization for Functional Output

A Bayesian optimization procedure is presented for calibrating a multimechanism micromechanical model for creep to experimental data of F82H steel. Reduced activation ferritic martensitic (RAFM) steels based on Fe(8–9)%Cr are the most promising candidates for some fusion reactor structures. Although there are indications that RAFM steel could be viable for fusion applications at temperatures up to 600°C, the maximum operating temperature will be determined by the creep properties of the structural material and the breeder material compatibility with the structural material. Due to the relative paucity of available creep data on F82H steel compared to other alloys such as Grade 91 steel, micromechanical models are sought for simulating creep based on relevant deformation mechanisms. As a point of departure, this work recalibrates a model form that was previously proposed for Grade 91 steel to match creep curves for F82H steel. Due to the large number of parameters (9) and cost of the nonlinear simulations, an automated approach for tuning the parameters is pursued using a recently developed Bayesian optimization for functional output (BOFO) framework (Huang et al., 2021, “Bayesian optimization of functional output in inverse problems,” Optim. Eng., 22, pp. 2553–2574). Incorporating extensions such as batch sequencing and weighted experimental load cases into BOFO, a reasonably small error between experimental and simulated creep curves at two load levels is achieved in a reasonable number of iterations. In conclusion, validation with an additional creep curve provides confidence in the fitted parameters obtained from the automated calibration procedure to describe the creep behavior of F82H steel.

42 ENGINEERING↗
DOE OSTI2018

Processing-induced strain glass states in a Ni 49.5 Ti 50.5 shape memory alloy

Shape memory alloys (SMAs) represent a revolutionary and innovative class of active materials which can provide potential solutions to many of today's engineering problems due to their compact form, high energy densities, and multifunctional capabilities. While many applications in the biomedical, aerospace, and automotive industries have already been investigated and realized for Nickel-Titanium (NiTi) based SMAs, the effects of restricting the ferroelastic transformation to nanosized domains is not well understood and the potential remains untapped. In binary NiTi, the martensitic transformation, which is characterized by long-range strain ordering (LRO), can be replaced with a strain glass transition, which consists of an LRO parent phase and a short-range strain ordered glassy phase. Such alloys have been named strain glass alloys (SGAs) due to the fact that they exhibit a glassy state which results from compositionally- or processing-induced strain. While SGAs do not exhibit a stress-free, temperature-induced macroscopic phase change, they still exhibit the strain recovery and actuation capabilities intrinsic to near equiatomic NiTi and other SMAs. It has been shown in the available literature that certain compositions, for example 51.5 at. % Nickel in binary NiTi, can create a strain glass; however, these compositionally-induced NiTi SGAs generally have transformation temperatures below 173 K and this will restrict their practical applications. In the present study, a new method for producing a strain glass phase in Ti-rich NiTi through sufficient plastic deformation via cold work is reported; the resulting SGA exhibits a temperature-induced ferroelastic recovery above room temperature. Additionally, the macroscopic actuation capabilities are improved when compared to both compositionally-induced SGAs and the base material due to the increased functional stresses of the SGA. To better understand the transition from an SMA to an SGA, Ni49.5Ti50.5 (at. %) rods were processed to several degrees of cold work and characterized via scanning and transmission electron microscopy, differential scanning calorimetry, thermomechanical testing, and synchrotron radiation x-ray diffraction. The experimental results indicate that twin size decreases with additional cold work and, around 45% thickness reduction, stress-free thermal cycling no longer results in a measurable phase transformation; however, mechanically-induced phase transformation is still possible, where fully recoverable strains in these SGAs were observed to be above 4.5% when loaded at room temperature and recovered at 150 °C.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Low-Lying Excited States of Linear All- Trans Polyenes: Insights from Analytic Gradient and Nonadiabatic Coupling Calculations Based on Multireference Configuration Interaction

Polyenes serve as a rigorous test for theoretical models and electronic structure methods, playing a key role in advancing computational and theoretical chemistry. Here, we present a high-level theoretical investigation of linear, all-trans polyenes using energy gradients and nonadiabatic coupling vectors based on an MR-CISD wave function to describe electronic transitions involving the ground state (1 1 A g – ) and three low-lying excited states (2 1 A g – , 1 1 B u + , and 2 1 B u – ) of hexatriene, octatetraene, and decapentaene. This approach enables accurate evaluation of both adiabatic and vertical excitation and emission energies, yielding results in excellent agreement with experiment, as well as locating minima on the crossing seam between adiabatic states. Our results show that vertical excitation energies to the 1 1 B u + state are blue-shifted by 0.2–0.3 eV relative to the experimental absorption maximum, whereas the vertical emission energy from the 2 1 A g – state is red-shifted by ∼0.2 eV relative to the experimental emission maximum. Upon relaxation from the Franck–Condon geometry, the 2 1 A g – state stabilizes by around 1 eV, compared to 0.2–0.3 eV for the 1 1 B u + state. An analysis of the S 1 /S 0 crossing seam in hexatriene shows that its minimum involves asymmetric backbone deformations and provides an efficient channel for ultrafast internal conversion to the ground state, consistent with the absence of detectable fluorescence in this molecule. These results demonstrate the power of analytic gradients and nonadiabatic coupling vectors based on an MR-CISD wave function for accurately characterizing the electronic structure and photophysics of polyenes.

Excited states↗
DOE OSTI2026

Observing neutron alignment at high angular momenta in heavy deformed actinides

An outstanding puzzle in the structure of very heavy nuclei lies in the rotation response of deformed transuranic isotopes. While the rotation alignment of pairs of 𝑖 13/2 protons is ubiquitous and well understood, observing the analogous alignment of the 𝑗 15/2 neutrons has been elusive. In this work we present a new quantitative analysis approach: from their moments of inertia as a function of frequency, we compare the alignments of rotational bands with both “low” and high-𝑗 intrinsic structures in the 𝑁 = 151 isotones of Pu, Cm, Cf, and No, thereby distinguishing between the alignments of neutron versus proton pairs through “blocking” effects. In this way, the alignment of 𝑗 15/2 neutrons can be clearly identified for the first time. The protons and neutrons begin aligning at a very similar frequency. The contribution of 𝑗 15/2 neutrons to the total alignment is found to grow with increasing 𝑍 compared to 𝑖 13/2 protons. To explain these observations, it is necessary to adopt a self-consistent set of single-particle energies corresponding to the nuclear shape that minimizes the mass. The shape description needs to include high-order deformations up to at least 𝛽 6 .

A ≥ 220↗
DOE OSTI2023

Comparison of Shape Optimization Methods for Heat Exchanger Fins Using Computational Fluid Dynamics

Inverse design techniques are one way to leverage advances in 3D printing, artificial intelligence, and computational resources to achieve increased performance of heat exchangers. Two optimization techniques (genetic algorithm and particle swarm) and three geometry representations (binary level set, composite Bézier, and free form deformation) are used to increase the heat transfer and reduce the pressure drop of a heat exchanger fin. After running 210,810 OpenFOAM simulations, results indicate that a significant performance increase of the fin can be realized in less than 48 hrs, allowing for such a process to be integrated in traditional design processes. The best design increased the performance of the objective function, compared to the baseline rectangular geometry by 75%. A custom distributed infrastructure was built, allowing for all methods to reach 95% of the final objective values in a little over 4 hrs, handling 1674 OpenFoam simulations per hour.

42 ENGINEERING↗
DOE OSTI2020

Deformation behavior of annealed Cu64Zr36 metallic glass via molecular dynamics simulations

In the present work, we investigate the isothermal annealing process of Cu 64 Zr 36 metallic glass (MG) by means of molecular dynamics (MD) simulations, and characterize the compression properties under different temperatures and compressive rates. The pair distribution functions, the bond-orientational order, and the cluster-type index method are modeled and analyzed to characterize changes in the structure. Results of the modeling and analysis reveal that the MgZn 2 -type Laves, MgCu 2 -type Laves and five-fold local symmetry structures can be formed during the isothermal holding at 950 K. Study on the compression properties shows that at a fixed strain rate, the yield strength and yield drop all increase with decreasing temperature. At a fixed temperature, however, the increase of the strain rate leads to a noticeable increase in the yield strength at 600 K, but has little effect on yield strength at 10 K. Moreover, the modeling and analysis of the structure at a temperature of 10 K and strain-rate of 5 × 10 7 s -1 demonstrate that the order-disorder transformation initiates the shear band.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Reciprocal swimming at intermediate Reynolds number

In Stokes flow, Purcell's scallop theorem forbids objects with time-reversible (reciprocal) swimming strokes from moving. In the presence of inertia, this restriction is eased and reciprocally deforming bodies can swim. A number of recent works have investigated dimer models that swim reciprocally at intermediate Reynolds numbers Re≈1–1000. These show interesting results (e.g. switches of the swim direction as a function of inertia) but the results vary and seem to be case specific. Here, in this work, we introduce a general model and investigate the behaviour of an asymmetric spherical dimer of oscillating length for small-amplitude motion at intermediate Re. In our analysis we make the important distinction between particle and fluid inertia, both of which need to be considered separately. We asymptotically expand the Navier–Stokes equations in the small-amplitude limit to obtain a system of linear partial differential equations. Using a combination of numerical (finite element) and analytical (reciprocal theorem, method of reflections) methods we solve the system to obtain the dimer's swim speed and show that there are two mechanisms that give rise to motion: boundary conditions (an effective slip velocity) and Reynolds stresses. Each mechanism is driven by two classes of sphere–sphere interactions, between one sphere's motion and (1) the oscillating background flow induced by the other's motion, and (2) a geometric asymmetry induced by the other's presence. We can thus unify and explain behaviours observed in other works. Our results show how sensitive, counterintuitive and rich motility is in the parameter space of finite inertia of particles and fluid.

42 ENGINEERING↗
DOE OSTI2022

Curated Collection of More than 20,000 Experimentally Reported One-Dimensional Metal–Organic Frameworks

Herein, a collection of more than 20,000 experimentally derived crystal structures for metal–organic frameworks (MOFs) that do not have two- or three-dimensional covalently bonded networks has been developed from the materials available at the Cambridge Crystallographic Data Centre. Of these 20,000 1D MOFs, more than 12,000 structures have been verified to be solvent-free and in exact agreement with the stoichiometry of the synthesized materials. More than 10% of the complete data set comprise materials including two or more distinct metals. The band gaps of more than 12,000 1D MOFs have been computed at the density functional theory-generalized gradient approximation level, finding more than 2000 materials that have a zero band gap. Molecular simulations of CH4 adsorption in a small number of 1D MOFs indicated that adsorbate-induced deformation plays a significant role in determining adsorption isotherms in these materials. As a result, methods that have been used previously for high-throughput predictions of molecular adsorption in 3D MOFs are not suitable for 1D MOFs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Thermomechanical conversion in metals: dislocation plasticity model evaluation of the Taylor-Quinney coefficient

Using a partitioned-energy thermodynamic framework which assigns energy to that of atomic configurational stored energy of cold work and kinetic-vibrational, in this study we derive an important constraint on the Taylor-Quinney coefficient, which quantifies the fraction of plastic work that is converted into heat during plastic deformation. Associated with the two energy contributions are two separate temperatures – the ordinary temperature for the thermal energy and the effective temperature for the configurational energy. We show that the Taylor-Quinney coefficient is a function of the thermodynamically defined effective temperature that measures the atomic configurational disorder in the material. Finite-element analysis of recently published experiments on the aluminum alloy 6016-T4 [1], using the thermodynamic dislocation theory (TDT), shows good agreement between theory and experiment for both stress-strain behavior and temporal evolution of the temperature. The simulations include both conductive and convective thermal energy loss during the experiments, and significant thermal gradients exist within the simulation results. Computed values of the differential Taylor-Quinney coefficient are also presented and suggest a value which differs between materials and increases with increasing strain.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Silver nanowires-templated metal oxide for broadband Schottky photodetector

Silver nanowires (AgNWs)-templated transparent metal oxide layer was applied for Si Schottky junction device, which remarked the record fastest photoresponse of 3.4 μs. Self-operating AgNWs-templated Schottky photodetector showed broad wavelength photodetection with high responsivity (42.4 A W{sup −1}) and detectivity (2.75 × 10{sup 15} Jones). AgNWs-templated indium-tin-oxide (ITO) showed band-to-band excitation due to the internal photoemission, resulting in significant carrier collection performances. Functional metal oxide layer was formed by AgNWs-templated from ITO structure. The grown ITO above AgNWs has a cylindrical shape and acts as a thermal protector of AgNWs for high temperature environment without any deformation. We developed thermal stable AgNWs-templated transparent oxide devices and demonstrated the working mechanism of AgNWs-templated Schottky devices. We may propose the high potential of hybrid transparent layer design for various photoelectric applications, including solar cells.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Electrostatically driven pattern formation in mixed charged–neutral multicomponent elastic membranes

Multicomponent crystalline and amorphous elastic shells exhibit heterogeneous surface patterns that provide distinctive functionalities in cellular environments. Such patterning typically arises from the competition between short-range attractive and long-range repulsive interactions in membranes. Here, we demonstrate that the intrinsic competition between electrostatic repulsion and elastic deformation is sufficient to drive spontaneous surface patterning in elastic shells, requiring no additional attractive interactions. Using numerical simulations, we demonstrate pattern formation in mechanically homogeneous membranes with heterogeneous surface charge composition across different topologies, including spheres, discs, and flat periodic membranes. We also examine patterns in crystalline and amorphous shells of coassembled charged and neutral components with different bending rigidities. At low charge fraction, discrete charged surface domains form. At intermediate charge fraction, the competition between electrostatics and elasticity leads to elongated domains (rods) of the charged component, which results in lamellar patterns at nearly equal fraction of the charged and neutral components. At high charge fraction, nanodomains of the neutral component form. Amorphous shells exhibit similar progressions but with disordered structures rather than ordered lamellar patterns. These pattern morphologies are observed in both the closed shells and flat membranes. As salt concentration increases, all patterns coarsen due to the screening of electrostatic interactions.

charge heterogeneity↗
DOE OSTI2022

Magnetic field induced hair structure in the charmonium gluon dissociation

We study the electromagnetic field effect on charmonium gluon dissociation in quark-gluon plasma. With the effective Hamiltonian derived from QCD multipole expansion under an external electromagnetic field, we first solve the two-body Schrödinger equation for a pair of charm quarks with mean field potentials for color and electromagnetic interactions and obtain the charmonium binding energies and wave functions, and then calculate the gluon dissociation cross section and decay width by taking the color electric and magnetic dipole interactions as perturbations above the mean field and employing Fermi’s golden rule. Considering the charmonium deformation in a magnetic field, the discrete Landau energy levels make the dissociation cross section grow hair, and the electric dipole channel is significantly changed, especially for the P -wave states χ c 0 and χ c ± . From our numerical calculation, the magnetic field strength e B = 5 m π 2 already changes the gluon dissociation strongly, which may indicate measurable effects in high-energy nuclear collisions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Optimizing thermal stability and mechanical behavior in segregation-engineered nanocrystalline Al–Ni–Ce alloys: A combinatorial study

The deliberate use of solute enrichment at grain boundaries, otherwise known as segregation engineering, is a promising approach to tailor the properties of interface-dominated materials such as nanocrystalline alloys. Here, the ensuing chemical and structural evolution at grain boundaries can give rise to thermal stability and excellent mechanical properties, but the interplay between enrichment, phase decomposition, grain growth, and mechanical behavior exists in a vast composition and processing space. In this study, a combinatorial synthesis and rapid characterization approach was applied to segregation-engineered nanocrystalline Al–Ni–Ce alloys to assess the evolution of microstructure and resulting mechanical behavior as functions of alloying content and annealing conditions. In addition to the identification of alloys and processing conditions that give rise to exceptional thermal stability, strength retention, and homogeneous plastic flow, we construct combined thermal stability and deformation mechanism maps that demarcate several important regimes of behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The energetics and dynamics of confinement in flexible frameworks and molecular confinement

Porous frameworks form the chemical and structural basis for critical technologies in separations, catalysis, nuclear waste containment and biomedical applications. Hundreds of zeolites and metal organic frameworks (MOFs) have been synthesized, and their ability to separate and store hydrogen, methane and carbon dioxide has been investigated both experimentally and theoretically. Nevertheless, a fundamental and systematic molecular-level understanding of the thermodynamic and structural factors governing the stability and guest-host interactions in these materials lags behind focused studies of specific systems. Because the guest molecules interact with each other and with the host framework, molecular confinement is a finely balanced and complex phenomenon. The ability of the guest molecules to bind and diffuse through the pores is determined by the nature of the host framework which, in turn, responds to the nature and concentration of guest molecules and to pressure and temperature. The work explores how framework flexibility, tailored by structure, composition, temperature and pressure, is a general phenomenon, similar in nature but variable in extent, in both zeolites and MOFs and is part of a free energy landscape in which framework-guest interactions, pressure, and temperature result in changes in framework geometry and, in some cases, phase transitions. These subtle and/or pronounced changes in lattice geometry, energetics, and dynamics can play a decisive role in confinement and in differentiating the binding of molecules of similar size. The free energy landscape created by these structural changes links polymorphism, amorphization, “breathing,” “gate opening” and confinement. Specifically, the generality of such behavior arises from commonalities in lattice dynamics and energetics of frameworks containing a combination of strong rigid bonds and weaker more flexible deformation modes. Identifying and describing these common and collective phenomena is the focus of the research on a selected group of zeolites and MOFs. Structural studies using X-ray and neutron diffraction explore the mechanical functionality of these important materials, specifically how framework materials respond to changes in temperature, pressure and guest loading. These structural studies are combined with calorimetric measurements, using techniques uniquely developed in the participating laboratories, of heats of formation, heat capacities and entropies, and guest-host interactions. The experimental thermodynamic studies are complemented by inelastic neutron scattering studies of the lattice dynamics related both to framework vibrations and to guest-host interactions. The report below summarizes the work done at UC Davis, which is complemented by work done at BYU (Woodfield) and Virginia Tech (Ross).

36 MATERIALS SCIENCE↗
DOE OSTI2022

The energetics and dynamics of confinement in flexible frameworks and molecular confinement

Porous frameworks form the chemical and structural basis for critical technologies in separations, catalysis, nuclear waste containment and biomedical applications. Hundreds of zeolites and metal organic frameworks (MOFs) have been synthesized, and their ability to separate and store hydrogen, methane and carbon dioxide has been investigated both experimentally and theoretically. Nevertheless, a fundamental and systematic molecular-level understanding of the thermodynamic and structural factors governing the stability and guest-host interactions in these materials lags behind focused studies of specific systems. Because the guest molecules interact with each other and with the host framework, molecular confinement is a finely balanced and complex phenomenon. The ability of the guest molecules to bind and diffuse through the pores is determined by the nature of the host framework which, in turn, responds to the nature and concentration of guest molecules and to pressure and temperature. The work explores how framework flexibility, tailored by structure, composition, temperature and pressure, is a general phenomenon, similar in nature but variable in extent, in both zeolites and MOFs and is part of a free energy landscape in which framework-guest interactions, pressure, and temperature result in changes in framework geometry and, in some cases, phase transitions. These subtle and/or pronounced changes in lattice geometry, energetics, and dynamics can play a decisive role in confinement and in differentiating the binding of molecules of similar size. The free energy landscape created by these structural changes links polymorphism, amorphization, “breathing,” “gate opening” and confinement. Specifically, the generality of such behavior arises from commonalities in lattice dynamics and energetics of frameworks containing a combination of strong rigid bonds and weaker more flexible deformation modes. Identifying and describing these common and collective phenomena is the focus of the research on a selected group of zeolites and MOFs. Structural studies using X-ray and neutron diffraction explore the mechanical functionality of these important materials, specifically how framework materials respond to changes in temperature, pressure and guest loading. These structural studies are combined with calorimetric measurements, using techniques uniquely developed in the participating laboratories, of heats of formation, heat capacities and entropies, and guest-host interactions. The experimental thermodynamic studies are complemented by inelastic neutron scattering studies of the lattice dynamics related both to framework vibrations and to guest-host interactions. This research was a collaboration between UC Davis (A. Navrotsky) BYU (B. Woodfield) and Virginia Tech (N. Ross).

36 MATERIALS SCIENCE↗
DOE OSTI2018

Transient Rheological Behavior of Semisolid SEED-Processed 7075 Aluminum Alloys in Rapid Compression

The transient rheological behavior and microstructure evolution of semisolid SEED-processed 7075 aluminum alloys were studied using the rapid compression tests. The effects of the TiB{sub 2} grain refinement on the grain morphology and size of semisolid slurries were investigated. Results indicated that the grain refiner could reduce the grain size and improve the globularity of α-Al grains. The grain-refined alloy can be easily deformed at a wide range of solid contents (0.42 to 0.53 Fs), in which the deformation level appears to be independent from the solid content. Under the transient state, the apparent viscosity decreased with increasing shear rate to a minimum value and followed by an increase as the shear rate decreased. The apparent viscosity of the base alloy exhibited a dependency on the solid content, while the apparent viscosity of the grain-refined alloy in the decreasing or increasing shear rate periods was not substantially influenced by the solid content. The viscosity as a function of applied shear rate can be described using the power law viscosity model. The differences in the flow behavior index (n) and the consistency index (k) for two alloys were discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

In situ study of microstructure evolution and α → ω phase transition in annealed and pre-deformed Zr under hydrostatic loading

The detailed study of the effect of the initial microstructure on its evolution under hydrostatic compression before, during, and after the irreversible α → ω phase transformation and during pressure release in Zr using in situ x-ray diffraction is presented. Two samples were studied: one is plastically pre-deformed Zr with saturated hardness and the other is annealed. Phase transformation α → ω initiates at lower pressure for a pre-deformed sample but for a volume fraction of ω Zr, c > 0.7, a larger volume fraction is observed for the annealed sample. This implies that the proportionality between the athermal resistance to the transformation and the yield strength in the continuum phase transformation theory is invalid; an advanced version of the theory is outlined. Phenomenological plasticity theory under hydrostatic loading is outlined in terms of microstructural parameters, and plastic strain is estimated. During transformation, the first rule is suggested, i.e., the average domain size, microstrain, and dislocation density in ω Zr for c < 0.8 are functions of the volume fraction, c of ω Zr only, which are independent of the plastic strain tensor prior to transformation and pressure. The microstructure is not inherited during phase transformation. Surprisingly, for the annealed sample, the final dislocation density and the average microstrain after pressure release in the ω phase are larger than for the severely pre-deformed sample. The results suggest that an extended experimental basis is required for the predictive models for the combined pressure-induced phase transformations and microstructure evolutions.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Strength anisotropy of shales deformed under uppermost crustal conditions

Conventional triaxial tests were performed on three sets of samples of Tournemire shale along different orientations relative to bedding (0°, 45°, and 90°). Experiments were carried out up to failure at increasing confining pressures ranging from 2.5 to 160 MPa, at strain rates ranging between 3 × 10 -7 s -1 and 3 × 10 -5 s -1 . This allowed us to determine the entire anisotropic elastic compliance matrix as a function of confining pressure. Results show that the orientation of principal stress relative to bedding plays an important role on the brittle strength, with 45° orientation being the weakest. We fit our results with a wing crack micromechanical model and an anisotropic fracture toughness. We found low values of internal friction coefficient and apparent friction coefficient in agreement with friction coefficient of clay minerals (between 0.2 and 0.3) and values of K /c comparable to that already published in the literature. We also showed that strain rate has a strong impact on peak stress and that dilatancy appears right before failure and hence highlighting the importance of plasticity mechanisms. Although brittle failure was systematically observed, stress drops and associated slips were slow and deformation always remained aseismic (no acoustic emission were detected). Lastly, this confirms that shales are good lithological candidates for shallow crust aseismic creep and slow slip events.

58 GEOSCIENCES↗