Low frequency bulk modulus of partially molten peridotite, KLB-1
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
This article describes a method to directly measure the bulk modulus of rubber using a piston-cylinder type apparatus. The function of the testing apparatus is to reliably transfer load to the water inside the vessel and allow for the accurate and continuous measurement of volumetric displacement as the contents are compressed. During the test, load is applied to a piston-cylinder containing the rubber test sample and water using a conventional mechanical testing load frame. Load and the travel of the piston are recorded by a load cell and extensometer, respectfully. A brass or stainless steel plug of the same dimensions as the test sample is tested as a control sample to isolate the compliance of the system, including the compressibility of water and compliance of the piston-cylinder. The results from this control test are used to disincorporate the compliance of the system from the sample test results. Load-linear displacement data are mathematically converted to pressure-volumetric displacement data, from which bulk modulus is calculated. The test apparatus and method are validated using the known bulk modulus of water. Finally, the bulk modulus for two similar rubbers is evaluated. The validity of the results and limitations of the device are discussed.
Plutonium metal exhibits an anomalously large softening of its bulk modulus at elevated temperatures that is made all the more extraordinary by the finding that it occurs irrespective of whether the thermal expansion coefficient is positive, negative, or zero—representing an extreme departure from conventional Grüneisen scaling. We show here that the cause of this softening is the compressibility of plutonium’s thermally excited electronic configurations, which has thus far not been considered in thermodynamic models. We show that when compressible electronic configurations are thermally activated, they invariably give rise to a softening of the bulk modulus regardless of the sign of their contribution to the thermal expansion. The electronically driven softening of the bulk modulus is shown to be in good agreement with elastic moduli measurements performed on the gallium-stabilized δ phase of plutonium over a range of temperatures and compositions and is shown to grow rapidly at small concentrations of gallium and at high temperatures, where it becomes extremely sensitive to hydrostatic pressure.
We report an experimental protocol for measuring the frequency dependence of the bulk modulus carried out in a synchrotron x-ray facility based multi-anvil high-pressure apparatus. An oscillating pressure perturbation characterized by x-ray diffraction produces a volume strain measured by imaging. Together, these yield the bulk modulus of the sample. In this work, we report data at 3 mHz as an example of the possibility of providing these data for the frequency range of 1 mHz–100 mHz.
Additive friction stir deposition (AFSD) is a rapidly developing solid-state additive manufacturing (AM) technology that enables the fabrication of customized products on a relatively large scale. The unique forge-liked process offers high quality products with refined equiaxed microstructure and exceptional mechanical performance. However, in AFSD, the spatially distributed mechanical properties within a deposited layer have not been established so far, due to the lack of a suitable large-scale mapping method. In the present work, an ultrasonic dynamic-elasticity mapping method was employed to map the distribution of physical properties in a AFSD 6061 Al. Strong spatial asymmetries in elastographies were mainly due to the residual stresses. The origin of the asymmetries was investigated through thermo-mechanical simulations of the AFSD process. It was observed that majority of the asymmetries were associated with the mechanical residual stresses instead of the symmetric thermal residual stresses during linear deposition of the AFSD process. However, the thermal residual stresses can be manipulated to become asymmetric with a 90° turn in the linear deposition path, due to the difference between the levels of cooling along the edges of the inner and outer corners. Hence, the mechanical and thermo-mechanical residual stresses can be additive or subtractive depending on the direction of the turns. Here, such numerically predicted behaviors were also experimentally observed on the AFSD fabricated samples by ultrasonic elastography. Furthermore, additional preliminary studies on multi-layer AFSD fabrications show the asymmetric mechanical stresses generated during the deposition process on different layers generate constructive or destructive interference which results to spatially variable dynamic bulk modulus depending on the deposition orientations of layers.
Efficient computation and visualization of elastic and mechanical properties are crucial in the selection of materials and the design of new materials. Here, the ElasTool v3.0 toolkit marks a significant advancement in the computational analysis and visualization of elastic and mechanical properties of materials, essential in material selection and design. This enhanced version extends beyond standard calculations like elastic tensor, Young's modulus, bulk modulus, and Poisson's ratio. It introduces capabilities for computing minimum thermal conductivity, linear compressibility, rendering the Christoffel equation, and elastic energy density. Notably, it integrates advanced visualization tools, including compatibility with Plotly and Elate web platforms for interactive web-based property exploration. A key feature of ElasTool v3.0 is the implementation of second-order elastic constants (SOECs) for tubular 2D-based nanostructures and nanotubes. Leveraging high-efficiency strain-matrix sets (OHESS), the toolkit now facilitates efficient computation of elastic constants and mechanical properties at both zero and finite temperatures for 1D, 2D, and 3D dimensions. ElasTool is openly accessible on GitHub: https://github.com/gmp007/elastool.
On the one hand, multi-principal element alloys (MPEAs) have created a paradigm shift in alloy design due to large compositional space, whereas on the other, they have presented enormous computational challenges for theory-based materials design, especially density functional theory (DFT), which is inherently computationally expensive even for traditional dilute alloys. In this paper, we present a machine learning framework, namely PREDICT (PRedict properties from Existing Database In Complex alloys Territory), that opens a pathway to predict elastic constants in large compositional space with little computational expense. The framework only relies on the DFT database of binary alloys and predicts Voigt–Reuss–Hill Young’s modulus, shear modulus, bulk modulus, elastic constants, and Poisson’s ratio in MPEAs. We show that the key descriptors of elastic constants are the A–B bond length and cohesive energy. The framework can predict elastic constants in hypothetical compositions as long as the constituent elements are present in the database, thereby enabling property exploration in multi-compositional systems. We illustrate predictions in a FCC Ni-Cu-Au-Pd-Pt system.
On the one hand, multi-principal element alloys (MPEAs) have created a paradigm shift in alloy design due to large compositional space, whereas on the other, they have presented enormous computational challenges for theory-based materials design, especially density functional theory (DFT), which is inherently computationally expensive even for traditional dilute alloys. In this project, we developed a machine learning framework, namely PREDICT ( PR edict properties from E xisting D atabase I n C omplex alloys T erritory), that opens a pathway to predict elastic constants in large compositional space with little computational expense. The framework only relies on the DFT database of binary alloys and predicts Voigt–Reuss–Hill Young’s modulus, shear modulus, bulk modulus, elastic constants, and Poisson’s ratio in MPEAs. We show that the key descriptors of elastic constants are the A–B bond length and cohesive energy. The framework can predict elastic constants in hypothetical compositions as long as the constituent elements are present in the database, thereby enabling property exploration in multi-compositional systems. We illustrate predictions in a FCC Ni-Cu-Au-Pd-Pt system.
Lab data on permeability, grain bulk modulus, drained bulk modulus, Skempton's pore pressure coefficient, and Biot's effective stress coefficient for Utah Forge Well 58-32 and Well 78-32. Work done as part of project number 5-2615 by the University of Oklahoma.
Explore the source record for details and available documents.
Through an extensive data analysis of powder X-ray diffraction data obtained at pressures up to at least 78 GPa, we report the experimental equations of state for all lanthanide nitrides between NdN and LuN, excluding the radioactive Pm. By fitting the obtained unit cell volumes as a function of pressure with the third order Birch–Murnaghan equation of state, we find that the bulk modulus increases with an increasing lanthanide number from K 0 = 146(12) GPa for NdN to 182(7) GPa in EuN. Hereafter, the bulk modulus reaches a plateau for the rest of the series except for TmN which has a lower bulk modulus. We find that the first derivative of the bulk modulus is around 4 for all compounds except TbN, which displays a significantly different compression behavior. In addition, we find a B1 to B10 pressure-induced phase transition in NdN, SmN, EuN, and GdN at increasingly higher pressures. In fact, we observe that the onset pressure of the phase transition increases linearly with Ln atomic number. From TbN and onwards, we do not observe any sign of a B1 to B10 transition indicating that the transition pressure exceeds the studied pressure range. Therefore, we believe that, for the heavier lanthanides, the linear relationship between the onset pressure and the lanthanide number does not hold and even higher pressures are needed to observe the transition. This coherent study of the series of lanthanide nitrides offers a unique opportunity for benchmark studies of computational methods applied to compounds with 4f electrons.
Microbial conversion of coal to methane is a promising technology during transition of coal-based energy from conventional coal mining to natural gas recovery from coal. Significant research advances have been made towards engineering ideal microbial communities and nutrients for bio-stimulation of coal. However, actual field applications require geomechanical and flow behavior characterization of coal during the bioconversion process as well as gas production over the life of the created “biogenic gas reservoirs”. This work presents the results of an experimental investigation to estimate gas/methane production using bioconversion of coal and analyses of variations in bulk modulus, strain and permeability. Bulk modulus of coal demonstrated time-dependent behavior with continued bioconversion of solid coal. The phenomenon was modeled using mass balancing in a closed environment and the logistic equation based biogenic gas production with time. The modeled gas production and changes in modulus of coal showed excellent agreement with the corresponding experimental results. Next, numerical simulation of biogenic conversion of coal in a constant stress in-situ condition, replicating the Huff ‘n Puff method under field conditions, was carried out. The results showed that bioconversion can lead to decrease in effective stress, increase in permeability and decrease in the modulus of coal with time. Repeated recharge of coal with nutrients and recovering the produced gas, that is, sequential cycles of Huff n’ Puff, showed further increase in permeability and decrease in coal strength, potentially leading to coal failure in-situ, further increasing the permeability, thus enhancing the prospect of field application of the technology.
Montmorillonite clay (MMT) has been widely used in engineering and environmental applications as a landfill barrier and toxic waste repository due to its unique property as an expandable clay mineral that can absorb water easily. This absorption process rendered MMT to be highly exothermic due to electrostatic interactions among molecules and hydrogen bonds between surface atoms. A detailed study of a large supercell model of structural clay enables us to predict long-term nuclear waste storage. Herein, a large solvent MMT model with 4071 atoms is studied using ab initio density functional theory. The DFT calculation and analysis clarify the important issues, such as bond strength, solvation effect, elasticity, and seismic wave velocities. These results are compared to our previous study on crystalline MMT (dry). The solvated MMT has reduced shear modulus (G), bulk modulus (K), and Young’s modulus (E). We observe that the conduction band (CB) in the density of states (DOS) of solvated MMT model has a single, conspicuous peak at –8.5 eV. Moreover, the atom-resolved partial density of states (PDOS) summarizes the roles played by each atom in the DOS. These findings illuminate numerous potential sophisticated applications of MMT clay.
Seismic tomography studies have revealed a high-velocity anomaly at depths between 100 km and 300 km in the Farallon plate. However, the reasons for the high-velocity anomaly continue to be debated. An analysis of the mineral proportions of eclogites exhumed from the Farallon plate shows that the average amount of zoisite in eclogite is ∼16.0 vol%. Therefore, the presence of zoisite eclogite needs to be considered to explain the high-velocity anomaly of the Farallon plate. However, the thermal equation of state and stability of zoisite have not been fully investigated under high pressure−temperature (P-T) conditions. We investigated the high-pressure and high-temperature behavior of natural zoisite utilizing synchrotron single crystal-X-ray diffraction (XRD). The results indicate that zoisite is metastable up to 24.8 GPa and 700 K. We obtained the ambient unit-cell volume V0 = 901.26(3) Å3 by synchrotron single crystal-XRD measurement. We also fitted the pressure-volume-temperature data to a high-temperature Birch-Murnaghan equation of state and obtained the zero-pressure bulk modulus K0 = 134.7(8) GPa, the temperature derivative of the bulk modulus (∂K/∂T)P = −0.011(4) GPa/K, and the thermal expansion coefficients α0V = 1.9(7) × 10−5 K−1 and α1V = 3(2) × 10−8 K−2 when the pressure derivative of bulk modulus (K′0) is fixed at 4. By incorporating the results from previous studies, we calculated the density and bulk sound velocity profiles of zoisite eclogite along the Farallon plate geotherm. Finally, we infer that zoisite could carry water to depths of ∼300 km within cold subducting slabs and that the bulk sound velocity of typical zoisite eclogite with ∼61.0 vol% omphacite, ∼23.0 vol% garnet, and ∼16.0 vol% zoisite could cause the high-velocity anomalies at depths of 100−300 km in the Farallon plate.
We present an integrated experimental and theoretical study of the structural behavior of U 3 Si 5 at highpressure conditions using angle-dispersive synchrotron X-ray diffraction (XRD) in a diamond anvil cell (DAC) and density functional theory (DFT) calculations. On increasing pressure, the ambient hexagonal structure of U 3 Si 5 with space group P6/mmm remains stable up to 16.7 GPa, the maximum pressure tested with DAC. The bulk modulus and the a- and c-axial moduli of U 3 Si 5 were experimentally determined to be 126 ± 4 GPa, 173 ± 8 GPa and 79.7 ± 4.3 GPa, respectively. Thus an anisotropy in the axial compressibility of U 3 Si 5 is observed with its c-axis being more compressible than the a-axis. Our DFT calculation results are in general agreement with the experimental values, including reproducing the compressibility anisotropy. A comparison of the bulk modulus of U 3 Si 5 to those of other U-Si compounds reveals a general trend that the bulk modulus of U-Si decreases with increasing U/(U+Si) ratio.
Here, we developed a PyTorch-based architecture called HydraGNN that implements graph convolutional neural networks (GCNNs) to predict the formation energy and the bulk modulus for models of solid solution alloys for various atomic crystal structures and relaxed volumes. We trained the GCNN surrogate model on a dataset for nickel–niobium (NiNb) generated by the embedded atom model (EAM) empirical interatomic potential for demonstration purposes. The dataset was generated by calculating the formation energy and the bulk modulus as a prototypical elastic property for optimized geometries starting from initial body-centered cubic (BCC), face-centered cubic (FCC), and hexagonal compact packed (HCP) crystal structures, with configurations spanning the possible compositional range for each of the three types of initial crystal structures. Numerical results show that the GCNN model effectively predicts both the formation energy and the bulk modulus as function of the optimized crystal structure, relaxed volume, and configurational entropy of the model structures for solid solution alloys.
Water (H 2 O), in all forms, is an important constituent in planetary bodies, controlling habitability and influencing geological activity. Under conditions found in the interior of many planets, as the pressure increases, the H-bonds in water gradually weaken and are replaced by ionic bonds. Recent experimental measurements of the water equation of state (EOS) showed both a new phase of H-bonded water ice, ice-VIIt, and a relatively low transition pressure just above 30 GPa to ionic bonded ice-X, which has a bulk modulus 2.5 times larger. The higher bulk modulus of ice-X produces larger planets for a given mass, thereby either reducing the atmospheric contribution to the volume of many exoplanets or limiting their water content. We investigate the impact of the new EOS measurements on the planetary mass–radius relation and interior structure for water-rich planets. We find that the change in the planet mass–radius relation caused by the systematic differences between previous and new experimental EOS measurements is comparable to the observational uncertainties in some planet sizes – an issue that will become more important as observations continue to improve.
In situ synchrotron x-ray diffraction was conducted on polycrystalline DyPO 4 to elucidate the details of the pressure-induced transition from the xenotime polymorph to the monazite polymorph. We used three different pressure-transmitting media (neon, a 16:3:1 methanol-ethanol-water mixture, and potassium chloride) to investigate the effect of hydrostaticity on the phase behavior. Specifically, our data clearly show a hydrostatic onset pressure of the xenotime-monazite transition of 9.1 GPa, considerably lower than the 15.3 GPa previously determined by Raman spectroscopy. Based on (quasi)hydrostatic data taken in a neon environment, third-order Birch-Murnaghan equation-of-state fits give a xenotime bulk modulus of 144 GPa and a monazite bulk modulus of 180 GPa (both with pressure derivatives of 4.0). Structural data and axial compressibilities show that DyPO 4 is sensitive to shear and has an anisotropic response to pressure. More highly deviatoric conditions cause the onset of the transition to shift to pressures at least as low as 7.0 GPa. We attribute early transition to shear-induced distortion of the PO4 tetrahedra. Finally, our characterization of the high-pressure behavior of DyPO 4 under variable hydrostaticity is critical for advancing rare earth orthophosphate fiber coating applications in ceramic matrix composites and may inform future tailoring of phase composition for controlled shear and pressure applications.