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.
Reduced-thermal-expansion, reduced-bulk-modulus silicone rubber for use in mold tooling made by incorporating silica powder into silicone rubber. Pressure exerted by thermal expansion reduced even further by allowing air bubbles to remain in silicone rubber instead of deaerating it. Bubbles reduce bulk modulus of material.
With increasing interest in molten salt reactors, there becomes a demand for investigations into thermophysical properties of salt systems. The LiF-NaF-KF (FLiNaK) salt system is a primary candidate for use in these reactors. However, the thermophysical properties of compositions outside the eutectic composition are still largely unknown. In this article, properties of ten unique compositions, including four ternary compositions, are investigated using ab initio molecular dynamics simulations across five temperatures between 900 K and 1300 K. The properties of interest are the density, thermal expansion, bulk modulus, compressibility, heat capacity, and enthalpy of mixing. In general, the results were found to be in good agreement with other literature and experimental results. The density and heat capacity had a tendency to be slightly underpredicted. No conclusions could be drawn about the bulk modulus and compressibility in terms of compositional dependence. The thermal expansion had a negative trend with respect to the LiF concentration and no trends were observed for the NaF or KF concentration. The enthalpy of mixing shows minima for the ternary compositions, with the near-equiatomic composition exhibiting the lowest values. Here this work shows the potential for compositional tailoring in the FLiNaK system to optimize thermophysical properties.
We report the P-V-T equation of state measurements of B 4 C to 50 GPa and approximately 2500 K in laser-heated diamond anvil cells. We obtain an ambient temperature, third-order Birch–Murnaghan fit to the P-V data that yields a bulk modulus K 0 of 221(2) GPa and derivative, (dK/dP) 0 of 3.3(1). These were used in fits with both a Mie–Grüneisen–Debye model and a temperature-dependent, Birch–Murnaghan equation of state that includes thermal pressure estimated by thermal expansion (α) and a temperature-dependent bulk modulus (dK 0 /dT). The ambient pressure thermal expansion coefficient (α 0 + α 1 T), Grüneisen γ(V) = γ 0 (V/V 0 ) q and volume-dependent Debye temperature, were used as input parameters for these fits and found to be sufficient to describe the data in the whole P-T range of this study.
Recently it has been argued based on theoretical calculations and experimental data that there is a universal form for the equation of state of solids. This observation was restricted to the range of temperatures and pressures such that there are no phase transitions. The use of this universal relation to estimate pressure-volume relations (i.e., isotherms) required three input parameters at each fixed temperature. It is shown that for many solids the input data needed to predict high temperature thermodynamical properties can be dramatically reduced. In particular, only four numbers are needed: (1) the zero pressure (P=0) isothermal bulk modulus; (2)it P=0 pressure derivative; (3) the P=0 volume; and (4) the P=0 thermal expansion; all evaluated at a single (reference) temperature. Explicit predictions are made for the high temperature isotherms, the thermal expansion as a function of temperature, and the temperature variation of the isothermal bulk modulus and its pressure derivative. These predictions are tested using experimental data for three representative solids: gold, sodium chloride, and xenon. Good agreement between theory and experiment is found.
Recently it has been argued based on theoretical calculations and experimental data that there is a universal form for the equation of state of solids. This observation was restricted to the range of temperatures and pressures such that there are no phase transitions. The use of this universal relation to estimate pressure-volume relations (i.e., isotherms) required three input parameters at each fixed temperature. It is shown that for many solids the input data needed to predict high temperature thermodynamical properties can be dramatically reduced. In particular, only four numbers are needed: (1) the zero pressure (P = 0) isothermal bulk modulus; (2) its P = 0 pressure derivative; (3) the P = 0 volume; and (4) the P = 0 thermal expansion; all evaluated at a single (reference) temperature. Explicit predictions are made for the high temperature isotherms, the thermal expansion as a function of temperature, and the temperature variation of the isothermal bulk modulus and its pressure derivative. These predictions are tested using experimental data for three representative solids: gold, sodium chloride, and xenon. Good agreement between theory and experiment is found.
The shock wave (Hugoniot) data on single crystal and porous anorthite (CaAl2Si208) to pressures of 120 GPa are presented. These data are inverted to yield high pressure values of the Grueneisen parameter, adiabatic bulk modulus, and coefficient of thermal expansion over a broad range of pressures and temperatures which in turn are used to reduce the raw Hugoniot data and construct an experimentally based, high pressure thermal equation of state for anorthite. The hypothesis that higher order anharmonic contributions to the thermal properties decrease more rapidly upon compression than the lowest order anharmonicities is supported. The properties of anorthite corrected to lower mantle conditions show that although the density of anorthite is comparable to that of the lower most mantle, its bulk modulus is considerably less, hence making enrichment in the mantle implausible except perhaps near its base.
This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called • BCC_opt • FCC_opt • HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string “A_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters “A_B_Câ€, three files are provided: • A_B_C_opt.xyz: The optimized geometries in xyz format • A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. • A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. • A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.
Differential scanning calorimetry data for samples of a 52 year old plutonium alloy with 3.3 at. % Ga that were heated beyond the melting point is analyzed using transition state theory to find activation energies for the δ to ε and ε to liquid phase transitions. A Bayesian statistical method involving a Gaussian process model is used to find mean values and confidence intervals for the activation energies. The activation energy for the δ to ε phase transition increases by 3.3 ± 3.8% per decade, relative to the case when all age related plutonium lattice point defects have been removed through annealing. The corresponding increase in activation energy for the ε to liquid transition is shown to be 7.1 ± 1.8% per decade. It is postulated that the change in activation energy with age for both phase transitions is caused, in part, by the accumulation of the same type of lattice point defects associated with the observed increase in elastic bulk modulus over time.
Bounds for plane strain bulk modulus and shear moduli of multiphase fiber reinforced materials of arbitrary transverse phase geometry
The application of NASTRAN (level 16.0.1) to the static analysis of a sonar dome rubber window (SDRW) was demonstrated. The assessment of the conventional model (neglecting the enclosed fluid) for the stress analysis of the SDRW was made by comparing its results to those based on a sophisticated model (including the enclosed fluid). The fluid was modeled with isoparametric linear hexahedron elements with approximate material properties whose shear modulus was much smaller than its bulk modulus. The effect of the chosen material property for the fluid is discussed.
Multi Principal Elemental Alloys (MPEA) have emerged as promising materials for next-generation nuclear reactors due to their exceptional irradiation resistance. Eight equiatomic MPEA based on the MoNbTi system, comprising of elements with low thermal neutron absorption cross-sections were explored using a combined approach employing empirical parameter estimations, and CALPHAD simulations by which the phases and elemental segregation observed in all the alloys in the as-cast state were predicted. Solution heat treatment at 1500°C transformed five alloys into single-phase matrix materials, enhancing homogeneity and reducing hardness. The densities of the alloys ranged between 6.47 to 7.68 g/cm3, hardness between 472 and 656 VHN, Young’s modulus between 142 GPa to 169 GPa, shear modulus between 54 GPa and 62 GPa, bulk modulus between 117 GPa to 194 GPa and Poisson’s ratio between 0.3 to 0.35. The in-situ high temperature Xray diffraction results, differential scanning calorimetry and dilatometry results up to 1000°C suggested the high temperature phase stability of the MPEA. Subsequent ageing heat treatment at 800 and 1000 oC for 96 hours revealed significant secondary phase precipitation in MoNbTiZr, MoNbTiZrV, and MoNbTiCrA. Oxidation studies at 800 oC for 24 hours in air revealed superior oxidation resistance and cubic rate law dependence in Cr containing MPEA especially in MoNbTiCrAl, while severe mass gain resulting in total disintegration and exfoliation in MoNbTiZr and MoNbTiZrV. This comprehensive study underscores the potential of novel MPEA as promising materials for advanced nuclear reactor applications, shedding light on their microstructural control, mechanical properties, thermal stability, and oxidation resistance.
The high-entropy alloy with composition AlCoCrFeNi 2.1 , additively manufactured with the laser powder-bed fusion technique, has a far-from-equilibrium BCC/FCC eutectic nanolamellar structure. We studied the high-pressure response of this alloy under both static compression and high-strain rate shock compression. The response to static compression using a diamond anvil cell was studied at pressures up to 302 GPa with synchrotron x-ray diffraction at the advanced photon source. The high-pressure FCC-only phase of the EHEA previously observed by Pope et al. [AIP Adv. 13, 035124 (2023)] is found to be stable up to the highest pressure achieved in this study with a volume compression of V/V 0 = 0.587 at ambient temperature. The shock experiments were performed by using GEKKO XII lasers at the Institute of Laser Engineering, Osaka University. The principal Hugoniot equation-of-state of the EHEA was measured up to a pressure of 515 GPa and a compression of V/V 0 = 0.613. Additionally, the thermal equation of state of the EHEA was measured up to 6.2 GPa and 1623 K using a large-volume Paris–Edinburgh cell to obtain the temperature dependence of bulk modulus and thermal expansion coefficients. The melting temperature for EHEA AlCoCrFeNi 2.1 at a pressure of 5.6 GPa was measured to be 1648 ± 25 K. These results can be used to refine stochastic (or special) quasi-random structure (SQS) models for high-pressure high-temperature behavior of high-entropy alloys.
Ternary I–III–VI 2 semiconductors, such as CuInSe 2 , exhibit diverse polymorphs with unique structural characteristics and optoelectronic properties. This study investigates the pressure-induced phase transitions of metastable wurtzite-like CuInSe 2 nanocrystals. Using a combination of synchrotron X-ray diffraction, pair distribution function analysis, and density functional theory calculations, we reveal a transition from cation-ordered wurtzite-like (Pmc2 1 ) to cation-disordered NaCl-like (Fm3̅m) structures at 7.7 GPa. The cation-disordered NaCl-like phase persists upon decompression. Bulk modulus calculations highlight size-dependent deviations from bulk material behavior. These findings deepen our understanding of phase stability in colloidal I–III–VI 2 semiconductor nanocrystals, with implications for tailoring functional materials under extreme conditions.
Protein crystals are so fragile that they often can not be handled by tweezers. Indeed, measurements of the Young modulus, E, of lysozyme crystals resulted in E approx. equals 0.1 - 1 GPa, the lower figures, 0.1 - 0.5 GPa, being obtained from triple point bending of as-grown and not cross-linked crystals sitting in solution. The bending strength was found to be approx.10(exp -2) E. On the other hand, ultrasound speed and Mandelstam-Raman-Brilloin light scattering experiments led to much higher figures, E approx. equals 2.7 GPa. The lower figures for E were found from static or low frequency crystal deformations measurements, while the higher moduli are based on high frequency lattice vibrations, 10(exp 7) - 10(exp 10) 1/s. The physical reason for the about an order of magnitude discrepancy is in different behavior of water filling space between protein molecules. At slow lattice deformation, the not-bound intermolecular water has enough time to flow from the compressed to expanded regions of the deformed crystal. At high deformation frequencies in the ultra- and hypersound waves, the water is confined in the intermolecular space and, on that scale, behaves like a solid, thus contributing to the elastic crystal moduli. In this case, the reciprocal crystal modulus is expected to be an average of the water protein and water compressibilities (reciprocal compressibilities): the bulk modulus for lysozyme is 26 GPa, for water it is 7 GPa. Anisotropy of the crystal moduli comes from intermolecular contacts within the lattice while the high frequency hardness comes from the bulk of protein molecules and water bulk moduli. These conclusions are based on the analysis of liquid flow in porous medium to be presented.
Boron suboxide (B 6 O) is a boron-rich compound derived from the α-rhombohedral boron lattice with extreme hardness and unusual semiconducting properties. In this work, density functional theory (DFT) was used to show that unit cell volume, mechanical strength, band gaps, and thermodynamic stabilities of B 6 O were influenced by the interstitial elements and point defects at the icosahedral sites. While the hexagonal unit cell volume (HUCV) varies with interstitial occupancy, it is the icosahedral defect that weakens the intrinsic bulk modulus of B 6 O. Using the hybrid HSE functional, we confirmed that the perfect B 6 O bulk is a p-type semiconductor with a direct band gap of 2.8 eV. Furthermore, by screening α-boron compounds systematically, we found that a simple octet rule may offer a consistent explanation for the variations in the computed electronic structures. The formation free energies calculated over a wide range of temperatures (0–2500 K) and pressures (0–80 GPa) predict that formations of interstitial defects become favorable only at higher temperatures (ca. 1800 K) in bulk B 6 O lattices. The nudged elastic band (NEB) method was employed to identify the minimum energy pathways for the diffusions of dislocated B and O atoms. The diffusion of icosahedral B atoms has an energy barrier of 0.16 eV. More complex B diffusion paths involving the reorganization of icosahedral boron atoms incur higher barriers (>1 eV). In contrast, the diffusion of interstitial O atoms is facile with a barrier of 0.4 eV. Lastly, successive O insertions into the α-B lattice were performed using DFT to generate a basic understanding of the oxidation process. Furthermore, these calculations provide fundamental atomistic insights into the growth of B 6 O crystals and control of their point defects.
A number of new glass compositions have been prepared with increased emphasis on compositions without beryllia. Glass preparations have been much more broadly based and have included the eutectic glass fields, and the mullite-rare earth glass systems. Of the new glasses, the best non-toxic composition is UARL 472 with a bulk modulus of only 18.20 million psi. A second experimental glass, UARL 417, was chosen for research in making large quantities of fiber in monofilament form. Tests of these UARL 417 epoxy resin samples in comparison to similar composites made with the DuPont organic fiber, PRD-49-1, show that the UARL composites have a compressive strength 41/2 times higher and a specific compressive strength at least 21/2 times greater. Much of the research effort attempted to answer the question of why a given glass should have an impact strength superior to other glasses. No definitive answer to the question was found.
Molecular simulations demonstrate that the enthalpic softening of elastomeric nanocomposites upon heating can arise naturally from a Poisson’s ratio mismatch between elastomer and nanoparticle networks, providing a more parsimonious explanation for this phenomenon than the widely accepted interpretation based on glassy interparticle bridging. Despite a century of use, the mechanism of nanoparticle-driven mechanical reinforcement of elastomers is unresolved. Here, a major hypothesis attributes it to glassy interparticle bridges, supported by an observed inversion of the variation of the modulus E(T) on heating – from entropic stiffening in elastomers to enthalpic softening in nanocomposites. Here, molecular simulations reveal that elastomer enthalpic softening can instead emerge from a competition over the preferred volumes between elastomer and nanoparticulate networks. A theory for this competition accounting for softening of the bulk modulus on heating predicts the simulated E(T) inversion, suggesting that reinforcement is driven by a volume-competition mechanism unique to cocontinuous systems of soft and rigid networks.
A simple, inexpensive, yet accurate method for measuring the dynamic complex modulus of elasticity is described. Using a 'free-free' bar selectively excited in three independent vibrational modes, the shear modulus is obtained by measuring the frequency of the torsional resonant mode and the Young's modulus is determined from measurement of either the longitudinal or flexural mode. The damping properties are obtained by measuring the quality factor (Q) for each mode. The Q is inversely proportional to the loss tangent. The viscoelastic behavior of the sample can be obtained by tracking a particular resonant mode (and thus a particular modulus) using a phase locked loop (PLL) and by changing the temperature of the sample. The change in the damping properties is obtained by measuring the in-phase amplitude of the PLL which is proportional to the Q of the material. The real and imaginary parts or the complex modulus can be obtained continuously as a function of parameters such as temperature, pressure, or humidity. For homogeneous and isotropic samples only two independent moduli are needed in order to characterize the complete set of elastic constants, thus, values can be obtained for the dynamic Poisson's ratio, bulk modulus, Lame constants, etc.