Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion process”

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 127 records · Page 7

Harnessing nanoscale densification for controlling gas selectivity in flexible zeolitic imidazolate frameworks

Flexible metal-organic frameworks (MOFs) are promising for gas separation due to their molecular recognition capabilities. However, gate-opening behavior in certain flexible MOFs can reduce gas selectivity, enabling simultaneous adsorption of multiple gases. We introduce a binder-free, pressure-free densification process to enhance gas selectivity in flexible zeolitic imidazolate frameworks (ZIF-7 and ZIF-9) by optimizing nanoparticle packing. Unlike powders, the densified ZIF-7 and ZIF-9 retain rhombohedral crystal structures after solvent removal, confirmed via synchrotron powder X-ray diffraction. This process modifies adsorption behavior by introducing a diffusion barrier, enabling diffusivity-driven kinetic separation and enhanced CO2/hydrocarbon selectivity. At 0.5 bar, ZIF-7 mono shows 2×, 4×, and 9.5× higher CO 2 selectivity over C 2 H 2 , C 2 H 4 , and C 2 H 6 , respectively, versus ZIF-7 pwd . ZIF-9 mono similarly shows 6.5×, 3×, and 39× improvements over ZIF-9 pwd . Both monoliths outperform powders in ternary gas mixtures. This study provides a novel strategy to control gate opening and enhance gas selectivity in flexible MOFs for practical separation applications.

CO2/hydrocarbon separation↗

Correlated dynamic disorder, octahedral tilts, and acoustic phonon softening in CsSnBr 3 and CsPbBr 3

Metal halide perovskites (MHPs) have emerged as highly promising materials for optoelectronic applications, with all-inorganic MHPs presenting enhanced stability compared to their hybrid counterparts. Here, in this study, we investigate the atomic dynamics and structural fluctuations in single crystals of CsSnBr⁢ 3 and CsPbBr⁢ 3 through systematic inelastic neutron scattering (INS) measurements as a function of temperature. Our experiments are compared with first-principle simulations, augmented with large-scale molecular dynamics modeling, based on machine-learned neural network potentials. Through both INS and simulations, we find quasi-elastic diffuse rods in reciprocal space in both compounds, originating from fluctuating planar domains featuring correlated tilts of Br octahedron. The diffuse rods exhibit a slow, overdamped dynamic process, modulated across 𝑸 space, reflecting the strong lattice anharmonicity of the inorganic framework. We do not find evidence for dynamic off-centering of the Sn 2+ ions besides phonon vibrations at the center of the Br octahedron. These results offer valuable insights into the unusual anharmonic atomic dynamics and intricate correlated structural distortions in MHPs, which will be critical for rationalizing and further tailoring their thermal and optoelectronic properties.

36 MATERIALS SCIENCE↗

Time-resolved tracking of cellulose biosynthesis and assembly during cell wall regeneration in live Arabidopsis protoplasts

Cellulose, the most abundant polysaccharide on earth composing plant cell walls, is synthesized by coordinated action of multiple enzymes in cellulose synthase complexes embedded within the plasma membrane. Multiple chains of cellulose fibrils form intertwined extracellular matrix networks. It remains largely unknown how newly synthesized cellulose is assembled into an intricate fibril network on cell surfaces. Here, we have established an in vivo time-resolved imaging platform to continuously visualize cellulose biosynthesis and fibril network assembly onArabidopsis thalianaprotoplast surfaces as the primary cell wall regenerates. Our observations provide the basis for a model of cellulose fibril network development in protoplasts driven by an interplay of multiscale dynamics that includes rapid diffusion and coalescence of nascent cellulose fibrils, processive elongation of single fibrils, and cellulose fibrillar network rearrangement during maturation. This study provides fresh insights into the dynamic and mechanistic aspects of cell wall synthesis at the single-cell level.

Science & Technology - Other Topics↗

Morphological Evolution and Dealloying During Corrosion of Ni20Cr (wt.%) in Molten FLiNaK Salts

The dealloying corrosion behavior of the FCC Ni20Cr (wt%) in molten LiF-NaF-KF (FLiNaK) salts at 600 °C under varying applied potentials was investigated. Using in-operando electrochemical techniques and a multi-modal suite of characterization methods, we connect electrochemical potential, thermodynamic stability, and electro-dissolution kinetics to the corrosion morphologies. Notably, under certain potential regimes, a micron-scale bicontinuous structure, characterized by a network of interconnected pores and ligaments riched with the composition of the more noble (MN) element, becomes prominent. At other potentials both MN and less noble (LN) elements dealloy but at different rates. The dealloying process consists of lattice and grain boundary diffusion of Cr to the metal/salt interface, interphase Cr oxidation, accompanied by surface diffusion of Ni to form interconnected ligaments. At higher potentials, the bicontinuous porous structure undergoes further surface coarsening. Concurrently, Cr(II), Cr(III), and Ni(II) begin to dissolve, with the dissolution of Ni occurring at a significantly slower rate. When solid-state transport of Cr is exceeded by the interfacial rates, dealloying depths are limited.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamic Metal–Support Interaction Dictates Cu Nanoparticle Sintering on Al 2 O 3 Surfaces

Nanoparticle sintering remains a critical challenge in heterogeneous catalysis. In this work, we present a unified deep potential (DP) model based on the Perdew–Burke–Ernzerhof approximation of density functional theory for Cu nanoparticles on three Al 2 O 3 surfaces (γ-Al 2 O 3 (100), γ-Al 2 O 3 (110), and α-Al 2 O 3 (0001)). Using DP-accelerated simulations, we reveal that the nanoparticle size-mobility relationship strongly depends on the supporting surface. The diffusion of nanoparticles on the two γ-Al 2 O 3 surfaces is almost independent of the size of the nanoparticle, while the diffusion on α-Al 2 O 3 (0001) decreases rapidly with increasing size. Interestingly, nanoparticles with fewer than 55 atoms diffuse several times faster on α-Al 2 O 3 (0001) than on γ-Al 2 O 3 (100) at 800 K while expected to be more sluggish based on their larger binding energy at 0 K. The diffusion on α-Al 2 O 3 (0001) is facilitated by dynamic metal–support interaction (MSI), where Al atoms move out of the surface plane to optimize contact with the nanoparticle and relax back to the plane as the nanoparticle moves away. In contrast, the MSI on γ-Al 2 O 3 (100) and on γ-Al 2 O 3 (110) is dominated by more stable and directional Cu–O bonds, consistent with the limited diffusion observed on these surfaces. Our extended MD simulations provide insight into the sintering processes, showing that the dispersity of the nanoparticles strongly influences the coalescence driven by nanoparticle diffusion. We observed that the coalescence of Cu 13 nanoparticles on α-Al 2 O 3 (0001) can occur in a short time (10 ns) at 800 K even with an initial internanoparticle distance increased to 3 nm, while the coalescence on the two γ-Al 2 O 3 surfaces are inhibited significantly by increasing the initial internanoparticle distance. These findings demonstrate that the dynamics of the supporting surface is crucial to understanding the sintering mechanism and offer guidance for designing sinter-resistant catalysts by engineering the support morphology.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Water Sorption Study of Various Polymers and Polymeric Composites

In the formulations of many composite materials, polymeric materials are often used as binder materials to enhance processability, functionality, and safety during manufacturing. For instance, poly(ester urethane), such as Estane® 5703, is used in the formulation of PBX 9501, and vinyl copolymer elastomer (VCE) is used in the formulation of VCE/filler compos ites. The stability of these binders directly impacts the overall stability and performance of their composites. Estane and VCE are prone to hydrolytic degradation, with water sorption behavior playing a significant role in their long-term stability and aging behavior. Hence, it is critical to study their water sorption behavior under various humidity and tempera ture conditions. To improve the accuracy and reduce labor-intensive work involved in water sorption experiments, an automated Dynamic Vapor Sorption (DVS) system was employed, which can generate large datasets. To efficiently process and analyze these datasets, Python scripts were developed. These scripts not only streamlined data processing, but also accu rately derived the Arrhenius parameters related to water diffusion and sorption properties for the tested materials. This approach offers a robust framework for evaluating water trans port and sorption processes in polymeric materials and adds valuable capabilities to future water sorption testing efforts.

36 MATERIALS SCIENCE↗

A two-dimensional numerical study of the magneto-Rayleigh–Taylor instability with FLASH: Application to the staged Z-pinch concept

Magnetically driven implosions involving a liner collapsing onto a target are inherently vulnerable to the magneto-Rayleigh–Taylor instability (MRTI). Among the various approaches proposed to achieve fusion conditions within the target, the staged Z-pinch (SZP) concept employs a high-Z liner, the advantages of which remain an active area of investigation. Consequently, ongoing design optimization efforts are essential, while critical physical processes such as magnetic field (B-field) diffusion and radiation transport increase the complexity of required simulations. In this study, we utilize the new capabilities of the FLASH code to simulate in 2D, for the first time using this code, staged Z-pinch configurations (designated SZP1 and SZP1*), focusing on their stability. First, a comparison of simulation results with theoretical predictions of MRTI growth provide new insights into the often-overlooked influences of high mesh resolution and initial perturbation seeding on instability dynamics. These findings then lead to a discussion on potential improvements for future SZP experiments. These include the use of an axial B-field stabilization and the optimization of radiation transport processes during the implosion. The results presented herein establish the framework for simulating multi-dimensional Z pinches using the FLASH code and pave the way for the development of innovative experimental configurations leveraging its advanced simulation capabilities.

Adaptive mesh refinement↗

Curvature effects in the reaction-diffusion system governing frontal polymerization

Frontal polymerization is a fast and energy-efficient curing process of thermoset polymers involving the self-sustained propagation of a sharp front that transforms a liquid monomer or gel into a polymer. The temperature, speed, and geometry of the polymerization front are closely intertwined quantities that can affect the curing process and the mechanical properties of the polymer. This work investigates the effect of the front curvature on the front propagation speed based on a reaction-diffusion model. The proposed analytical relations between front speed and curvature are verified by finite element simulations and validated through experiments. The models show that a concave front has a higher front temperature and speed than its convex counterpart. When the curvature of a convex front reaches a critical value, the reaction can no longer be sustained, and the front is quenched. The relation between the speed of a steadily propagating front and its geometry is derived, showing good agreement with experimental observations. This study also provides analytical estimates of the thickness of the polymerization front and insight into the relationship between the front curvature and the kinetics of the exothermic reaction.

42 ENGINEERING↗

Space-Time Finite Element Tensor Network Approach for the Time-Dependent Convection–Diffusion–Reaction Equation with Variable Coefficients

In this paper, we present a new space-time Galerkin-like method, where we treat the discretization of spatial and temporal domains simultaneously. This method utilizes a mixed formulation of the tensor-train (TT) and quantized tensor-train (QTT) (please see Section Tensor-Train Decomposition), designed for the finite element discretization (Q1-FEM) of the time-dependent convection–diffusion–reaction (CDR) equation. We reformulate the assembly process of the finite element discretized CDR to enhance its compatibility with tensor operations and introduce a low-rank tensor structure for the finite element operators. Recognizing the banded structure inherent in the finite element framework’s discrete operators, we further exploit the QTT format of the CDR to achieve greater speed and compression. Additionally, we present a comprehensive approach for integrating variable coefficients of CDR into the global discrete operators within the TT/QTT framework. The effectiveness of the proposed method, in terms of memory efficiency and computational complexity, is demonstrated through a series of numerical experiments, including a semi-linear example.

convection–diffusion–reaction equation↗

Emergent facilitation and glassy dynamics in supercooled liquids

In supercooled liquids, dynamical facilitation refers to a phenomenon where microscopic motion begets further motion nearby, resulting in spatially heterogeneous dynamics. This is central to the glassy relaxation dynamics of such liquids, which show super-Arrhenius growth of relaxation timescales with decreasing temperature. Despite the importance of dynamical facilitation, there is no theoretical understanding of how facilitation emerges and impacts relaxation dynamics. Here, we present a theory that explains the microscopic origins of dynamical facilitation. We show that dynamics proceeds by localized bond-exchange events, also known as excitations, resulting in the accumulation of elastic stresses with which new excitations can interact. At low temperatures, these elastic interactions dominate and facilitate the creation of new excitations near prior excitations. Using the theory of linear elasticity and Markov processes, we simulate a model, which reproduces multiple aspects of glassy dynamics observed in experiments and molecular simulations, including the stretched exponential decay of relaxation functions, the super-Arrhenius behavior of relaxation timescales as well as their two-dimensional finite-size effects. The model also predicts the subdiffusive behavior of the mean squared displacement (MSD) on short, intermediate timescales. Furthermore, we derive the phonon contributions to diffusion and relaxation, which when combined with the excitation contributions produce the two-step relaxation processes, and the ballistic–subdiffusive–diffusive crossover MSD behaviors commonly found in supercooled liquids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Decoupling Li out-diffusion and surface diffusion in the lithiation-assisted epitaxial growth of lithium tungstate

Lithiation-assisted epitaxy offers a flexible and robust approach for synthesizing high-quality Li-containing materials and interfaces with precise control. Here, in this study, we use lithium tungstate (Li x WO 3+x/2 , where x = 0 to 2) as a model system to investigate the intertwined effects of Li out-diffusion-induced compositional changes and surface-diffusion-induced morphological changes. By systematically varying synthesis and processing conditions, we uncover their impact on lithium tungstate film formation. Comprehensive characterizations, including X-ray diffraction, atomic force microscopy, X-ray photoemission spectroscopy and time-of-flight secondary ion mass spectrometry, reveal that low-temperature growth (< 300 °C) followed by high-temperature annealing yields continuous lithium tungstate films with significantly reduced surface roughness. In contrast, high-temperature deposition (≥ 300 °C) accelerates surface diffusion and Li out-diffusion, leading to island formation. Furthermore, in situ scanning transmission electron microscopy demonstrates the beam sensitivity of Li 2 WO 4 and reveals a phase transition from Li 2 WO 4 to LiWO 3.5 under prolonged electron beam exposure. These findings deepen our understanding of how to control composition and morphology of Li-containing films, providing valuable insights for the design and integration of energy materials.

Shi, Jueli [Pacific Northwest National Laboratory ↗

Modification of Ni-20Cr corrosion dealloying behavior in molten fluorides via cold work induced plastic deformation

The corrosion dealloying behavior of cold-worked (CW) Ni20Cr alloy (wt%) was studied in molten LiF-NaF-KF (or FLiNaK) salts at 600 °C, equal to a homologous temperature (TH) of 0.52. Alloys were cold-rolled to achieve reductions of thickness of 10%, 30%, and 50% introducing plastic deformation and a high density of dislocations. Potentiostatic holds (Eapplied) were applied in two different electrode potential regimes. At 1.75VK+/K$${\text{V}}_{{\text{K}}^{+}/{\text{K}}}$$, Cr dealloying to Cr(II) and Cr(III) is predominant, while at 1.90 VK+/K$${\text{V}}_{{\text{K}}^{+}/{\text{K}}}$$, both Ni and Cr are oxidized in molten FLiNaK at 600 °C. In these potential regimes, dealloyed Ni20Cr displayed bicontinuous porosity within the grain interior and at grain boundaries, driven by the high driving force for Cr dissolution and sustained by defect mediated outward solid state diffusion of Cr in parallel with surface diffusion of Ni. The bicontinuous porous structure developed was observed to undergo further coarsening and densification of the Ni-rich ligaments at a higher electrode potential. The main effect of CW observed is the introduction of plastic deformation and dislocation substructures that serve as short-circuit paths for Cr solid state diffusion to surfaces exposed to FLiNaK. This modified the evolution of the bicontinuous porous structure which increased with CW substantially. Kinetic analysis reveals that the Cr dealloying at +1.75 VK+/K$${\text{V}}_{{\text{K}}^{+}/{\text{K}}}$$ and 1.90 VK+/K$${\text{V}}_{{\text{K}}^{+}/{\text{K}}}$$ is initially charge transfer controlled, except for the 50% CW condition at 1.90 VK+/K$${\text{V}}_{{\text{K}}^{+}/{\text{K}}}$$, where the process becomes limited by slow Cr defect mediated bulk diffusion. The rate determining factors are explored and compared to experimental results.

Chan, Ho Lun↗

Enhanced accuracy through ensembling of randomly initialized auto-regressive models for dynamical systems

Computational mechanics simulations using traditional finite element methods (FEM) require prohibitively expensive computational resources for real-time engineering applications, design optimization, and digital twin implementations. While machine learning (ML) surrogate models offer significant computational speedups, autoregressive ML models for time-dependent mechanical systems suffer from error accumulation that compromises long-term prediction reliability - a critical concern for engineering applications where accuracy over extended time horizons is essential for safety and performance assessments. Here, we propose a deep ensemble framework specifically designed to address this challenge in computational mechanics applications, where multiple ML surrogate models with random weight initializations are trained in parallel and their predictions aggregated during inference. This approach leverages statistical diversity to maximize information gain from a fixed set of training data and to mitigate error propagation, while maintaining the computational efficiency that makes ML surrogates attractive for engineering practice. We validate the framework on three representative problems spanning critical areas of computational mechanics: stress field evolution in heterogeneous microstructures under complex loading (relevant to advanced materials design and composite analysis), planetary-scale shallow water dynamics (applicable to environmental and geotechnical engineering), and Gray-Scott reaction-diffusion systems (relevant to mass transport and chemical process engineering). Across all test cases, the ensemble approach demonstrates consistent error reduction of 15-33% compared to individual models. The codes for this work are available on GitHub (https://github.com/Graham-Brady-Research-Group/AutoregressiveEnsemble_SpatioTemporal_Evolution).

autoregressive prediction↗

Electrode Modification by Laser Ablation to Overcome Mechanical and Scalability Limitations of Solid-State Prealkylation of Electrode Architectures

In this work, laser ablation of Si electrode surfaces is utilized to provide localized strain relief during the Solid State Prealkylation of Electrode Architectures (SPEAR) process, in addition to volumetric expansion relief in the form of free volume channels across the electrode. A hexagonal mesh patterning was employed via laser ablation, where the hexagonal diameter (pitch) was varied from 250 to 1000 µm. It was demonstrated that these channels serve as expansion voids, accommodating the increase in electrode thickness observed during SPEAR and preventing cracking at the electrode surface. Optical microscopy and Raman mapping confirm that Li diffusion is the rate-limiting solid-state process, requiring a 24 h rest period to achieve homogeneous LixSi alloying through the electrode bulk. Electrochemical evaluations in full cells show that a 1000 μm pitch pattern yields a high initial discharge capacity of 2485 mA h g–1 and 86.4% capacity retention after 80 cycles, whereas smaller pitches (<500 μm) result in lithium plating at 4.1 V in full cells paired with NMC 811 due to excessive current collector exposure from the ablation process. XPS analysis reveals that the laser treatment modifies the electrode surface by increasing the Si–O and C–O species, changing the SEI composition observed via an increase in fluorophosphates upon cycling.

Musgrove, Amanda [ORNL] (ORCID:0009000220910389)↗

Interfacial Reactivity-Triggered Oscillatory Lattice Strains of Nanoalloys

Understanding the structure evolution of nanoalloys under reaction conditions is vital to the design of active and durable catalysts. Herein, we report an operando measurement of the dynamic lattice strains of dual-noble-metal alloyed with an earth-abundant metal as a model electrocatalyst in a working proton-exchange membrane fuel cell using synchrotron high-energy X-ray diffraction coupled with pair distribution function analysis. Additionally, the results reveal an interfacial reaction-triggered oscillatory lattice strain in the alloy nanoparticles upon surface dealloying. Analysis of the lattice strains with an apparent oscillatory irregularity in terms of frequency and amplitude using time-frequency domain transformation and theoretical calculation reveals its origin from a metal atom vacancy diffusion pathway to facilitate realloying upon dealloying. This process, coupled with surface metal partial oxidation, constitutes a key factor for the nanoalloy's durability under the electrocatalytic oxygen reduction reaction condition, which serves as a new guiding principle for engineering durable or self-healable electrocatalysts for sustainable fuel cell energy conversion.

Wu, Zhi-Peng [State University of New York (SUNY),↗

Emergent hydrodynamic mode on SU(2) plaquette chains and quantum simulation

We search for emergent hydrodynamic modes in real-time Hamiltonian dynamics of 2+1-dimensional SU(2) lattice gauge theory on a quasi-one-dimensional plaquette chain, by numerically computing symmetric correlation functions of energy densities on lattice sizes of about 20 with the local Hilbert space truncated at 𝑗 max = $\frac{1}{2}$. Because of the Umklapp processes, we only find a mode for energy diffusion. The symmetric correlator exhibits transport peak near zero frequency with a width approximately proportional to momentum squared at small momentum, when the system is fully quantum ergodic, as indicated by the eigenenergy level statistics. This transport peak leads to a power-law 𝑡 −$\frac{1}{2}$ decay of the symmetric correlator at late time, also known as the long-time tail, as well as diffusionlike spreading in position space. We also introduce a quantum algorithm for computing the symmetric correlator on a quantum computer and find it gives results consistent with exact diagonalization when tested on the IBM emulator. Finally we discuss the future prospect of searching for the sound modes.

Hamiltonian systems↗

Modeling Microwave-Enhanced Chemical Vapor Infiltration Process for Preventing Premature Pore Closure

The chemical vapor infiltration (CVI) process involves infiltrating a porous preform with reacting gases that undergo chemical transformation at high temperatures to deposit the ceramic phase within the pores, ultimately leading to a dense composite. The conventional CVI process in composite manufacturing needs to follow an isothermal approach to minimize temperature differences between the external and internal surfaces of the preform, ensuring that reactive gases infiltrate internal pores before external surfaces seal. Here, this study addresses the challenge of premature pore closure in CVI processes through microwave heating. A frequency-domain microwave solver is developed in OpenFOAM to investigate volumetric heating mechanisms within the preform. Through numerical studies, we demonstrate the capability of microwave heating of creating an inside-out temperature inversion. This inversion accelerates reactions proximal to the preform center, effectively mitigating the risk of premature external pore closure and ensuring uniform densification. The results reveal a significant enhancement in temperature inversion when high-permittivity reflectors are incorporated to generate resonant waves. This microwave heating strategy is then coupled with high-fidelity direct numerical simulation (DNS) of reacting flow, enabling the analysis of resulting densification processes. The DNS includes detailed chemistry and realistic diffusion coefficients. The numerical results can be used to estimate the impact of microwave-induced temperature inversion on densification in productions.

42 ENGINEERING↗

Biochar-compost-based controlled-release nitrogen fertilizer intended for an active microbial community

Nitrogen (N) fertilizers in agriculture suffer losses by volatilization of N to the air, surface runoff and leaching into the soil, resulting in low N use efficiency (NUE) (< 50%) and raising severe environmental pollutions. Controlled- release nitrogen fertilizers (CRNFs) can control the release of N nutrients to NUE in crop production. Different methods were used to develop new CRNFs. However, different CRNF technologies are still underdeveloped due to inadequate controlling on N releasing time and/or unsustainable diffusion. The study on the influences of CRNF processing parameters on microbial conditions are lacking when the CRNFs composed of various bio-ingredients such as biochar, composts, and biowaste. The complexity of processing methods, material biodegradability, and other physical properties make current CRNFs of questionable value in agricultural production. This research aims to develop a novel biochar-compost-based controlled-release urea fertilizer (BCRUF) to preserve microbial properties carried by the compost. The BCRUF was synthesized by pelletizing the 50:50 (dry, wt/wt) mixture of biochar and compost. BCRUF was loaded with urea and then spray-coated with polylactic acid (PLA). The releasing time of two types of BCRUFs, coated and uncoated with PLA, for 80% of N release in water was up to 6 h at three different temperatures (4, 23, and 40 °C), compared to conventional urea fertilizer and commercial environmentally smart N (ESN) fertilizer. The releasing time of coated BCRUF for 80% N release in soil was up to 192 h (8 days). Fourier-transform infrared spectroscopy (FTIR) analysis revealed that no new functional groups were found in the release solution, indicating no new chemical hazards generated. The differential scanning calorimetry (DSC) tests also verified that its thermal stability could be up to 160 °C. The microbe populations in the BCRUF pellets were reduced after the pelleting and drying processes in BCRUF fabrication, but a few bacteria can endure in the air-drying process. BCRUF pellets soaked in water for 4 days retained some bacteria. The BCRUF showed very promising characteristics to improve NUE and sustainability in agricultural production.

58 GEOSCIENCES↗