Search NASA⌕ Search

SEARCH · Search NASA

Results for “STRESS DISTRIBUTION”

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

Direct observation of strain-enhanced hydrogen segregation and failure at high-angle grain boundaries in nickel

Understanding the mechanisms underlying hydrogen embrittlement remains difficult, even in single-element metals. Both microstructure and stress state influence hydrogen distribution in metals and alloys, which impacts deformation and failure. Here, in this work, we use in-situ Kelvin probe force microscopy (KPFM) to monitor the hydrogen distribution in pure nickel over time at 1.1 % and 3.5 % strain. The sample strained to 3.5 % results in preferential hydrogen segregation to high-angle grain boundaries whereas the sample strained to 1.1 % does not exhibit preferential hydrogen segregation. Optical digital image correlation (DIC) shows that hydrogen charging results in both localized and reduced strains during tensile testing of a notched sample. Later stages of deformation and failure (e.g. microcracking) are studied by using in-situ transmission X-ray microscopy (TXM). TXM reveals nanoscale structural changes to a propagating crack in a hydrogen environment. Localized void growth and secondary cracking occur at grain boundaries near the primary crack front. Correlative electron backscattering diffraction (EBSD) is used to relate the cracking at grain boundaries to the hydrogen segregation observed in KPFM. These findings are unified in a proposed hydrogen embrittlement mechanism that describes the interaction of hydrogen with grain boundaries, and the role of grain boundaries in hydrogen embrittlement.

Grain boundary fracture↗
DOE OSTI2024

Investigation of inter-subject variation in ultrafine particle deposition across human nasal airways: A study involving children, adults, and the elderly

Occupational and environmental exposure to toxic nanoparticles, driven by the rapid expansion of nanotechnology, raises significant respiratory health concern. Numerous studies have explored airflow and particle dynamics in adult nasal airways, but understanding the impact of age-related anatomical changes in children and the elderly remains limited. This study systematically investigates age-related anatomical variations and associated influence on nasal airflow dynamics and ultrafine particle deposition characteristics. Using Computational Fluid-Particle Dynamics (CFPD) method, simulation was conducted under diverse inhalation conditions spanning a wide age range, including: two children (5 years old), two young adults (in their twenties), and two elderly (over 77 years old). Our results reveal distinctive variations across age groups in anatomical dimensions, which affect distribution of wall shear stress where the elderly and children display unique patterns distinct from the young adults. While total deposition efficiency differs significantly between children and adults, filtration efficiency in the subregion with most deposition, main respiratory, remains consistent. However, inter-subject differences are observed in the vestibular and olfactory regions,emphasizing nuanced impact of age-related anatomical variations. Overall and subregional empirical equations for deposition efficiency were developed by incorporating the combined diffusion parameter, Sc a Δ b , corroborating the use of geometrical characteristic parameters for each specific subject in predicting nasal deposition efficiency across age groups. Our findings contribute to predictive nanoparticle exposure analysis in nasal airways across different age groups, thereby enhancing respiratory healthcare for individuals across the life span.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Orientation-Dependent Nanomechanical Behavior of Pentaerythritol Tetranitrate as Probed by Multiple Nanoindentation Tip Geometries

Nanoindentation can be leveraged to aid in the high fidelity modeling of dislocation mediated plasticity in pentaerythritol tetranitrate (PETN), an anisotropic energetic molecular crystal. Moreover, nanoindentation tip parameters such as tip geometry, size, and degree of acuity can be utilized to target anisotropic behavior. In this work, nanoindentation was conducted across a range of orientations on the (110) face of PETN to characterize resultant yield behavior, mechanical property measurements, and resultant slip behavior and fracture initiation. Three different indentation tips were utilized: a 3-sided pyramidal Berkovich tip, a 4-sided high aspect ratio Knoop tip, and a 90° conical tip. Ultimately, indenter tip radius was documented to impact yield behavior, whereas tip geometry affected larger scale processes such as slip, and tip acuity was the dominating factor that led to fracture. The axisymmetric conical tip, serving as a baseline, showed the least amount of variation in mechanical property measurements but also the largest distribution of maximum shear stress at which initial yielding occurred. Its high degree of acuity, however, was more prone to induce fracture at higher loads. The Knoop tip was shown to be suitable for average measurements, but also for elucidation of certain anisotropic features. A distinctly higher perceived hardness at 45° was measured with the Knoop tip, indicating less dislocation motion in that direction also observed in this work via scanning probe microscopy. Lastly, the commonly used Berkovich tip was a good compromise whereby it provided a representative volume element describing the average behavior of the material. These results can be utilized to target desired anisotropic behavior in a wider range of molecular crystals, as well as to inform theoretical considerations for dislocation mediated plasticity in PETN.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Multiscale Mechanisms of Twisted Carbon Nanotube Yarns Probed In Situ by Soft X-rays during Tensile Loading

Piecing together carbon nanotubes (CNTs) into assemblies has so far failed to achieve the same elite strength performance metrics as individual CNTs, highlighting a critical deficiency in understanding the effects that the processing of individual nanostructures have on the performance of their derived macroscale assemblies, thereby hindering the development of a process-structure-performance map for these materials. Here, in this work, we propose a method to decouple the distribution orientation of nanoscale tortuosity and the microscale twist of CNT dry-spun yarns under applied loads via in situ soft X-ray probing at high energy (1200 eV) and low energy (280 eV), respectively. With this decoupling enabled by in situ soft X-ray scattering, we acquired a deeper understanding of the deformation mechanisms of these yarns. We found that for untreated yarns, the twist angle of collective CNT bundles at the macroscale is more sensitive to applied stress than the nanoscale alignment distribution. We also found that increasing nominal twist densities of yarns as well as increased strengthening via plasma treatments and polymer infiltration act to decrease the yarns’ sensitivity to realignment at the nanoscale and prevent failure by the slip mechanism.

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

3D TRISO particle-explicit compact meshing

The TRI-structural ISOtropic (TRISO) layered fuel particle is a robust nuclear fuel form offering enhanced safety and performance for advanced reactor concepts, including high-temperature gas-cooled reactors and other Generation IV designs. These poppy-seed-sized particles are embedded in a graphite matrix to form fuel elements that must withstand elevated temperatures and high burn-up levels. The heterogeneous nature of these fuel elements — comprising thousands of randomly distributed TRISO particles — produces complex stress fields and thermal gradients that one- and two-dimensional models cannot accurately capture. While three-dimensional modeling has improved predictions of dimensional changes, internal pressure buildup, and fission product transport under irradiation, current approaches rely on homogenized material properties that are known to have considerable divergence from experimental observations. This work presents a methodology for optimized random packing of TRISO fuel compacts and full three-dimensional mesh generation within the BISON fuel performance code, with each particle coating layer individually discretized. The resulting mesh was demonstrated through heat conduction simulations under representative in-reactor operating conditions, showing strong agreement with expected behavior. This capability enables detailed analysis of particle-to-particle interactions, matrix cracking mechanisms, and the statistical distribution of coating layer failures — all of which directly govern fuel performance and safety margins.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Plant genotype and rhizobia strain combinations strongly influence the transcriptome under heavy metal stress conditions in Medicago truncatula

Heavy metals such as cadmium (Cd) and mercury (Hg) pose significant threats to plant health and food safety as they are absorbed from the environment. Legumes are generally considered sensitive to heavy metals but possess standing genetic variation for accumulation and tolerance to toxic ions. We conducted a transcriptomic analysis on hydroponically and soil grown Medicago truncatula plants to investigate gene expression responses to Cd and Hg exposure in roots, leaves, and nodules. By using plant genotypes with varying metal tolerance or accumulation levels, we observed distinct clustering of gene ontologies, indicating tissue-specific, genotype-specific, and metal-specific gene expression patterns. Considering the symbiotic relationship between legumes and nitrogen-fixing bacteria, we further examined plant phenotypes and transcriptomes of plant genotypes with contrasting Hg accumulation levels and inoculated them with high or low Hg-tolerant Sinorhizobium medicae strains that have presence-absence variation for a mercury reductase (Mer) operon. Host plants inoculated with the Hg-tolerant rhizobia strain possessing a Mer operon exhibited less reduction in nodule number and plant biomass. A smaller reduction in iron (Fe) distribution in nodules after Hg stress was measured using X-ray Fluorescence (XRF) imaging. Dual transcriptome (host plant and bacteria) analysis of nodules revealed a remarkable decrease in the number of differentially expressed genes (DEGs) and clustering of gene ontologies in plants inoculated with the Hg-tolerant rhizobia strain, including symbiosis related genes. This finding suggests that the Hg-tolerant rhizobia strain has the potential to mitigate Hg stress in host plants. Furthermore, we observed genotype by-genotype interactions between the high Hg accumulating plant genotype and the Hg-tolerant rhizobia strain. These findings provide insights into enhancing plant resilience in contaminated environments through optimizing legume-rhizobia interactions for heavy metal tolerance.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Statistical evaluation of microscale stress conditions leading to void nucleation in the weak shock regime

Here, we investigate the heterogeneity of the stress state driven by anisotropic deformation response at the single crystal level through five statistical volume element (SVE) calculations of polycrystalline BCC tantalum. This work focuses on grain boundaries as a prominent material defect type prone to void nucleation based upon experimental observations of predominantly intergranular void nucleation in this material. The SVEs are constructed to be statistically representative of larger volumes of material and are meshed such that mean and standard deviation of grain size and orientation information is reconstructed. The computational meshes feature hexahedral (brick) elements and smooth conformal grain boundaries where significant stress concentration is known to occur, a tail effect of interest in the extreme events process of dynamic ductile damage. An existing micromechanical crystallographic plasticity model shown to capture the single crystal behavior of BCC tantalum well is used to perform the polycrystal calculations. The model includes representation of the non-Schmid effect of non-planar screw dislocation kinetics in tantalum. A three-dimensional stress state time profile predicted by damage modeling of a flyer plate impact experiment is applied as boundary conditions to each SVE. Resulting grain boundary stress state statistics are strongly non-Gaussian. Significant structural evolution is observed within the compressive hold before unloading into tension in the stress profile. Strong angular dependence of grain boundary traction magnitude with shock direction is observed. Non-Schmid effects continue to suggest their influence on propensity of microstructural defect types to nucleate voids. A general void nucleation criterion is proposed using probability theory. The general framework is specified to polycrystalline BCC tantalum in the weak shock regime to include the SVE calculations and literature molecular dynamics calculations of grain boundary void nucleation strength. Probability density functions (PDFs) are used to describe the interaction between the local stress state heterogeneity and the distributed grain boundary void nucleation strength state. A causation entropy maximization procedure removes the requirement for ad hoc selection of a PDF functional form and provides a rigorous procedure for data-based PDF determination. The resulting physically informed PDF describes the spatial appearance frequency of nucleated voids as a function of applied macroscale pressure. Lower length scale physics are thus packaged in a precise and computationally efficient way to provide computational plasticity insight to macroscale dynamic ductile damage models.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Photovoltaic Cable Jackets: A Comparison of Representative Products Using Combined-Accelerated Stress Testing [Slides]

Photovoltaic (PV) cables facilitate the distribution of electricity collected from modules to the energy grid. Durable cabling enables continuous operation of PV installations, whereas cables with a lifetime less than the modules must be replaced - reducing electricity generation and adding to the operating expense. This study primarily focusses on the aging of the key cable types using the combined-accelerated stress testing (C-AST) protocol. Representative cables for utility, building, and control/auxiliary applications were examined. Cable jacket materials examined include: polyolefin, polyethylene, polyamide, poly(vinyl chloride), chlorinated polyethylene, thermoplastic elastomer, and ethylene propylene diene monomer rubber. Specimen characterizations applied include: optical microscopy, mechanical profilometry, instrumented indentation, scanning electron microscopy (SEM), energy-dispersive X-ray spectroscopy (EDS), Fourier-transform infrared spectroscopy (FTIR), thermogravimetric analysis (TGA), and differential scanning calorimetry (DSC). A variety of performance and durability characteristics were observed, depending on the base material, polymer formulation, and jacket color. The results from C-AST are analyzed and discussed relative to a recent industry survey on electronic balance of system components in addition to a recent study where similar cables were aged using steady state ultraviolet weathering (International Electrotechnical Commission Technical Specification 62788-7-2). Recommendations are made for the screening, industry qualification, and service life prediction of PV cable jackets.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

A Gaussian Process-Based extended Goldak heat source model for finite element simulation of laser powder bed fusion additive manufacturing process

In this study, laser powder bed fusion (L-PBF) additive manufacturing (AM) is a key enabling technology to manufacture highly complex and integrated metallic structures. In L-PBF AM process, the melting of the metal powders and the layers underneath can be governed by either “conduction mode” or “keyhole mode”, with the keyhole mode reportedly leading to porosity and decreased strength and ductility by many studies. In part scale simulations, finite element (FE) model is often used to study the temperature distribution during printing and to predict the residual stress, where a volumetric heat flux with a Gaussian or a double ellipsoidal (Goldak) distribution is often applied as the laser heat source. However, the above heat source models can only capture the melt pool shape in the conduction mode, and fail to capture the transition to keyhole melting mode when the process parameters change. To overcome this inaccuracy, an extended Goldak heat source model is proposed by introducing a laser penetration term as a function of laser parameters obtained from a Gaussian-Process (GP) model. The model is validated by “2D pad” AlSi10Mg L-PBF experiments under a wide range of laser power, scan speed, and laser focus offset, and the results show the model successfully captures the measured melt pool shape in all conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Evaluation of the Simplified Assessment Peak Equivalent Stress Design Limit Probability of Failure

The ASME Boiler and Pressure Vessel Code Section III Division 5 Subsection HH Subpart A (HHA) outlines requirements and qualification methods for graphite core components in high temperature reactor environments. One such structural qualification method is the Simplified Assessment outlined in HHA-3220, which represents a comparison of service level specific equivalent stresses to allowable stress values derived from target probability of failures (POFs). The target POF is based on the structural reliability class (SRC) for the components and service condition of interest, and is used to establish design limits for combined membrane and peak equivalent stresses. Combined membrane stress limits are derived from two parameter Weibull distribution fits of the graphite grade specific tensile strength data and the target POF. The combined membrane stress limits are then scaled based on the ratio of flexural strength to tensile strength (referred hereinafter as Rtf) to establish the corresponding peak equivalent stress limits. It is desired that the target POF is at least maintained for both combined membrane and peak equivalent design limits, especially with the use of Rtf. The assessment presented herein examines tensile and flexural data for medium and fine grade graphites relative to the calculated combined membrane and peak equivalent stress limits per HHA-3220 and illustrates the simplified assessment method results in anticipated, or at least conservative, peak equivalent design stress limits.

36 - MATERIALS SCIENCE↗
DOE OSTI2024

4D Electrical Resistivity Imaging of Stress Perturbations Induced During High–Pressure Shear Stimulation Tests

Fluid flow through fractured media is typically governed by the distribution of fracture apertures, which are in turn governed by stress. Consequently, understanding subsurface stress is critical for understanding and predicting subsurface fluid flow. Although laboratory–scale studies have established a sensitive relationship between effective stress and bulk electrical conductivity in crystalline rock, that relationship has not been extensively leveraged to monitor stress evolution at the field scale using electrical or electromagnetic geophysical monitoring approaches. In this paper we demonstrate the use time–lapse 3–dimensional (4D) electrical resistivity tomography to image perturbations in the stress field generated by pressurized borehole packers deployed during shear–stimulation attempts in a 1.25 km deep metamorphic crystalline rock formation.

58 GEOSCIENCES↗
DOE OSTI2025

Time-series elemental imaging reveals CAX-dependent redistribution patterns for anoxia recovery

Flooding-induced oxygen deprivation (anoxia) is a challenge to plant survival, necessitating adaptive mechanisms for recovery. This study investigated elemental redistribution during anoxia recovery using time-series elemental imaging to show changes in nutrient distribution. Focusing on the role of Cation/H + Exchangers (CAXs) in Arabidopsis thaliana, we show how mutants deficient in specific CAX transporters (cax1 and the cax1-4 quadruple mutant) respond to anoxia and metal stress. Mutants showed reduced lipid peroxidation and increased expression of flood-tolerance proteins during recovery. X-ray fluorescence microscopy and laser ablation–inductively coupled plasma mass spectrometry were used to show elemental redistribution over time. In wild-type plants (Col-0), post-anoxia elemental distribution resembled the elemental distribution of CAX mutants under normoxic conditions, suggesting that CAX-mediated elemental distribution before anoxia enables faster recovery post-anoxia, rather than affecting remobilization post-anoxia. Although CAX mutants had altered tolerance to excess manganese and copper, leaf metal distribution during metal stress was not altered. Here, these findings introduce the potential utility of time-series elemental imaging to show stress-response phenotypes and the importance of elemental distribution to recovery after anoxia. The novelty of this work lies in resolving spatial distribution patterns in a non-static system to gain insight into mechanisms of stress resilience in plants.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Factors controlling injection-induced rupture of intersecting faults during geological sequestration of CO 2

This study addresses coupled multiphase fluid flow and geomechanics effects on potential fault activation associated with subsurface CO 2 injection around intersecting faults. An enhanced fault-representation model is used to capture geomechanical responses of two intersecting faults with finite length during CO 2 injection. The faults are embedded in a strike-slip stress regime of a caprock-reservoir-basement system with the faults represented by zero-thickness interfaces with adjacent finite-thickness damage zones. A sensitivity analysis is conducted to study the effect of fault permeability, slip-weakening behavior, well location relative to the orientation of faults, and well placement (the number and location of injection wells). Five metrics (pressure, CO 2 plume, shear state on the fault, as well as shear displacement and stress path at selected fault monitoring points) are selected to assess CO 2 migration and reactivation of intersecting faults. The results show that induced ruptures are favored by low permeability faults due to high pressure buildup and by slip-weakening behavior resulting from fault strength reduction. The location of one injection well relative to fault orientation determines the magnitude of changes in effective normal stress and shear stress, affecting the location of induced ruptures. Well placement (two injection wells used in the paper) dominates pressure diffusion around the intersection and tips of faults. This redistributes changes in effective normal stress caused by each injection well, influencing the spatial distribution of ruptures along faults. A larger injection volume induces far-field ruptures that are controlled by stress transfer within the injection layer. The findings presented here can provide valuable insights into engineering operations for a long-term, safe, and reliable geologic CO 2 storage.

Fault permeability↗
DOE OSTI2024

Verification of Adaptive Protection in Hardware in the Loop for Coordination with Solar Variability

As inverter-based resources continue to be installed at all levels of the electric grid, the fixed protection schemes used at the distribution level will continue to be stressed until they no longer ensure the protection of the grid. Adaptive protection has been proposed as a solution with the ability to update the protection schemes in near real-time to ensure reliability and increase the resilience of the grid. However, weather variability poses a significant challenge to the ability of these methods to keep the selectivity and reliability of these schemes coordinated. If the ramp rates, due to solar variability, of the inverter-based resources change faster than the adaptive protection can issue new settings, the protection system could be uncoordinated, with the wrong device responding to a system fault. The proposed adaptive protection method ensures that due to solar variability, communication, and protection calculation latency, it can issu e coordinated protection settings promptly, in one minute or less. The hardware-in-the-loop results show the protection settings being issued and maintaining system coordination in under a minute.

Summers, Adam↗
DOE OSTI2024

Tailoring Piezoelectricity of 3D Printing PVDF-MoS 2 Nanocomposite via In Situ Induced Shear Stress

3D printing of unique structures with tunable properties offers significant advantages in the fabrication of complex and customized electronic devices. Here, this study introduces a process-microstructure–property-guided manufacturing route to fabricate PVDF-2D MoS 2 piezoelectric nanocomposites with tunable piezoelectric properties without having a postprocess. We control PVDF’s microstructure through direct ink writing (DIW) 3D printing while tuning PVDF-MoS 2 interfacial strain by controlling rheology and 3D printing parameters, such as nozzle size and printing speed. Our approach demonstrates tunable piezoelectricity in PVDF-MoS 2 , achieving a 15-fold increase in the piezoelectric coefficient (d 33 ) at a printing-induced shear stress of 6685 Pa. This enhancement arises from the electrostatic interactions between PVDF and MoS 2 and the filler distribution and alignment caused by the in situ shear stress in 3D printing, as confirmed by XPS and Raman mapping analyses. Our findings advance the understanding of piezoelectric mechanisms in PVDF-based nanocomposites, laying the foundation for 3D printing of piezoelectric sensors in wearable device applications with enhanced performance and customization capabilities.

2D MoS2↗
DOE OSTI2026

Flaw Tolerance Assessment for DOE Standard SNF Dry Storage Canisters - 26550

The U.S. DOE has designed four spent nuclear fuel (SNF) dry storage canisters for storing DOE standardized SNFs. The DOE standard canisters are cylindrical shells with a diameter of 24 inches (610 m) or 18 inches (457 m), a wall thickness of 0.5 inches (12.7 m) or 0.375 inches (9.53 m), and a length of 15 feet (4.57 m) or 10 feet (3.05 m). These DOE canister geometries are completely different from commercial canisters. The latter may experience chloride-induced stress cracking corrosion (CI-SCC) because they are stored near coastal regions. The former may not experience CI-SCC but face different challenges because they are stored in the SNF storage facilities. Because of large residual stresses, mechanical flaws may occur in the DOE canisters during long-distance transportation or lifting handling. To date, only limited structural integrity analyses were carried out through drop tests on the DOE canisters, but a more general flaw tolerance assessment has not been performed. Therefore, the failure assessment diagram (FAD)-based fracture mechanics method, as codified by the latest API 579-1/ASME FFS-1-2021 Edition, is adopted in this work to assess surface flaw tolerance for DOE canisters under operation loading and welding residual stresses (WRS), where the new code-recommended WRS distributions are used. To more adequately consider the transverse distribution of WRS, an equivalent WRS distribution is proposed to account for the WRS reduction with distance from the weld centerline. Moreover, the closed-form solutions of stress intensity factor K, which serves as the crack driving force during subcritical crack growth, are developed from the tabular data of the K factors provided in API 579-1/ASME FFS-1 and used to determine more accurate flaw sizes at flaw instability. Subsequently, the Level 2 assessment procedures with 12 assessment steps, as codified and detailed in API 579-1 and ASME FFS-1, are followed to assess the flaw tolerance for the surface flaws in the DOE standard canisters with consideration of normal or accident operation loads combined with WRS. The assessment results show that the four designs of DOE standard canisters can tolerate all surface flaws that meet the code permitted maximum sizes of a flaw length of 8 inches (i.e., 200 mm) and a flaw depth of 80% wall thickness. This demonstrates that all designs of DOE standard canisters are robust and reliable.

DOE standard canister↗
DOE OSTI2024

Stressful crystal histories recorded around melt inclusions in volcanic quartz

Abstract Magma ascent and eruption are driven by a set of internally and externally generated stresses that act upon the magma. We present microstructural maps around melt inclusions in quartz crystals from six large rhyolitic eruptions using synchrotron Laue X-ray microdiffraction to quantify elastic residual strain and stress. We measure plastic strain using average diffraction peak width and lattice misorientation, highlighting dislocations and subgrain boundaries. Quartz crystals across studied magma systems preserve similar and relatively small magnitudes of elastic residual stress (mean 53–135 MPa, median 46–116 MPa) in comparison to the strength of quartz (~ 10 GPa). However, the distribution of strain in the lattice around inclusions varies between samples. We hypothesize that dislocation and twin systems may be established during compaction of crystal-rich magma, which affects the magnitude and distribution of preserved elastic strains. Given the lack of stress-free haloes around faceted inclusions, we conclude that most residual strain and stress was imparted after inclusion faceting. Fragmentation may be one of the final strain events that superimposes stresses of ~ 100 MPa across all studied crystals. Overall, volcanic quartz crystals preserve complex, overprinted deformation textures indicating that quartz crystals have prolonged deformation histories throughout storage, fragmentation, and eruption.

58 GEOSCIENCES↗