Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interior structure”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A high internal heat flux and large core in a warm Neptune exoplanet

Interactions between exoplanetary atmospheres and internal properties have long been proposed to be drivers of the inflation mechanisms of gaseous planets and apparent atmospheric chemical disequilibrium conditions. However, transmission spectra of exoplanets have been limited in their ability to observationally confirm these theories owing to the limited wavelength coverage of the Hubble Space Telescope (HST) and inferences of single molecules, mostly H 2 O. Here, in this work, we present the panchromatic transmission spectrum of the approximately 750 K, low-density, Neptune-sized exoplanet WASP-107b using a combination of HST Wide Field Camera 3 (WFC3) and JWST Near-Infrared Camera (NIRCam) and Mid-Infrared Instrument (MIRI). From this spectrum, we detect spectroscopic features resulting from H 2 O (21σ), CH 4 (5σ), CO (7σ), CO 2 (29σ), SO 2 (9σ) and NH 3 (6σ). The presence of these molecules enables constraints on the atmospheric metal enrichment (M/H is 10–18× solar), vertical mixing strength (log 10 K zz = 8.4–9.0 cm 2 s -1 ) and internal temperature (>345 K). The high internal temperature is suggestive of tidally driven inflation4 acting on a Neptune-like internal structure, which can naturally explain the large radius and low density of the planet. These findings suggest that eccentricity-driven tidal heating is a critical process governing atmospheric chemistry and interior-structure inferences for most of the cool (<1,000 K) super-Earth-to-Saturn-mass exoplanet population.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

The formation and structure of iron-dominated planetesimals

Metal-rich asteroids and iron meteorites are considered core remnants of differentiated planetesimals and/or products of oxygen-depleted accretion. Investigating the origins of iron-rich planetesimals could provide key insights into planet formation mechanisms. Using differentiation models, we evaluate the interior structure and composition of representative-sized planetesimals (~200 km diameter), while varying oxygen fugacity and initial bulk meteoritic composition. Under the oxygen-poor conditions that likely existed early in the inner regions of the Solar System and other protoplanetary disks, core fractions remain relatively consistent across a range of bulk compositions (CI, H, EH, and CBa). Some of these cores could incorporate significant amounts of silicon (10–30 weight%) and explain the metal fractions of Fe-rich bodies in the absence of mantle stripping. Conversely, planetesimals forming under more oxidizing conditions, such as beyond snow lines, could exhibit smaller cores, enriched in carbon, sulfur (>1 wt%), and oxides. Sulfur-rich cores, like those formed from EH and H bulk compositions, could remain partly molten, sustain dynamos, and even drive sulfur-rich volcanism. Additionally, bodies with high carbon contents, such as CI compositions, can form graphitic outer layers. These variations highlight the importance of initial formation conditions in shaping planetesimal structures. Future missions, such as NASA’s Psyche mission, offer an opportunity to measure the relative abundances of key elements (Fe, Ni, Si, and S) necessary to distinguish among formation scenarios and structure models for Fe-rich and reduced planetesimals.

meteorites↗
DOE OSTI2025

A small core in Vesta inferred from Dawn’s observations

Vesta’s large-scale interior structure had previously been constrained primarily using the gravity and shape data from the Dawn mission. However, these data alone still allow a wide range of possibilities for the differentiation state of the body. The moment of inertia is arguably the most diagnostic parameter related to the radial density distribution of a planetary body, making it crucial for assessing the body’s state of internal differentiation. Determining the moment of inertia requires additional measurements of the amplitudes of small rotational motions, such as precession and nutation. Here we report an updated estimate of the moment of inertia of Vesta inferred from Dawn’s Doppler tracking via the Deep Space Network and onboard imaging data. The recovered value for Vesta’s normalized polar moment of inertia is $\overline{C}$/MR 2 = 0.4208 ± 0.0047 (where M is the mass of Vesta and R is the reference radius), which is only 6.6% lower than the homogeneous value of 0.4505. This value, combined with the gravity field and global shape, suggests that Vesta’s interior has limited density stratification beneath its howardite–eucrite–diogenite-dominated crust. We propose two possible origin scenarios that are consistent with the observed constraints. In the first scenario, Vesta’s interior did not undergo full differentiation due to late accretion. In the second scenario, Vesta originated as an impact remnant of a larger differentiated body re-accreted with non-chondritic bulk composition produced from a catastrophic impact. Vesta did not experience complete differentiation in either scenario, suggesting that its current state reflects a complex interplay between its accretion timing, thermal evolution, redistribution of 26 Al bearing melt and/or impact processes.

CNEOS 2014-01-08 bolide↗
DOE OSTI2025

Cation Disorder of ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$ in Super‐Earth Mantles

Understanding the mineralogy of exoplanets is essential for unraveling their interior structures, dynamics, and evolution. For large super-Earths, the post-post spinel ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$, one of the major mantle phases, may undergo the order-disorder transition (ODT) at high temperatures. However, the ODT phase boundary of ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$ has not been rigorously constrained. Additionally, fundamental thermodynamic properties of the disordered ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$ remain poorly investigated. Here, we develop a unified machine learning potential (MLP) for ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$ of ab initio accuracy under super-Earth mantle conditions. With the efficient MLP, we extensively calculate the free energy of post-post spinel ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$ via the thermodynamic integration method. The results are used to constrain the ODT phase boundary. Furthermore, we report the P-V-T equation of state and Grüneisen parameters for post-post spinel ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$ across various degrees of disorder. These thermodynamic properties are further applied to update the adiabatic thermal profiles and the mass-radius relation of super-Earths.

36 MATERIALS SCIENCE↗
DOE OSTI2025

X-ray and Neutron Radiography for Quantitative Material Reconstructions

Radiography is a powerful tool to determine the interior structure of objects. X-ray radiography is widely used and provides high-resolution images, though X-rays have limited transmission through materials of high atomic number (Z) and density. In contrast, neutrons can penetrate many materials that are heavily attenuating to X-rays, such as metals, providing contrast in the inner layers of highly attenuating items. Past work has shown the value in using both X-ray and neutron radiography for estimating material thicknesses, though that work was limited to simulated data. Here, we demonstrate quantitative material reconstructions using experimental X-ray and neutron radiography data from lab-based systems, accurately modeling radiography system responses to within a few percent to enable quantitative measures of material thickness. We demonstrate the utility of neutron radiography and X-ray radiography for these quantitative reconstructions and introduce methods for using their complementarity to improve image quality and optimize experimental design.

Gilbert, Andrew J.↗
DOE OSTI2025

High-Resolution Full-Field Structural Microscopy of the Voltage-Induced Filament Formation in VO 2 -Based Neuromorphic Devices

In order to make neuromorphic functions in memristive devices more efficient, information about the structural properties of filament formation at the micro- and mesoscopic scales is necessary. Despite extensive research on VO 2 , a key material due to its filament formation, local operando structural measurements remain challenging and often involve destructive specimen preparation and long rastering times, greatly limiting the scope of experimental studies. Utilizing dark-field X-ray microscopy (DFXM), a fullfield imaging modality, structural signatures of the filament formation process operando are revealed in VO 2 devices. DFXM experiments illustrate that rutile filaments contain isolated monoclinic clusters, indicating structural nonuniformity interior to the filament. The formation of the rutile phase beneath device electrodes was shown to precede filament development, followed by the formation of filament paths guided by nucleation sites within the device. Finally, a medium-term (<30 min) memory mechanism is observed in VO 2 , mediated by sites within the device gap that tend to switch at significantly lower voltages after electrical cycling, a tendency that persists through a brief thermal reset. High spatial resolution, large field-of-view, structure selectivity, and fast signal acquisition of DFXM provided insight into structural features of the filamentary channel and surrounding regions during voltage cycling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Thermally anisotropic composites for thermal management in building environments

An improved system for thermal management is provided. The system includes thermally anisotropic composites coupled with a thermal loop to re-direct, reduce, and shape heat flows through a building envelope, having the potential to (1) significantly reduce envelope-generated heating and cooling loads and (2) provide grid services such as decreasing peak loads and shaping energy use. In one embodiment, the thermal management system includes an anisotropic composite that consists of alternating layers of thermal insulation and thermally conductive materials that are immediately adjacent to each other, including polyisocyanurate foam boards and aluminum sheets. The thermal management system also includes a thermal loop along the long edge or the entire the perimeter of the anisotropic composite, the thermal loop having dynamically controlled or floating temperature that is maintained at lower than an outdoor ambient temperature (for cooling). An interior wall structure is inwardly adjacent to the anisotropic composite.

Shrestha, Som S.↗
DOE OSTI2025

Lightning Induced Interior Fields And Voltage Bounds For Coaxial Topologies

We assemble bounding formulas for the interior fields and pin voltages inside a cylindrical coaxial Faraday cage which has been struck by lightning. Approximate formulas for penetrations through a circumferential door slot with subsequent coupling to the interior center conductor structure. Fields at the opposite open end are estimated and used to drive a capped connector and estimate interior pin voltages. Finally, penetrations through small circular holes and direct diffusion through the barrier are also addressed.

42 ENGINEERING↗
DOE OSTI2025

Convective shells in the interior of Cepheid variable stars: Overshooting models based on hydrodynamic simulations

Context. Because Cepheid variable stars have long been used as a cosmic benchmark for scaling distances in our Galaxy and beyond, the accuracy of stellar evolution models for Cepheids have wide-reaching effects. However, our understanding of the dynamics in the interiors of these physically complex stars is limited. Aims. Our goal is to provide a detailed multi-dimensional picture of hydrodynamic convection and convective boundary mixing in the interior of Cepheids. Methods. Using the Modules for Experiments in Stellar Astrophysics (MESA), we studied the structure of intermediate-mass stars that cross the instability strip. Then, we performed two-dimensional hydrodynamic simulations of six stars with the fully compressible Multidimensional Stellar Implicit Code (MUSIC). Our simulations did not model the radial pulsations but focused on the interior structure of this family of stars. We developed and applied a new statistical analysis to examine convection and convective boundary mixing in the interior of these stellar simulations. Results. Based on a grid of MESA models, we demonstrated that a common structure for intermediate mass Cepheids includes an interior convective shell as well as a thin outer convective envelope. Using the extreme value theory approach to analyze our MUSIC simulation data, we found that overshooting above the convective shell fills the space between these convectively unstable layers. We developed a new statistical analysis that provides a clearer picture of how overshooting fills this layer; it also allowed us to formulate a detailed comparison between overshooting above and below the convective shell. Our analysis effectively decomposes the overshooting layer into two layers: a weak overshooting layer and a strong overshooting layer. Statistically, this is accomplished by decomposing the strongly non-Gaussian probability density function into a mixture of gamma distributions. Using our mixture model, we showed that the ratio of overshooting lengths above and below the convective shell depends directly on the radial extent of the convective shell as well as its depth in the star. We proposed a new form for the diffusion coefficient that addresses the need for overlapping overshooting layers between convective shells. We introduced the idea of a “super-mixing layer” where overshooting from both the convective shell and the convective envelope results in efficient mixing and could be viewed as merging the two adjacent convective zones.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Topological Characterization and Uncertainty Visualization of Atmospheric Rivers

Atmospheric rivers (ARs) are long, narrow regions of water vapor in the Earth's atmosphere that transport heat and moisture from the tropics to the mid-latitudes. ARs are often associated with extreme weather events in North America and contribute significantly to water supply and flood risk. However, characterizing ARs has been a major challenge due to the lack of a universal definition and their structural variations. Existing AR detection tools (ARDTs) produce distinct AR boundaries for the same event, making the risk assessment of ARs a difficult task. Understanding these uncertainties is crucial to improving the predictability of AR impacts, including their landfall areas and associated precipitation, which could cause catastrophic flooding and landslides over the coastal regions. In this work, we develop an uncertainty visualization framework that captures boundary and interior uncertainties, i.e., structural variations, of an ensemble of ARs that arise from a set of ARDTs. We first provide a statistical overview of the AR boundaries using the contour boxplots of Whitaker et al. that highlight the structural variations of AR boundaries based on their nesting relationships. We then introduce the topological skeletons of ARs based on Morse complexes that characterize the interior variation of an ensemble of ARs. We propose an uncertainty visualization of these topological skeletons, inspired by MetroSets of Jacobson et al. that emphasizes the agreements and disagreements across the ensemble members. Through case studies and expert feedback, we demonstrate that the two approaches complement each other, and together they could facilitate an effective comparative analysis process and provide a more confident outlook on an AR's shape, area, and onshore impact.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Constraining the Dynamo Layers in Jupiter and Saturn with Observations and Scaling Laws

The dipole-dominated magnetic fields of Jupiter and Saturn provide evidence for active dynamos operating within their deep interiors, yet the depth of the convecting dynamo layers remains poorly constrained. While magnetic field observations, gravity data, and interior models each provide partial insight, they have not been combined into a single, self-consistent picture of the internal structure. Here, we develop a framework that links observed magnetic field strength with intrinsic heat flux and gravity-constrained interior structure using energy-based dynamo scaling laws. By relating the axial magnetic field strength to the convective power, we infer the radial thickness of the dynamo-generating region for both Jupiter and Saturn. The constants of proportionality in the scaling relations are derived using independent constraints from Earth observations, Jupiter observations, and numerical dynamo simulations. Applied to Jupiter, this framework shows how the inferred dynamo layer thickness is coupled to the outer boundary of the dynamo region. Thinner dynamo layers are predicted when the outer boundary shifts to shallower depths, and no solutions are possible when the outer boundary is less than 73% of Jupiter’s radius. These results constrain plausible geometries for future numerical dynamo simulations. Extending the analysis to Saturn, we find a thick, deep-seated dynamo layer with an outer radius at 42% of the radius to be most plausible. An alternative solution with an inner radius of the dynamo region at 60% of the planetary radius, as suggested by ring seismology models, requires a very thin dynamo layer, occupying only 2%–3% of the total radius.

Geosciences↗
DOE OSTI2024

Structural evolution of liquid silicates under conditions in Super-Earth interiors

Molten silicates at depth are crucial for planetary evolution, yet their local structure and physical properties under extreme conditions remain elusive due to experimental challenges. In this study, we utilize in situ X-ray diffraction (XRD) at the Matter in Extreme Conditions (MEC) end-station of the Linear Coherent Linac Source (LCLS) at SLAC National Accelerator Laboratory to investigate liquid silicates. Using an ultrabright X-ray source and a high-power optical laser, we probed the local atomic arrangement of shock-compressed liquid (Mg,Fe)SiO 3 with varying Fe content, at pressures from 81(9) to 385(40) GPa. We compared these findings to ab initio molecular dynamics simulations under similar conditions. Results indicate continuous densification of theO-O and Mg-Si networks beyond Earth’s interior pressure range, potentially altering melt properties at extreme conditions. This could have significant implications for early planetary evolution, leading to notable differences in differentiation processes between smaller rocky planets, such as Earth and Venus, and super-Earths, which are exoplanets withmasses nearly three times that of Earth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

2D nanoconfinement distorts the solvation structure of hydroxide but not of hydronium

Understanding ion-specific behavior in nanoconfined water is essential for controlling charge transport and selectivity in two-dimensional membranes. Motivated by recent experiments revealing anomalous dielectric and transport behavior of water confined between hBN sheets, we use machine-learning-accelerated first-principles molecular dynamics to investigate the interfacial propensities of hydronium and hydroxide ions under similar confinement. We find that hydronium remains interfacial across all confinement regimes, whereas hydroxide shifts toward the interior as the environment becomes more bulk-like. This contrasting behavior reflects the combined influence of hydrogen bonding, interfacial water layering, and the polarization of hBN, which collectively stabilize hydronium at the surface while making hydroxide slightly more favorable within the structured interior. These findings expose an asymmetry in ion-surface coupling and establish a microscopic origin for hydronium’s enhanced interfacial affinity. The results provide mechanistic insight into ion partitioning in two-dimensional channels and highlight the collective structuring of confined water as a key determinant of interfacial ion thermodynamics.

organic↗
DOE OSTI2026

A mild calcium carbonate treatment improves moisture resistance and mechanical/interfacial properties of jute and hemp fibers

Natural fibers are attractive, sustainable reinforcements, but often exhibit high moisture uptake and weak bonding with polymer matrices. In this study, jute and hemp were treated by immersion in aqueous calcium carbonate (CaCO3) suspensions (10–40?wt.%, 1.5?h, room temperature) and were characterized by scanning electron microscopy (SEM), X-ray diffraction (XRD), moisture sorption testing, single-fiber tensile testing, and bundle pull-out testing. SEM showed granular CaCO3 deposition on fiber surfaces, with particulates bridging surface voids and microcracks. At 40 wt.% CaCO3, equilibrium moisture uptake decreased to ~5% (hemp) and ~3% (jute), while tensile strength increased to ~904 MPa (hemp) and ~960 MPa (jute). Apparent interfacial shear strength from pull-out testing increased up to ~1.35 MPa (hemp) and ~0.56 MPa (jute). These results indicate that CaCO3 deposition is a mild, scalable surface treatment that improves moisture resistance, fiber tensile strength, and interfacial load transfer relevant to semi-structural composite applications, including interior and secondary automotive components.

Calcium carbonate treatment↗
DOE OSTI2021

Nonlinear nanoelectrodynamics of a Weyl metal

Chiral Weyl fermions with linear energy-momentum dispersion in the bulk accompanied by Fermi-arc states on the surfaces prompt a host of enticing optical effects. While new Weyl semimetal materials keep emerging, the available optical probes are limited. In particular, isolating bulk and surface electrodynamics in Weyl conductors remains a challenge. We devised an approach to the problem based on near-field photocurrent imaging at the nanoscale and applied this technique to a prototypical Weyl semimetal TaIrTe 4 . As a first step, we visualized nano-photocurrent patterns in real space and demonstrated their connection to bulk nonlinear conductivity tensors through extensive modeling augmented with density functional theory calculations. Notably, our nanoscale probe gives access to not only the in-plane but also the out-of-plane electric fields so that it is feasible to interrogate all allowed nonlinear tensors including those that remained dormant in conventional far-field optics. Surface- and bulk-related nonlinear contributions are distinguished through their “symmetry fingerprints” in the photocurrent maps. Robust photocurrents also appear at mirror-symmetry breaking edges of TaIrTe 4 single crystals that we assign to nonlinear conductivity tensors forbidden in the bulk. Here, nano-photocurrent spectroscopy at the boundary reveals a strong resonance structure absent in the interior of the sample, providing evidence for elusive surface states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Experiments to validate Thermodynamics and Transport models of Strongly Coupled Dusty Plasma Matter (Final Technical Report for DE-SC0023416)

The goal of this two-year grant is to provide access to the PI to dusty plasma experimental facilities at the DOE-funded Collaborative Research Facility Magnetized Plasma Research Laboratory, Auburn University to become a user of that facility, to obtain experimental data to support another ongoing grant DE-SC0021146 (an Early Career Award to the PI that is focused on modeling of dusty plasma thermodynamics and transport processes), to generate experimental data for funding proposals, and to provide exposure to University of Memphis students to advanced experimental techniques. The following technical accomplishments were made: 1. Development of a novel Bidirectional Electrode Control Arms Assembly (BECAA) for producing perfect 2D grain layers for complex plasma experimentation. BECAA uses movable electrode arms to tilt or move the electrode in a RF discharge from outside the chamber, allowing for the manipulation of grain clouds without needing to change the plasma parameters or gas pressure. This work addresses a longstanding gap in the literature for a method to produce clusters of selectable number of grains and that are perfectly two dimensional as opposed to being only quasi-2D. 2. Experimental investigation of the structural properties of finite-N clusters with N=2 to 50. Individual particle behavior in clusters could vary from grain to grain and this study measured systematically produced clusters for two different grain sizes. Analysis (funded by another grant DE-SC0021146) is currently underway to quantify the differences between grains that are found on the surface vs. the interior of clusters, the shell structure, and the decay of correlations in position, velocity, and kinetic energy. 3. An experimental method to measure the structural entropy of clusters was developed by observing the self-induced structural transitions between various possible arrangements. In a series of heating and cooling cycles, the number of times each possible arrangement was attained was experimentally observed and used to compute the probability of existence of that arrangement, and subsequently the configurational entropy of the cluster. Analysis (funded by another grant DE-SC0021146) is currently underway to produce the entropy of clusters as a function of the number of grains and use the same to compute thermodynamic state variables for 2D complex plasma/grain clusters. 4. A preliminary experimental study of multibody collisions between N grains (N=2 – 10) was conducted. The clusters were produced using the BECAA technique and velocities were imparted to the grains using manipulation laser pulses. Analysis (funded by another grant DE-SC0021146) is currently underway to develop a theoretical framework to describe multibody collisions analogous to classical two-body interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Electron ptychography reveals a ferroelectricity dominated by anion displacements

Sodium niobate, a lead-free ferroic material, hosts delicately balanced, competing order parameters, including ferroelectric states that can be stabilized by epitaxial strain. Here we show that the resulting macroscopic ferroelectricity exhibits an unconventional microscopic structure using multislice electron ptychography. This technique overcomes multiple scattering artefacts limiting conventional electron microscopy, enabling both lateral spatial resolution beyond the diffraction limit and recovery of three-dimensional structural information. These imaging capabilities allow us to separate the ferroelectric interior of the sample from the relaxed surface structure and identify the soft phonon mode and related structural distortions with picometre precision. Unlike conventional ferroelectric perovskites, we find that the polar distortion in this material involves minimal distortions of the cation sublattices and is instead dominated by anion displacements relative to the niobium sublattice. We establish limits on film thickness for interfacial octahedral rotation engineering and directly visualize a random octahedral rotation pattern, arising from the flat dispersion of the associated phonon mode.

ferroelectrics and multiferroics↗
DOE OSTI2024

Different structural behavior of MgSiO 3 and CaSiO 3 glasses at high pressures

Knowledge of the structural behavior of silicate melts and/or glasses at high pressures provides fundamental information for discussing the nature and properties of silicate magmas in the Earth’s interior. The behavior of Si-O structures under high-pressure conditions has been widely studied, while the effect of cation atoms on the high-pressure structural behavior of silicate melts or glasses has not been well investigated. Here, in this study, we investigated the structures of MgSiO 3 and CaSiO 3 glasses up to 5.4 GPa by in situ X-ray pair distribution function measurements to understand the effect of different cations (Mg 2+ and Ca 2+ ) on high-pressure structural behavior of silicate glasses. We found that the structural behavior of MgSiO 3 and CaSiO 3 glasses are different at high pressures. The structure of MgSiO 3 glass changes by shrinking of Si-O-Si angle with increasing pressures, which is consistent with previous studies for SiO 2 and MgSiO 3 glasses. On the other hand, CaSiO 3 glass shows almost no change in Si-Si distance at high pressures, while the intensities of two peaks at ~3.0 and ~3.5 Å change with increasing pressure. The structural change in CaSiO 3 glass at high pressure is interpreted as the change in the fraction of the edge-shared and corner-shared CaO 6 -SiO 4 structures. The different high-pressure structural behavior observed in MgSiO 3 and CaSiO 3 glasses may be the origin of differences in properties, such as viscosity between MgSiO 3 and CaSiO 3 melts at high pressures. This signifies the importance of different structural behaviors due to different cations in investigations of the nature and properties of silicate magmas in Earth’s interior.

36 MATERIALS SCIENCE↗