Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational inference”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Gravity and topography

The paper summarizes the fundamental gravity field constants for Mars and a brief historical review of early determinations and current-day accurate estimates. These include the planetary gravitational constant, global figure, dynamical oblateness, mean density, and rotational period. Topographic results from data acquired from the 1967 opposition to the most recent, 1988, opposition are presented. Both global and selected local topographic variations and features are discussed. The inertia tensor and the nonhydrostatic component of Mars are examined in detail. The dimensionless moment of inertia about the rotational axis is 0.4 for a body of uniform density and 0.37621 if Mars were in hydrostatic equilibrium. By comparing models of both gravity and topography, inferences are made about the degree and depth of compensation in the interior and stresses in the lithosphere.

Esposito, P. B.↗
NASA NTRS1994

A three-dimensional ring current decay model

This work is an extension of a previous ring current decay model. In the previous work, a two-dimensional kinetic model was constructed to study the temporal variations of the equatorially mirroring ring current ions, considering charge exchange and Coulomb drag losses along drift paths in a magnetic dipole field. In this work, particles with arbitrary pitch angle are considered. By bounce averaging the kinetic equation of the phase space density, information along magnetic field lines can be inferred from the equator. The three-dimensional model is used to simulate the recovery phase of a model great magnetic storm, similar to that which occurred in early February 1986. The initial distribution of ring current ions (at the minimum Dst) is extrapolated to all local times from AMPTE/CCE spacecraft observations on the dawn and dusk sides of the inner magnetosphere spanning the L value range L = 2.25 to 6.75. Observations by AMPTE/CCE of ring current distributions over subsequent orbits during the storm recovery phase are compared to model outputs. In general, the calculated ion fluxes are consistent with observations, except for H+ fluxes at tens of keV, which are always over-estimated. A newly-invented visualization idea, designated as a chromogram, is used to display the spatial and energy dependence of the ring current ion differential flux. Important features of storm-time ring current, such as day-night asymmetry during injection and drift hole on the dayside at low energies (less than 10 keV), are manifested in the chromogram representation. The pitch angle distribution is well fit by the function, j(sub o)(1+Ay(exp n)), where y is sine of the equatorial pitch angle. The evolution of the index n is a combined effect of charge exchange loss and particle drift. At low energies (less than 30 keV), both drift dispersion and charge exchange are important in determining n.

Fok, Mei-Ching↗
NASA NTRS1995

Three-dimensional ring current decay model

This work is an extension of a previous ring current decay model. In the previous work, a two-dimensional kinetic model was constructed to study the temporal variations of the equatorially mirroring ring current ions, considering charge exchange and Coulomb drag losses along drift paths in a magnetic dipole field. In this work, particles with arbitrary pitch angle are considered. By bounce averaging the kinetic equation of the phase space density, information along magnetic field lines can be inferred from the equator. The three-dimensional model is used to simulate the recovery phase of a model great magnetic storm, similar to that which occurred in early February 1986. The initial distribution of ring current ions (at the minimum Dst) is extrapolated to all local times from AMPTE/CCE spacecraft observations on the dawnside and duskside of the inner magnetosphere spanning the L value range L = 2.25 to 6.75. Observations by AMPTE/CCE of ring current distributions over subsequent orbits during the storm recovery phase are compared to model outputs. In general, the calculated ion fluxes are consistent with observations, except for H(+) fluxes at tens of keV, which are always overestimated. A newly invented visualization idea, designated as a chromogram, is used to display the spatial and energy dependence of the ring current ion differential flux. Important features of storm time ring current, such as day-night asymmetry during injection and drift hole on the dayside at low energies (less than 10 keV), are manifested in the chromogram representation. The pitch angle distribution is well fit by the function, J(sub o)(1 + Ay(sup n)), where y is sine of the equatorial pitch angle. The evolution of the index n is a combined effect of charge exchange loss and particle drift. At low energies (less than 30 keV), both drift dispersion and charge exchange are important in determining n.

Fok, Mei-Ching↗
DOE OSTI2022

Ca-isotopes as a robust tracer of magmatic differentiation

The large mass difference (~10%) between the two most abundant isotopes of calcium, 40 Ca and 44 Ca, gives Ca great potential in tracking mass-dependent fractionation during magmatic processes. Resolvable Ca-isotope fractionation during fractional crystallization of magma, particularly by feldspar in evolved melts, has been theoretically inferred but not robustly tested in nature. To further explore the effects of magmatic differentiation on Ca-isotope systematics, we studied the late-Permian alkaline igneous suite of the Øyangen Caldera, Oslo Rift, Norway, consisting of volcanic and intrusive units ranging from basanitic to rhyolitic compositions. Major and trace element variations and modeling demonstrate that the main series of samples (N = 21), including basanites, ring-dyke syenites, and central-dome syenites, likely documents a co-genetic and closed-system fractional crystallization sequence. Our data show minimal δ 44/40 Ca variation (< 0.05 ‰) in the intermediate magma and a marked increase in δ 44/40 Ca in the felsic magma of the Øyangen Caldera (from 0.62 ± 0.02 ‰ to 1.15 ± 0.03 ‰ relative to Ca standard, SRM915a). The systematic increase is best explained by equilibrium isotopic fractionation dominated by alkali feldspar in the fractionating mineral assemblage. This is further supported by strong correlations between δ 44/40 Ca, CaO, and Eu/Eu* in the main-series samples. Implementing a Monte Carlo approach, isotopic modeling of the liquid line of descent using Rayleigh fractionation is highly consistent with the observed Ca-isotope evolution. For the first time, we confirm prominent Ca stable isotope fractionation in felsic-stage differentiation of alkaline magma and constrain the isotope fractionation factors of plagioclase and K-feldspar. Integrated with extant estimations on mineral fractionation factors from the literature, our results suggest increasing fractionation effects of rock-forming minerals with decreasing Ca content. Finally, the affirmation of significant Ca-isotope fractionation in alkaline magma by feldspar empowers the application of Ca as a versatile tracer of crustal evolution, allowing further tests in other magmatic conditions across various planetary objects.

58 GEOSCIENCES↗
DOE OSTI2023

Laser Wakefield Accelerator modelling with Variational Neural Networks

A machine learning model was created to predict the electron spectrum generated by a GeV-class laser wakefield accelerator. The model was constructed from variational convolutional neural networks, which mapped the results of secondary laser and plasma diagnostics to the generated electron spectrum. An ensemble of trained networks was used to predict the electron spectrum and to provide an estimation of the uncertainty of that prediction. It is anticipated that this approach will be useful for inferring the electron spectrum prior to undergoing any process that can alter or destroy the beam. In addition, the model provides insight into the scaling of electron beam properties due to stochastic fluctuations in the laser energy and plasma electron density.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Disentangling top-down drivers of mortality underlying diel population dynamics of Prochlorococcus in the North Pacific Subtropical Gyre

Photosynthesis fuels primary production at the base of marine food webs. Yet, in many surface ocean ecosystems, diel-driven primary production is tightly coupled to daily loss. This tight coupling raises the question: which top-down drivers predominate in maintaining persistently stable picocyanobacterial populations over longer time scales? Motivated by high-frequency surface water measurements taken in the North Pacific Subtropical Gyre (NPSG), we developed multitrophic models to investigate bottom-up and top-down mechanisms underlying the balanced control of Prochlorococcus populations. Here, we find that incorporating photosynthetic growth with viral- and predator-induced mortality is sufficient to recapitulate daily oscillations of Prochlorococcus abundances with baseline community abundances. In doing so, we infer that grazers in this environment function as the predominant top-down factor despite high standing viral particle densities. The model-data fits also reveal the ecological relevance of light-dependent viral traits and non-canonical factors to cellular loss. Finally, we leverage sensitivity analyses to demonstrate how variation in life history traits across distinct oceanic contexts, including variation in viral adsorption and grazer clearance rates, can transform the quantitative and even qualitative importance of top-down controls in shaping Prochlorococcus population dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Fiducial-cosmology-dependent systematics for the DESI 2024 full-shape analysis

We assess the impact of the fiducial cosmology choice on cosmological inference from full-shape (FS) fits of the galaxy power spectrum in the DESI 2024 Data Release 1 (DR1). Using a suite of AbacusSummit DR1 mock catalogues based on the Planck 2018 best-fit cosmology, we quantify potential systematic shifts introduced by analysing the data under five secondary cosmologies — featuring variations in matter density, thawing dark energy, higher effective number of neutrino species, reduced clustering amplitude, and the DESI DR1 BAO best-fit w 0 w a CDM cosmology — relative to DESI's baseline Planck 2018 cosmology. We investigate two complementary FS analysis approaches: full-modelling (FM) and ShapeFit (SF), each with distinct sensitivities to the assumed fiducial model. Across all tracers, we find for FM that systematic shifts induced by fiducial cosmology mismatches remain well below the DESI DR1 statistical uncertainties, with maximum deviations of 0.22σ DR1 in ΛCDM scenarios and 0.12σ DR1+SN when including SN Ia mock data in extended w 0 w a CDM fits. For SF, the shifts in the compressed parameters remain below 0.45σ DR1 for all tracers and cosmologies.

dark energy experiments↗
DOE OSTI2022

The half-mass radius of MaNGA galaxies: effect of IMF gradients

Gradients in the stellar populations (SP) of galaxies – e.g. in age, metallicity, stellar initial mass function (IMF) – can result in gradients in the stellar-mass-to-light ratio, M * / L . Such gradients imply that the distribution of the stellar mass and light is different. For old SPs, e.g. in early-type galaxies at z ∼ 0, the M * / L gradients are weak if driven by variations in age and metallicity, but significantly larger if driven by the IMF. A gradient which has larger M * / L in the centre increases the estimated total stellar mass ( M * ) and reduces the scale which contains half this mass ( R e,* ), compared to when the gradient is ignored. For the IMF gradients inferred from fitting MILES simple SP models to the H β, 〈Fe〉, [MgFe], and TiO 2SDSS absorption lines measured in spatially resolved spectra of early-type galaxies in the MaNGA survey, the fractional change in R e,* can be significantly larger than that in M * , especially when the light is more centrally concentrated. The R e,* – M * correlation which results from accounting for IMF gradients is offset to smaller sizes by 0.3 dex compared to when these gradients are ignored. Comparisons with ‘quiescent’ galaxies at higher z must account for evolution in SP gradients (especially age and IMF) and in the light profile before drawing conclusions about how R e,* and M * evolve. Here, the implied merging between higher z and the present is less contrived if R e,* / R e at z ∼ 0 is closer to our IMF-driven gradient calibration than to unity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Sequential Monte Carlo for Cut-Bayesian Posterior Computation

We propose a sequential Monte Carlo (SMC) method to efficiently and accurately compute cut-Bayesian posterior quantities of interest, variations of standard Bayesian approaches constructed primarily to account for model misspecification. We prove finite sample concentration bounds for estimators derived from the proposed method along with a linear tempering extension and apply these results to a realistic setting where a computer model is misspecified. We then illustrate the SMC method for inference in a modular chemical reactor example that includes submodels for reaction kinetics, turbulence, mass transfer, and diffusion. The samples obtained are commensurate with a direct-sampling approach that consists of running multiple Markov chains, with computational efficiency gains using the SMC method. Overall, the SMC method presented yields a novel, rigorous approach to computing with cut-Bayesian posterior distributions.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1998

Evidence for Solar-Cycle Forcing and Secular Variation in the Armagh Observatory Temperature Record

A prominent feature of previous long-term temperature studies has been the appearance of warming since the 1880s, this often being taken as evidence for anthropogenic-induced global warming. In this investigation, the long-term, annual, mean temperature record (1844-1992) of the Armagh Observatory (Armagh, North Ireland), a set of temperature data based on maximum and minimum thermometers that predates the 1880s and correlates well with northern hemispheric and global standards, is examined for evidence of systematic variation, in particular, as related to solar-cycle forcing and secular variation. Indeed, both appear to be embedded within the Armagh data. Removal of these effects, each contributing about 8% to the overall reduction in variance, yields residuals that are randomly distributed. Application of the 10-year moving average to the residuals, furthermore, strongly suggests that the behavior of the residuals is episodic, inferring that (for extended periods of time) temperatures at Armagh sometimes were warmer or cooler (than expected), while at other times they were stable. Comparison of cyclic averages of annual mean temperatures against the lengths of the associated Hale cycles (i.e., the length of two, sequentially numbered, even-odd sunspot cycle pairs) strongly suggests that the temperatures correlate inversely (r = -0.886 at less than 2% level of significance) against the length of the associated Hale cycle. Because sunspot cycle 22 ended in 1996, the present Hale cycle probably will be shorter than average, implying that temperatures at Armagh over this Hale cycle will be warmer (about 9.31 q 0.23 C at the 90% confidence level) than average (= 9.00 C).

Wilson, Robert M.↗
NASA NTRS1998

Water Isotopes in Precipitation: Data/Model Comparison for Present-Day and Past Climates

Variations of HDO and H2O-18 concentrations are observed in precipitation both on a geographical and on a temporal basis. These variations, resulting from successive isotopic fractionation processes at each phase change of water during its atmospheric cycle, are well documented through the IAEA/WMO network and other sources. Isotope concentrations are, in middle and high latitudes, linearly related to the annual mean temperature at the precipitation site. Paleoclimatologists have used this relationship to infer paleotemperatures from isotope paleodata extractable from ice cores, deep groundwater and other such sources. For this application to be valid, however, the spatial relationship must also hold in time at a given location as the location undergoes a series of climatic changes. Progress in water isotope modeling aimed at examining and evaluating this assumption has been recently reviewed with a focus on polar regions and, more specifically, on Greenland. This article was largely based on the results obtained using the isotopic version of the NASA/GISS Atmospheric General Circulation Model (AGCM) fitted with isotope tracer diagnostics. We extend this review in comparing the results of two different isotopic AGCMs (NASA/GISS and ECHAM) and in examining, with a more global perspective, the validity of the above assumption, i.e. the equivalence of the spatial and temporal isotope-temperature relationship. We also examine recent progress made in modeling the relationship between the conditions prevailing in moisture source regions for precipitation and the deuterium-excess of that precipitation.

Jouzel, J.↗
NASA NTRS2000

Feeling Gravity's Pull: Gravity Modeling. The Gravity Field of Mars

Most people take the constant presence of gravitys pull for granted. However, the Earth's gravitational strength actually varies from location to location. This variation occurs because mass, which influences an object's gravitational pull, is not evenly distributed within the planet. Changes in topography, such as glacial movement, an earthquake, or a rise in the ocean level, can subtly affect the gravity field. An accurate measurement of the Earth's gravity field helps us understand the distribution of mass beneath the surface. This insight can assist us in locating petroleum, mineral deposits, ground water, and other valuable substances. Gravity mapping can also help notice or verify changes in sea surface height and other ocean characteristics. Such changes may indicate climate change from polar ice melting and other phenomena. In addition, gravity mapping can indicate how land moves under the surface after earthquakes and other plate tectonic processes. Finally, changes in the Earth's gravity field might indicate a shift in water distribution that could affect agriculture, water supplies for population centers, and long-term weather prediction. Scientists can map out the Earth's gravity field by watching satellite orbits. When a satellite shifts in vertical position, it might be passing over an area where gravity changes in strength. Gravity is only one factor that may shape a satellite's orbital path. To derive a gravity measurement from satellite movement, scientists must remove other factors that might affect a satellite's position: 1. Drag from atmospheric friction. 2. Pressure from solar radiation as it heads toward Earth and. as it is reflected off the surface of the Earth 3. Gravitational pull from the Sun, the Moon, and other planets in the Solar System. 4. The effect of tides. 5. Relativistic effects. Scientists must also correct for the satellite tracking process. For example, the tracking signal must be corrected for refraction through the atmosphere of the Earth. Supercomputers can calculate the effect of gravity for specific locations in space following a mathematical process known as spherical harmonics, which quantifies the gravity field of a planetary body. The process is based on Laplace's fundamental differential equation of gravity. The accuracy of a spherical harmonic solution is rated by its degree and order. Minute variations in gravity are measured against the geoid, a surface of constant gravity acceleration at mean sea level. The geoid reference gravity model strength includes the central body gravitational attraction (9.8 m/sq s) and a geopotential variation in latitude partially caused by the rotation of the Earth. The rotational effect modifies the shape of the geoid to be more like an ellipsoid, rather than a perfect, circle. Variations of gravity strength from the ellipsoidal reference model are measured in units called milli-Galileos (mGals). One mGal equals 10(exp -5) m/sq s. Research projects have also measured the gravity fields of other planetary bodies, as noted in the user profile that follows. From this information, we may make inferences about our own planet's internal structure and evolution. Moreover, mapping the gravity fields of other planets can help scientists plot the most fuel-efficient course for spacecraft expeditions to those planets.

Lemoine, Frank↗
NASA NTRS2013

Intercomparison of Satellite Dust Retrieval Products over the West African Sahara During the Fennec Campaign in June 2011

Dust retrievals over the Sahara Desert during June 2011 from the IASI, MISR, MODIS, and SEVIRI satellite instruments are compared against each other in order to understand the strengths and weaknesses of each retrieval approach. Particular attention is paid to the effects of meteorological conditions, land surface properties, and the magnitude of the dust loading. The period of study corresponds to the time of the first Fennec intensive measurement campaign, which provides new ground-based and aircraft measurements of the dust characteristics and loading. Validation using ground-based AERONET sunphotometer data indicate that of the satellite instruments, SEVIRI is most able to retrieve dust during optically thick dust events, whereas IASI and MODIS perform better at low dust loadings. This may significantly affect observations of dust emission and the mean dust climatology. MISR and MODIS are least sensitive to variations in meteorological conditions, while SEVIRI tends to overestimate the aerosol optical depth (AOD) under moist conditions (with a bias against AERONET of 0.31), especially at low dust loadings where the AOD<1. Further comparisons are made with airborne LIDAR measurements taken during the Fennec campaign, which provide further evidence for the inferences made from the AERONET comparisons. The effect of surface properties on the retrievals is also investigated. Over elevated surfaces IASI retrieves AODs which are most consistent with AERONET observations, while the AODs retrieved by MODIS tend to be biased low. In contrast, over the least emissive surfaces IASI significantly underestimates the AOD (with a bias of -0.41), while MISR and SEVIRI show closest agreement.

Fennec Campaign↗
NASA NTRS1980

Spectrophotometry of Comet West

Postperihelion observations of Comet West (1975n = 1976 VI) have been made with a Fourier transform spectrometer at heliocentric distances from 0.57 to 1.68 AU. Measurements were made of the emission bands of C2, CN, C3, CH, and NH2, as well as the emission lines of Na D and forbidden (O I), and the flux in the continuum in nine different bandpasses. Several ratios of the band strengths of CN have been used to determine the two free parameters in the fluorescence equilibrium model of CN of Danks and Arpigny (1973). From the values of the parameters it is inferred that the vibrational transition probability for the ground electronic state is between 0.025 and 0.075 per sec and that the ratio of oscillator strengths between the (0-0) bands of the violet and red systems is between 25 and 30. When corrected for field-of-view effects, NH2 shows no systematic variation in abundance relative to C2 while CH shows a small increase. The cometary continuum is found to be slightly redder than the solar continuum, consistent with results for other bright, dusty comets. The equivalent width of the Delta u = 0 sequence of C2 shows a marked decrease at r(H) = 1.2 AU.

Ahearn, M. F.↗
NASA NTRS1984

Mass transfer from a circular cylinder: Effects of flow unsteadiness and slight nonuniformities

Experiments were performed to determine the effect of periodic variations in the angle of the flow incident to a turbine blade on its leading edge heat load. To model this situation, measurements were made on a circular cylinder oscillating rotationally in a uniform steady flow. A naphthalene mass transfer technique was developed and used in the experiments and heat transfer rates are inferred from the results. The investigation consisted of two parts. In the first, a stationary cylinder was used and the transfer rate was measured for Re = 75,000 to 110,000 and turbulence levels from .34 percent to 4.9 percent. Comparisons with both theory and the results of others demonstrate that the accuracy and repeatability of the developed mass transfer technique is about + or - 2 percent, a large improvement over similar methods. In the second part identical flow conditions were used but the cylinder was oscillated. A Strouhal number range from .0071 to .1406 was covered. Comparisons of the unsteady and steady results indicate that the magnitude of the effect of oscillation is small and dependent on the incident turbulence conditions.

Marziale, M. L.↗
NASA NTRS1991

Magellan - Electrical and physical properties of Venus' surface

Magellan probes Venus'surface by 12.6-cm-wavelength vertical and oblique radar scattering and measures microwave thermal emission. Emissivity and root-mean-square slope maps between 330 deg and 30 deg E and 90 deg N and 80 deg S are dissimilar, although some local features are exceptions. Inferred surface emissivities typically are 0.85, but vary from 0.35 at Maxwell to 0.95 northeast of Gula Mons and other locations. Lowest emissivities appear in topographically high areas; this relation suggests that a phase change or differences in chemical weathering occur at about 6055-kilometer radius. Initial results indicate that there are significant variations in the surface scattering function.

Tyler, G. Leonard↗
NASA NTRSDate not supplied

Remote Sensing of Clouds, Water Vapor, and Near-Surface Winds Under Various Aerosol Conditions During CAMP2Ex

The Cloud, Aerosol and Monsoon Processes Philippines Experiment (CAMP2Ex) field campaign took place from August – October 2019. CAMP2Ex science objectives include understanding aerosol-cloud interaction under various Meteorological conditions, and how these Meteorological phenomena affect aerosol distribution in the atmosphere. Numerous instruments were flown on the National Aeronautics and Space Administration (NASA) P-3 aircraft during 19 CAMP2Ex science flights, including the Advanced Microwave Precipitation Radiometer (AMPR) the Airborne Third Generation Precipitation Radar (APR-3), and the High Spectral Resolution Lidar 2 (HSRL2). Geophysical retrieval equations for cloud liquid water, columnar water vapor, and near-surface wind speed using AMPR brightness temperatures have been derived and tested using CAMP2Ex data. Water vapor and wind speed retrievals have been validated using in situ observations, such as dropsondes and shipborne instruments. This presentation will demonstrate these AMPR-based retrievals during CAMP2Ex, including discussions of retrieval uncertainties. Specific CAMP2Ex flights will be presented in greater detail, including inferences of other Meteorological phenomena in clear air and around precipitating and non-precipitating clouds, such as background wind speeds and potential identification of gust fronts. These analyses will be supplemented with APR-3 reflectivity and Doppler velocity data to infer microphysical processes occurring within and around selected clouds of interest. In cases where multiple remote-sensing measurements were obtained for the same cloud, a time-evolution analysis of these features may be performed. Throughout these analyses, HSRL2 data will be used to quantify aerosol properties beneath the P-3, and variations in the AMPR/APR-3 observations with respect to various aerosol conditions will be discussed.

Corey Amiot↗
DOE OSTI2023

MHz free electron laser x-ray diffraction and modeling of pulsed laser heated diamond anvil cell

A new diamond anvil cell experimental approach has been implemented at the European x-ray Free Electron Laser, combining pulsed laser heating with MHz x-ray diffraction. In this report we use this setup to determine liquidus temperatures under extreme conditions, based on the determination of time-resolved crystallization. The focus is on a Fe-Si-O ternary system, relevant for planetary cores. This time-resolved diagnostic is complemented by a finite-element model, reproducing temporal temperature profiles measured experimentally using streaked optical pyrometry. This model calculates the temperature and strain fields by including (i) pressure and temperature dependencies of material properties, and (ii) the heat-induced thermal stress, including feedback effect on material parameter variations. Making our model more realistic, these improvements are critical as they give 7000 K temperature differences compared to previous models. Laser intensities are determined by seeking minimal deviation between measured and modeled temperatures. Combining models and streak optical pyrometry data extends temperature determination below detection limit. The presented approach can be used to infer the liquidus temperature by the appearance of SiO 2 diffraction spots. In addition, temperatures obtained by the model agree with crystallization temperatures reported for Fe–Si alloys. Our model reproduces the planetary relevant experimental conditions, providing temperature, pressure, and volume conditions. Those predictions are then used to determine liquidus temperatures at experimental timescales where chemical migration is limited. This synergy of novel time-resolved experiments and finite-element modeling pushes further the interpretation capabilities in diamond anvil cell experiments.

36 MATERIALS SCIENCE↗