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 757 records · Page 42

Ultrastrong colloidal crystal metamaterials engineered with DNA

Lattice-based constructs, often made by additive manufacturing, are attractive for many applications. Typically, such constructs are made from microscale or larger elements; however, smaller nanoscale components can lead to more unusual properties, including greater strength, lighter weight, and unprecedented resiliencies. Here, solid and hollow nanoparticles (nanoframes and nanocages; frame size: ~15 nanometers) were assembled into colloidal crystals using DNA, and their mechanical strengths were studied. Nanosolid, nanocage, and nanoframe lattices with identical crystal symmetries exhibit markedly different specific stiffnesses and strengths. Unexpectedly, the nanoframe lattice is approximately six times stronger than the nanosolid lattice. Nanomechanical experiments, electron microscopy, and finite element analysis show that this property results from the buckling, densification, and size-dependent strain hardening of nanoframe lattices. Last, these unusual open architectures show that lattices with structural elements as small as 15 nanometers can retain a high degree of strength, and as such, they represent target components for making and exploring a variety of miniaturized devices.

42 ENGINEERING↗

Tailored ultrasound propagation in microscale metamaterials via inertia design

The quasi-static properties of micro-architected (meta)materials have been extensively studied over the past decade, but their dynamic responses, especially in acoustic metamaterials with engineered wave propagation behavior, represent a new frontier. However, challenges in miniaturizing and characterizing acoustic metamaterials in high-frequency (megahertz) regimes have hindered progress toward experimentally implementing ultrasonic-wave control. Here, we present an inertia design framework based on positioning microspheres to tune responses of 3D microscale metamaterials. We demonstrate tunable quasi-static stiffness by up to 75% and dynamic longitudinal-wave velocities by up to 25% while maintaining identical material density. Using noncontact laser-based dynamic experiments of tunable elastodynamic properties and numerical demonstrations of spatio-temporal ultrasound wave propagation, we explore the tunable static and elastodynamic property relation. This design framework expands the quasi-static and dynamic metamaterial property space through simple geometric changes, enabling facile design and fabrication of metamaterials for applications in medical ultrasound and analog computing.

42 ENGINEERING↗

Molecular origin of anisotropic shear elastoplasticity in chitin

Chitin nano- and mesoscale structures present in the exoskeleton of crustaceans exhibit exceptional longitudinal stiffness and toughness, rivaling or even exceeding that of many synthetic polymer architectures. Here, we reveal the origin of the asymmetric shear response in chitin multiscale architectures, marked by pronounced anisotropy in deformation. Under shear aligned with the molecular axis, chitin accommodates strain through coherent atomic rearrangements that enable elastic recovery. In contrast, shear applied perpendicular to this axis induces liquid-like plasticity via localized sliding. These results demonstrate the intrinsic mechanical anisotropy of chitin, underpinning a dual function in resisting repetitive loading while dissipating internal stress during high-strain events. Our findings establish the molecular basis of shear elastoplasticity in multiscale chitin structures, wherein axial elasticity supports energy storage in load-bearing regions, whereas transverse plasticity enables controlled energy dissipation. These atomic-scale insights lay a foundation for the predictive design of chitin-based materials with tunable strength-toughness profiles.

Wan, Zhangmin [University of British Columbia, Van↗

Molecular insights into the complex mechanics of plant epidermal cell walls

Computational analysis of cell walls Layers of intertwined fibers make up plant cell walls. The various types of fibers respond differently to deformation. Cellulose microfibrils, for example, can stretch or curve, changing their end-to-end length, and can also slide past each other, reorient relative directions, and bundle with neighboring microfibrils. Zhang et al. developed a computational model based on observations of onion skin epidermis that describes how these complex changes in space govern cell wall mechanics. The results inform ways to engineer multifunctional fibrous materials. Science , this issue p. 706

Science & Technology - Other Topics↗

Asymptotically Compatible Reproducing Kernel Collocation and Meshfree Integration for Nonlocal Diffusion

Reproducing kernel (RK) approximations are meshfree methods that construct shape functions from sets of scattered data. We present an asymptotically compatible (AC) RK collocation method for nonlocal diffusion models with Dirichlet boundary condition. The numerical scheme is shown to be convergent to both nonlocal diffusion and its corresponding local limit as nonlocal interaction vanishes. The analysis is carried out on a special family of rectilinear Cartesian grids for a linear RK method with designed kernel support. The key idea for the stability of the RK collocation scheme is to compare the collocation scheme with the standard Galerkin scheme, which is stable. In addition, assembling the stiffness matrix of the nonlocal problem requires costly computational resources because high-order Gaussian quadrature is necessary to evaluate the integral. We thus provide a remedy to the problem by introducing a quasi-discrete nonlocal diffusion operator for which no numerical quadrature is further needed after applying the RK collocation scheme. The quasi-discrete nonlocal diffusion operator combined with RK collocation is shown to be convergent to the correct local diffusion problem by taking the limits of nonlocal interaction and spatial resolution simultaneously. The theoretical results are then validated with numerical experiments. We additionally illustrate a connection between the proposed technique and an existing optimization based approach based on generalized moving least squares.

97 MATHEMATICS AND COMPUTING↗

Spectral Equivalence of Low-Order Discretizations for High-Order H(curl) and H(div) Spaces

In this paper, we present spectral equivalence results for high-order tensor product edge- and face-based finite elements for the H(curl) and H(div) function spaces. Specifically, we show for certain choices of shape functions that the mass and stiffness matrices of the high-order elements are spectrally equivalent to those for an assembly of low-order elements on the associated Gauss--Lobatto--Legendre mesh. Based on this equivalence, efficient preconditioners can be designed with favorable computational complexity. Numerical results are presented which confirm the theory and demonstrate the benefits of the equivalence results for overlapping Schwarz preconditioners.

97 MATHEMATICS AND COMPUTING↗

High Order Strong Stability Preserving MultiDerivative Implicit and IMEX Runge--Kutta Methods with Asymptotic Preserving Properties

In this article we present a class of high order unconditionally strong stability preserving (SSP) implicit two-derivative Runge--Kutta schemes and SSP implicit-explicit (IMEX) multi-derivative Runge--Kutta schemes where the time-step restriction is independent of the stiff term. The unconditional SSP property for a method of order $p>2$ is unique among SSP methods and depends on a backward-in-time assumption on the derivative of the operator. We show that this backward derivative condition is satisfied in many relevant cases where SSP IMEX schemes are desired. We devise unconditionally SSP implicit Runge--Kutta schemes of order up to $p=4$ and IMEX Runge--Kutta schemes of order up to $p=3$. For the multiderivative IMEX schemes, we also derive and present the order conditions, which have not appeared previously. The unconditional SSP condition ensures that these methods are positivity preserving, and we present sufficient conditions under which such methods are also asymptotic preserving when applied to a range of problems, including a hyperbolic relaxation system, the Broadwell model, and the Bhatnagar--Gross--Krook kinetic equation. We present numerical results to support the theoretical results on a variety of problems.

97 MATHEMATICS AND COMPUTING↗

Impact of symmetry energy on sound speed and spinodal decomposition in dense neutron-rich matter

Using a meta model for nuclear Equation of State (EOS) with its parameters constrained by astrophysical observations and terrestrial nuclear experiments, we examine effects of nuclear EOS especially its symmetry energy $E$ sym ($ρ$) term on the speed of sound squared $C$ $^2_s$($ρ$) and the critical density ρt where $C$ $^2_s$($ρ_t$) vanishes (indicating the onset of spinodal decomposition) in both dense neutron-rich nucleonic matter relevant for relativistic heavy-ion collisions and the cold $n$ + $p$ + $e$ + $μ$ matter in neutron stars at $β$- equilibrium. Unlike in nucleonic matter with fixed values of the isospin asymmetry $δ$, in neutron stars with a density dependent isospin profile $δ$($ρ$) determined consistently by the $β$-equilibrium and charge neutrality conditions, the $C$ $^2_s$($ρ$) almost always show a peak and then vanishes at $ρ_t$. Furthermore, the latter strongly depends on the high-density behavior of $E$ sym ($ρ$) if the skewness parameter $J$ 0 characterizing the stiffness of high-density symmetric nuclear matter (SNM) EOS is not too far above its currently known most probable value of about -190 MeV inferred from recent Bayesian analyses of neutron star observables. Moreover, in the case of having a super-soft $E$ sym ($ρ$) that is decreasing with increasing density above about twice the saturation density of nuclear matter, the $ρ$ $t$ is significantly lower than the density where the $E$ sym ($ρ$) vanishes (indicating the onset of isospin-separation instability and pure neutron matter formation) in neutron star cores.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Role of humidity and surface roughness on direct wafer bonding

Bodies made from elastically stiff material usually bind very weakly unless the surfaces are flat and extremely smooth. In direct wafer bonding flat surfaces bind by capillary bridges and by the van der Waals interaction, which act between all solid objects. Here we study the dependency of the work of adhesion on the humidity and surface roughness in hydrophilic direct wafer bonding. We show that the long-wavelength roughness (usually denoted waviness) has a negligible influence on the strength of wafer bonding (the work of adhesion) from the menisci that form from capillary condensation of water vapor.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Mechano-Electrochemical Interaction in Solid-State Lithium Batteries

A solid-state lithium (Li) battery primarily consists of Li metal anode, solid electrolyte separator, and cathode. The asymmetric volume changes, originating from ion transport and interfacial Li growth during plating, lead to significant stresses in the layered architecture. In this study, we develop a coupled mechanics-electrochemistry formalism for polymer electrolyte based solid-state batteries, in particular, focusing on the stress effect on electrochemical performance. By means of a coupling coefficient, it is found that stress-assisted ion transport in the electrolyte results in a delayed Sand’s time and increased critical current density of unstable electrodeposition, and consequently, alleviates the propensity of dendrite formation. Stress at the Li metal-electrolyte interface affects the electrochemical reaction kinetics, and the influences from the deviatoric stress and hydrostatic pressure vary with Li plating time. In addition, a low restraint stiffness to the layered structure could elastically buffer the volumetric changes and thus reduce the stress during Li plating. This fundamental study provides guidance for the design of solid-state batteries, aimed at stable electrodeposition and mechanical integrity.

25 ENERGY STORAGE↗

Orthotropic Viscoelastic Modeling of Polymeric Battery Separator

The separator is one of the most important components for battery safety since its failure can cause short-circuit that may lead to thermal runaway. As such, understanding the mechanical behaviors of separators is critical to improve battery safety. Nevertheless, modeling the mechanical behaviors of separators is challenging since they usually depend on multiple factors such as orientation, temperature and strain rate due to the microstructure of separators. Here we develop an orthotropic linear viscoelastic model for separators, which takes both anisotropy and strain rate effects into account. This model is the first step of a future model that considers more factors such as temperature and large deformation. The model has been implemented as a user material model in a finite element package, where a discretization algorithm is developed to evaluate the stiffness-based hereditary integral with a kernel of Prony series. The user model has been applied to a polypropylene (PP) separator Celgard ® 2400. The simulation results agree well with experimental data, which include the rate dependent uniaxial tensile behavior of the separator in the two in-plane material directions and in shear, and the strain contour patterns of every strain component during biaxial tension.

25 ENERGY STORAGE↗

Electrochemical-Mechanical Coupling Strongly Affects the Performance of Nanopore, Thin-Film, and Solid-State Batteries

Here, the present study focuses on the electrochemical-mechanical (ECM) coupling effects of a thin-film, solid-state battery with only stiff, ceramic materials, in contrast to prior investigations that focus on individual active material particles or aggregations of particles. We model the impacts of ECM couplings including stress-transport and stress-equilibrium potential on the full-cell performance and potential mechanical failure modes of a thin-film battery conformally deposited in a nanopore scaffold, which is an experimentally achievable device. Model results indicate electrode volume changes due to lithium insertion or removal, along with mechanical boundary conditions, result in stress gradients that alter the lithium-ion flux, reduce lithium concentration gradients, and improve cell rate capability. However, the high stress levels in the cell can also lead to mechanics-related failure such as the separation of cell layers. For the parameter set in this work, stress-transport coupling has a much greater influence on rate capability than stress-potential coupling. Optimization of thin-film batteries to harness the benefits of ECM coupling effects requires leveraging geometric design and material selection. The current work underscores the need for further theoretical and experimental investigation into ECM coupling effects in thin-film batteries to enhance their understanding and design optimization.

25 ENERGY STORAGE↗

Dynamic Response of a Semiactive Suspension System with Hysteretic Nonlinear Energy Sink Based on Random Excitation by means of Computer Simulation

This paper aims to investigate the property and behavior of the hysteretic nonlinear energy sink (HNES) coupled to a half vehicle system which is a nine-degree-of-freedom, nonlinear, and semiactive suspension system in order to improve the ride comfort and increase the stability in shock mitigation by using the computer simulation method. The HNES model is a semiactive suspension device, which comprises the famous Bouc–Wen (B-W) model employed to describe the force produced by both the purely hysteretic spring and linear elastic spring of potentially negative stiffness connected in parallel, for the half vehicle system. Nine nonlinear motion equations of the half vehicle system are derived in terms of the seven displacements and the two dimensionless hysteretic variables, which are integrated numerically by employing the direct time integration method for studying both the variables of vertical displacements, velocities, accelerations, chassis pitch angle, and the ride comfort and driver safety, respectively, based on the bump and random road inputs of the pseudoexcitation method as excitation signal. Simulation results show that, compared with the HNES model and the magnetorheological (MR) model coupled to the half vehicle system, the ride comfort and stability have been evidently improved. A successful validation process has been performed, which indicated that both the ride comfort and driver safety properties of the HNES model coupled to half vehicle significantly improved.

Chen, Hui↗

MoorPy (Quasi-Static Mooring Analysis in Python)

MoorPy is a quasi-static mooring model and a suite of associated functions for mooring system analysis. The core model supports quasi-static analysis of moored floating systems including any arrangement of mooring lines and floating platforms. It solves the distributed position and tension of each mooring line segment using standard catenary equations. Floating platforms can be represented with linear hydrostatic characteristics. MoorPy automatically computes a floating system's equilibrium state and can be queried to identify a mooring system's nonlinear force-displacement relationships. Linearized stiffness matrices are efficiently computed using semi-analytic Jacobians. MoorPy also includes plotting functions and a library of mooring component property and cost coefficients. MoorPy can be used directly from Python scripts to perform mooring design and analysis tasks, or it can be coupled with other tools to compute quasi-static mooring reactions as part of a larger simulation.

Hall, Mathew↗

Experimental investigation of the crush performance of prepreg platelet molding compound tubes

This work presents an experimental study characterizing the crush performance of hollow cylindrical tubes made with prepreg platelet molding compound (PPMC). PPMC is a composite material system that uses platelets from chopped and slit unidirectional prepreg as the basis for a molding compound. This material system has a higher fiber volume fraction and better mechanical properties than traditional short fiber systems and can be molded into complex geometries unlike continuous fiber systems. As such, this material system shows promise for use with complex structural members in vehicles. The failure morphology and specific energy absorption of the material are evaluated with different thickness-to-diameter geometries and test speeds. In addition, this work investigates how PPMC components compare to traditional continuous fiber components and finds that PPMC performs as well as continuous fiber layups with similar effective laminate stiffnesses.

Materials Science↗

Analyzing photovoltaic module mechanics using composite plate theories and finite element solutions

Deflection and stress calculated from an experimentally validated, high-fidelity finite element model (FEM) of a photovoltaic module experiencing mechanical load was compared to results from a simplified FEM treating the module laminate as a homogenized composite using a rule of mixtures approach, and further compared to analytical calculations treating the module as a Kirchoff-Love flat plate. The goal of this study was to determine the error incurred by analyzing module mechanics with varying levels of simplification, since resolving the aspect ratios of a module is computationally expensive. Homogenized FEMs were found to underpredict peak deflection under a 1.0 kPa load by between 13 and 19% for lower and upper bound application of the rule of mixtures. However, module shape was captured, implying that a useful replication of a resolved model could be achieved with a reduced, calibrated material stiffness. Homogenized stress results captured glass layer tensile stress components to within 46 to 52% at a sample location of interest, though agreement was poor through the remainder of the laminate due to the lack of material resolution. For plate theory, deflection was overpredicted by 45 to 67% for upper and lower bound homogenizations, and frame-adjacent module shapes were not adequately replicated. Stress results mirrored FEM trends but magnitudes were not well correlated to resolved model values. These results support the use of homogenized laminate models for module shape derivation, though resolved models remain necessary for stress analyses. The accuracy of plate theory was found to be inadequate for most applications.

14 SOLAR ENERGY↗

A GPU-based compressible combustion solver for applications exhibiting disparate space and time scales

High-speed chemically active flows pose significant computational challenges due to their disparate space and time scales, with stiff chemistry often dominating simulation time. While modern scientific computing programs achieve exascale performance by leveraging graphics processing units (GPUs), existing GPU-based compressible combustion solvers face critical limitations in memory management, load balancing, and handling the highly localized nature of chemical reactions. To this end, we present a high-performance compressible reacting flow solver built on the AMReX framework and optimized for multi-GPU settings. Here, our approach addresses three GPU performance bottlenecks: memory access patterns through column-major storage optimization, computational workload variability via a bulk-sparse integration strategy for chemical kinetics, and multi-GPU load distribution for adaptive mesh refinement applications. The solver adapts existing matrix-based chemical kinetics formulations to multi-grid contexts. Using representative combustion applications, including 2D and 3D detonations and a 3D jet-in-crossflow configuration, we demonstrate 1.4–5× performance improvements over initial implementations on an in-house cluster of NVIDIA H100 GPUs, and near-ideal weak scaling on the Frontier supercomputer (Oak Ridge Leadership Computing Facility) with up to 1024 AMD Instinct MI250X GPUs. Roofline analysis reveals substantial improvements in arithmetic intensity for both convection (∼ 10 ×) and chemistry (∼ 4 ×) routines, confirming efficient utilization of GPU memory bandwidth and computational resources.

42 ENGINEERING↗

Characterization of micro-sandwich structures via direct ink writing epoxy based cores

Sandwich structured (SS) composites demonstrate considerable flexural stiffness and high strength-to-weight ratios and can be tailored as functional materials. Historically they have been constrained to specific material types and geometry due to limitations in manufacturing methods. However, employing additive manufacturing (AM), specifically direct ink writing (DIW), can provide an alternative method for making SS composites with complex and controllable micro and mesostructures with multifunctionality targeted at desired mechanical, thermal, and electrical properties. DIW, an extrusion-based AM technique, uses a viscous and thixotropic ink with desired components that, once printed, is cured to obtain the final complex net shape parts. In this paper, a novel hybrid AM technique is employed to manufacture SS composite materials containing bisphenol A-based epoxy core and carbon fiber reinforced polymer (CFRP) face sheets that are fabricated via DIW and vacuum infusion process (VIP), respectfully. We demonstrate that the fabrication of these SS composites can be tailored from a thermosetting material, from which additives and/or various lattice structures can be manufactured to achieve enhanced and desirable mechanical integrity with functional properties. Surface topology and mechanical testing techniques are used to characterize the fabricated hybrid SS composites to study and assess mechanical stability. A rheo-kinetic cure model was developed for the core material to allow for additive manufacturing process requirements while ensuring complete cross-linking for the thermoset-based core material. Because of the ability to obtain relatively small core-thickness and controlled architecture, this method now allows for fabricating layered micro-sandwich structures for realizing further light-weighting in relevant applications.

36 MATERIALS SCIENCE↗