SEARCH · Search NASA

Results for “Mass diffusivity”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2025

The dark matter diffused supernova neutrino background

We consider neutrinos scattering off Milky Way dark matter and the impact of this scattering on supernovae neutrinos. This can take the form of attenuation on the initial flux of neutrinos and a time-delayed flux of scattered neutrinos. Considering dark matter masses above 100 MeV and past Milky Way supernovae, we find this time-delayed flux is nearly constant in time. We call this flux the Dark Matter Diffused Supernova Neutrino Background (DMDSNB), and use Super-K limits on the Diffuse Supernova Neutrino Background (DSNB) flux to set limits on the dark matter-neutrino scattering cross section. We find σ DM-ν /m DM ≲ 2.4 × 10 -24 cm 2 /GeV for m DM ≳ 1 GeV, which is the strongest bound to date on dark matter-neutrino scatterings at MeV energies, and stronger than bounds set from SN1987A neutrino attenuation by an order of magnitude. We end by discussing how the DMDSNB could be distinguished from the DSNB.

dark matter theory↗

DOE OSTI2024

Weak Gravitational Lensing around Low Surface Brightness Galaxies in the DES Year 3 Data

We present galaxy-galaxy lensing measurements using a sample of low surface brightness galaxies (LSBGs) drawn from the Dark Energy Survey Year 3 (Y3) data as lenses. LSBGs are diffuse galaxies with a surface brightness dimmer than the ambient night sky. These dark-matter-dominated objects are intriguing due to potentially unusual formation channels that lead to their diffuse stellar component. Given the faintness of LSBGs, using standard observational techniques to characterize their total masses proves challenging. Weak gravitational lensing, which is less sensitive to the stellar component of galaxies, could be a promising avenue to estimate the masses of LSBGs. Our LSBG sample consists of 23,790 galaxies separated into red and blue color types at g - i ≥ 0.60 and g - i < 0.60 , respectively. Combined with the DES Y3 shear catalog, we measure the tangential shear around these LSBGs and find signal-to-noise ratios of 6.67 for the red sample, 2.17 for the blue sample, and 5.30 for the full sample. We use the clustering redshifts method to obtain redshift distributions for the red and blue LSBG samples. Assuming all red LSBGs are satellites, we fit a simple model to the measurements and estimate the host halo mass of these LSBGs to be log(M host /M ⊙ ) = $12.98^{+0.10}_{-0.11}$. We place a 95% upper bound on the subhalo mass at log(M sub /M ⊙ ) < 11.51. By contrast, we assume the blue LSBGs are centrals, and place a 95% upper bound on the halo mass at log(M host /M ⊙ ) < 11.84. We find that the stellar-to-halo mass ratio of the LSBG samples is consistent with that of the general galaxy population. This work illustrates the viability of using weak gravitational lensing to constrain the halo masses of LSBGs.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Dust-void evolution driven by turbulent dust flux can induce runaway migration of Earth-mass planets

Torques from asymmetric dust structures (so-called dust-void and filamentary structures) formed around low-mass planets embedded in a nonturbulent dust-gas disk can exceed the torques produced by the gas disk component and then go on to dominate the planet’s orbital dynamics. Here, we investigate how these structures (hence the dust torque) change when the effect of turbulent dust diffusion and dust feedback are included, along with the direct implications on the migration of Earth-like planets. Using the FARGO3D code, we performed 2D and 3D multifluid hydrodynamic simulations, focusing on a non-migrating planet with a mass of M p = 1.5 M ⊕ in 2D and on migrating planets with M p ∈ [1.5, 12] M ⊕ in 3D. We varied the δ-dimensionless diffusivity parameter in the range [0, 3 × 10 −3 ] and considered three different Stokes numbers, St = {0.04, 0.26, 0.55}, which are representative of the gas-dominated, the transitional, and the gravity-dominated regimes, respectively. In our 2D models, we find that turbulent diffusion of dust prevents the formation of the dust-void and filamentary structures when δ > 3 × 10 −4 . Otherwise, dust structures survive turbulent diffusion flow. However, dust and total torques become positive only in transitional and gravity-dominated regimes. In our 3D models, we find that the dust-void is drastically modified and the high-density ring-shaped barrier delineating the dust-void disappears if δ ≳ 10 −4 , due to the effect of dust turbulent diffusion along with the back-reaction of the dust. For all values of δ, the filament in front of the planet is replaced by a low-density trench. Remarkably, as we have allowed the planets to migrate, the evolving dust-void can drive either runaway migration or outward (inward) oscillatory-torque migration. Our study thus suggests that low-mass Earth-like planets can undergo runaway migration in dusty disks.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Tandem bulk oxygen diffusion and surface reactions in reducible metal oxides control redox cycle dynamics

The interplay between bulk oxygen diffusion and surface reactions in reducible metal oxides is key in heterogeneous catalysts, but direct measurements of oxygen mobility, transient kinetics, and in situ spectroscopies have been lacking. Here, we reveal complex dynamic behavior of ceria-zirconia by H 2 using transient kinetics via mass spectrometry and in situ Raman and near-ambient pressure x-ray photoelectron spectroscopies. Molecular dynamics simulations with a machine learning potential delineate competitive oxygen diffusion mechanisms, with an optimal mobility at intermediate reductions. We expose a compensation between vacancy availability and lattice distortion at intermediate to high reductions and Frenkel defects at low reductions, underscoring a potential deficiency of 16 O/ 18 O exchange experiments in deducing oxygen mobility. Vacancies in proximity require electron localization on Ce atoms further away. The continuous replenishment of surface oxygen results in a varying reduction rate, with H 2 dissociation being the rate-limiting step. Multiscale transient simulations, consistent with experiments, indicate catalysts of potentially spatially varying oxidation states. The approach is broadly applicable to reducible oxide materials.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Exploring nuclear structure with multiparticle azimuthal correlations at the LHC

Details of the nuclear structure of Xe, such as the quadrupole deformation and the nuclear diffuseness, are studied by extensive measurements of anisotropic-flow-related observables in Xe–Xe collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{NN}}$ = 5.44 TeV with the ALICE detector at the LHC. The results are compared with those from Pb–Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for a baseline, given that the 208 Pb nucleus exhibits a very weak deformation. Furthermore, comprehensive comparisons are performed with a state-of-the-art hybrid model using IP-Glasma+MUSIC+UrQMD. It is found that among various IP-Glasma+MUSIC+UrQMD calculations with different values of nuclear parameters, the one using a nuclear diffuseness parameter of α 0 = 0.492 and a nuclear quadrupole deformation parameter of β 2 = 0.207 provides a better description of the presented flow measurements. These studies represent the first systematic exploration of nuclear structure at TeV energies, utilizing a comprehensive set of anisotropic flow observables. The measurements serve as a critical experimental benchmark for rigorously testing the interplay between nuclear structure inputs and heavy-ion theoretical models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

High resolution identification and quantification of diffuse deep groundwater discharge in mountain rivers using continuous boat-mounted helium measurements

Discharge of deeply sourced groundwater to streams is difficult to locate and quantify, particularly where both discrete and diffuse discharge points exist, but diffuse discharge is one of the primary controls on solute budgets in mountainous watersheds. The noble gas helium is a unique identifier of deep groundwater discharge because groundwater with long residence times is commonly enriched in helium. In this study, a portable mass spectrometer was used to measure longitudinal variation in dissolved helium concentrations in two mountainous rivers at high spatial resolution not feasible with traditional sampling techniques. Helium profiles were then simulated using a mass-balance model to quantify longitudinal variation in groundwater discharge to the receiving rivers. Results indicate helium concentrations were enriched by multiple orders of magnitude above atmospheric equilibrium in both rivers and that this persisted for up to 18 km below observed pulse inputs in the Colorado River. Helium mass-balance models match observed longitudinal patterns with the exception of sharp initial increases in helium observed in the rivers. Increased longitudinal groundwater discharge rates correspond to mapped geologic structures in both watersheds that likely transport deep geothermal water. Models show variable sensitivity to spatial assignment of input variables representing the groundwater source, illustrating the importance of collecting data from discrete groundwater discharges where possible. The methodology shows promise for field experiments designed to assess air–water exchange rates and to quantify total groundwater discharge from a combination of discrete and diffuse sources.

Deep groundwater↗

DOE OSTI2025

Lithium‐Containing Hybrid SEI Layer Enabling High Mass Loading and Anode‐Less Sodium Metal Batteries

Abstract The continuous rupturing and rebuilding of unstable solid electrolyte interphase (SEI) layer during cycling would block Na + diffusion and induce Na dendrite formation, ultimately limiting the practical application of high‐energy‐density sodium metal batteries. Herein, a hybrid SEI layer containing Li‐species is dexterously constructed on the surface of sodium metal anode. Li‐containing inorganic components (Li 3 N, LiF, and Li 2 CO 3 ) are introduced to stabilize the Na/electrolyte interface and enhance the mechanical and diffusion kinetic properties of the SEI layer, which can reduce the side reactions and gas generation, regulate Na + flux during cycling and promote rapid Na + migration for uniform dendrite‐free Na deposition. As a result, the constructed Na symmetric cells achieve low overpotential and long cycle life of 5900, 1800, and 500 h at current densities of 3, 10, and 30 mA cm −2 , respectively. Furthermore, the full cells paired with the Na₃V₂(PO₄)₃ cathode demonstrate high specific capacity and excellent cycle stability, even at an ultra‐high cathode loading of 39.3 mg cm −2 and a low N/P ratio (negative/positive electrode capacity ratio of 1.21).

Xia, Li [State Key Laboratory of Physical Chemistr↗

DOE OSTI2026

Thermodynamic Control of Interface Directs MnO 2 Nucleation Chemistry for Dense and Conformal Electrodeposition

Manganese dioxide (MnO 2 ) is widely recognized as a promising material for high-energy-density energy storage systems due to its broad applicability and facile electrodeposition. However, achieving uniform, thin, and high-mass-loading MnO 2 coatings on high-surface-area electrodes remains a significant challenge. Conventional electrodeposition methods typically yield nonuniform, thick layers with poor conductivity and limited material utilization, restricting their practical use. Here, we uncover a thermodynamically engineered vanadyl/pervanadyl (VO 2+ /VO 2 + ) interface that fundamentally reshapes MnO 2 electrodeposition chemistry, enabling highly uniform and dense coatings. Here, combining in situ AFM measurement, Classical Nucleation Theory, and Johnson–Mehl–Avrami–Kolmogorov modeling, we show that this interface reduces early-stage detectable MnO 2 island size by 35-fold and shifts the MnO 2 growth from diffusion-limited to reaction-limited progressive nucleation. This thermodynamically controlled interface yields highly dense and conformal MnO 2 films with record-high mass loading of 241 mg cm –2 (1607 mg cm –3 ) on 3D-printed graphene aerogels, without compromising porosity or inducing thickness gradient. As a prototype demonstration, the resulting MnO 2 electrodes deliver record-setting volumetric performance in both capacitors (106 F cm –3 ) and Zn//MnO 2 pouch cells (162 mAh cm –3 ). Beyond energy storage, our findings demonstrate the significance of thermodynamic interface control in MnO 2 nucleation chemistry for achieving dense and uniform coatings on various substrates, with implications for electrocatalysis, semiconductor processing, and advanced materials manufacturing.

Batteries↗

DOE OSTI2025

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 ↗

DOE OSTI2026

Observation of the Jet Diffusion Wake Using Dijets in Heavy-Ion Collisions

Energetic quarks and gluons traversing a hot and dense quark-gluon plasma deposit energy and momentum into the medium before hadronizing to collimated sprays of particles, known as jets. This energy-momentum deposition is expected to produce medium responses, collectively known as jet wakes, with “diffusion wake” denoting a depletion of particles in the direction opposite to the propagating jet. The diffusion wake is studied by comparing dijet-hadron correlations measured in lead-lead (PbPb) and proton-proton (pp) collisions. The analysis uses PbPb and pp data recorded at a nucleon-nucleon (NN) center-of-mass energy s NN = 5.02 TeV with the CMS detector at the CERN LHC. By exploring how the dijet-hadron correlation distributions differ for various pseudorapidity separations of the two jets in the dijet, the presence of a jet diffusion wake is firmly established. The diffusion wake signal has a significance greater than 5 standard deviations for charged particles in the transverse momentum range 1 < p T < 2 GeV . The measurements are compared with model predictions with and without jet wake effects, providing new insights into quark-gluon plasma properties and the formation of jet-induced wakes.

Hayrapetyan, Aram [Yerevan Phys. Inst.]↗

DOE OSTI2025

Microstructural features and deuterium diffusion in lithium penta-aluminate pellets under He + and D + ion irradiation

Lithium (Li) penta-aluminate (LiAl 5 O 8 ) is investigated as a potential tritium (T) breeding material, with a focus on microstructural response to ion irradiation and deuterium (D) diffusion behavior. Under high-fluence ion irradiation (2 x 10 17 (He + +D + )/cm 2 ) at 773 K, LiAl 5 O 8 exhibits significant disorder on the Li sublattice, as revealed by atomic-resolution scanning transmission electron microscopy, while the Al and O sublattices remain stable, demonstrating strong resistance to structural amorphization. Irradiation induces the formation of platelet-shaped antiphase boundaries (APBs), which may serve as effective D trapping sites. Atom probe tomography suggests the presence of 6 LiD clusters in the mass spectra, though definite conclusions regarding APB composition are hindered by signal overlap and limited data statistics. Time-of-flight secondary ion mass spectrometry reveals that D retention approaches to saturation at 3 x 10 17 (He + +D + )/cm 2 . Isothermal and isochronal annealing studies determine an average diffusivity of 1.6 x 10 -13 at 773 K and an effective activation energy of 0.8 ± 0.1 eV for D migration. Compared to γ-LiAlO 2 , LiAl 5 O 8 demonstrates superior irradiation resistance, minimal Li loss, and enhanced D retention, underscoring its potential as a durable breeder material for T production. In conclusion, these findings provide key insights into the microstructural evolution, defect dynamics, and D retention mechanisms in LiAl 5 O 8 under reactor-relevant conditions.

42 ENGINEERING↗

DOE OSTI2024

Oxide dissolution and oxygen diffusion scenarios in niobium and implications on the Bean–Livingston barrier in superconducting cavities

We generalize a native Nb2O5 dissolution model [G. Ciovati, Appl. Phys. Lett. 89, 022507 (2006)] to sequential overlayer dissolutions, multilayer dissolution, and realistic temperature profiles, which may be applicable to other materials. The model is applied to secondary ion mass spectrometry depth profile measurements for varying temperature profiles and two-step oxide dissolution in Nb and found to agree well. In the context of the Meissner screening response due to impurity profiles on the length scale of the London penetration depth, the shallow diffusion of O impurities results in a substantial decrease in the peak supercurrent density near the surface. In this framework, oxide dissolution and oxygen diffusion can account for a rise in peak supportable magnetic field in SRF cavities with baking time and a suppression after the optimal baking time is reached, in good agreement with peak-field baking temperatures and times as well as recent quench field measurements.

Lechner, Eric↗

DOE OSTI2026

Substrate-Mediated Evaporation and Stochastic Evolution of Supported Au Nanoparticles

Here, we use in situ transmission electron microscopy with automated tracking to study supported gold nanoparticles (NPs) during high-temperature vacuum annealing. The average mass loss per NP is governed by a flat, nearly size-independent substrate-mediated evaporation profile. On top of this mean shrinkage, individual NPs show significant fluctuations in apparent growth or shrinkage, and NP volume follows a random-walk-like trajectory. To rationalize both the ensemble-mean behavior and the particle-resolved variability, we develop a self-consistent theory that couples substrate-mediated evaporation to collective 2D Ostwald-type mass exchange through a shared adatom field, described in terms of a renormalized screening length and background concentration. In the experimentally relevant regime, the theory predicts an approximately size-independent mean shrinkage rate and clarifies how net mass loss suppresses classical coarsening. Superimposed on this deterministic drift, we quantify stochastic volume trajectories and capture their fluctuation spectrum with a minimal Langevin description consistent with intermittent adatom attachment and detachment events. In addition, we characterize the lateral diffusive motion of NPs, which is responsible for their coalescence. Altogether, our results highlight that stochasticity is intrinsic at the nanoscale and that predicting the evolution of supported NPs at early and intermediate times requires a unified framework combining substrate-mediated evaporation, collective mass exchange, and stochastic fluctuations.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2026

Iron Oxide Reduction Rate Affects Iron Product Morphology

Magnetite concentrates from different parts of the Mesabi Range, exhibiting varying tendencies to encapsulate, were reduced using thermogravimetric analysis (TGA) and a high-temperature confocal scanning laser microscope (CSLM). Gas-phase mass-transfer conditions were intentionally varied by employing two reactors with different geometries: the TGA produced significantly slower reduction, whereas the CSLM setup yielded reduction rates up to 20 times faster under mass-transfer-controlled conditions. A preliminary mechanism is proposed to explain the transition from a porous iron product to a dense encapsulating layer. Encapsulation is interpreted as the outcome of competition between pore creation, driven by the reduction rate, and pore elimination, driven by surface diffusion. When pore formation is insufficient to counteract pore coarsening, a dense iron layer develops and restricts further reduction.

encapsulation↗

DOE OSTI2025

Removal of Selenium Oxyanions from Aqueous Solutions by Ion Exchange: Equilibrium, Kinetics, and Mechanistic Modeling

Selenium (Se) is an essential micronutrient but toxic at high concentrations, posing challenges for water treatment. This study investigated the removal of selenate (SeO 4 2– ) and selenite (SeO 3 2– ) using the strong-base anion-exchange resin IRA-900, particularly in the presence of competing sulfate (SO 4 2– ). The performance of the commercially available resin IRA-900 was systematically investigated. The batch equilibrium behavior was studied in both single- and binary-component systems, and the kinetic behavior was investigated in single-component systems. Results confirmed a selectivity order of SeO 4 2– > SO 4 2– > SeO 3 2– , indicating preferential SeO 4 2– removal over competing SO 4 2– but lower affinity for SeO 3 2– . The maximum total exchange capacity was determined to be 2.04 mequiv/g. Furthermore, SeO 3 2– uptake was found to be pH-dependent, whereas SeO 4 2– uptake remained stable across a broad pH range. From a modeling perspective, the Law of Mass Action model effectively described equilibrium data, and a transport–reaction modeling framework captured removal kinetics of oxyanions including film and intraparticle diffusion. Finally, X-ray photoelectron spectroscopy confirmed ion exchange between chloride and Se oxyanions as the primary removal mechanism. These findings provide fundamental insights into the removal of Se oxyanions from aqueous solutions by ion exchange.

IRA-900↗

DOE OSTI2026

Comparing Tandem Cell Designs for Electrochemical CO 2 Reduction to Ethylene

Electrochemical carbon dioxide reduction (CO 2 R) is a promising approach for the decentralized production of fuels such as ethylene (C 2 H 4 ). However, the use of Cu, the most efficient metal CO 2 R catalyst for the generation of C 2 H 4 known to date, generally yields a product stream with poor selectivity. In an effort to increase selectivity, the reaction from CO 2 to C 2 H 4 can be broken down into two steps using tandem CO 2 R electrolyzers: formation of CO from CO 2 and subsequent reduction of CO to C 2 H 4 . Here, in this study, we present two novel tandem electrolyzer architectures that closely integrate two cathodes, one for CO generation and one for conversion to C 2 H 4 , while still enabling independent electrical control of the cathodic surfaces. Cathode segmentation in each of these designs also permits the controlled sequencing of mass flow of chemical intermediates in the order of Au to Cu cathode catalysts, in contrast to earlier work relying on uncontrolled, passive diffusion to facilitate the flow of chemical intermediates between catalysts. When comparing the performance of the newly developed electrolyzer cell designs with a dual electrolyzer system, we found that the dual electrolyzer system yields the highest C 2 H 4 faradaic efficiencies (FEs) of 31% and C 2 H 4 concentrations (∼8 mol %). However, a single Cu-containing electrolyzer outperformed all three tandem systems in terms of C 2 H 4 FE (34%). Our findings, enabled by independent control of the two tandem cathode surfaces, indicate that tandem CO 2 R systems need to be evaluated carefully by testing them at various relevant current densities.

C2H4↗

DOE OSTI2026

Boosting CO2R Performance of Ag Electrocatalysts by Sulfur-Doped Carbon Support

We find that S-doped carbon support can boost the CO2 reduction (CO2R) performance of Ag electrocatalysts. Firstly, surface science enabled electrocatalysis showed that Ag supported on S-doped highly oriented pyrolytic graphite (HOPG), a model electrocatalyst, demonstrated 100% higher CO turnover frequency (TOFCO = 3.6 ± 0.2 CO/atomAg/s) than that supported on S-free HOPG (TOFCO = 1.8 ± 0.2 CO/atomAg/s). Computational modeling based on density functional theory (DFT) revealed a more stabilized *COOH intermediate on Ag supported on S-doped carbon and thus a more favorable energetic pathway of CO2-to-CO, consistent with experimental results from the model electrocatalysts studies. Finally, this proof of concept was translated to the synthesis of powder electrocatalyst with 2 wt% Ag supported on S-doped carbon black, demonstrating > 40-fold high CO mass activity than a commercial Ag cathode with steady FECO ~ 96% at 100 mA/cm2 for 50 hours of continuous operation in a gas diffusion electrode (GDE) electrolyzer. For comparison, 2 wt% Ag supported on carbon black without S- doping showed a maximum FECO ~ 70% at 100 mA/cm2. This work demonstrates a successful bottom-up design of CO2R electrocatalysts guided by surface science enabled electrocatalysis.

CO2 conversion↗

DOE OSTI2025

Effects of Temperature Fluctuations on Surface Mobility of Atomic Steps and Oxidation Dynamics in High-Temperature Alloys

In contrast to the traditional perspective that thermal fluctuations are insignificant in surface dynamics, here we report their influence on surface reaction dynamics. Using real-time low-energy electron microscopy imaging of NiAl(100) under both vacuum and O 2 atmospheres, we demonstrate that transient temperature variations substantially alter the direction of atom diffusion between the surface and bulk, leading to markedly different oxidation outcomes. During heating, substantial outward diffusion of atoms from the bulk to the surface results in step growth. Conversely, cooling induces considerable inward diffusion of adatoms, producing a distinct oxide morphology. In both scenarios, initially formed oxide islands impede local atomic step mobility, thereby increasing step length due to mass transfer between the surface and bulk, with atomic steps acting as adatom sinks during heating and sources during cooling. Furthermore, we show that this pinning effect on atomic step mobility can be mitigated by applying persistent temperature fluctuations. As a result, understanding these nuances is vital for accurately predicting and dynamically manipulating the performance of active materials in various chemical processes under transient thermal conditions.

36 MATERIALS SCIENCE↗