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.
We place metamaterials in the context of underpinning physical phenomena, including negative refraction, bandgaps, wave focusing, and negative Poisson’s ratio. The designs, mechanisms, governing equations, and effective parameters are discussed.
Catalysts with high carrier mobility, high activity, and an active basal plane have been highly sought for the hydrogen evolution reaction (HER). However, combining these advantages into one single material is a grand challenge. Herein, using first principles computations, we predicted that a two-dimensional (2D) Dirac nodal loop semimetal, namely the RuB 4 monolayer, is promising as a superior catalyst for the HER. Our systematic computations showed that the single layer RuB 4 is thermodynamically, dynamically, mechanically, and thermally stable and presents multiple Ru and B sites for the HER on the basal plane. The estimated Gibbs free energy for hydrogen adsorption at a Ru site is approaching zero (–8.8 meV), suggesting its excellent HER performance. The RuB 4 monolayer is a Dirac nodal loop semimetal with high Fermi velocities, which can accelerate charge transfer between catalysts and reaction intermediates. The RuB 4 monolayer is an auxetic material with an out-of-plane negative Poisson's ratio, implying its novel mechanical properties. Finally, this work provides an example of using a Dirac nodal loop semimetal for high-performance HER catalysts, which is a promising alternative to the known catalysts with trivial metallic properties.
A nanoscale metastructure with tunable auxetic properties is constructed with DNA origami. This DNA nanostar structure can reconfigure in two orientations that can lead to distinct negative Poisson's ratios.
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.
Compressional (P) and shear (S) wave velocities of tungsten single crystals along the [100] and [110] directions were measured using ultrasonic interferometry at room temperature up to 11.3 GPa. Least-squares fitting of $V$$^{[100]}_{P}$, $V$$^{[100]}_{S}$, $V$$^{[110]}_{P}$, and pressure to finite strain (FS) equations yields the elastic constants: C 11 = 523.5(5) GPa, C 12 = 205.1(32) GPa, and C 44 = 160.8(4) GPa, along with their respective pressure derivatives: $C$$^{′}_{11}$ = 6.20(2), $C$$^{′}_{12}$ = 3.35(1), and $C$$^{′}_{44}$ = 1.65(6). Using the Voigt–Reuss–Hill approximation, the elastic moduli were derived as K S0 = 311.2(22) GPa and G 0 = 160.2(7)GPa, along with their respective pressure derivatives: $K$$^{′}_{S0}$ = 4.30(4) and $G$$^{′}_{0}$ = 1.56(1). The Debye temperature of tungsten was determined to be 380.7(8) K, showing good agreement with previous calorimetric measurements. The elastic anisotropy increases slightly from 1.01(1) at ambient pressure to 1.03(1) at 11.3 GPa, suggesting that tungsten remains nearly isotropic under compression. Poisson's ratio slightly increased from 0.281(3) to 0.288(3) with pressure. Additionally, Pugh's ratio decreased from 0.512(4) to 0.494(4), while Pettifor's ratio increased from 0.148(10) to 0.175(10) as pressure reached 11.3 GPa. These results suggest that tungsten is weakly ductile compared to the critical thresholds for ductile behavior (i.e., Pugh's ratio <0.6 and Pettifor's ratio >0, indicating ductile behavior), but that its ductility increases marginally under pressure.
We modify the traditional linear-elastic field of an edge dislocation. The modifications stem from symmetry and energy requirements imposed in terms of embedded (deformed) coordinates. These requirements are satisfied if a line force is added to the dislocation field. The modification to the stress field is expressed by coefficients that are a function only of Poisson’s ratio. Qualitatively, the field of the dislocation is increased in the glide direction and decreased in the climb direction. There are small changes to the field of a screw dislocation that do not entail the addition of a line force. Effects of nonlinearity, anisotropic elasticity, and image forces are briefly considered.
Isostructural transitions have been reported to occur in several pure elements; however, most of these have been disproven. Zr is one of those elements which remains controversial with two reports of an isostructural bcc-to-bcc transition at ~56–60 GPa at an ambient temperature and a more recent report proposing a second bcc-to-bcc transition at 110 GPa. Here, we report a detailed experimental and theoretical study of zirconium in the pressure region of the proposed isostructural transition(s). We conducted three room temperature static compression experiments using Ne as a pressure transmitting medium up to 80 GPa. We see no evidence for an isostructural transition finding a smooth compression curve for β-Zr. We determined the theoretical volume, elastic parameters, and Poisson’s ratio for β-Zr up to 717 GPa and no obvious anomalies are identified. Furthermore, our new results join the growing body of evidence that does not find evidence of isostructural transitions in β-Zr.
Bi₂Sr₂CaCu₂O8−x (Bi-2212) multi-filament round wire is a high-temperature superconductor (HTS) capable of carrying high transport currents, which makes it suitable for high-field magnet applications. However, its weak Ag–Mg sheath leaves it vulnerable to mechanical stress, posing challenges for high-field magnet design. To better understand and improve mechanical stress management in Bi-2212 winding packs, we conducted an experimental study evaluating the axial stress–strain behavior of five winding pack configurations with varying insulation materials, reinforcement strategies, and construction quality. Using uniaxial tensile testing at 77 K, we measured Young’s modulus and Poisson’s ratio for each composition. Our results show that pure alumina braid insulation and co-wind reinforcements significantly enhance stiffness compared to aluminosilicate braids, with more than 2.5 times increased winding pack Young’s modulus. Rule of mixtures analysis further quantified the contribution of non-wire composite components to overall stiffness. These findings highlight the critical role of insulation material selection and reinforcement design in optimizing Bi-2212 coil performance under stress, providing a foundation for improved mechanical models and more reliable high-field HTS magnet designs.
SUMMARY Analysis of tectonic and earthquake-cycle associated deformation of the crust can provide valuable insights into the underlying deformation processes including fault slip. How those processes are expressed at the surface depends on the lateral and depth variations of rock properties. The effect of such variations is often tested by forward models based on a priori geological or geophysical information. Here, we first develop a novel technique based on an open-source finite-element computational framework to invert geodetic constraints directly for heterogeneous media properties. We focus on the elastic, coseismic problem and seek to constrain variations in shear modulus and Poisson’s ratio, proxies for the effects of lithology and/or temperature and porous flow, respectively. The corresponding nonlinear inversion is implemented using adjoint-based optimization that efficiently reduces the cost function that includes the misfit between the calculated and observed displacements and a penalty term. We then extend our theoretical and numerical framework to simultaneously infer both heterogeneous Earth’s structure and fault slip from surface deformation. Based on a range of 2-D synthetic cases, we find that both model parameters can be satisfactorily estimated for the megathrust setting-inspired test problems considered. Within limits, this is the case even in the presence of noise and if the fault geometry is not perfectly known. Our method lays the foundation for a future reassessment of the information contained in increasingly data-rich settings, for example, geodetic GNSS constraints for large earthquakes such as the 2011 Tohoku-oki M9 event, or distributed deformation along plate boundaries as constrained from InSAR.
Founded on a novel analytical formulation that led to a rigorous yet compact path-integral representation of the time-domain elastodynamic half-space Green's function, a unified analysis of the possible occurrence of different singular wavefront behaviour in the response under arbitrary impulsive internal or surface point loads at arbitrary source-receiver locations is presented. With the decomposition of the general solution into distinct initiating and reflected wave group integrals that share a common factored format and simple contour definitions, the mathematical framework is shown to allow a straightforward identification of the specific conditions and the particular wave groups that are responsible for the singular wavefront phenomena without resorting to advanced analytic function theories or asymptotic methods. Analytic characterizations of the nature, strength and direction of all intrinsic singular wavefront behaviours of the three-dimensional Green's function in three canonical cases of source-receiver configurations are given in a dual integral-closed form format to facilitate their theoretical understanding as well as computational applications. Graphical illustrations of their variation with the source-receiver configuration and the medium's Poisson's ratio together with relevant comparison and clarifications of some classical treatments are included.
The high-temperature and high-pressure response of magnesium has been investigated through shock-release experiments performed up to shock melting. The longitudinal and bulk sound speeds of Mg are reported along the Hugoniot from 25 to 56 GPa and used to determine the elastic properties and Grüneisen parameter. The hexagonal close-packed (hcp)–body-centered cubic (bcc) phase transition is marked by a reduction in the determined shear wave speed. Thermal softening is observed to begin between 40 and 44 GPa, with incipient melt at 55.5 GPa, in close agreement with previous diffraction measurements under laser shock loading. Examination of the release profiles showed that two different responses were observed, depending on the peak stress. When plotted in a stress-energy phase diagram, the two responses are observed to form separate lines that intersect the static hcp-bcc phase transition and incipient melt. The results indicate that the hcp-bcc transition occurs on the Hugoniot at 28.4 GPa and place the hcp-bcc-liquid triple point ~20 GPa. In conclusion, a multiphase equation of state is developed which places the melt boundary below the previously reported static measurements.
A bonded particle model is used to explore how variations in the material properties of brittle, isotropic solids affect critical behavior in fragmentation. To control material properties, a model is proposed which includes breakable two- and three-body particle interactions to calibrate elastic moduli and mode I and mode II fracture toughnesses. In the quasistatic limit, fragmentation leads to a power-law distribution of grain sizes which is truncated at a maximum grain mass that grows as a nontrivial power of system size. In the high-rate limit, truncation occurs at a mass that decreases as a power of increasing rate. A scaling description is used to characterize this behavior by collapsing the mean-square grain mass across rates and system sizes. Consistent scaling persists across all material properties studied, although there are differences in the evolution of grain size distributions with strain as the initial number of grains at fracture and their subsequent rate of production depend on Poisson's ratio. Importantly, this evolving granular structure is found to induce a unique rheology where the ratio of the shear stress to pressure, an internal friction coefficient, decays approximately as the logarithm of increasing strain rate. The stress ratio also decreases at all rates with increasing strain as fragmentation progresses and depends on elastic properties of the solid.
Ultrathin films and multilayers, with controlled thickness down to single atomic layers, are critical for advanced technologies ranging from nanoelectronics to spintronics to quantum devices. However, for thicknesses less than 10 nm, surfaces and dopants contribute significantly to the film properties, which can differ dramatically from that of bulk materials. Therefore, for amorphous films being developed as low dielectric constant interfaces for nanoelectronics, the presence of surfaces or dopants can soften films and degrade their mechanical performance. Here we use coherent short-wavelength light to fully and nondestructively characterize the mechanical properties of individual films as thin as 5 nm within a bilayer. In general, we find that the mechanical properties depend both on the amount of doping and the presence of surfaces. In very thin (5-nm) silicon carbide bilayers with low hydrogen doping, surface effects induce a substantial softening - by almost an order of magnitude - compared with the same doping in thicker (46-nm) bilayers. These findings are important for informed design of ultrathin films for a host of nano- and quantum technologies, and for improving the switching speed and efficiency of next-generation electronics.
The interplay between metallic glasses (MGs) mechanical properties, fracture energy (G), and glass-forming ability (GFA) and their dependence on alloy composition remain poorly understood. Here, we perform molecular dynamics simulations to investigate the intrinsic composition-dependence of G in Cu x Zr 100-x MGs (x = 20, 30, 40, 50, 64). The results indicate that the value of G increases with Cu content. In addition, it is revealed that MGs with higher G values display higher Poisson’s ratio (ν) and GFA, suggesting a close correlation between fracture toughness, mechanical properties, and GFA. This correlation between G, ν and GFA can be understood based on the fragility (m) of supercooled liquids, which is directly related to the structural heterogeneity in MGs. Larger m values are related to dynamic slowdown and supercooled liquid stabilization, which enhance GFA and the formation of pronounced structural heterogeneity, comprised of loosely packed regions that favor β-relaxation and the activation of shear transformation zones. Those concurrently promote the expansion of the plastic zone at the crack tip, enhancing the observed value of G. Furthermore, these simulation results shed light on the intrinsic relationship between fracture toughness, mechanical properties, and alloy composition in MGs.
Disordered solids often change their elastic response as they slowly age. Using experiments and simulations, we study how aging disordered planar networks under an applied stress affects their nonlinear elastic response. We are able to modify dramatically the elastic properties of our systems in the nonlinear regime. Using simulations, we study two models for the microscopic evolution of properties of such a material; the first considers changes in the material strength, while the second considers distortions in the microscopic geometry. Both models capture different aspects of the experiments including the encoding of memories of the aging history of the system and the dramatic effects on the material’s nonlinear elastic properties. Our results demonstrate how aging can be used to create complex elastic behavior in the nonlinear regime.
Understanding epitaxial strain relaxation is one of the key challenges in functional thin films with strong structure–property relations. Herein, we employ an emerging data analytics approach to quantitatively evaluate the underlying relationships between critical thickness ($h_c$) of strain relaxation and various physical and chemical features, despite the sporadic experimental data points available. First, we have collected and refined the reported $h_c$ of the perovskite oxide thin film/substrate system to construct a consistent sub-dataset which captures a common trend among the varying experimental details. Then, we employ correlation analyses and feature engineering to find the most relevant feature set which includes Poisson's ratio and lattice mismatch. With the insight offered by correlation analyses and feature engineering, machine learning (ML) models have been trained to deduce a decent accuracy, which has been further validated experimentally. In this work, the demonstrated framework is expected to be efficiently extended to the other classes of thin films in understanding $h_c$.
With the mining depth of coal resources increasing, the thickness of traditional lining for deep mines becomes large. The bearing capacity of the outer lining cannot be fully utilized, so a new type of functionally graded lining structure (FGL) with radial Young’s modulus varying in gradient is proposed. In this study, through theoretical analysis and numerical simulation, the mechanical properties of the functionally graded lining were studied, including the characteristics of the elastic working state and the ultimate bearing state. The influences of the structural parameters of the functionally graded lining, including the inner radius, the thickness, Young’s modulus, the compressive strength of the concrete, Poisson’s ratio, and the number of layers on the mechanical properties were analyzed. The calculation formula of the ultimate bearing capacity of the lining and the calculation formula of the maximum tangential strain at the time of lining failure were put forward. The accuracy of the formula was verified by comparing with the numerical simulation results. The research results provide the basis for the construction of the design theory of a functionally graded lining structure and have great engineering significance for the construction of the kilometer-deep mines in the future.
Fractionating whole-plant corn (WPC) in a single-pass harvesting system requires studies on the WPC-to-equipment interaction for improved property control, as well as mechanical and air-driven separation processes compared to the traditional multi-pass grain and stover harvesting system. The discrete element method (DEM) technique has the potential to simulate WPC mechanical fractionation and support simulation-based design of WPC separation processes. In this study, methods to develop DEM particle models of WPC (kernel, cob, stalk, and husk) and their material properties for simulating mass fractionation using the ASABE standard mechanical shaker were proposed. Measurement was done on the axial dimensions (major, intermediate, and minor) and mass of each WPC type (mean sample size is 56), sampled from single-pass harvesting. Applying gaussian multivariate regression and bootstrapping re-sampling techniques, a DEM particle approximate to each WPC was developed. Sensitivity analysis of the DEM Young‘s modulus, Poisson‘s ratio, and interaction parameters of coefficient of restitution, coefficient of rolling friction, and coefficient of static friction on mass fraction was performed after 156 ASABE sieve-shaking DEM simulation runs, generated using Latin Hypercube Design (LHD) design of experiment (DOE) from 19 DEM material parameters. DEM simulation using Hertz-Mindlin with flexible bond contact laws and DOE optimized material properties successfully reproduced the mass fractions retained in ASABE sieves at 9.8% mean relative error and a coefficient of determination of R2 = 0.87. Here, the DEM methodology developed for mechanical WPC mass fractionation could be deployed to perform virtual design of feedstock handling equipment and performance analysis of mechanical fraction systems.