Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reacting”

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 955 records · Page 53

Design and Analyses of the NSTX-U PF1A Poloidal Field Coil Support

New poloidal field (PF) coils and supports have been designed for the National Spherical Torus Experiment Upgrade (NSTX-U) Recovery project. The PF1A coil structural requirements demand a large preload to pre-compress the coil insulation. Coil support design accounts for all operational conditions and loads of NSTX-U including bakeout, seismic events, and electromagnetic loads. Background magnetic fields produce coil Lorentz forces that react on the structural supports. These Lorentz forces include interactions from the coil current and all plasma events: normal operations, disruptions, and vertical displacement events. This paper will detail the design and design loads and highlight some of the analyses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Century: Zap Energy’s 100-kW-Scale Repetitive Sheared-Flow-Stabilized Z -Pinch System with Liquid Metal Cooling

Zap Energy is developing the sheared-flow-stabilized (SFS) Z-pinch concept for commercial applications. The SFS Z pinch relies on plasma self-organization, in the sense that plasma dynamics play a critical role in confinement. Using plasma axial current for confinement and compression eliminates the need for external confinement or heating technologies. This compact magnetic confinement technology could, in turn, provide the basis for a cost-effective deuterium-tritium fusion power plant. In addition to a robust experimental program pushing plasma performance towards breakeven conditions, Zap Energy has parallel programs developing power handling systems suitable for future power plants. Technologies under development include high average-power repetitive pulsed power, high duty-cycle cathodes, and liquid metal wall systems. Century is the name of Zap Energy’s first effort to integrate these three components into an operational system capable of firing non-reacting hydrogen SFS Z-pinch plasmas into a liquid-metal-lined container at sustained repetition rates on the order of 0.1 Hz. Here, the pulsed power driver and liquid metal heat exchanger are both designed to sustain input powers of 100 kW. Construction and initial operations with an interim ~10 kW liquid metal heat exchanger are described.

Century↗

Quantitative x-ray scattering of free molecules

Advances in x-ray free electron lasers have made ultrafast scattering a powerful method for investigating molecular reaction kinetics and dynamics. Accurate measurement of the ground-state, static scattering signals of the reacting molecules is pivotal for these pump-probe x-ray scattering experiments as they are the cornerstone for interpreting the observed structural dynamics. Here, this article presents a data calibration procedure, designed for gas-phase x-ray scattering experiments conducted at the Linac Coherent Light Source x-ray Free-Electron Laser at SLAC National Accelerator Laboratory, that makes it possible to derive a quantitative dependence of the scattering signal on the scattering vector. A self-calibration algorithm that optimizes the detector position without reference to a computed pattern is introduced. Angle-of-scattering corrections that account for several small experimental non-idealities are reported. Their implementation leads to near quantitative agreement with theoretical scattering patterns calculated with ab-initio methods as illustrated for two x-ray photon energies and several molecular test systems.

74 ATOMIC AND MOLECULAR PHYSICS↗

Efficient passivation of n-type and p-type silicon surface defects by hydrogen sulfide gas reaction

An efficient surface defect passivation is observed by reacting clean Si in a dilute hydrogen sulfide-argon gas mixture (< 5% H2S in Ar) for both n-type and p-type Si wafers with planar and textured surfaces. Surface recombination velocities of 1.5 and 8 cm/s are achieved on n-type and p-type Si wafers, respectively, at an optimum reaction temperature of 550oC that are comparable to the best surface passivation quality used in high efficiency Si solar cells. Surface chemical analysis using x-ray photoelectron spectroscopy shows that sulfur is primarily bonded in a sulfide environment, and synchrotron-based soft x-ray emission spectroscopy of the adsorbed sulfur atoms suggests the formation of S-Si bonds. Furthermore, the sulfur surface passivation layer is unstable in air, attributed to surface oxide formation and a simultaneous decrease of sulfide bonds. However, the passivation can be stabilized by a low-temperature (300oC) deposited amorphous silicon nitride (a-Si:NX:H) capping layer.

14 SOLAR ENERGY↗

Novel inverse oxide/metal catalysts for methanol synthesis: impact of oxide–metal interactions and reversible morphological changes

Inverse oxide/metal catalysts have proved to be excellent systems for the generation of methanol by CO 2 hydrogenation or the partial oxidation of methane. These systems can exhibit unique structural and chemical properties due to the nano size of the oxide component and strong oxide–metal interactions. Recent studies for ZnO/Cu, CeO 2 /Cu, TiO 2/ Cu, MgO/Cu, In 2 O 3 /Cu, and In 2 O 3 /Au catalysts have shown large variations in the composition and morphology of the oxide overlayer as a function of temperature and chemical environment. These oxide–metal interfaces are able to react with CO 2 and CH 4 at room temperature, and both reactants have a strong influence on the physical and chemical properties of the catalysts. Under reaction conditions, switches between oxide–metal and metal–metal interfaces can take in the catalyst surface. A dynamic behavior that can be linked to a high selectivity for methanol production over systems such as ZnO/Cu, CeO 2 /Cu and In 2 O 3 /Au. Finally, this type of correlation deserves additional systematic studies since it could be a powerful tool for designing highly efficient catalysts for methanol synthesis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Robust two-colour pyrometry uncertainty analysis to acquire spatially-resolved measurements

Two-colour pyrometry (2CP) has been used over several decades to study engine-relevant combustion processes, but results are generally regarded as qualitative or semi-quantitative. In many current 2CP systems, large measurement errors are introduced by parallax because the two measured wavelengths are not from the same line of sight. Here this work presents a spatially-resolved 2CP system with the objective of quantifying and reducing measurement uncertainty. An optical setup that eliminates parallax in 2CP is used together with pixel-by-pixel calibration of the camera sensor to increase measurement accuracy. Primary uncertainty terms are identified, and an error propagation analysis is performed to compute uncertainties in the final results of soot temperature, soot concentration parameter, KL, and soot mass. These methodologies are applied to investigate an auto-igniting fuel spray in a constant pressure flow rig at diesel-like conditions of high ambient pressure and temperature. Results show bias uncertainty of around 200 K (≈10%) for temperature and about 40%–60% for KL. High uncertainty was found to occur on the diffusion flame front where both optical thickness and soot concentrations are small. However, these uncertain measurement zones with relatively low soot concentrations contribute minimally to the total soot mass present in the reacting jet during the temporal evolution of the flame.

42 ENGINEERING↗

Ammonia generation in Ns pulse and Ns pulse/RF discharges over a catalytic surface

Plasma-catalytic ammonia synthesis in a ns pulse discharge and a 'hybrid' ns pulse/RF discharge in plane-to-plane geometry is studied by Fourier Transform infrared absorption spectroscopy. The data are taken in a preheated H 2 –N 2 mixture, with and without Ni/γ-Al 2 O 3 or Co/γ-Al 2 O 3 catalyst placed in the discharge section. The measurement results are taken using two different approaches. The first is a 'single-stage' process, where a ns pulse discharge in the H 2 –N 2 mixture is sustained continuously. In this case, the ammonia yield increases slowly, over a period of tens of minutes. The second is a 'two-stage' process, where the catalyst is first activated by the ns pulse discharge sustained in pure nitrogen, and then the activated catalyst is exposed to the H 2 –N 2 flow, with or without the discharge. In this case, a strong overshoot of the NH 3 number density at the reactor exit is detected, by over a factor of two compared to the single-stage process. After the initial overshoot, the ammonia yield gradually decreases to the 'single stage' value (with the discharge on), or to near zero (with the discharge off). The results demonstrate that the ammonia yield in the plasma-catalytic reactor is controlled by the N atom accumulation on the catalyst surface, which reacts with H atoms thermally dissociated on the catalyst or generated in the plasma. The results also show that the plasma-catalytic ammonia yield is significantly higher compared to that in the ns pulse discharge without the catalyst. The accumulation of H atoms on the catalyst, with their subsequent reactions with N atoms generated in the plasma, is of relatively minor importance at the present conditions. An additional series of measurements was made with a sub-breakdown RF waveform overlapped with the ns pulse discharge train, to enhance the vibrational excitation of nitrogen. The ammonia yield measured with the RF waveform added is approximately 20% higher compared to that at the baseline ns pulse discharge conditions, both with and without the catalyst. This effect is weaker compared to that of the catalyst activation by N atoms. Additional data are necessary to isolate the possible effect of the vibrationally excited N 2 molecules on the ammonia synthesis in the plasma catalytic reactions.

surface reactions↗

Plasma catalysis: separating plasma and surface contributions for an Ar/N 2 /O 2 atmospheric discharge interacting with a Pt catalyst

Atmospheric pressure non-equilibrium plasmas can form nitrogen oxide (NO x ) compounds directly from nitrogen and oxygen without a catalyst, and at lower catalyst temperatures than would be possible without plasma. In this work, the oxidation of plasma-produced NO from an Ar/N 2 /O 2 non-equilibrium atmospheric-pressure plasma-jet (APPJ) over a platinum-on-alumina powder catalyst was investigated with in-situ infrared spectroscopy. Products downstream of the catalyst bed were analyzed along with catalyst surface species. The catalyst was exposed to plasma at both constant temperature and a cyclic temperature ramp in order to study long-lasting and transient surface changes. Primary incident reactive species to the catalyst were assessed to be NO and O 3 . Pt-Al 2 O 3 at 350 °C increased oxidation of NO relative to Al 2 O 3 or an empty chamber. The surface state of Pt-Al 2 O 3 evolves during plasma-effluent exposure and requires upwards of 20 min exposure for stabilization compared to Al 2 O 3 . Once stable surface conditions are achieved, thermal cycling reveals a repeatable hysteresis pattern in downstream products. At low temperature, oxygen and NO x accumulate on the catalyst surface and react at elevated temperatures to form NO 2 . Increasing plasma power and O 2 :N 2 ratio increases the hysteresis of the heating relative to the cooling curves in the pattern of NO 2 formation. The limitation on NO oxidation at high temperatures was assessed to be Pt-O which is depleted as the catalyst is heated. Once stored species have been depleted, NO oxidation rates are determined by incoming reactants. Two overlapping NO oxidation patterns are identified, one determined by surface reactants formed at low temperature, and the other by reactants arriving at the surface at high temperature. As a result, the plasma is responsible for providing the reactants to the catalyst surface, while the catalyst enables reaction at high temperature or storage at low temperature for subsequent reaction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Availability and reactivity of N 2 ( v ) for NH 3 synthesis by plasma catalysis

Production of vibrationally excited N 2 (N 2 (v)) in atmospheric pressure nonthermal plasma and loss of N 2 (v) by gas-phase reactions and reactions on catalytic surfaces are analyzed to examine the role of N 2 (v) in NH 3 formation by plasma catalysis. Vibrational state-to-state kinetic models complemented with molecular beam mass spectrometry (MBMS) measurements demonstrate that N 2 (v> 0) is produced with densities 100× greater than the density of N radicals by a radiofrequency atmospheric pressure plasma jet. The experimentally measured loss of N 2 (v) corresponds with a state-to-state kinetic model that describes loss of N 2 (v) by surface-mediated vibrational relaxation without consideration of reactions that convert N 2 (v) to NH 3 over the catalyst surface. Rate constants for vibrational relaxation of N 2 (v) on catalyst surfaces exceed upper bounds on proposed rate constants for NH 3 formation reactions from N 2 (v) over Fe when v < 9, Ni when v < 18, and Ag when v < 39, which indicates that only higher vibrational levels can possibly contribute to catalytic NH 3 formation faster than they undergo vibrational relaxation on the surface. Densities of N 2 (v> 8), vibrational levels that can possibly react over Fe to form NH 3 faster than they undergo vibrational relaxation, are less than or similar to N densities at the inlet of the catalyst bed and measured NH 3 formation for the investigated conditions in this work, while densities of N 2 (v> 17) and N 2 (v> 38) are orders of magnitude below the N density at the inlet of the catalyst bed and the measured NH 3 formation. The loss of N 2 (v) by vibrational relaxation on the surface limits the ability of N 2 (v) to contribute to catalytic NH 3 formation and explains why N 2 (v) does not produce NH 3 in quantities that are comparable to NH 3 formation from N even though N 2 (v > 0) is more abundantly produced by the plasma.

state-to-state kinetic modeling↗

Rapid Nb 3 Sn film growth by sputtering Nb on hot bronze

Nb deposited by magnetron sputtering onto hot Cu-15 wt.%Sn bronze substrates at temperatures above 700 °C achieved Nb 3 Sn film growth at a rate of 33 nm min -1 , which was an order of magnitude faster than that achieved for deposition of Nb on bronze at low temperature followed by in situ post reaction at the same high temperatures. Tin content in the Nb 3 Sn films made on hot bronze was ~26.3%, which is significantly higher than ~24.5% obtained by post-reaction as well as for typical bulk reactions between Nb and α -bronze. The Nb 3 Sn lattice parameter was consistent with measured tin content and predicted elastic strain for both routes. Critical temperatures of 14 K–16 K, instead of 18 K, were consistent with elastic strain due to coefficient of thermal expansion mismatch between the Nb 3 Sn and bronze substrate and, for the hot-bronze samples, stress related to the growth mechanism. Films were fully coalesced and had surface roughness values <20 nm over a 100 μ m 2 scan. Grain structure of the Nb 3 Sn films produced by Nb sputtering on hot bronze resembles zone 2 in the Thornton structure-zone diagram, in contrast to the equiaxed grain structure reminiscent of microstructure observed in reacted Nb 3 Sn wires exhibited by the post-reaction route.

Physics↗

Significant reduction in the low-field magnetization of Nb 3 Sn superconducting strands using the internal oxidation APC approach

Nb 3 Sn superconductors are promising for building accelerator magnets for future energy-frontier circular colliders. A critical factor for this application is the low-field persistent-current magnetization because it leads to several critical issues: e.g. low-field instability (including flux jumps), hysteresis loss, and field errors in magnet bores. Suppression of low-field magnetization requires reduction of low-field critical current density (J c ) or effective subelement size (d eff ). However, reduction of deff of state-of-the-art Nb 3 Sn conductors—the restacked-rod-process (RRP ® ) type—below 40–50 μm without a pronounced decrease in high-field J c is difficult. On the other hand, the internal oxidation method which forms artificial pinning centers (APC) in Nb 3 Sn offers an alternative approach to reducing the low-field magnetization. Compared with a conventional Nb 3 Sn conductor whose flux pinning force versus field (F p –B) curve peaks at ~20% of its irreversibility field (B irr ), the F p –B curve peaks of APC conductors shift to higher fields due to the point pinning effect, leading to flattening of the J c –B curves. The goal of this paper is to quantitatively study how much the APC approach can reduce the low-field magnetization. Here we measured the J c –B curves of an RRP ® conductor and two APC conductors (reacted at 700 °C) from zero field to B irr using a high-field vibrating sample magnetometer. Here the results showed that the APC conductors have higher non-Cu J c at high fields (e.g. 32%–41% higher at 16 T) and simultaneously lower non-Cu Jc at low fields (e.g. 28%–34% lower at 1 T) compared with the RRP ® . This effect is due to a competition between their Nb3Sn layer fraction ratios and layer F p ratios. Suppose they reach the same 16 T non-Cu J c , then the 1 T non-Cu J c and magnetization of the APC conductors are only half or even less compared with the RRP ® conductor.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Constraints on ion velocity distributions from fusion product spectroscopy

Abstract Recent inertial confinement fusion experiments have shown primary fusion spectral moments which are incompatible with a Maxwellian velocity distribution description. These results show that an ion kinetic description of the reacting ions is necessary. We develop a theoretical classification of non-Maxwellian ion velocity distributions using the spectral moments. At the mesoscopic level, a monoenergetic decomposition of the velocity distribution reveals there are constraints on the space of spectral moments accessible by isotropic distributions. General expressions for the directionally dependent spectral moments of anisotropic distributions are derived. At the macroscopic level, a distribution of fluid element velocities modifies the spectral moments in a constrained manner. Experimental observations can be compared to these constraints to identify the character and isotropy of the underlying reactant ion velocity distribution and determine if the plasma is hydrodynamic or kinetic.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Feasibility of fusion plasma burn control via real-time, sub-divertor neutral gas isotopic and compositional analysis

The ability to provide fusion burn control without requiring physical access through the first wall and fuel breeding blankets, would be vital for any future, magnetically confined fusion power reactor. A multi-sensor, fusion fuel cycle exhaust, neutral gas analysis system on JET, capable of delivering real time data, and accessing only the sub-divertor region, provides an excellent example of such capability. Optimized for and operated during the deuterium–tritium experimental campaigns 2 and 3 (DTE2, DTE3), it is proving valuable for planning to explore fusion reactor burn control in ITER with a comparable diagnostic system called the Diagnostic Residual Gas Analyzer (DRGA). This paper aims to show feasibility of developing model-based controllers for ITER and next generation, reactor-relevant devices, by building both on the empirical experience in JET-DTE2, and on the already emerging experience on developing such models specifically for ITER. The paper begins with a specific use-case from JET-DTE2, pertaining to the observed sensitivity of the fusion neutron yield on the concentration of isotopic helium-3 ( 3 He), with data from one of the high-performance DT shots exhibited with emphasis on the 3 He measurement via the sub-divertor. Then, a first model is developed and then explored with simulations that aim to discover how well the controllers in the model react to either insufficient levels of 3 He or excessive levels of 3 He. The simulations then explore potential impact from a delay in the measurement (or the response) that would be comparable to the ∼1 s, conductance limited response for the ITER DRGA system, currently in its final design. The simulations show that control is feasible, and that its effectiveness is not significantly impacted by such delay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental demonstration of real-time electron temperature profile control in DIII-D

Future tokamak reactor operation will require the ability to maintain a given plasma scenario for extended periods of time. This will necessitate the capability to react to changes in the plasma state and return the plasma to the target scenario; the principal method to achieve this is through feedback control. Thus, it is necessary to develop and test feedback controllers for the plasma profiles that define a target scenario. In this work, a feedback controller for the electron temperature (Te) profile is tested experimentally in DIII-D. This experiment relied on the ability to ascertain the electron temperature profile in real time, which was achieved using an observer algorithm. The observer relies on both diagnostic data and a predictive model of the electron temperature profile evolution; this predictive model includes contributions from neural network surrogate models. Because of these dependencies, a number of capabilities needed to be added to the real-time PCS for DIII-D in order to support the Te profile control experiment. The neural network surrogates needed to be integrated into the PCS to be called in real time. An observer algorithm for the Te profile needed to be added and connected to the Thomson scattering system to allow access to the current state of the profile in real time. When tested, the observer was shown to produce Te profiles that are consistent with the shape of the Thomson scattering data while rejecting much of the noise in the diagnostic data. Finally, the controller itself was tested in real time. This experiment showed that the controller is capable of tracking the electron temperature target at locations across the spatial profile.

Morosohk, Shira [Oak Ridge Associated Universities↗

On the attribution of the impacts of extreme weather events to anthropogenic climate change

Investigations into the role of anthropogenic climate change in extreme weather events are now starting to extend into analysis of anthropogenic impacts on non-climate (e.g. socio-economic) systems. However, care needs to be taken when making this extension, because methodological choices regarding extreme weather attribution can become crucial when considering the events’ impacts. The fraction of attributable risk (FAR) method, useful in extreme weather attribution research, has a very specific interpretation concerning a class of events, and there is potential to misinterpret results from weather event analyses as being applicable to specific events and their impact outcomes. Using two case studies of meteorological extremes and their impacts, we argue that FAR is not generally appropriate when estimating the magnitude of the anthropogenic signal behind a specific impact. Attribution assessments on impacts should always be carried out in addition to assessment of the associated meteorological event, since it cannot be assumed that the anthropogenic signal behind the weather is equivalent to the signal behind the impact because of lags and nonlinearities in the processes through which the impact system reacts to weather. Whilst there are situations where employing FAR to understand the climate change signal behind a class of impacts is useful (e.g. ‘system breaking’ events), more useful results will generally be produced if attribution questions on specific impacts are reframed to focus on changes in the impact return value and magnitude across large samples of factual and counterfactual climate model and impact simulations. We advocate for constant interdisciplinary collaboration as essential for effective and robust impact attribution assessments.

54 ENVIRONMENTAL SCIENCES↗

Statistical properties of filaments in the cosmic web

ABSTRACT In the context of the cosmological and constrained Exploring the Local Universe with the reConstructed Initial Density field (ELUCID) simulation, this study explores the statistical characteristics of filaments within the cosmic web, focussing on aspects such as the distribution of filament lengths and their radial density profiles. Using the classification of the cosmic web environment through the Hessian matrix of the density field, our primary focus is on how cosmic structures react to the two variables $R_{\rm s}$ and $\lambda _{\rm th}$. The findings show that the volume fractions of knots, filaments, sheets, and voids are highly influenced by the threshold parameter $\lambda _{\rm th}$, with only a slight influence from the smoothing length $R_{\rm s}$. The central axis of the cylindrical filament is pinpointed using the medial-axis thinning algorithm of the COsmic Web Skeleton (COWS) method. It is observed that median filament lengths tend to increase as the smoothing lengths increase. Analysis of filament length functions at different values of $R_{\rm s}$ indicates a reduction in shorter filaments and an increase in longer filaments as $R_{\rm s}$ increases, peaking around $2.5R_{\rm s}$. The study also shows that the radial density profiles of filaments are markedly affected by the parameters $R_{\rm s}$ and $\lambda _{\rm th}$, showing a valley at approximately $2R_{\rm s}$, with increases in the threshold leading to higher amplitudes of the density profile. Moreover, shorter filaments tend to have denser profiles than their longer counterparts.

Zhang, Youcai (ORCID:0000000319674091)↗

Structural basis for the inability of chloramphenicol to inhibit peptide bond formation in the presence of A-site glycine

Ribosome serves as a universal molecular machine capable of synthesis of all the proteins in a cell. Small-molecule inhibitors, such as ribosome-targeting antibiotics, can compromise the catalytic versatility of the ribosome in a context-dependent fashion, preventing transpeptidation only between particular combinations of substrates. Classic peptidyl transferase center inhibitor chloramphenicol (CHL) fails to inhibit transpeptidation reaction when the incoming A site acceptor substrate is glycine, and the molecular basis for this phenomenon is unknown. Here, we present a set of high-resolution X-ray crystal structures that explain why CHL is unable to inhibit peptide bond formation between the incoming glycyl-tRNA and a nascent peptide that otherwise is conducive to the drug action. Our structures reveal that fully accommodated glycine residue can co-exist in the A site with the ribosome-bound CHL. Moreover, binding of CHL to a ribosome complex carrying glycyl-tRNA does not affect the positions of the reacting substrates, leaving the peptide bond formation reaction unperturbed. These data exemplify how small-molecule inhibitors can reshape the A-site amino acid binding pocket rendering it permissive only for specific amino acid residues and rejective for the other substrates extending our detailed understanding of the modes of action of ribosomal antibiotics.

59 BASIC BIOLOGICAL SCIENCES↗

Mineralogy of the deep lower mantle in the presence of H 2 O

Understanding the mineralogy of the Earth's interior is a prerequisite for unravelling the evolution and dynamics of our planet. Here, we conducted high pressure-temperature experiments mimicking the conditions of the deep lower mantle (DLM, 1800–2890 km in depth) and observed surprising mineralogical transformations in the presence of water. Ferropericlase, (Mg, Fe)O, which is the most abundant oxide mineral in Earth, reacts with H 2 O to form a previously unknown (Mg, Fe)O 2 H x (x ≤ 1) phase. The (Mg, Fe)O 2 H x has a pyrite structure and it coexists with the dominant silicate phases, bridgmanite and post-perovskite. Depending on Mg content and geotherm temperatures, the transformation may occur at 1800 km for (Mg 0.6 Fe 0.4 )O or beyond 2300 km for (Mg 0.7 Fe 0.3 )O. The (Mg, Fe)O 2 H x is an oxygen excess phase that stores an excessive amount of oxygen beyond the charge balance of maximum cation valences (Mg 2+ , Fe 3+ and H + ). This important phase has a number of far-reaching implications including extreme redox inhomogeneity, deep-oxygen reservoirs in the DLM and an internal source for modulating oxygen in the atmosphere.

58 GEOSCIENCES↗