Engineering Papers⌕ Search

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.

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

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Adaptive Stiffness Power Take-Off System for a Hydraulic Hose-Pump Point Absorber Wave Energy Converter

Two crucial challenges faced in the design of ocean wave energy converters (WECs) are maintaining high capacity factors and remaining robust in everchanging, harsh ocean environments. This study introduces a novel mechanism for achieving reactive WEC control with a soft hydraulic power take-off, thereby addressing both challenges. The device of interest involves a hydraulic pump, constructed from a fluidic flexible matrix composite (F2MC), which is cyclically stretched by driving waves and pressurizes working fluid through a turbine-generator. This class of fiber-reinforced stretch hoses is established for ocean use as robust moorings and shock absorbers. The dynamic properties of the F2MC pump depend on internal pressure, which can be controlled by a continuously variable transmission between the turbine and generator, allowing the system to be adapted for better performance over a wide range of wave conditions. This work presents a dynamic model for operation and adaptive behavior of this power take-off system, experimentally validates components of the model, and performs dynamic simulations of the system that show increases in power production and capacity factor compared to non-adaptive operation for a variety of driving wave conditions. This work develops the basis for future large-scale experimental testing and controller design of this system.

Wave energy conversion, fluidic flexible matrix co↗
DOE OSTI2026

Mechanical Characterization of the NIF Ignition Target TMPSA Bonding Flexure

The Thermo-Mechanical Package Sub-Assembly (TMPSA) provides a critical mechanical, thermal, and electrical interface between a silicon arm and a TMP aluminum can in a NIF ignition target. During assembly, sixteen silicon pads are bonded to the aluminum using a fixture that positions the components and applies a repeatable prescribed displacement through a compliant flexure. The flexure converts fixture interference into displacement and reaction force. Because bondline thickness variation must be maintained within +/-4 µm, consistent flexure behavior is important to the assembly process. With a recent string of TMPSAs exhibiting low bond strength, the flexures were inspected. Despite being manufactured to the same specifications, flexures were found to exhibit variation in measured stiffness. Additionally, the measured stiffness values did not always follow the presumed beam mechanics model. This work addresses two questions: whether the bonding fixture is working as intended, and whether the previously made stiffness measurements are accurate representations of the use case. The flexures are analyzed using Euler-Bernoulli beam theory, measured beam dimensions, finite element analysis, and tolerance stack-up calculations. The analysis shows that the fixed-guided beam mechanics model appropriately represents the flexure during TMPSA bonding, but the chisel-tip stiffness measurement method introduces a deformation to the inner ring of the flexure that is not represented during use. The measured stiffness values should therefore be interpreted as test-condition stiffness values rather than direct measurements of operational flexure stiffness. The discrepancy is therefore attributed primarily to the measurement boundary condition rather than to failure of the fixed-guided beam model. Recommendations are provided for GD&T, dimensional inspection, and a representative stiffness testing method to better control bondline variation.

42 ENGINEERING↗
DOE OSTI2026

Additively Reinforced Thermoformable PETG Composite Sheets for Improved Structural Efficiency

Thermoforming of short-fiber reinforced thermoplastic sheets offers a viable pathway for producing lightweight composite components; however, inherent anisotropy in fiber-reinforced sheets can limit structural performance under multidirectional loading. In this work, short carbon fiber, glass fiber, and hybrid fiber–reinforced PETG sheets were evaluated as candidate feedstock materials for thermoforming, with flexural and tensile testing performed both along the primary fiber direction and in the off-axis orientation to establish baseline stiffness, strength, and anisotropy. As expected, short carbon fiber PETG exhibited the highest stiffness and strength in the primary fiber direction, while all systems showed reduced performance in the off-axis direction. This off-axis performance reduction provides clear justification for the use of additive reinforcement when such thermoformed sheets are intended for structural applications. The intended manufacturing sequence involves thermoforming the reinforced sheet first, followed by the application of additively manufactured lattice reinforcement; therefore, the reinforcement strategy does not impose limitations on sheet formability during thermoforming. Post-forming lattice reinforcement significantly reduced load-normalized displacement by approximately 95–99% relative to non-reinforced sheets and improved weight-normalized stiffness by ~70%. These findings demonstrate that geometry-driven additive reinforcement can effectively compensate for off-axis property reductions in thermoformed PETG composites, enabling enhanced multidirectional structural performance without compromising manufacturability.

Talabi, Isaac [ORNL] (ORCID:0000000340215594)↗
DOE OSTI2025

Regulating hydrogel mechanical properties with an electric field

Stimuli-responsive polymeric materials have attracted significant attention due to their ability to change properties in response to various external stimuli. Using an electric field as the stimulus is of particular interest as it possesses the potential for seamless integration of materials with electronic systems. While many materials with electric field responsive actuation have an associated mechanical property change, it is beneficial to develop materials that exhibit mechanical property changes without accompanying significant shape deformation. To address this challenge, here we designed a semi-interpenetrating polymer network (semi-IPN) hydrogel system containing both polyelectrolytes and salt ions, which enables electric field induced changes in mechanical properties while minimizing actuation. We first successfully verified the viability of our design by removing salt ions through a diffusion-only method where we witnessed the stiffness increased to 4.5 times the initial value while still being highly deformable. After this, we applied an electric field to transport the salt ions out of the hydrogel, as shown by both Raman spectroscopy and scanning electron microscopy. We were able to show a time-dependent stiffness increase, the maximum of which was 5 times the original stiffness. We quantified ion transport and water-splitting in the hydrogel by both experiments and simulations. Following this, we showed functional system reversibility by reversing the direction of the current to reinject salt ions into the semi-IPN hydrogel and reducing its stiffness to below the initial value. It's worth noting that our simulations enable us to understand the governing mechanisms behind ion generation and salt transport that leads to mechanical property changes. Finally, we were able to fabricate a spatially variable stiffness haptic interface with our hydrogel, with demonstrated reversibility and cyclability. This research can possibly find applications in soft robotics and haptics and also inspire the development of bio-compatible electronics related devices.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effect of inversion asymmetry on the superconducting and exciton condensates of bilayer graphene

Inversion asymmetry in bilayer graphene can be tuned by the displacement field. As a result, the band dispersion in biased bilayer graphene acquires flatband regions near the Dirac points along with a nontrivial band geometry. We analyze the effect of inversion asymmetry on the critical temperature and superfluid stiffness of the superconducting state of AB-stacked graphene bilayer and the exciton condensate in double layers formed by two AB-stacked graphene bilayers. We find that the geometric superfluid stiffness in bilayer graphene superconductors is negligible due to the small superconducting gap. Furthermore, since the geometric superfluid stiffness is maximized for a constant order parameter, it can be neglected in biased bilayer graphene superconductors with any pairing symmetry. In contrast, the displacement field enhances the geometric superfluid stiffness in exciton condensates. It is most prominent at low densities and high displacement fields. Here, a consequence of the geometric superfluid stiffness is a modest enhancement of the Berezinskii-Kosterlitz-Thouless transition temperature in bilayer graphene’s exciton condensate.

BKT transition↗
DOE OSTI2026

Mechanically graded granular scaffolds for osteochondral tissue engineering

Engineered scaffolds designed to approximate the mechanical microenvironment of the osteochondral unit often address this complexity using discrete, two-phase architectures that introduce mechanical discontinuities and interfacial stress concentrations rather than a contiguous stiffness transition. To address this challenge, we created a photoannealed polyethylene glycol (PEG) granular scaffold with a spatially controlled stiffness gradient within a cell-permissive, macroporous architecture. Stiffness was dictated by photoannealing microgels using a photomask. We tuned void volume and available surface area by varying microgel diameter and tested how mesenchymal stromal cells (MSCs) interpret local mechanical environments. MSCs exhibited position-dependent differences in morphology, cytoskeletal structure, matrix deposition, and lineage-specific gene expression within the gradient scaffolds. Softer regions supported rounded cell morphology and deposition of a glycosaminoglycan-rich matrix, whereas stiffer regions promoted cell elongation, increased cytoskeletal tension, and expression of mineral-associated markers. Gradients formed from smaller microgels magnified these spatial responses by increasing cellular confinement and adhesion site availability. Disruption of actomyosin contractility eliminated these regional differences, demonstrating that MSCs rely on tension-dependent mechanotransduction to interpret the gradient. These findings reveal that coupling microgel architecture with continuous stiffness transitions provides a tractable platform to study multiscale mechanobiologic regulation and spatially guide osteochondral tissue formation.

Biological and medical sciences↗
DOE OSTI2025

Toward engineering lattice structures with the material point method (MPM)

This study examines the potential of two variants of the material point method—the generalized interpolation material point (GIMP) and dual domain material point (DDMP) methods—in developing a robust computational framework for engineering lattice structures under different loading conditions. The study begins with assessing the ability of the two methods in predicting elastic buckling phenomena using column geometries with and without initial geometric imperfections. The results indicate that both methods effectively capture buckling phenomena when initial geometric imperfections are introduced. After this verification step, we create several models of tetrahedral lattice structures with varying strut diameter and orientation and subject them to quasi-static loading. We then validate the numerical results using laboratory test results. The results show that, while both methods accurately predict load–displacement curves in the pre-buckling regime, their predictive capabilities diminish in the post-buckling regime. Through visual comparison between the numerical and experimental deformed shapes, it appears that the discrepancies between model and experimental results are attributed to initial geometric imperfections in the lattices that occurred during 3D printing. We then establish a second set of lattice models where different types of initial geometric imperfections are considered. The results from these models show that imperfections have a negligible influence in the pre-buckling regime but affect the behavior considerably in the post-buckling regime. As a final step in this work, we subject the lattice models to impact loading and employ hypothetical soft and stiff materials. These results show that the lattice stiffness, which depends on material stiffness, strut diameter, and orientation, significantly influences the ability of a lattice structure to resist impact. In particular, we find that a stiffer lattice (i.e., one made with a stiff material and thicker struts) is capable of absorbing more energy than a softer one during impact. Although material nonlinearities, inelasticity, and detailed contact formulations are not considered in this study, the findings obtained herein lay the groundwork for engineering lattice structures under extreme loading conditions through a simulation-driven framework based on particle-based methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Nonlinear viscoelastic response of silicone additively manufactured direct ink write (DIW) foams under repetitive compression

To investigate dynamic fatigue behavior of foam in military protective applications, such as helmets, additively manufactured (AM) foams were compressively strained into the plateau region using a reduced design of experiments. A simple power law was found to govern the decline in dynamic stiffness (complex modulus) as the foams underwent the purchase order requirement of 10,000 cycles of small deformation in the plateau region. This rate of decline was newly found to correlate with the degree of nonlinearity in the material’s deformation, quantified using total harmonic distortion. Materials with low nonlinearity exhibited relatively stable stiffness across cycles, while those with high nonlinearity experienced greater losses. The observed nonlinearity depended on both applied stress and strain rate. A strong linear correlation (R 2 = 0.78) was identified between second-order nonlinearity and the time-dependent stiffness response. Two lattice structures were examined: face-centered tetragonal (FCT) and simple cubic (SC). The SC material exhibited higher total harmonic distortion (5%) and lower stiffness retention than the FCT (2%). In conclusion, these results suggest that for cyclic compression applications in a wide variety of industries such as packaging, personal protective equipment, or aerospace, selecting materials with lower stress and greater structural uniformity can enhance the stability of dynamic performance.

42 ENGINEERING↗
DOE OSTI2025

Correlations between maximum mass of neutron stars and the nuclear matter properties and the constraints from PSR J0740+6620 and GW190814

The correlations between the maximum mass of neutron stars (NSs) with the properties of nuclear matter at saturation density ρ 0 , i.e., the incompressibility K 0 , the symmetry energy E sym (ρ 0 ), and its slope L(ρ 0 ), the isoscalar effective mass m$^{*}_{s0}$, and the parameter f I , have been investigated using several sets of Skyrme interactions based on the momentum-dependent (MD) Skyrme interactions SAMi-J27 (soft), SAMi-J31 (stiff), and SAMi-J35 (super-stiff), as well as the momentum-independent (MI) Skyrme interaction SkT5. Only one physical quantity is varied at a time for a given set of Skyrme interactions. It is found that the maximum mass of NSs has (1) a strong positive linear correlation with K 0 for stiff and super-stiff equation of state (EOS) and no correlation with K 0 for soft EOS; (2) a negative logarithm correlation with the symmetry energy E sym (ρ 0 ) for all the EOSs; (3) a positive logarithm correlation with the slope L(ρ 0 ) for all the EOSs; (4) a positive logarithm correlation with the parameter f I for MD EOS; (5) a negative power-law correlation with m$^{*}_{s0}$ for MD EOS. Based on these established correlations, the constraints to the nuclear matter properties from the mass of PSR J0740+6620 and GW190814 are investigated. Furthermore, we find that while rotation enhances the maximum mass of NSs, it does not effect the underlying correlations between the maximum mass and the aforementioned nuclear matter parameters.

Zhou, M. [Shaanxi Normal University, Xi’an (China)↗
DOE OSTI2026

Influence of Numerical Modeling Approaches on Damped Behavior of Flexible Beams: Preprint

Composites structures are widely used in aerospace and wind energy applications for their excellent stiffness and strength-to-weight properties. In these structures, structural damping is critical to predict vibration amplitudes, performance, and reliability. Structural damping is of particular interest for slender wings, rotorcraft blades, and wind turbine blades that can exhibit complex vibration phenomena and are frequently modeled with geometrically exact beam theory (GEBT). Standard approaches of stiffness proportional or modal damping merely assign user defined values and cannot predict damping behavior. This work compares stiffness proportional damping to two more advanced damping approaches: modal strain energy and Prony series. The modal strain energy approach uses a sectional analysis tool to calculate the beam stiffness and postprocess internal stresses from GEBT simulations. The internal stresses are then used to calculate modal damping factors. The Prony series is implemented within GEBT to directly model viscoelastic behavior of the composites. These approaches are compared by modeling the evolution of the damping factors of a realistic flexible wind turbine blade with varying rotational speed. Discrepancies between the approaches suggest areas for future modeling development, but differences in nonlinear damping values are less than current uncertainties about the magnitude of structural damping.

17 WIND ENERGY↗
DOE OSTI2025

Dynamic response of a freely rotating butterfly valve in the advanced test reactor − dynamic coefficients modeling

Here, in evaluating the water hammer issue pertaining to the primary-coolant-regulating butterfly valve in the Advanced Test Reactor, the dynamic fluid body interaction (DFBI) approach was implemented in the analysis covered in Part I. Although DFBI modeling accurately and simultaneously solved the dynamic motion of the valve’s disk along with the flow field of the surrounding fluid, it shed little light on the reason behind such motion. For Part II, the reacting torque of the fluid on the disk was decomposed into representations of the dynamic coefficients in terms of stiffness, damping, and added mass. These were evaluated via simulations with steady-state static (stiffness), constant angular speed (damping), and variable angular speed (added mass) disks. Substituting the dynamic coefficients into Newton’s second law enabled the response trajectories to be obtained. Stable (by average) and unstable equilibrium positions and thrust tendencies of the valve were determined based on the stiffness coefficient (or static torque), the response amplitude was dampened or enlarged by the damping coefficient (minorly affected by added mass), and the response frequency was altered by the damping and added mass coefficients. Although the dynamic coefficient approach renders slightly different trajectories, due to the averaging effect of the torque in comparison to the DFBI method, the overall trend of the response aligns with the DFBI simulation, thus confirming the conclusion in Part I that a fix to the current butterfly valve is necessary.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

A review of high order strong stability preserving two-derivative explicit, implicit, and IMEX methods

High order strong stability preserving time discretizations ensure the nonlinear non-inner-product strong stability properties of spatial discretizations suited for the stable simulation of hyperbolic PDEs in a wide variety of application areas including fluid dynamics, magnetohydrodynamics, semiconductor devices, electromagnetics, and astrophysics. Over the past decade multiderivative time-stepping have been increasingly used for the time-evolution hyperbolic PDEs, so that the strong stability properties of these methods have become important. In this work we review sufficient conditions for a two-derivative multistage method to preserve the strong stability properties of spatial discretizations in a forward Euler and different conditions on the second derivative. In particular we present the strong stability preserving theory for explicit and implicit two-derivative Runge–Kutta schemes, including a special condition on the second derivative under which these implicit methods may be unconditionally strong stability preserving. This special condition is natural for the stiff component of wide range of plasma physics problems, and can be useful in the context of strong stability preserving implicit-explicit multi-derivative Runge–Kutta schemes, where the time-step restriction is then independent of the stiff term. Lastly, we present the strong stability preserving theory for implicit-explicit multi-derivative general linear methods, and some novel second and third order methods where the time-step restriction is independent of the stiff term.

97 MATHEMATICS AND COMPUTING↗
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 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↗