Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “modulus”

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 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Large-scale atomistic studies of sliding friction in polycrystalline aluminum interfaces

In this work, we discuss large-scale non-equilibrium molecular dynamics (NEMD) simulations of ductile metal sliding comprising up to 1.8 × 10 9 atoms over time scales of 100 ns. The results of these simulations have identified a variety of physical mechanisms that are important in determining the steady-state frictional force for a wide range of velocities at compressed metal–metal interfaces. These include grain growth and refinement, the evolution of large plastic strains and strain rates, material mixing, and melting. These phenomena can be included in a strain, strain rate, and grain size model that gives good agreement with the NEMD simulations and can be applied to macroscopic continua.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

An investigation of shock-induced phase transition in soda-lime glass

There exists a large body of evidence from experiments and molecular dynamics simulations to suggest the occurrence of phase transitions in soda-lime glass (SLG) and other silica glasses subject to shock compression to pressures above 3 GPa. In light of these findings, the current work investigated the existence of phase transition in SLG using shock and release experiments. The experiments employed symmetric SLG–SLG impact to achieve complete unloading to zero stress after shock compression to stresses in the range of 3–7 GPa. The stress–strain response and the Lagrangian release wave speed behavior of SLG obtained from these experiments are seen to reveal a mismatch between the loading and unloading paths of the pressure–strain curve for the material, which serves as compelling evidence for the occurrence of a shock-induced phase transition in the material at relatively low pressures. Furthermore, the release wave speed vs strain data obtained from experiments were used to construct a methodology for modeling the shock and release behavior of SLG. Lastly, this scheme implemented in numerical simulations was able to capture the release behavior of shock compressed SLG, for which a robust and satisfactory model was previously unavailable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

The effect of polymer stiffness on magnetization reversal of magnetorheological elastomers

Ultrasoft magnetorheological elastomers (MREs) offer convenient real-time magnetic field control of mechanical properties that provides a means to mimic mechanical cues and regulators of cells in vitro. Here, we systematically investigate the effect of polymer stiffness on magnetization reversal of MREs using a combination of magnetometry measurements and computational modeling. Poly-dimethylsiloxane-based MREs with Young’s moduli that range over two orders of magnitude were synthesized using commercial polymers Sylgard™ 527, Sylgard 184, and carbonyl iron powder. The magnetic hysteresis loops of the softer MREs exhibit a characteristic pinched loop shape with almost zero remanence and loop widening at intermediate fields that monotonically decreases with increasing polymer stiffness. A simple two-dipole model that incorporates magneto-mechanical coupling not only confirms that micrometer-scale particle motion along the applied magnetic field direction plays a defining role in the magnetic hysteresis of ultrasoft MREs but also reproduces the observed loop shapes and widening trends for MREs with varying polymer stiffnesses.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Resonant ultrasound spectroscopy for crystalline samples containing initial strain

Resonant ultrasound spectroscopy (RUS) is a mature and robust technique for the nondestructive characterization of the elastic properties of solids capable of providing the elastic constants of anisotropic crystalline solids. The traditional method is based on assuming that the solid is linear elastic and, therefore, obeys a linearized Hookean constitutive relationship (Hooke’s law). In this work, an alternative constitutive law is provided that allows for an initial strain or prestrain to be present stemming from residual stress. Then, the constitutive relationship is integrated into the RUS framework. The model is demonstrated using a realistic prestrain field obtained by simulating shot-peening processing of a polycrystalline Cu sample. The sensitivity of the resonances to the developed prestrain is established and discussed. This work allows researchers to consider the influence of initial strain or residual stress in their samples and the potential influence on accurate elastic constant estimates. The model also supports the potential of RUS for the nondestructive characterization of prestrain in materials.

42 ENGINEERING↗
DOE OSTI2022

An Eulerian crystal plasticity framework for modeling large anisotropic deformations in energetic materials under shocks

Here, this paper demonstrates a novel Eulerian computational framework for modeling anisotropic elastoplastic deformations of organic crystalline energetic materials (EM) under shocks. While Eulerian formulations are advantageous for handling large deformations, constitutive laws in such formulations have been limited to isotropic elastoplastic models, which may not fully capture the shock response of crystalline EM. The present Eulerian framework for high-strain rates, large deformation material dynamics of EM incorporates anisotropic isochoric elasticity via a hypo-elastic constitutive law and visco-plastic single-crystal models. The calculations are validated against atomistic calculations and experimental data and benchmarked against Lagrangian (finite element) crystal plasticity computations for shock-propagation in a monoclinic organic crystal, octahydro-1,3,5,7-tetranitro-1,3,5,7 tetrazocine (β-HMX). The Cauchy stress components and the resolved shear stresses calculated using the present Eulerian approach are shown to be in good agreement with the Lagrangian computations for different crystal orientations. The Eulerian framework is then used for computations of shock-induced inert void collapse in β-HMX to study the effects of crystal orientations on hotspot formation under different loading intensities. The computations show that the hotspot temperature distributions and the collapse profiles are sensitive to the crystal orientations at lower impact velocities (viz., 500 m/s); when the impact velocity is increased to 1000 m/s, the collapse is predominantly hydrodynamic and the role of anisotropy is modest. The present methodology will be useful to simulate energy localization in shocked porous energetic material microstructures and other situations where large deformations of single and polycrystals govern the thermomechanical response.

42 ENGINEERING↗
DOE OSTI2022

Yield strength of CeO 2 measured from static compression in a radial diamond anvil cell

Cerium oxide (ceria, CeO 2 ) is frequently used as a standard in applications such as synchrotron and x-ray free electron lasers for calibrating x-ray wavelengths and offers the potential for understanding the high pressure properties and deformation mechanisms in a wide range of similar face centered cubic (fcc) materials. In this study, the pressure dependence of the strength of ceria was investigated up to 38 GPa using angle dispersive x-ray diffraction in a radial geometry in a diamond anvil cell. In this experiment, the difference in the stress along the axis of compression and perpendicular to the direction of compression can be determined, giving a quantity known as the differential stress. It was found that the differential stress (t), a measure of the lower bound for yield strength, initially increases rapidly from 0.35 ± 0.06 GPa to 2.2 ± 0.4 GPa at pressures of 1.8 and 3.8 GPa, respectively. Above 4 GPa, t increases more slowly to 13.8 ± 2.6 GPa at a pressure of 38 GPa. The changes in the preferred orientation (texture) of CeO 2 with pressure were also measured, allowing for the determination of active deformation mechanisms using an elasto-viscoplastic self-consistent model (EVPSC). It was found that as pressure increased, the [001] direction had a slight preferred orientation along the axis of compression. Our EVPSC model of experimental fiber (cylindrically symmetric) textures and lattice strains were most consistent with dominant slip activity along ${111}$ $\langle$$1\bar10$$\rangle$.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Drag reduction and the Vogel exponent of a flexible beam in transient shear flows

Interactions between a flexible beam and a fluid in a channel are of great relevance to biological hairy surfaces, aquatic vegetation, marine life (e.g., fish gills), and many industrial systems alike. While steady state response of a beam to such flows is fairly well-explored, their behavior in the transient regime is not fully understood. Here, a series of numerical simulations are performed to study the laminar Couette flow of an incompressible viscous fluid past an elastic beam in a two-dimensional channel. The flexible beam is perpendicular to the direction of flow, and its base is fixed to the stationary bottom of the channel. We measure the evolution of the Vogel exponent, drag reduction, and reconfiguration number during the transient and steady-state response of the fluid–structure system for different geometrical and physical properties. Our benchmark shows a good agreement between numerical and experimental observations. Our results show that the system's steady-state response at different bulk-fluid velocities can be reproduced by investigating the shear flow response during the transient regime. We define a new variable that characterizes the evolution of the local velocity profile in the proximity of the free end of the beam and use that to characterize the transient-regime response. The analysis yields insight into the competing effects of elasticity of the beam and non-linear flow response.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Contactless excitation of acoustic resonance in insulating wafers

Contactless excitation and detection of high harmonic acoustic overtones in a thin insulator single crystal are described using radio frequency spectroscopy techniques. Single crystal [001] silicon wafer samples were investigated, one side covered with a Nb thin film, the common starting point for the fabrication of quantum devices. Additionally, the coupling between electromagnetic signals and mechanical oscillation is achieved from the Lorentz force generated by an external magnetic field. This method is suitable for any sample with a metallic surface or covered with a thin metal film. High resolution measurements of the temperature dependence of the sound velocity and elastic constants of silicon are reported and compared with known results.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Exploring the hardness and high-pressure behavior of osmium and ruthenium-doped rhenium diboride solid solutions

Rhenium diboride (ReB 2 ) exhibits high differential strain due to its puckered boron sheets that impede shear deformation. Here, we demonstrate the use of solid solution formation to enhance the Vickers hardness and differential strain of ReB 2 . ReB 2 -structured solid solutions (Re 0.98 Os 0.02 B 2 and Re 0.98 Ru 0.02 B 2 , noted as “ReOsB 2 ” and “ReRuB 2 ”) were synthesized via arc-melting from the pure elements. In-situ high-pressure radial x-ray diffraction was performed in the diamond anvil cell to study the incompressibility and lattice strain of ReOsB 2 and ReRuB 2 up to ~56 GPa. Both solid solutions exhibit higher incompressibility and differential strain than pure ReB 2 . However, while all lattice planes are strengthened by doping osmium (Os) into the ReB 2 structure, only the weakest ReB 2 lattice plane is enhanced with ruthenium (Ru). These results are in agreement with the Vickers hardness measurements of the two systems, where higher hardness was observed in ReOsB 2 . The combination of high-pressure studies with experimentally observed hardness data provides lattice specific information about the strengthening mechanisms behind the intrinsic hardness enhancement of the ReB 2 system.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Nanolamellar phase transition in an additively manufactured eutectic high-entropy alloy under high pressures

Much is unknown about how phase transitions link to micro-/nano-structures in high-entropy systems, especially under extreme pressure and temperature conditions. This work studies the evolution of dual-phase nanolamellar eutectic high-entropy alloy phases of AlCoCrFeNi 2.1 generated by laser powder-bed fusion (L-PBF) for pressures up to 42 GPa. We compare quasi-hydrostatic high pressure synchrotron x-ray diffraction studies on L-PBF printed cylindrical samples up to 5.5 GPa (large-volume Paris–Edinburgh cell) to those carried out on an L-PBF printed foil in a diamond anvil cell where the pressure reached 42 GPa. Our results show that the initially alternating face-centered cubic (FCC) and body-centered cubic (BCC) nanolamellar structure of AlCoCrFeNi 2.1 transformed into single-phase FCC nanolamellae under high pressure with BCC–FCC phase transformation completion at 21 ± 3 GPa. Our results indicate a diffusionless BCC–FCC transformation in this additively manufactured far-from-equilibrium microstructure and demonstrate that the FCC phase is stable up to very high pressures. The measured equation of state for the FCC phase of AlCoCrFeNi 2.1 is presented up to 42 GPa and shows excellent agreement between the data obtained in large-volume press and diamond anvil cell experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Using deep machine learning to interpret proton radiography data from a pulsed power experiment

Deep machine learning is used to analyze a proton radiograph from a tin pulsed power experiment and determine density values for each pixel in the image. Two promising convolutional neural network architectures that have proven to be effective for image analysis in other applications are applied to analyze a proton radiograph and find density values. The process of creating a suitable training dataset is described, involving the Lagrangian hydrodynamic model used for simulations of the experiment, the proton radiography forward model to make synthetic images for training, and the manner in which data augmentation is used to expand the resulting image dataset. It is shown that machine learning not only produces a reasonable density field but is also able to predict features in the density field that are suggested by the proton radiograph but not captured by simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Symmetry and scaling in one-dimensional compressible two-phase flow

Investigations of shock compression of heterogeneous materials often focus on the shock front width and overall profile. The number of experiments required to fully characterize the dynamic response of a material often belie the structure–property relationships governing these aspects of a shock wave. Recent observations measured a pronounced shock-front width on the order of 10 s of ns in particulate composites. We focus on particulate composites with disparate densities and investigate whether the mechanical interactions between the phases are adequate to describe this emergent behavior. The analysis proceeds with a general Mie–Grüneisen equation of state for the matrix material, a general drag force law with general power-law scaling for the particle-matrix coupling of the phases, and a volume fraction-dependent viscosity. Lie group analysis is applied to one-dimensional hydrodynamic flow equations for the self-consistent interaction of particles embedded in a matrix material. The particle phase is characterized by a particle size and volume fraction. The Lie group analysis results in self-similar solutions reflecting the symmetries of the flow. The symmetries lead to well-defined scaling laws, which may be used to characterize the propagation of shock waves in particle composites. An example of the derived scaling laws for shock attenuation and rise time is shown for experimental data on shock-driven tungsten-loaded polymers. A key result of the Lie analysis is that there is a relationship between the exponents characterizing the form of the drag force and the exponent characterizing the shock velocity and its attenuation in a particulate composite. Comparison to recent experiments results in a single exponent that corresponds to a conventional drag force.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Computational synthesis of a new generation of 2D-based perovskite quantum materials

Perovskite-based optoelectronic devices have emerged as a promising energy source due to their potential for scalable production. This study introduces “perovskene,” a novel class of 2D materials derived from the ABC3-like perovskites, synthesized via a data-driven, high-throughput computational strategy. We harness machine learning and multitarget deep neural networks to systematically investigate the structure–property relations, paving the way for targeted material design and optimization in fields such as renewable energy, electronics, and catalysis. The characterization of over 1500 synthesized structures shows that more than 500 structures are stable, revealing properties such as ultra-low work function and large magnetic moment, underscoring the potential for advanced technological applications.

2D materials↗
DOE OSTI2024

Demonstrating resonant ultrasound spectroscopy as a viable technique to characterize thermally conditioned high explosive materials

We present results of resonant ultrasound spectroscopy (RUS) measurements applied to granular high explosive materials at different bulk pressing densities and degree of thermal conditioning. The material chosen in this study is a ubiquitously used explosive material known as pentaerythritol tetranitrate (PETN), which is used commercially in civil and defense applications both as a binderized plastic bonded explosive material and an unbinderized neat material. However, changes in granular PETN bulk elastic properties due to thermal conditioning, which could have implications for better understanding environmental aging-related effects, have not been well studied even though it is believed that elasticity may play an important role in explosive material initiation mechanisms. Furthermore, monitoring elastic property changes in granular explosive pressings has not yet been demonstrated using RUS, which is an appealing non-destructive characterization tool that requires only dry point contact with the explosive material. To this end, we report the first study using RUS to quantify the elastic properties of binderized and neat PETN pressings as well as to quantify changes in elastic properties as a function of both thermal conditioning and bulk pressing density. Elastic stiffness coefficients, sometimes more commonly referred to as elastic constants, calculated from the RUS measurements on the different PETN-based materials show a significant increase for the post-conditioned samples compared to the pre-conditioned samples. This trend of increasing elastic properties with thermal conditioning was consistent for different density pressings, different thermal exposure conditions, and even different neat PETN pressings of differing average crystal sizes and/or specific surface areas.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Thermally frustrated phase transition at high pressure in B 2-ordered FeV

X-ray diffraction measurements of equiatomic B2-ordered FeV were performed in a diamond-anvil cell at room temperature at several pressure points up to 80 GPa that showed the cubic phase to be stable with no indication of structural phase transitions. Density functional theory at 0 K predicts Fermi surface nesting, an electronic topological transition, and a phonon dynamical instability within the experimentally investigated pressure range. Nevertheless, the instability is absent in phonon dispersion curves extracted from ab initio molecular dynamics simulations below the critical volume at temperatures as low as 10 K, indicating that thermal atomic displacements can frustrate the phase transition by renormalizing the phonon dispersion curves. Ferrimagnetism is critical for the stability of the cubic phase at low temperature, but thermal atomic displacements are enough to support the structure at and above the Néel temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Enhancing high-fidelity neural network potentials through low-fidelity sampling

The efficacy of neural network potentials (NNPs) critically depends on the quality of the configurational datasets used for training. Prior research using empirical potentials has shown that well-selected liquid–solid transitional configurations of a metallic system can be translated to other metallic systems. This study demonstrates that such validated configurations can be relabeled using density functional theory (DFT) calculations, thereby enhancing the development of high-fidelity NNPs. Training strategies and sampling approaches are efficiently assessed using empirical potentials and subsequently relabeled via DFT in a highly parallelized fashion for high-fidelity NNP training. Our results reveal that relying solely on energy and force for NNP training is inadequate to prevent overfitting, highlighting the necessity of incorporating stress terms into the loss functions. To optimize training involving force and stress terms, we propose employing transfer learning to fine-tune the weights, ensuring that the potential surface is smooth for these quantities composed of energy derivatives. This approach markedly improves the accuracy of elastic constants derived from simulations in both empirical potential-based NNPs and relabeled DFT-based NNPs. Overall, this study offers significant insights into leveraging empirical potentials to expedite the development of reliable and robust NNPs at the DFT level.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Proton radiation effects in indium oxide using cascade molecular dynamics simulations

Metal oxide (MO) semiconductors, characterized by their wide band gaps and notable charge transport properties, are promising candidates for electronic applications in extreme environments, including near-Earth space. However, atomistic simulations of radiation–matter interactions in MOs remain challenging due to the limitations of existing interatomic potentials, which often fail to capture both the short-range repulsive forces essential for radiation damage modeling and the long-range electrostatic effects governing defect evolution. In this work, we develop a customized interatomic potential tailored for radiation damage simulations in indium oxide (In 2 O 3 ) as a model system, a representative MO material. Our potential integrates the Ziegler-Biersack-Littmark potential to accurately describe short-range interactions with Buckingham and Coulombic potentials to account for long-range forces. We perform molecular dynamics simulations of low-energy proton irradiation using this custom potential. We employ the primary knock-on atom (PKA) cascade method to study atomic displacements and primary defect formation. Simulations were conducted for 1 keV proton irradiation in four randomly chosen directions, and PKA-driven defect analyses at 5, 10, and 15 keV to examine the effects of direction and energy level on damage generation. Our results provide insight into the impact of irradiation direction and energy level on the cascade evolution and defect formation mechanisms.

Atomistic simulations↗
DOE OSTI2025

Raman scattering of rhenium for secondary pressure calibration

With the increasing number of 100 s GPa experiments in the diamond anvil cell (DAC), improved accuracy in secondary pressure calibrations to extreme pressures is essential. The rhenium equation of state has been proposed as a pressure calibrant via x-ray diffraction with potentially broad applications as it is commonly used as a gasket material in DAC experiments. In this work, we conducted Raman spectroscopy experiments on rhenium in the DAC and report the pressure shift of the E2g mode, a refined high-pressure C44 and mode-Grüneisen parameter above 200 GPa. We used flat, beveled, and toroidal diamond anvils under quasi-hydrostatic and non-hydrostatic conditions. By measuring the E2g mode from the culet edge to the center, we analyzed pressure distribution based on culet type and distance from the anvil center. The shift in the E2g mode can be expressed as a function of pressure, and diamond edge measurements appear reliable across all anvil types. Comparing the center and edge pressures reveals anvil cupping, offering insights into predicting or preventing anvil failure during materials properties measurements at extreme conditions.

Diamond anvil cells↗