Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “observation geometry”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Scale-Dependent Bedform Migration and Deformation in the Physical and Spectral Domains

Multi-scale bathymetries observed in laboratory channels and natural rivers have posed several challenges to the description of bedform geometry and kinematics. To better quantify scale-dependent migration velocity, high resolution spatio-temporal bed evolution data are analyzed using three independent methods: (i) a bedform tracking method that identifies individual bedforms in longitudinal bed elevation profiles, (ii) a cross-correlation based, image matching technique in consecutive bathymetry scans similar to Particle Image Velocimetry, and (iii) two-dimensional frequency-wavenumber spectra marking a dispersion relation between bedform period and length, allowing estimation of Fourier-based, scale-dependent migration velocities. The comparative studies show that the spectral and bedform tracking methods present small differences in the migration velocity of medium-large bedforms. However, significant deviations occur in the small bedforms, mostly due to bedform sheltering or amplified exposure to the flow, which also contribute to their enhanced scale-dependent deformation. The analysis also shows that some form of scale separation is needed to extract the large and slowly migrating bedforms as they are masked by secondary features when several bedform orders coexist. Bathymetric Image Velocimetry reasonably estimates the averaged migration velocity and exhibits potential to study relations between spanwise bed surface velocity and local bedform kinematics.

58 GEOSCIENCES↗
DOE OSTI2026

Morphological descriptors of nanoparticles: The link between atomistic structures and x-ray absorption spectra

Understanding and quantifying the morphology of nanoparticles are essential for linking their atomic structure to diverse applications and verifying theoretical models. While experimental information on the structure of nanoparticles in the size range below ∼5 nm can be extracted from x-ray absorption spectroscopy using a small number of descriptors—most commonly coordination numbers—developing an understanding of morphology descriptors from experimental data remains a challenge. Here, in this study, we introduce NanoGene, a genetic algorithm-based method for generating structurally diverse nanoparticle models guided by user-defined descriptors. We establish correlations among structural, size-related, and morphological descriptors and demonstrate how experimentally accessible parameters, such as coordination numbers, can be leveraged to infer otherwise inaccessible ones, such as the generalized coordination number or particle oblateness. Principal component and clustering analyses reveal the relative importance of descriptors, with the number of atoms emerging as a key discriminant of the nanoparticle structure. By providing both the methodology and an extensive dataset of nanoparticle geometries, this work offers a practical foundation for descriptor-based analysis and interpretation of experimental observations, bridging the gap between local atomic coordinates and global morphological characterization.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Enhanced orbital anisotropy through the proximity to a SrTiO 3 layer in the perovskite iridate superlattices

We have used angle-dependent soft x-ray absorption spectroscopy (XAS) at the O K edge and first-principles calculations to investigate the electronic structures of iridate-based superlattices ( SrIr O 3 ) m / ( SrTi O 3 ) ( m = 1 , 2, 3, and ∞). We focus on the pre-edge Ir 5 d t 2 g – O 2 p orbital hybridization feature in the XAS spectra. By varying the measurement geometry relative to the incident photon polarization, we are able to extract the dichroic contrast and observe the systematic increase in the anisotropy of Ir 5 d orbitals as m decreases. First-principles calculations elucidate the orbital anisotropy coming mainly from the enhanced out-of-plane compression of Ir O 6 octahedra in the SrIr O 3 layers that are adjacent to the inserted SrTi O 3 layers. As m decreases, the increased volume fraction of these interfacial SrIr O 3 layers and their contact with the SrTi O 3 layers within the ( SrIr O 3 ) m / ( SrTi O 3 ) supercell lead to enhanced orbital anisotropy. Furthermore, the tilt and rotation of Ir O 6 octahedra are shown to be essential to understand the subtle orbital anisotropy in these superlattices, and constraining these degrees of freedom will give an incorrect trend. Overall, our results demonstrate that the structural constraint from the inserted SrTi O 3 layers, in addition to other electronic means such as polar interface and charge transfer, can serve as a knob to control the orbital degree of freedom in these iridate-based superlattices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Piezoelectric Actuation of Graphene-Coated Polar Structures

Ferroelectric materials based on lead zirconate titanate (PZT) are widely used as sensors and actuators because of their strong piezoelectric activity. However, their application is limited because of the high processing temperature, brittleness, lack of conformal deposition, and a limited possibility to be integrated with the microelectromechanical systems (MEMS). Recent studies on the piezoelectricity in the 2-D materials have demonstrated their potential in these applications, essentially due to their flexibility and integrability with the MEMS. In this work, we deposited a few layer graphene (FLG) on the amorphous oxidized Si 3 N 4 membranes and studied their piezoelectric response by sensitive laser interferometry and rigorous finite-element modeling (FEM) analysis. Modal analysis by FEM and comparison with the experimental results show that the driving force for the piezoelectric-like response can be a polar interface layer formed between the residual oxygen in Si 3 N 4 and the FLG. The response was about 14 nm/V at resonance and could be further enhanced by adjusting the geometry of the device. These phenomena are fully consistent with the earlier piezoresponse force microscopy (PFM) observations of the piezoelectricity of the graphene on SiO 2 and open up an avenue for using graphene-coated structures in the MEMS.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Model Form Error Correction for a Black-Box Thermal Battery Heat Transfer Simulation

Thermal batteries are crucial for supplying power to high-consequence engineering applications such as rockets. Computational simulations have been developed to predict thermal battery behavior, but these simulations often suffer from modeling errors, including model form uncertainty. Addressing this uncertainty can be achieved by quantifying either the model discrepancy in the output or the model form error (MFE) in the governing equation. MFE is particularly valuable as it can be better extrapolated beyond observed outputs, which is essential for predictions involving changes in external system loading, system configuration and geometry, or output quantities. This paper employs a state estimation approach to estimate MFE using experimental data and then utilizes machine learning (ML) to model its relationship with state variables. A nonintrusive technique is used to estimate MFE in a black-box thermal battery heat transfer simulation. The trained machine learning model for MFE is then applied to correct simulation predictions under extrapolated initial conditions and battery configurations. In conclusion, the methodology's performance is evaluated using additional experimental data, demonstrating its effectiveness in improving prediction accuracy.

Batteries↗
DOE OSTI2020

Widely tunable, heterogeneously integrated quantum-dot O-band lasers on silicon

Heterogeneously integrated lasers in the O-band are a key component in realizing low-power optical interconnects for data centers and high-performance computing. Quantum-dot-based materials have been particularly appealing for light generation due to their ultralow lasing thresholds, small linewidth enhancement factor, and low sensitivity to reflections. Here, we present widely tunable quantum-dot lasers heterogeneously integrated on silicon-on-insulator substrate. The tuning mechanism is based on Vernier dual-ring geometry, and a 47 nm tuning range with 52 dB side-mode suppression ratio is observed. These parameters show an increase to 52 nm and 58 dB, respectively, when an additional wavelength filter in the form of a Mach–Zehnder interferometer is added to the cavity. The Lorentzian linewidth of the lasers is measured as low as 5.3 kHz.

Malik, Aditya (ORCID:0000000334487577)↗
DOE OSTI2024

Probabilistic Analysis of Long-Term Degradation of Microwave Cavity Flow Sensor

We are investigating a microwave resonant cavity transducer for flow sensing in the vessel of a high temperature fluid advanced reactor (AR), such as a molten salt cooled reactor (MSCR) or a sodium fast reactor (SFR). This transducer is a hollow metallic cylindrical cavity, with the flat wall of the cylinder flexible enough to undergo microscopic deflection due to dynamic fluid pressure. Membrane deflection leads to a shift in the resonant frequency, which can be detected with a spectrum analyzer. We have performed a proof-of-concept experiment of flow sensing with the transducer in liquid sodium at 340°C in impinging liquid jet geometry. The transducer remained in liquid sodium for 70 days. After removal, no structural damage was observed, and the expected transducer response was verified in a water test. Because long-term (multi-year) experimental tests of transducer resilience to harsh environment are not practical, we have developed a probabilistic model of creep to estimate transducer resilience to the harsh environment. The probabilistic model considers diffusion creep under the condition of high temperature and low stress, where the stress and temperature are allowed to be random variables with Gaussian distributions. Using the probabilistic model, we estimate inelastic membrane deflections due to creep for several temperature ranges. We conclude that for temperatures less than 650°C, creep has negligible long-term effect on the transducer performance. Since a yellowish residue was observed on the transducer surface after 70 days of immersion in liquid sodium, we have investigated possible evidence of corrosion. Chromium depletion is a typical indicator of the corrosion process in stainless steel. Scraping off a residue from the transducer and performing scanning electron microscopy (SEM) with energy dispersive analysis (EDS) did not find any chromium in the residue. Approximately 60% of the residue consisted of copper, which can be attributed to contamination of sodium due to powder residue from machining of copper and brass components of the transducer.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Precursory Off-Fault Deformation in Restraining and Releasing Step Overs: Insights From Discrete Element Method Models

Accelerating geophysical activity is detected preceding some, but not all, large earthquakes. This observation may indicate that no precursors occur before some earthquakes, or that the instrumentation lacks the required sensitivity. Here, to aid crustal monitoring efforts, we use discrete element method models to identify the locations and styles of deformation that may provide useful information about approaching fault reactivation. We model the reactivation of two healed rough faults in a variety of step over configurations, embedded in a host rock with varying amounts of damage subject to shear velocity loading parallel to the faults. Both the fault geometry and ratio of fault to host rock strength control the amount of off-fault deformation. Consistent with field observations, models with larger steps and more preexisting host rock damage produce higher amounts of off-fault deformation. We assess the size of the continuous regions of high velocities and strains to compare the value of the precursory information of each velocity and strain component. Comparing the three components of the velocity vector suggests that the fault-parallel velocity produces the largest and most temporally continuous regions of elevated velocity. The size of these regions increases toward failure, indicating the usefulness of tracking this component. Comparing the volumetric and shear components of the three-dimensional strain tensor suggests that during most of the interseismic period, the shear strain provides more information about approaching fault slip than the volumetric strain. However, in the days and months preceding fault reactivation, both the shear and volumetric strains provide similarly valuable information.

58 GEOSCIENCES↗
DOE OSTI2021

Arctic observations and numerical simulations of surface wind effects on Multi-Angle Snowflake Camera measurements

Abstract. Ground-based measurements of frozen precipitation are heavily influenced by interactions of surface winds with gauge-shield geometry. The Multi-Angle Snowflake Camera (MASC), which photographs hydrometeors in free-fall from three different angles while simultaneously measuring their fall speed, has been used in the field at multiple midlatitude and polar locations both with and without wind shielding. Here, we present an analysis of Arctic field observations – with and without a Belfort double Alter shield – and compare the results to computational fluid dynamics (CFD) simulations of the airflow and corresponding particle trajectories around the unshielded MASC. MASC-measured fall speeds compare well with Ka-band Atmospheric Radiation Measurement (ARM) Zenith Radar (KAZR) mean Doppler velocities only when winds are light (≤5ms-1) and the MASC is shielded. MASC-measured fall speeds that do not match KAZR-measured velocities tend to fall below a threshold value that increases approximately linearly with wind speed but is generally <0.5ms-1. For those events with wind speeds ≤1.5ms-1, hydrometeors fall with an orientation angle mode of 12∘ from the horizontal plane, and large, low-density aggregates are as much as 5 times more likely to be observed. Simulations in the absence of a wind shield show a separation of flow at the upstream side of the instrument, with an upward velocity component just above the aperture, which decreases the mean particle fall speed by 55 % (74 %) for a wind speed of 5 m s−1 (10 m s−1). We conclude that accurate MASC observations of the microphysical, orientation, and fall speed characteristics of snow particles require shielding by a double wind fence and restriction of analysis to events where winds are light (≤5ms-1). Hydrometeors do not generally fall in still air, so adjustments to these properties' distributions within natural turbulence remain to be determined.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Non-premixed counterflow methane flames in DC/AC/NS electric fields

The response of counterflow diffusion flames to sub-breakdown DC and AC electric fields, as well as their superposition with ns pulse discharge waveforms, is studied in the plane-to-plane electrode geometry. Sub-breakdown DC and low-frequency AC electric fields cause the flame displacement toward the grounded electrode, in the direction of the applied field, indicating that the body force on the positive ions exceeds that on the electrons and negative ions. As the AC frequency increases, the flame response becomes less pronounced, due to the reduction of the electrohydrodynamic (EHD) body force impulse over the AC half-period. The electric field in the electrode gap is determined by ps Electric Field Induced Second Harmonic (E-FISH) generation, with absolute calibration using sub-breakdown ns pulses overlapped with the measured electric field waveform. The results show that the electric field distribution across the flame in the current saturation regime follows the Laplacian field. This indicates that the space charge density in the gap is too low to distort the applied DC or AC field, consistent with the kinetic modeling predictions. Combining a nanosecond pulse discharge with a sub-breakdown DC field generates a diffuse plasma across the entire gap. Time-resolved and spatially resolved measurements of the electric field in the discharge indicate the ionization wave propagation between the electrodes. The present results do not exhibit a detectable flame displacement enhancement by ns discharge pulses combined with a sub-breakdown field, observed previously. Kinetic modeling calculations show that the absence of this effect in the plane-to-plane geometry is due to the rapid plasma self-shielding. This indicates that alternative electrode geometries limiting the self-shielding would be more effective for the plasma / electric field enhanced flameholding and flame stabilization applications.

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

Binding of Li + to Negatively Charged and Neutral Ligands in Polymer Electrolytes

Conceptually, single-ion polymer electrolytes (SIPE) with the anion bound to the polymer could solve major issues in Li-ion batteries, but their conductivity is too low. Experimentally, weakly interacting anionic groups have the best conductivity. To provide a theoretical basis for this result, density functional theory calculations of the optimized geometries and energies are performed for charged ligands used in SIPE. Comparison is made to neutral ligands found in dual-ion conductors, which demonstrate higher conductivity. Further, the free energy differences between adding and subtracting a ligand are small enough for the neutral ligands to have the conductivity seen experimentally. However, charged ligands have large barriers, implying that lithium transport will coincide with the slow polymer diffusion, as observed in experiments. Overall, SIPE will require additional solvent to achieve a sufficiently high conductivity. Additionally, the binding of mono- and bidentate geometries varies, providing a simple and clear reason that polarizable force fields are required for detailed interactions.

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

Measurement of prompt ${\mathrm{D^0}}$ and ${\mathrm{\overline{D}}{}^0}$ meson azimuthal anisotropy and search for strong electric fields in PbPb collisions at $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV

The strong Coulomb field created in ultrarelativistic heavy ion collisions is expected to produce a rapidity-dependent difference ( Δv 2 ) in the second Fourier coefficient of the azimuthal distribution (elliptic flow, v 2 ) between D 0 ( $\bar{u}c$ ) and ${\mathrm{\overline{D}}{}^0}$ ( $u\bar{c}$ ) mesons. Motivated by the search for evidence of this field, the CMS detector at the LHC is used to perform the first measurement of Δv 2 . The rapidity-averaged value is found to be &lang;Δv2&rang;=0.001±0.001(stat)±0.003(syst) in PbPb collisions at $\sqrt{s_{NN}}$=5.02 TeV . In addition, the influence of the collision geometry is explored by measuring the D 0 and ${\mathrm{\overline{D}}{}^0}$ mesons v 2 and triangular flow coefficient ( v 3 ) as functions of rapidity, transverse momentum ( p T ), and event centrality (a measure of the overlap of the two Pb nuclei). A clear centrality dependence of prompt D 0 meson v 2 values is observed, while the v 3 is largely independent of centrality. These trends are consistent with expectations of flow driven by the initial-state geometry.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Interfacial and Confinement-Mediated Organization of Gas Hydrates, Water, Organic Fluids, and Nanoparticles for the Utilization of Subsurface Energy and Geological Resources

Harnessing the subsurface geologic environments in an efficient and environmentally sustainable manner is challenged by uncertainties associated with predicting the fate of fluids and sustaining porosity and permeability in subsurface geologic environments. Some of these uncertainties arise from confined and interfacially induced structures of fluids in subsurface geologic environments. The formation of gas hydrates, phase transitions of confined fluids, assembly and deposition of heavy hydrocarbons, and the agglomeration and fate of nanoparticles in confined environments are summarized in this review. Nanoscale confinement contributes to anisotropic structures and dynamics of fluids, which is the basis for anomalous phase transition thermodynamics, reactivity, transport, and geomechanical behavior. In this review, we discuss the structures of confined fluids and deviation in observed properties from bulk fluids. The factors influencing the structures of confined fluids can be generally divided into two groups: (a) pore characteristics including pore size, pore surface chemistry, and pore geometry and (b) confined fluid/solid characteristics such as molecular structure, concentrations, charges, pore filling, and presence of additives. Scientific advancements and knowledge gaps in our understanding of the structures of confined fluids and the associated differences in observed properties compared to bulk fluids are discussed. Here, the phenomena discussed in this review are of particular relevance to our efforts in harnessing the subsurface environments for a low carbon future by increasing the utilization of geothermal energy, using CO 2 as a working fluid, and storing CO 2 in subsurface geologic environments.

42 ENGINEERING↗
DOE OSTI2023

Multifidelity validation of digital surrogates using variable-density turbulent mixing models

In this study, ensembles of experimental data are presented and utilized to compare and validate two models used in the simulation of variable-density (Atwood = 0.22), compressible turbulent mixing. Though models of this kind (Reynolds averaged NavierStokes and large-eddy simulations) have been validated extensively with more canonical flows in previous studies, here the present approach offers novelty in the complexity of the geometry, the ensemble-based validation, and the uniformity of the computational framework on which the models are tested. Moreover, all experimental and computational tasks were completed by the authors which has led to a tightly coupled experimental configuration with its “digital twin.” The experimental divergent-shock-tube facility and its data acquisition methods are described and replicated in simulation space. A 2D Euler model which neglects the turbulent mixing at the interface is optimized to experimental data using a Gaussian process. This model then serves as the basis for both the 2D RANS and 3D LES studies that make comparisons to the mixing-layer data from the experiment. A relatively simple RANS model is shown to produce good agreement with experimental data only at late flow development times. The LES ensembles generally show good agreement with experimental data but display sensitivity to the characterization of initial conditions. Resolution-dependent behavior is also observed for certain higher-order statistics of interest. Overall, the LES model successfully captures the effects of divergent geometry, compressibility, and combined nonlinear instabilities inherent to the problem. The successful prediction of mixing width and its growth rate highlight the existence of three distinct regimes in the development of the instability, each with similarities to previously studied instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Anisotropic thermal-conductivity degradation in the phase-field method accounting for crack directionality

Dynamic loads applied to metals may lead to brittle or ductile fracture depending on the imposed strain rates, material properties, and specimen geometry. With an increase in the applied velocity, brittle to ductile failure mode transition may also be observed. Furthermore, at high strain rates, ductile fracture may be preceded by shear bands, which are narrow bands of intense plastic deformation, typically accompanied by a significant rise in temperature. Reliable models are needed to predict the response of metals subject to dynamic loads. In particular, capturing the interplay between heat conduction and crack propagation using the phase-field fracture method is still an open research field. To accurately capture the heat transfer physics across crack surfaces, damage models degrading thermal-conductivity are necessary. While isotropic thermal-conductivity degradation models were proposed, they may lead to errors as they are indifferent to crack directionality. In our current work, a new anisotropic approach is proposed in which thermal-conductivity, which depends on the phase-field gradient, is degraded solely across the crack. Here, it is shown that this approach improves the near-field approximation of temperature and heat flux compared with isotropic degradation, when taking the discontinuous crack solutions as reference. Additionally, the anisotropic degradation approach is implemented in a unified dynamic fracture model and a couple of examples demonstrate the viability of this approach.

42 ENGINEERING↗
DOE OSTI2024

2D kinetic-ion simulations of inverted corona fusion targets

Laser-driven “inverted corona” fusion targets have attracted interest as a low-convergence neutron source and platform for studying kinetic physics. The scheme consists of a hollow or gas-filled spherical shell made of deuterated plastic. The shell has one or more laser entrance holes (LEH), resembling a spherical hohlraum. The laser passes through the LEH’s and illuminates the interior surface of the shell, ablating a plasma that travels inward towards the target center. Long ion mean free paths in the converging plasma can lead to significant interpenetration, atomic mix, and other kinetic effects. Here, in this work we report on numerical simulations of inverted corona targets using the kinetic-ion, fluid–electron hybrid particle-in-cell (PIC) approach in 2D RZ geometry. 2D simulations suggest that shape effects do not have a significant impact on plasma evolution and observed yield trends are primarily the result of 1D kinetic mix mechanisms. Simulations are also compared against available experimental data recorded at the OMEGA laser facility. In particular, synthetic x-ray emission images show good qualitative agreement with experimental results, albeit with an apparent timing discrepancy for the two-sided vacuum target. More generally, we demonstrate the potential of hybrid-PIC simulations for full-system modeling and experimental design, including collisional absorption of laser energy, plasma evolution, mix, and fusion burn.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The UV Photoinduced Ring-Closing Reaction of Cyclopentadiene Probed with Ultrafast Electron Diffraction

Conjugated cyclic organic molecules are common across many fields such as pharmaceuticals, are naturally occurring in biological systems, and are used in synthetic materials. One particular area of interest from a photochemical point of view is the formation of highly strained cyclic organics. We investigate the photoinduced reaction of cyclopentadiene, a five-membered organic ring molecule which can form strained three and four carbon rings after photoexcitation with UV light, with the gas-phase ultrafast electron diffraction instrument at the SLAC MeV-UED facility. Electron diffraction offers a direct probe sensitive to the nuclear geometry during the reaction, allowing for the determination of the distribution of products formed following photoexcitation. We observe the simultaneous formation of the highly strained ring- closed bicyclo[2.1.0]pentene and vibrationally hot cyclopentadiene within the temporal resolution of the experiment, and determine the relative yield of all reaction products. Furthermore, the experimental results are in good agreement with the predictions of trajectory simulations.

Computational Chemistry↗
DOE OSTI2021

A Pore‐Scale Investigation of Mineral Precipitation Driven Diffusivity Change at the Column‐Scale

Abstract Mineral precipitation affects the pore structure and thus transport properties of porous media. In this study, we investigated the pore‐scale dynamics of precipitation in diffusion controlled systems and the resulting impacts on the effective diffusivity, using a micro‐continuum reactive transport model. Forty two‐dimensional pore structures representing both idealized and realistic geometries were simulated with consideration of different precipitation scenarios and rates. A homogeneous nucleation scenario reproduced patterns observed in previous experimental study showing mixing‐induced precipitation, and a surface growth scenario captured the pattern for mineral precipitation on a substrate with the same or similar mineral structures. In all cases, local precipitation resulted in the reduction in the average porosity of the domain ( Φ ) until the cessation of diffusive transport and the termination of precipitation. The minimum porosity reached was referred to as the critical porosity ( Φ c ). The effective diffusivity ( D eff ) decreased with Φ and dropped sharply to effectively zero, that is, the critical effective diffusivity ( D eff c ), as Φ c was reached. These pore‐scale dynamics can be captured by a revised relationship that explicitly considers the critical porosity and the corresponding effective diffusivity, and the pre‐exponential coefficient and the exponent of the relationship varied with initial pore structure and the precipitation kinetics. Overall, the homogeneous nucleation scenario results in systematically larger Φ c and coefficients that give rise to a sharper decrease in diffusivity as Φ c is approached, compared to the surface growth scenario. The revised relationship was also implemented at continuum scale and used to examine column scale diffusivity change and reactions.

58 GEOSCIENCES↗