SEARCH · Engineering Papers

Results for “Stiffness”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Design optimization of lightweight automotive seatback through additive manufacturing compression overmolding of metal polymer composites

With the growing demand for enhanced automotive fuel efficiency and environmental sustainability, there is a need for lightweighting automotive components through innovative design and manufacturing processes. Here, this study leverages a combination of numerical iterative design optimization and hybrid additive manufacturing–compression molding (AM-CM) technique for metal polymer composites to lightweight an automotive seatback. The AM-CM process enables robust mechanical interlocking between metals and composites, boasting high stiffness and strength with low overall density. Replacing metallic components with such metal polymer composites allows for comparable mechanical performance while significantly reducing the overall weight. First, the automotive seatback design space is reduced to critical load carrying regions using topology optimization and high stress concentration areas are identified using finite element analysis. Next, a lightweight metal polymer subcomponent is designed for a high stress concentration region. The full seatback frame with spatially heterogeneous material-specific design is then iteratively optimized to enable enhanced stiffness with minimal weight. Overall, the automotive seatback frame designed with location-specific metal, polymer, and metal polymer composite materials weighs 20% less than the metal-only design while exhibiting similar stiffness.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Testing and modeling of CLT-to-CLT joints made using hardwood dowels and screws with varying spacing

With the increasing popularity of mass timber, utilizing large-diameter wood screws and hardwood dowels can offer an easy and cost-efficient way to join structural elements. However, the existing research on hardwood dowels is limited to single fastener joints under monotonic loading. No research is available on the group effects of joints with multiple hardwood dowels. Here, this paper presents the results of 96 monotonic and cyclic single-shear plane CLT-to-CLT joint tests with hardwood dowels (Red Oak and Yellow Birch, 25.4 mm diameter) and self-tapping wood screws (8 mm diameter) in Grand-Fir CLT, where screws were inserted both perpendicular and at a 45-degree angle to the grain for comparison to hardwood dowels. The study evaluated the mechanical properties of the joints by changing fastener spacing. Mechanical properties such as yielding and peak strength, displacement, elastic and yielding stiffness, and ductility were evaluated. The results showed equal strength and stiffness properties for hardwood dowels comparable with large-diameter screws. Additionally, the hardwood dowel joints demonstrated moderate ductility properties. The strength and stiffness of the tested joints were compared with analytical equations found in the literature by considering the group action factor of the fasteners. Lastly, a nonlinear force-displacement model was presented and compared to the experimental results and analytical calculations, which gives the possibility to predict and optimize the mechanical behavior of hardwood dowel and screw joints with varying numbers and spacing between fasteners.

36 MATERIALS SCIENCE↗

DOE OSTI2024

One-sweep moment-based semi-implicit-explicit integration for gray thermal radiation transport

Thermal radiation transport (TRT) is a time dependent, high dimensional partial integro-differential equation. In practical applications such as inertial confinement fusion, TRT is coupled to other physics such as hydrodynamics, plasmas, etc., and the timescales one is interested in capturing are often much slower than the radiation timescale. As a result, TRT is treated implicitly, and due to its stiffness and high dimensionality, is often a dominant computational cost in multiphysics simulations. Here we develop a new approach for implicit-explicit (IMEX) integration of gray TRT in the deterministic SN setting, which requires only one sweep per stage, with the simplest first-order method requiring only one sweep per time step. The partitioning of equations is done via a moment-based high-order low-order formulation of TRT, where the streaming operator and first two moments are used to capture the asymptotic stiff regimes of the streaming limit and diffusion limit. Absorption-reemission is treated explicitly, and although stiff, is sufficiently damped by the implicit solve that we achieve stable accurate time integration without incorporating the coupling of the high order and low order equations implicitly. Due to nonlinear coupling of the high-order and low-order equations through temperature-dependent opacities, to facilitate IMEX partitioning and higher-order methods, we use a semi-implicit integration approach amenable to nonlinear partitions. In conclusion, results are demonstrated on thick Marshak and crooked pipe benchmark problems, demonstrating orders of magnitude improvement in accuracy and wallclock compared with the standard first-order implicit integration typically used.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Hourglass control in staggered-grid hydrodynamics using virtual element stabilization techniques

Numerical simulations using the staggered-grid hydrodynamics (SGH) discretization suffer from hourglass instabilities. In this work, we develop a stabilization method to suppress the hourglass instabilities using techniques from the virtual element method (VEM). The stiffness matrix of the VEM consists of two terms: the consistency matrix which is rank deficient and the stability matrix. Here, we first show that in two dimensions and on general polygons, the stiffness matrix of the SGH is identical to the consistency matrix of the linear VEM for both the diffusion equation and the linear elasticity equation. These analyses explain the origin of the hourglass instabilities of the SGH discretization method, and establish a theoretical foundation for our proposed stabilization method by augmenting the stiffness matrix of the SGH discretization using the VEM stability matrix. Then, we present numerical examples using Lagrangian SGH simulations. The numerical experiments demonstrate that the proposed VEM stabilization method is effective at eliminating hourglass modes in the SGH discretization.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

A parametric study of slow dynamic nonlinear elasticity with comparisons to models

Several phenomenological models that aspire to quantitative description of anomalous nonlinear mesoscopic elasticity are reviewed and compared with laboratory measurements. This class of nonlinearity, best known perhaps for slow dynamics and aging, is seen widely in imperfectly consolidated granular solids but is not well understood. Typical slow dynamic tests show that a modest conditioning oscillatory "pump" strain depresses material stiffness, which then recovers like the logarithm of time after conditioning ceases. Several phenomenological models based on physical arguments have been proposed that predict the material stiffness response to arbitrary pump strain histories during conditioning and recovery. Approximate closed form and numerical solutions to the models are presented that predict the quantitative influence of three key pump parameters: the pump's strain amplitude, the pump's strain rate, and the pump’s duration. Laboratory measurements on Berea sandstone, concrete and a confined single aluminum bead find that slow dynamic responses are linear in pump strain and independent of pump frequency. Measurements also show that, after pump-off, stiffness recovers over times far longer than the pump duration. These observations and others are compared to model predictions. One of the considered models, based on a picture of fast brittle damage and slow healing, successfully matches all these behaviors.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Graph neural networks for mechanical property prediction of 2D fiber composites

This work investigates the ability of graph neural networks (GNNs) to homogenize 2D fiber composite microstructures. We use different inhomogeneity and anisotropy indices to motivate and show that the Volume Elements (VEs) used in ML methods should ideally be far from their Representative Volume Element (RVE) size limit and, consequently, are notably anisotropic. Hence, training only the isotropic limit properties may not be acceptable. Another aspect is the need to normalize elastic stiffness values for ML, especially when high elastic contrast ratios are encountered between composite phases or in the material set. We introduce a normalization technique based on the mean-field method (MFM) to handle such high contrast ratios and train for the entire stiffness tensor. We show that the proposed GNN approaches exhibit high accuracy and efficiency compared to traditional methods and convolutional neural networks, utilizing unstructured graphs constructed from microstructure topology. Our model successfully predicts the stiffness tensor, peak strength under bulk damage, and brittle fracture initiation strength across diverse microstructure configurations while maintaining high accuracy even for extreme material contrasts and volume fractions. We also present a method to improve prediction accuracy for small dataset sizes using Voronoi partitioning.

Brittle strength↗

DOE OSTI2025

Computational alchemy clarifies origins of alloy strengthening

Solid solution strengthening (SSS) is widely used to enhance mechanical properties of metals. Originally developed for dilute alloys, classical SSS theories are presently challenged by the rise of complex concentrated alloys (CCA) with nearly equiatomic compositions. Here, we propose and develop a method of “computational alchemy” in which interatomic interactions are modified to systematically vary two key physical parameters defining SSS - atomic size misfit and elastic stiffness misfit - over a maximally wide range of two misfits. The resulting alchemical alloys are subjected to massive (~10 8 atoms) molecular dynamics (MD) simulations reproducing full complexity of plastic strength response. At variance with prevailing views, stiffness misfit is observed to contribute to SSS on par if not more than size misfit. Furthermore, depending on exactly how two misfits are combined, they result in synergistic (amplification) or antagonistic (compensation) effect on alloy strengthening. Unlike real CCAs in which each component element comes with its own specific size and stiffness, our alchemical model alloys span the space of two misfits continuously revealing trends in alloy strengthening unrecognized so far. Our study demonstrates unique value of intentionally unrealistic models for gaining deep physical insights into material behaviors that are difficult to reveal otherwise.

36 MATERIALS SCIENCE↗

DOE OSTI2025

An eigenvalue-based method for computing the relaxed pressure in compressible multiphase flow with N phases

The modeling of compressible multiphase flows is a decades-old area of study with many applications across various fields. Many of these application areas use stiff pressure relaxation. This process involves the solution of a nonlinear system with N + 1 equations and N + 1 unknowns, where N is the number of phases. The resolution of this system with general equations of state (EOSs) is difficult. Furthermore, nonlinear systems can admit multiple solutions, and current solution methods do not address this possibility. Very recently, a thermodynamic relaxation method was introduced, which effectively maps a relatively simple predictor equation of state onto a more complex target equation of state. In this context, the target EOSs are the chosen EOSs for the thermodynamic model. Furthermore, this thermodynamic relaxation has the benefit of simplifying the stiff pressure relaxation system of equations. In this article, we show this system reduces to a polynomial of degree N, which can be recast as an eigenvalue problem through the use of the associated companion matrix. We show that although this eigenvalue method is generally less efficient than Newton–Raphson iteration, it does not suffer from convergence issues and finds all N roots of the polynomial. Hence, the method provides a fail-safe for root-finding iterative methods and a way to address the issue of multiple solutions to the nonlinear system of equations in stiff pressure relaxation.

Eigenvalue algorithm↗

DOE OSTI2025

Mineralized sclerites in the gorgonian coral Leptogorgia chilensis as a natural jamming system

The soft corals (Cnidaria, Octocorallia), a diverse group of colonial marine invertebrates, can reversibly tune their body stiffness in response to external stimuli. This capability is attributed to their dynamic skeletal systems, which consist of thousands of mineralized skeletal elements, called sclerites, embedded within a gel-like matrix that swells/deswells and unjams/jams the sclerites, thus modulating skeletal stiffness. While sclerite morphology is widely used for species identification, its role in the mechanical performance of a soft coral’s skeletal system is largely unknown. Here, we investigated structure-jamming relationships in sclerite-based skeletal architectures using the red gorgonian octocoral Leptogorgia chilensis as a model system. The sclerites of L. chilensis exhibit a shaft-like geometry with two axial branches and two sets of triradiate side branches, which are aligned with the crystallographic symmetry of the constituent magnesium-containing calcite. By combining multiscale three-dimensional (3D) structural characterization, parametric geometrical modeling, 3D printing, mechanical testing, and discrete element simulations, we demonstrate how sclerite geometry achieves a balanced jamming performance in terms of stiffness, weight, strength, and fracture resistance in comparison to alternative geometries parametrically modified from the native sclerites (e.g., changes in the length and number of side branches). Here, we also found that these performance metrics are achieved through the effective interlocking among side and axial branches, which is further enhanced by the fractal-like microscopic spikes on the branch tips. The findings in this natural jamming system offer insights for designing synthetic mechanotunable material architectures for a wide range of applications, from soft robotics to mechanical dampeners.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Impact of composition and symmetry energy on the temperature of quasiprojectiles simulated with antisymmetrized molecular dynamics

The equation of state describes the emergent physical properties of matter. Experimental data is needed to help constrain the equation of state for nuclear matter. These constraints can help distinguish between an “asy-stiff” and an “asy-soft” equation of state, which has astrophysical implications. One path to help constrain the models is to analyze the nuclear caloric curve; some experiments have shown dependence on neutron excess, and may thus be sensitive to the asymmetry. A difference in the caloric curve based on the asymmetry of the reconstructed quasiprojectile (QP) had been observed using 70 Zn on 70 Zn at 35 MeV/nucleon taken with the NIMROD array. Antisymmetrized molecular dynamics calculations were performed for the same system and deexcited with gemini++. Both Gogny (asy-soft) and Gogny-as (asy-stiff) data sets were generated. The particles were then filtered based on detector geometric acceptance and thresholds. From the accepted particles, the excitation energy and temperature were calculated in the same way as for experimental data. Additionally, filter effects on the observed nuclear caloric curves were investigated. A tendency for the asy-stiff nuclear caloric curves to have higher temperatures than their asy-soft counterparts was observed for a number of probes. In addition, some probes may show sensitivity to the reconstructed composition of the QP, but this is inconclusive due to high statistical fluctuations and a large dependence on the exact method of event selection.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Bayesian inference of fine features of the nuclear equation of state from future neutron star radius measurements to 0.1 km accuracy

To more precisely constrain the equation of state (EOS) of supradense neutron-rich nuclear matter, future high-precision x-ray and gravitational wave observatories are proposed to measure the radii of neutron stars (NSs) with an accuracy better than about 0.1 km. However, it remains unclear what particular aspects (other than the stiffness generally spoken of in the literature) of the EOS and to what precision they will be better constrained. In this work, within a Bayesian framework using a metamodel EOS for NSs, we infer the posterior probability distribution functions (PDFs) of incompressibility K 0 and skewness J 0 of symmetric nuclear matter (SNM) as well as the slope L, curvature K sym , and skewness J sym characterizing the density dependence of nuclear symmetry energy E sym ⁡(ρ), respectively, from mean values of NS radii consistent with existing observations and an expected accuracy Δ⁢R ranging from about 1.0 to 0.1 km. Here, we found that (1) the Δ⁢R has little effect on inferring the stiffness of SNM at suprasaturation densities, (2) smaller Δ⁢R reveals more accurately not only the PDFs but also pairwise correlations among parameters characterizing high-density E sym ⁡(ρ), (3) a double-peak feature of the PDF(K sym ) corresponding to the strong K sym – J sym and K sym – L anticorrelations is revealed when Δ⁢R is less than about 0.2 km, and the locations of the two peaks are sensitive to the maximum value of J sym reflecting the stiffness of E sym ⁡(ρ) above about 3 times the saturation density ρ 0 of SNM, and (4) the high-precision radius measurement for canonical NSs is more useful than that for massive ones for constraining the EOS of nucleonic matter around (2–3)⁢ρ 0 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Electronic density of states as the descriptor of elastic bond strength, ductility, and local lattice distortion in BCC refractory alloys

Although electronic density of states (DOS) is fundamental to materials properties, its general relationship to mechanical properties of alloys is not well established. In this paper, using density functional theory (DFT) calculations, we show that the electronic occupancy at the Fermi level, N(E f ), obtained from DOS is a key descriptor of alloy strength and ductility. Our comprehensive analysis of numerous body centered cubic (BCC) refractory high entropy alloys (RHEAs) shows an overwhelming correlation that low N(E f ) indicates strong bonds that have high stiffness resulting in high elastic constants. High bond stiffness indicates presence of covalent nature of bonds that are directional in nature resulting in resistance to deformation leading to high bulk (B) and shear (G) moduli. Consequently, N(E f ) provides a direct correlation to the tendency of alloy ductility evidenced in the Pugh ratio (G/B). As stiffer bonds result in lower local lattice distortion (LLD), N(E f ) are LLD are also found to be corelated which opens up a correlation to solid solution strengthening and yield strength. Thus, this work unveils fundamental correlations between N(E f ) and (1) elastic bond strength, (2) ductility, and (3) LLD. These correlations open opportunities for the design of high strength high ductile RHEAs.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Foundational insights into the mechanical and molecular evolution of porcine skin gelatin during gelation

Gelatin is a widely used material in biomedical fields, particularly in regenerative medicine and tissue engineering, due to its biocompatibility and versatile properties. While prior research has explored methods to enhance gelatin's mechanical strength and stability, fundamental studies on gelatin, specifically its curing process, mechanical stiffness, and chemical evolution during gelation, remain limited. This study uses ultrasonic testing and Fourier Transform Infrared Spectroscopy (FTIR) to examine gelatin's stiffness and molecular changes during gelation. Samples of 175 and 300 Porcine Skin Bloom Strength Gelatin at concentrations of 2% and 6% (w/v) were analyzed. Through transmission ultrasonic testing helped identify key transition points in gelation, with higher concentrations exhibiting delayed transitions. FTIR revealed that C-N bond formation peaks early while N-H bond deformation persists. A correlation emerged between sound speed and peak absorbance, suggesting that changes in molecular mobility may contribute to the observed sound speed behavior during periods of active bond formation. However, as gelation continues, fewer bonding components may be available, potentially decreasing molecular movement and contributing to the observed increase in sound speed. These findings provide insights into gelatin's mechanical and chemical evolution, offering a framework for improved control over its gelation kinetics. Swept-Frequency Acoustic Interferometry (SFAI) was performed at the end of the curing process to measure the sound speed, enabling the calculation of the bulk moduli of the gelatin samples. The combined use of ultrasonic and FTIR testing provides a non-destructive method for characterizing gelatin and other biomaterials. This approach advances understanding of gelatin curing behavior and supports the development of safer biomaterials with tailored mechanical properties for various applications such as tissue engineering and regenerative medicine.

Biomaterials↗

DOE OSTI2025

A Printed Microscopic Universal Gradient Interface for Super Stretchable Strain‐Insensitive Bioelectronics

Abstract Stretchable electronics capable of conforming to nonplanar and dynamic human body surfaces are central for creating implantable and on‐skin devices for high‐fidelity monitoring of diverse physiological signals. While various strategies have been developed to produce stretchable devices, the signals collected from such devices are often highly sensitive to local strain, resulting in inevitable convolution with surface strain‐induced motion artifacts that are difficult to distinguish from intrinsic physiological signals. Here all‐printed super stretchable strain‐insensitive bioelectronics using a unique universal gradient interface (UGI) are reported to bridge the gap between soft biomaterials and stiff electronic materials. Leveraging a versatile aerosol‐based multi‐materials printing technique that allows precise spatial control over the local stiffnesses with submicron resolution, the UGI enables strain‐insensitive electronic devices with negligible resistivity changes under a 180% uniaxial stretch ratio. Various stretchable devices are directly printed on the UGI for on‐skin health monitoring with high signal quality and near‐perfect immunity to motion artifacts, including semiconductor‐based photodetectors for sensing blood oxygen saturation levels and metal‐based temperature sensors. The concept in this work will significantly simplify the fabrication and accelerate the development of a broad range of wearable and implantable bioelectronics for real‐time health monitoring and personalized therapeutics.

Song, Kaidong [Department of Aerospace and Mechani↗

DOE OSTI2026

Physicochemical and Performance Characterization of Six Commercial Organic Solvent Nanofiltration Membranes

This work introduces a novel, gradient-free metamaterial design method based on Gaussian process regression to represent the density field of a unit cell. The dimension of the design space is determined by the covariance matrix dimension in the Gaussian process regression. We propose compressing this matrix using an autoencoder, enabling the decoder to generate the density field and effectively reduce the originally large design space to a lower-dimensional subspace. In this compressed space, we employ an active learning method, Bayesian Adaptive Direct Search (BADS), for efficient exploration of the design space. We demonstrate that for simple 2D designs aimed at maximizing unit cell stiffness, our method yields results comparable to those of standard topology optimization. Furthermore, we extend our approach to various mechanical problems, from linear elasticity to hyperelastic large deformation and elasto-plasticity under finite deformation, to 3D metamaterial design. This illustrates the method’s versatility and effectiveness across a range of applications.

Wu, Haoran↗

DOE OSTI2026

Extremized nonlinear and linearized responses in soft metamaterials enabled by gradient-based design and grayscale digital light processing

In this study, we develop a gradient-based design approach that exploits grayscale digital light processing (DLP) 3D printing for extremizing the nonlinear and linearized response of soft metamaterials — materials that harness engineered geometric instabilities to undergo large and programmable changes in configuration. Grayscale DLP approaches modulate local mechanical properties at the pixel scale by tuning the light intensity within a single grayscale image, unlocking an exceptionally large design space. To effectively navigate this space, we develop smooth mappings between local light intensity values and global quantities of interest that characterize the behavior of soft metamaterials. Enabling these smooth mappings are robust and differentiable nonlinear finite element simulations powered by a trust region solver. A PDE-constrained optimization problem is then solved to invert these mappings and produce light intensity distributions that endow the printed part with varying stiffness and flexibility in distinctive regions. It is shown that optimizing the distribution of soft and stiff phases throughout a metamaterial structure results in markedly different buckling and self-contact configurations to drive extremized nonlinear compression and linearized vibration responses. Optimized light intensity distributions are translated to grayscale images and directly used to print soft metamaterial samples, showing remarkable agreement between the buckling and self-contact response in simulated and measured deformed configurations.

Additive manufacturing↗

DOE OSTI2026

Ionic Interdiffusion at Cathode|Solid-Electrolyte Interface: A Machine Learning–Assisted Multiscale Investigation and Mitigation Strategies

Future lithium batteries are expected to use solid electrolytes to achieve higher energy density and fast charge capabilities. However, most solid electrolytes are thermodynamically unstable against layered oxide cathodes. In this study, the stability of LiCoO2 (LCO) cathode with Li10GeP2S12 (LGPS) solid electrolyte is investigated using ab initio molecular dynamics (AIMD) and machine learning molecular dynamics (MLMD). The propensity of ionic interdiffusion, formation of a passivating interphase layer, and corresponding decay in cell performance is addressed using a continuum model. Large-scale MLMD simulations confirm that the LCO|LGPS interface permits interdiffusion of cobalt (Co) and other ionic species, leading to the formation and growth of a resistive interphase and to dramatic capacity fade even in the first cycle. We examine the literature evidence that incorporating a thin layer of LiNb0.5Ta0.5O3 (LNTO) between LCO and LGPS prevents the interdiffusion of ions. Atomistic simulations suggest that substituting lithium (Li) in LNTO with Co is thermodynamically unfavorable, thereby inhibiting ionic interdiffusion. The stable Nb5+/Ta5+ states form a rigid metal-oxide framework, which consequently also prevents the substitution of niobium (Nb) or tantalum (Ta). However, continuum-level analysis suggests that the higher mechanical stiffness of LNTO can lead to interfacial delamination between the LCO and LNTO. This phenomenon reduces the effectiveness of the protective layer. This paper, therefore, highlights the need to develop novel interlayers that balance low ionic interdiffusion with low mechanical stiffness.

Ncube, Musawenkosi K.↗

DOE OSTI2025

Elucidating texture and grain morphology contributions to the micromechanical response of additively manufactured Inconel 625

Microstructural variation of additively manufactured (AM) metal components in comparison to wrought counterparts makes certification for critical applications a challenge. Microscale simulations leveraging modern computational tools may be used to supplement testing of AM microstructures, thus accelerating certification by reducing the number of experiments needed. However, as micromechanical response is closely tied to critical properties like fatigue-life and fracture, utilization of these simulations with macroscale experimental data alone is insufficient. One means to attain microscale experimental data is in situ diffraction data collected from synchrotron X-ray sources. In this work, such data were collected during in situ compression of AM Inconel 625 superalloy. Interpretation of experimental results was assisted by massive (8M element) complementary micromechanical simulations performed on sets of virtual microstructures generated using cellular automata. Together, micromechanical data from diffraction experiments and simulations were used to probe the effects of textured “track” microstructures generated during laser powder bed fusion and directional strength-to-stiffness on micromechanical response. Though fiber-averaged directional strength-to-stiffness ratios were expected to dominate given the high elastic anisotropy of the material, the combination of small variations in texture and specific grain configurations unique to AM microstructures lead to significant variability in micromechanical response after yield. The findings emphasize the importance of high-fidelity microstructural representation that captures key texture components and AM-specific morphology for property prediction of AM metals.

36 MATERIALS SCIENCE↗