Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “elastic materials”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Nonlinear response of graphite-epoxy wide columns subject to eccentric load

It is shown that including the effects of bend-twist coupling and material nonlinearity on the large displacement flexural response of graphite-epoxy wide columns is necessary to adequately model their response. Results from static tests on 0(30)T and (30/0/-30 sub 5)s laminated specimens under eccentric axial compression are compared to a nonlinear analysis. The analysis is based on a Kantrovich expansion of the displacements in the width coordinate and a finite element representation in the axial coordinate. An incremental, updated Lagrangian formulation and element corotational coordinates are employed. Nonlinear elastic material behavior was measured in a separate set of material characterization tests, and these data are used to compute tangent stillnesses in the updated configuration.

Griffin, O. H., Jr.↗
NASA NTRS1994

Effects of assumed tow architecture on the predicted moduli and stresses in woven composites

This study deals with the effect of assumed tow architecture on the elastic material properties and stress distributions of plain weave woven composites. Specifically, the examination of how a cross-section is assumed to sweep-out the tows of the composite is examined in great detail. The two methods studied are extrusion and translation. This effect is also examined to determine how sensitive this assumption is to changes in waviness ratio. 3D finite elements were used to study a T300/Epoxy plain weave composite with symmetrically stacked mats. 1/32nd of the unit cell is shown to be adequate for analysis of this type of configuration with the appropriate set of boundary conditions. At low waviness, results indicate that for prediction of elastic properties, either method is adequate. At high waviness, certain elastic properties become more sensitive to the method used. Stress distributions at high waviness ratio are shown to vary greatly depending on the type of loading applied. At low waviness, both methods produce similar results.

Chapman, Clinton Dane↗
DOE OSTI2023

Thermo-acoustoelastic effect of Rayleigh wave: Theory and experimental verification

Previous studies showed that the thermally-induced ultrasonic bulk wave velocity change could be used to measure acoustoelastic coefficients and third-order elastic constants of elastic materials. This method is naturally immune from the ambient temperature effect and has improved sensitivity and a simpler test setup than the conventional acoustoelastic test. However, Rayleigh wave is preferred for thick components or structures with only one accessible surface. In this work, the thermo-hyperelastic constitutive equation, along with acoustoelastic theory, is used to derive the expression of the thermo-acoustoelastic coefficient (TAEC) of Rayleigh wave. The numerical relationship between the TAEC of Rayleigh wave and Murnaghan constants ( l, m and n ) are given for common metals. Further, the TAEC expressions for Rayleigh wave and shear wave are similar, and both are dominated by the constant m . The TAEC of Rayleigh wave was measured on an aluminum 6061 specimen using the thermal modulation experiment in a temperature range of 22 ~ 35°C. The measured TAEC shows good agreement with the theoretical calculation. Then the third-order elastic constants were calculated based on TAECs of bulk waves and Rayleigh wave.

42 ENGINEERING↗
DOE OSTI2022

A fourth-order phase-field fracture model: Formulation and numerical solution using a continuous/discontinuous Galerkin method

Modeling crack initiation and propagation in brittle materials is of great importance to be able to predict sudden loss of load-carrying capacity and prevent catastrophic failure under severe dynamic loading conditions. Second-order phase-field fracture models have gained wide adoption given their ability to capture the formation of complex fracture patterns, e.g. via crack merging and branching, and their suitability for implementation within the context of the conventional finite element method. Higher-order phase-field models have also been proposed to increase the regularity of the exact solution and thus increase the spatial convergence rate of its numerical approximation. However, they require special numerical techniques to enforce the necessary continuity of the phase field solution. In this paper, we derive a fourth-order phase-field model of fracture in two independent ways; namely, from Hamilton’s principle and from a higher-order micromechanics-based approach. The latter approach is novel, and provides a physical interpretation of the higher-order terms in the model. In addition, we propose a continuous/discontinuous Galerkin (C/DG) method for use in computing the approximate phase-field solution. This method employs Lagrange polynomial shape functions to guarantee -continuity of the solution at inter-element boundaries, and enforces the required regularity with the aid of additional variational and interior penalty terms in the weak form. Finally, the phase-field equation is coupled with the momentum balance equation to model dynamic fracture problems in hyper-elastic materials. Two benchmark problems are presented to compare the numerical behavior of the C/DG method with mixed finite element methods.

42 ENGINEERING↗
DOE OSTI2020

Experimental and numerical studies of seismic fluid-structure interaction in a base-supported cylindrical vessel

Seismic design and qualification of advanced reactors will rely heavily on the use of verified and validated numerical models capable of capturing the interaction of the vessel, its contained fluid, and the internal equipment: fluid-structure interaction (FSI) analysis. Analytical solutions can be used for preliminary sizing and design of such vessels but their application is limited to simple geometries and boundary conditions, and small amplitude, translational (and rotational) inputs. To validate numerical models for seismic FSI analysis in finite element codes, a comprehensive set of experiments was performed on a liquid-filled cylindrical vessel, using a 6 degree-of-freedom earthquake simulator. Results in terms of sloshing frequency, damping ratio in sloshing modes, and hydrodynamic responses (wave height, hydrodynamic pressure, base shear, and base moment) for multi-directional earthquake simulator inputs are reported and compared with analytical solutions for liquid-filled vessels. The impact of seismic (base) isolation on hydrodynamic responses was studied using earthquake simulator inputs generated using a virtual isolation system. Data from the experiments are used to validate a numerical model of the fluid-filled vessel using the Arbitrary Lagrangian Eulerian (ALE) solver in the commercial finite element program LS-DYNA. Validation studies are presented for multi-directional seismic inputs, including rocking motions. Lagrangian modeling approaches using an elastic material formulation for the fluid are also investigated and their limitations and possible applications are identified. Here, the results are broadly applicable to the seismic response of base supported, liquid-filled vessels.

58 GEOSCIENCES↗
DOE OSTI2020

An anisotropic viscoplasticity model for shale based on layered microstructure homogenization

Viscoplastic deformation of shale is frequently observed in many subsurface applications. Many studies have suggested that this viscoplastic behavior is anisotropic—specifically, transversely isotropic—and closely linked to the layered composite structure at the microscale. In this work, we develop a two-scale constitutive model for shale in which anisotropic viscoplastic behavior naturally emerges from semianalytical homogenization of a bilayer microstructure. The microstructure is modeled as a composite of soft layers, representing a ductile matrix formed by clay and organics, and hard layers, corresponding to a brittle matrix composed of stiff minerals. This layered microstructure renders the macroscopic behavior anisotropic, even when the individual layers are modeled with isotropic constitutive laws. Using a common correlation between clay and organic content and magnitude of creep, we apply a viscoplastic modified Cam-Clay plasticity model to the soft layers, while treating the hard layers as a linear elastic material to minimize the number of calibration parameters. We then describe the implementation of the proposed model in a standard material update subroutine. The model is validated with laboratory creep data on samples from three gas shale formations. We also demonstrate the computational behavior of the proposed model through simulation of time-dependent borehole closure in a shale formation with different bedding plane directions.

58 GEOSCIENCES↗
DOE OSTI2026

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗
DOE OSTI2026

Propagation of elastic waves in a flexomagnetic solid

Flexomagnetism is the coupling between magnetism and strain gradients and is a technologically relevant phenomenon. We present a theory of elastic wave propagation in a linear elastic flexomagnetic material with microstructure and strain gradient elastic interactions. The expressions of frequency, phase velocity, and group velocity of longitudinal and transverse waves are derived and are shown to depend on the flexomagnetic coefficient and microstructure. We also show that the effect of flexomagnetism and microstructure can lead to some interesting phenomena in wave propagation, which are not observed in classical linear elasticity theory of waves. Specifically, in contrast to classical linear elastic materials, where wave propagation is non-dispersive, flexomagnetic materials with microstructure can exhibit both normal and abnormal dispersion. It is also noteworthy that, in flexomagnetic materials with gradient elasticity, the phase velocities of transverse waves can exceed those of longitudinal waves, which is atypical in classical elasticity. Furthermore, waves can also attenuate for a certain range of wavenumbers that depend on the flexomagnetic coefficient and microstructural parameters. Finally, we explore the possibility of waves exhibiting zero group velocity modes, where waves are non-propagating but have strong local energy confinement, negative group velocity modes, where the wave packet moves in the opposite direction to that of wave propagation, and the phenomenon of wave freezing, where a propagating wave stops in space without diffusing or spreading.

Ghosh, Swarnava [ORNL] (ORCID:0000000338005264)↗
DOE OSTI2020

Exploiting Nonlinear Elasticity for Anomalous Magnetoresponsive Stiffening

Here, a paradigm for enhanced magnetorheological elastic materials is introduced and experimentally established. We show that a nonlinearly stiffening polymer matrix can be exploited to achieve anomalous magneto-elastomer stiffening exceeding standard magneto-elastomer theory and experiment in terms of percentage stiffness change and sensitivity to applied magnetic flux. Using a model system of a semifexible fibrin network embedded with micron sized carbonyl iron particles, we demonstrate that even at a modest particle volume fraction (0.5-4%), a coupling between the magnetically-interacting dipoles and a strain-stiffening polymer mesh provides previously unexplored opportunities for material design. Our experiments indicate that confined particles within thefibrin network internally tension and stiffen the polymer mesh when an external field is applied, resulting in a field-dependent stiffening response from the polymer mesh that superposes with the magnetic interparticle interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Vector elastic deconvolution migration with dual wavefield decomposition

Here, we report elastic-wave imaging using multicomponent data can provide more useful subsurface information than acoustic-wave imaging, but it is usually algorithmically challenging. We have developed a vector elastic deconvolution migration method for high-resolution imaging of subsurface structures in isotropic and anisotropic elastic media. Our new method uses a vector deconvolution imaging condition based on dual wavefield decomposition, including an explicit directional wavefield separation using the Hilbert transform and a P/S vector wavefield decomposition using the low-rank decomposition method. Using three elastic models, we numerically determine that our new method produces notably higher resolution and more amplitude-balanced elastic images compared with a crosscorrelation-based vector elastic reverse time migration method.

58 GEOSCIENCES↗
DOE OSTI2025

Influence of gradation on failure mode of saturated sand at particle scale

Extensive research has been reported in the literature to characterize the failure mode of dry sand using various experimental techniques such as surface optical imaging, photo-elastic materials, three-dimensional (3D) computed tomography (CT), and 3D synchrotron micro-computed tomography (SMT). However, there is a limited literature about the behavior of saturated sand. This paper presents the results of axisymmetric triaxial compression (ATC) experiments that were conducted on saturated sand specimens. The behavior of specimens composed of a uniform sand with grain size between US sieves #40 and #50 is compared to specimens conducted on the same sand that has a wider gradation. 3D SMT technique was used to acquire 3D scans while shearing the specimens to probe localized events that are completely missed or misinterpreted when analyzing ATC measurements based on global standard measurements. The results show a higher effective principal stress ratio (EPSR) for the non-uniform specimen and a thicker shear band when compared to uniform specimen.

Elnur, Mohammed [University of Tennessee]↗
DOE OSTI2022

The Schwarz Alternating Method for the Seamless Coupling of Nonlinear Reduced Order Models and Full Order Models

Projection-based model order reduction allows for the parsimonious representation of full order models (FOMs), typically obtained through the discretization of a set of partial differential equations (PDEs) using conventional techniques (e.g., finite element, finite volume, finite difference methods) where the discretization may contain a very large number of degrees of freedom. As a result of this more compact representation, the resulting projection-based reduced order models (ROMs) can achieve considerable computational speedups, which are especially useful in real-time or multi-query analyses. One known deficiency of projection-based ROMs is that they can suffer from a lack of robustness, stability and accuracy, especially in the predictive regime, which ultimately limits their useful application. Another research gap that has prevented the widespread adoption of ROMs within the modeling and simulation community is the lack of theoretical and algorithmic foundations necessary for the “plug-and-play” integration of these models into existing multi-scale and multi-physics frameworks. This paper describes a new methodology that has the potential to address both of the aforementioned deficiencies by coupling projection-based ROMs with each other as well as with conventional FOMs by means of the Schwarz alternating method [41]. Leveraging recent work that adapted the Schwarz alternating method to enable consistent and concurrent multiscale coupling of finite element FOMs in solid mechanics [35, 36], we present a new extension of the Schwarz framework that enables FOM-ROM and ROM-ROM coupling, following a domain decomposition of the physical geometry on which a PDE is posed. In order to maintain efficiency and achieve computation speed-ups, we employ hyper-reduction via the Energy-Conserving Sampling and Weighting (ECSW) approach [13]. We evaluate the proposed coupling approach in the reproductive as well as in the predictive regime on a canonical test case that involves the dynamic propagation of a traveling wave in a nonlinear hyper-elastic material.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Ultrasonic waveguide reflectometry for quench detection (CRADA Final Report)

Early detection of thermal runaway (quench) is critical for safely operating high-temperature superconductor (HTS) magnets. Voltage-based techniques provide late detection of quenching in HTS, which can lead to magnet burnout and system damage. Various non-voltage quench-detection approaches are therefore investigated nowadays that are based on magnetic, acoustic, and optical sensing. LBNL has been actively developing the quench detection principle based on diffuse ultrasonic wave propagation since 2017 and has demonstrated the approach's viability with several peer-reviewed publications and a patent application. Over the last decade, Etegent, LBNL partner for this proposal and effort, has been developing a sensing technology that utilizes the propagation of stress waves through the solid metallic and nonmetallic waveguides (WG). The method is based on the dependence of materials' elastic properties and the speed of sound on temperature. Compared with optical analog, acoustic WG operates at much lower frequencies, eliminating the need for expensive, complex, highly sensitive receivers and data processing. Etegent proved this technology in many harsh environments, mainly concentrated high-temperature measurements. The initial contact between LBNL and Etegent started in 2020, aiming at combining our knowledge and effort toward developing a viable ultrasound-based quench detection technique for use with HTS magnets. A preliminary study conducted jointly in 2022 found no principle restrictions for expanding the WG technology toward cryogenic temperature operation. The present CRADA work aimed at designing, building, and testing ultrasonic waveguide temperature sensors at relevant cryogenic conditions, and demonstrating this technology's applicability to detecting heating and quenching in High Temperature Superconductor (HTS)-based cables and magnets. The technical objectives of the project included: (1) Design and fabricate ultrasonic WG sensors suitable for cryogenic operation, (2) Test ultrasonic WGs under cryogenic conditions, (3) Determining a suitable way of insulating WGs to prevent ultrasonic leakage while maintaining good thermal contact with the surrounding medium, (4) WG integration and cryogenic test with HTS conductors and HTS magnets, and (5) Design and basic prototyping of a standalone WG-based quench detection system. This project adds value to the present LBNL magnet diagnostics portfolio and enables a breakthrough in addressing the HTS magnet quench detection challenge.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Predicting Mechanical Properties from Microstructure Images in Fiber-Reinforced Polymers Using Convolutional Neural Networks

Evaluating the mechanical response of fiber-reinforced composites can be extremely time-consuming and expensive. Machine learning (ML) techniques offer a means for faster predictions via models trained on existing input–output pairs and have exhibited success in composite research. This paper explores a fully convolutional neural network modified from StressNet, which was originally used for linear elastic materials, and extended here for a non-linear finite element (FE) simulation to predict the stress field in 2D slices of segmented tomography images of a fiber-reinforced polymer specimen. The network was trained and evaluated on data generated from the FE simulations of the exact microstructure. The testing results show that the trained network accurately captures the characteristics of the stress distribution, especially on fibers, solely from the segmented microstructure images. The trained model can make predictions within seconds in a single forward pass on an ordinary laptop, given the input microstructure, compared to 92.5 h to run the full FE simulation on a high-performance computing cluster. These results show promise in using ML techniques to conduct fast structural analysis for fiber-reinforced composites and suggest a corollary that the trained model can be used to identify the location of potential damage sites in fiber-reinforced polymers.

Sun, Yixuan (ORCID:0000000311093380)↗
DOE OSTI2023

A guide to the design of the virtual element methods for second- and fourth-order partial differential equations

Here we discuss the design and implementation details of two conforming virtual element methods for the numerical approximation of two partial differential equations that emerge in phase-field modeling of fracture propagation in elastic material. The two partial differential equations are: (i) a linear hyperbolic equation describing the momentum balance and (ii) a fourth-order elliptic equation modeling the damage of the material. Inspired by, we develop a new conforming VEM for the discretization of the two equations, which is implementation-friendly, i.e., different terms can be implemented by exploiting a single projection operator. We use C 0 and C 1 virtual elements for the second-and fourth-order partial differential equation, respectively. For both equations, we review the formulation of the virtual element approximation and discuss the details pertaining the implementation.

42 ENGINEERING↗