Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Elasticity”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Mesoscale Linear Elastic Modeling and Homogenization of Marine Energy Composites

The design of fiber-reinforced composite (FRC)-based components for marine energy applications necessitates a fundamental understanding of material properties and the resulting geometry to predict long-term performance. In this work, we present a modeling workflow to predict linear elastic and diffusive bulk properties at the mesoscale for an idealized geometry based on knowledge of fiber and resin properties. A parametric study was performed to identify the key model input parameters that influence bulk properties. Furthermore, we demonstrate how bulk properties can be leveraged in high-fidelity image-based simulations, where imperfections in tow geometry and voids captured during X-ray computed tomography imaging are explicitly represented within the simulation. Bulk properties of interest include moduli, Poisson’s ratios, hygroscopic swelling, diffusivity, and moisture uptake, which are key parameters for characterizing FRC performance within marine environments. Modeling predictions agreed well with experimental data, except for estimating swelling coefficients, likely due to crack accumulation as a function of moisture uptake. The mesoscale modeling workflow ultimately highlights a versatile framework for understanding the influence of material and geometric properties, which can be leveraged to rapidly assess new FRC-based components.

computational mechanics↗
DOE OSTI2025

The Elastic Analysis Facility's (EAF's) Contribution to the Future of Analysis at Multi-Experiment Institutions and Future Colliders

The Elastic Analysis Facility (EAF) hosted at Fermi National Accelerator Laboratory (Fermilab) is a platform being developed with the goal of providing a fast and efficient facility for physics analysis. As high-energy physics moves towards collecting larger datasets, such as those from the High-Luminosity LHC, the EAF strives to provide a powerful and adaptable framework for future colliders and multi-experiment institutions. Currently, the EAF supports several experiments including CMS, NOvA, and DUNE as well as serving accelerator physicists and beam line operations through integrated software and secure connections to Fermilab's computing resources. In addition, the EAF was designed with a user-friendly interface, intended to be more intuitive for emerging generations of physicists, that is still accessible for established styles of analysis. The EAF can also achieve better analysis efficiency due to the modernization of software and tools that can better utilize Fermilab's computing power. Furthermore, its design incorporates industry standards whenever possible, enhancing its sustainability and making it a possible template for other national or international laboratories and research facilities. Overall, the EAF is a forward-looking solution that will meet the evolving needs of particle physics, ensuring readiness for future colliders and multi-experiment research institutions.

Chavez, Elise [Wisconsin U., Madison]↗
DOE OSTI2024

Refactoring the elastic–viscous–plastic solver from the sea ice model CICE v6.5.1 for improved performance

This study focuses on the performance of the elastic–viscous–plastic (EVP) dynamical solver within the sea ice model, CICE v6.5.1. The study has been conducted in two steps. First, the standard EVP solver was extracted from CICE for experiments with refactored versions, which are used for performance testing. Second, one refactored version was integrated and tested in the full CICE model to demonstrate that the new algorithms do not significantly impact the physical results. The study reveals two dominant bottlenecks, namely (1) the number of Message Parsing Interface (MPI) and Open Multi-Processing (OpenMP) synchronization points required for halo exchanges during each time step combined with the irregular domain of active sea ice points and (2) the lack of single-instruction, multiple-data (SIMD) code generation. The standard EVP solver has been refactored based on two generic patterns. The first pattern exposes how general finite differences on masked multi-dimensional arrays can be expressed in order to produce significantly better code generation by changing the memory access pattern from random access to direct access. The second pattern takes an alternative approach to handle static grid properties. The measured single-core performance improvement is more than a factor of 5 compared to the standard implementation. The refactored implementation of strong scales on the Intel® Xeon® Scalable Processors series node until the available bandwidth of the node is used. For the Intel® Xeon® CPU Max series, there is sufficient bandwidth to allow the strong scaling to continue for all the cores on the node, resulting in a single-node improvement factor of 35 over the standard implementation. This study also demonstrates improved performance on GPU processors.

58 GEOSCIENCES↗
DOE OSTI2026

Measurement of the Neutron Elastic Electromagnetic Form Factor Ratio at Large Momentum Transfer

Exploring nucleon structure is vital both for understanding its origin and existence as well as for the advancement of the sciences. It helps us answer key questions such as how quark and gluon dynamics create 99% of the nucleon mass. Electron- nucleon scattering has been widely used for precision studies of the nucleon and nuclear structure since the Nobel Prize winning investigations by Robert Hofstadter and collaborators in the 1950s. These studies provide information about the spatial charge and current densities of the nucleon in terms of the electromagnetic form factors. The form factors are functions of four momentum transfer squared (Q2). Extending the electromagnetic form factor measurements to higher Q2 plays a critical role in furthering the understanding of nucleon structure. This motivated the Super BigBite Spectrometer (SBS) program at Jefferson Lab. The open nature of the spectrometers and the direct line of sight from the target to the tracking detector locations in experimental setups such as SBS creates high levels of background at the detectors. This necessitates the use of tracking detectors with high rate capability and good position resolution. Gas Electron Multiplier (GEM) detectors are an excellent choice for tracking detectors in such experiments. Understanding the performance of the GEM detectors is important not just for SBS experiments but also for future high-luminosity experiments. This thesis reports the exploratory results from the measurement of the neutron elastic electromagnetic form factor ratio (Gn E/Gn M) at high momentum transfer. A longitudinally polarized electron beam was scattered off a polarized 3He target, used as an effective polarized neutron target. In this experiment, the polarized 3He target achieved a world record polarization weighted luminosity at a beam current of 45 µA. Double spin asymmetry of the scattered neutron events is used to extract the neutron form factor ratio. Measurements were taken at Q2 = 3.0, 6.8, 9.8 (GeV/c)2. The lowest Q2 measurement is in good agreement with the existing world data, and the higher-Q2 measurements extend the Q2 reach well beyond the existing world data and are expected to remain unmatched for a long time.

Gamage, Vimukthi Haththotuwa [Univ. of Virginia, C↗
DOE OSTI2026

Measurement of the $\bar ν_μ-$Hydrogen Charged-Current Quasi-Elastic Cross Section using the NOvA Near Detector

We report a measurement of the total cross section for muon antineutrino charged-current quasi-elastic scattering on hydrogen, $\bar ν_μ{\rm H} \to μ^+ n$, in the NOvA near detector using a $1.2\times10^{21}$ proton-on-target exposure in the NuMI beam. A selection based on topological and kinematic constraints yields 35,509 signal events in the hydrogen-rich ($10.8\%$) detector, providing the highest statistics of (anti)neutrino--hydrogen interactions measured to date. Backgrounds from (anti)neutrino interactions on heavier nuclei are constrained using dedicated data control samples, significantly reducing the related systematic uncertainties. We obtain a value $σ(\bar ν_μ{\rm H} \to μ^+ n) = 0.538 \pm 0.009 ({\rm stat}) \pm 0.010 ({\rm syst}) \pm0.055 ({\rm flux}) \times 10^{-38}$ cm$^2$ for the total cross section at an average energy of 1.9 GeV, the most precise total cross-section measurement of this process to date. The combined statistical and non-flux systematic uncertainty is more than four times smaller than the flux uncertainty, allowing a future use of this measurement to constrain the absolute $\bar ν_μ$ flux.

Abubakar, S. [Erciyes U.]↗
DOE OSTI2026

Measurement of the $\bar ν_μ-$Hydrogen Charged-Current Quasi-Elastic Cross Section using the NOvA Near Detector

We report a measurement of the total cross section for muon antineutrino charged-current quasi-elastic scattering on hydrogen, $\bar ν_μ{\rm H} \to μ^+ n$, in the NOvA near detector using a $1.2\times10^{21}$ proton-on-target exposure in the NuMI beam. A selection based on topological and kinematic constraints yields 35,509 signal events in the hydrogen-rich ($10.8\%$) detector, providing the highest statistics of (anti)neutrino--hydrogen interactions measured to date. Backgrounds from (anti)neutrino interactions on heavier nuclei are constrained using dedicated data control samples, significantly reducing the related systematic uncertainties. We obtain a value $σ(\bar ν_μ{\rm H} \to μ^+ n) = 0.538 \pm 0.009 ({\rm stat}) \pm 0.010 ({\rm syst}) \pm0.055 ({\rm flux}) \times 10^{-38}$ cm$^2$ for the total cross section at an average energy of 1.9 GeV, the most precise total cross-section measurement of this process to date. The combined statistical and non-flux systematic uncertainty is more than four times smaller than the flux uncertainty, allowing a future use of this measurement to constrain the absolute $\bar ν_μ$ flux.

Abubakar, S. [Erciyes U.]↗
DOE OSTI2025

Coupling Approaches with Non-matching Grids for Classical Linear Elasticity and Bond-based Peridynamic Models in 1D

Local-nonlocal coupling approaches provide a means to combine the computational efficiency of local models and the accuracy of nonlocal models. To facilitate the coupling of the two models, non-matching grids are often desirable as nonlocal grids usually require a finer resolution than local grids. In that case, it is often convenient to resort to interpolation operators so that models can exchange information in the overlap regions when nodes from the two grids do not coincide. This paper studies three existing coupling approaches, namely 1) a method that enforces matching displacements in an overlap region, 2) a variant that enforces a constraint on the stresses instead, and 3) a method that considers a variable horizon in the vicinity of the interfaces. Further, the effect of the interpolation order and of the grid ratio on the performance of the three coupling methods with non-matching grids is carefully studied on one-dimensional examples using polynomial manufactured solutions. The numerical results show that the degree of the interpolants should be chosen with care to avoid introducing additional modeling errors, or simply minimize these errors, in the coupling approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

High-pressure elasticity and equation of state of the fluoroelastomer Viton® A-500

Viton® A is a semi-crystalline copolymer of polyvinylidene fluoride and hexafluoropropylene used in various engineering applications due to its mechanical properties and chemical inertness. In situ ultrasonic spectroscopy and x-ray radiography measurements were performed in a Paris–Edinburgh press to measure the pressure dependence of the transverse and longitudinal acoustic velocity of the fluoroelastomer A-500 from 2.7 to 5.7 GPa at 296 K. In addition, we performed high-pressure Brillouin scattering measurements to obtain acoustic velocities from ambient pressure to 5.7 GPa to supplement the ultrasonic measurements, especially at low pressures. The acoustic velocities were then used to calculate a pressure–volume (P–V) equation of state, the bulk and shear moduli, and the Poisson's ratio. These quantities are compared with the reported pressure-dependent properties of related polymers over this range of pressures.

Acoustical properties↗
DOE OSTI2025

Multimaterial 3D Printing in Activating Bath Enables In Situ Polymerization of Thermosets with Intricate Geometries and Diverse Elastic Behaviors

Polydicyclopentadiene, p(DCPD), is a high‐performance thermoset valued for its exceptional toughness, strength, and stiffness. When copolymerized with 1,5‐cyclooctadiene (COD), its mechanical properties can be tuned from glassy to rubbery at room temperature. While frontal polymerization enables a rapid and energy‐efficient route to 3D print DCPD‐based materials, challenges such as ink shelf life and gravitational distortion, especially in direct ink writing of soft COD‐rich formulations, must be considered. Here, a complementary chemical strategy is presented, embedded 3D printing, that enables localized in situ polymerization of printed DCPD/COD inks within a reactive support matrix. The matrix provides both physical support and a reservoir of chemical activator, which diffuses into the ink, activates a latent bis(N‐heterocyclic carbene) Ru precatalyst, and initiates ring‐opening metathesis polymerization. Curing begins at the ink–matrix interface and propagates inward via diffusion, stabilizing the interface and preventing capillary‐driven deformation regardless of the matrix yield stress. This approach eliminates the need for cold storage, external curing, or photoinitiation, significantly expanding the processing window. Using this method, diverse thermosetting and elastomeric architectures are fabricated with features as small as 5 µm and aspect ratios of 100, including interlinked chains, shallow spherical shells exhibiting snap‐through buckling, and hair‐like fin arrays inaccessible through traditional techniques.

chemical activation↗
DOE OSTI2024

Data‐driven variational method for discrepancy modeling: Dynamics with small‐strain nonlinear elasticity and viscoelasticity

Abstract The effective inclusion of a priori knowledge when embedding known data in physics‐based models of dynamical systems can ensure that the reconstructed model respects physical principles, while simultaneously improving the accuracy of the solution in the previously unseen regions of state space. This paper presents a physics‐constrained data‐driven discrepancy modeling method that variationally embeds known data in the modeling framework. The hierarchical structure of the method yields fine scale variational equations that facilitate the derivation of residuals which are comprised of the first‐principles theory and sensor‐based data from the dynamical system. The embedding of the sensor data via residual terms leads to discrepancy‐informed closure models that yield a method which is driven not only by boundary and initial conditions, but also by measurements that are taken at only a few observation points in the target system. Specifically, the data‐embedding term serves as residual‐based least‐squares loss function, thus retaining variational consistency. Another important relation arises from the interpretation of the stabilization tensor as a kernel function, thereby incorporating a priori knowledge of the problem and adding computational intelligence to the modeling framework. Numerical test cases show that when known data is taken into account, the data driven variational (DDV) method can correctly predict the system response in the presence of several types of discrepancies. Specifically, the damped solution and correct energy time histories are recovered by including known data in the undamped situation. Morlet wavelet analyses reveal that the surrogate problem with embedded data recovers the fundamental frequency band of the target system. The enhanced stability and accuracy of the DDV method is manifested via reconstructed displacement and velocity fields that yield time histories of strain and kinetic energies which match the target systems. The proposed DDV method also serves as a procedure for restoring eigenvalues and eigenvectors of a deficient dynamical system when known data is taken into account, as shown in the numerical test cases presented here.

Masud, Arif↗
DOE OSTI2024

Temporal and Spatial Evolution of Non-Elastic Strain Accumulation in Stanstead Granite During Brittle Creep

Understanding the long-term behavior of brittle rocks requires fundamental consideration of time-dependent strain evolution and brittle creep processes. Previous studies have evaluated sub-critical crack growth during time-dependent deformation and damage evolution in brittle rocks; however, there is an incomplete knowledge of how damage evolves spatially and temporally within the body of intact rocks, where distributed regions of damage interact and coalesce during creep. This paper presents laboratory research focusing on evaluating brittle creep damage processes in Stanstead granite (SG) using 2-dimensional digital image correlation (2D-DIC). In the laboratory, the prismatic SG specimens were loaded beyond an estimated Crack Damage stress threshold (CD) level and then maintained a constant stress to initiate the creep process. DIC was used to characterize full-field spatiotemporal strain evolution, which was then interpreted in the context of local regions of “damage”, determined according to a strain-based criterion. Here, a method was proposed for identifying “existing” and “new” damage regions over specified intervals during the test, followed by spatial clustering of these regions to assess their spatiotemporal evolution. The clustering analysis results demonstrated the extension of existing damage regions was the main damage process during brittle creep, which is consistent with existing models of sub-critical crack growth. In addition, temporal analysis of tensile and shear strains on a point-by-point basis revealed both new damage formation and the strain concentration within existing damaged regions significantly contribute to overall specimen strain during primary creep. In contrast, during secondary creep, increases in specimen deformation are influenced by the accumulation of strains within already damaged regions.

58 GEOSCIENCES↗
DOE OSTI2025

Embedded 3D printing of isotropically enhanced lattice structure with programmable elasticity via three-dimensional fiber alignment

Engineered composites with fiber orientation in a complete three-dimensional (3D) context are highly desired for maximizing structural performance as they mimic the well-organized fiber arrangement in natural composites. However, most fabrication techniques lack the capability for through-plane (z-axis) alignment, limiting the possibility for structural design and performance optimization. In this study, we present an embedded 3D printing approach that enables complete 3D fiber alignment, including region-specific alignment control along the through-plane direction. Our method utilizes the extrusion and suspension of composite inks within a supporting matrix. By regionally adjusting the ink extrusion rate and nozzle translation rate, we locally manipulated the fiber arrangement within the composite 3D architecture. Here, we demonstrate that the 3D alignment of carbon fibers in a structural composite significantly enhances both mechanical properties and thermal conductivity in the in-plane and through-plane directions. Furthermore, by strategically arranging regional fiber orientations within micro-architected materials, we fabricated metamaterial engineered composites with programmable anisotropy ratio and directional moduli.

36 MATERIALS SCIENCE↗