Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear 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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A compute-bound formulation of Galerkin model reduction for linear time-invariant dynamical systems

This work aims to advance computational methods for projection-based reduced-order models (ROMs) of linear time-invariant (LTI) dynamical systems. For such systems, current practice relies on ROM formulations expressing the state as a rank-1 tensor (i.e., a vector), leading to computational kernels that are memory bandwidth bound and, therefore, ill-suited for scalable performance on modern architectures. This weakness can be particularly limiting when tackling many-query studies, where one needs to run a large number of simulations. This work introduces a reformulation, called rank-2 Galerkin, of the Galerkin ROM for LTI dynamical systems which converts the nature of the ROM problem from memory bandwidth to compute bound. We present the details of the formulation and its implementation, and demonstrate its utility through numerical experiments using, as a test case, the simulation of elastic seismic shear waves in an axisymmetric domain. We quantify and analyze performance and scaling results for varying numbers of threads and problem sizes. In conclusion, we present an end-to-end demonstration of using the rank-2 Galerkin ROM for a Monte Carlo sampling study. We show that the rank-2 Galerkin ROM is one order of magnitude more efficient than the rank-1 Galerkin ROM (the current practice) and about 970 times more efficient than the full-order model, while maintaining accuracy in both the mean and statistics of the field.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Universal Maximum Strength of Solid Metals and Alloys

As described here, interstitial electron density ρ o is offered as a direct metric for maximum strength in metals, arising from universal properties derived from an electron-gas – ρ o sets the exchange-correlation parameter r s in density-functional theory (DFT). It holds also for maximum shear strength τ max in polycrystals [Phys. Rev. Lett. 124, 125501]. Elastic moduli and τ max for polycrystalline (amorphous) metals are linear with ρ o and melting T m (glass-transition T g ) temperature. ρ o or r s , even with rule-of-mixture estimate, predicts relative strength for rapid, reliable selection of high-strength alloys with ductility, as confirmed for elements to steels to complex solid-solutions, and validated experimentally.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Crystal structure of calcium-ferrite type NaAlSiO4 up to 45 GPa

Abstract Alkali-rich aluminous high-pressure phases including calcium-ferrite (CF) type NaAlSiO4 are thought to constitute ~20% by volume of subducted mid-ocean ridge basalt (MORB) under lower mantle conditions. As a potentially significant host for incompatible elements in the deep mantle, knowledge of the crystal structure and physical properties of CF-type phases is therefore important to understanding the crystal chemistry of alkali storage and recycling in the Earth’s mantle. We determined the evolution of the crystal structure of pure CF-NaAlSiO4 and Fe-bearing CF-NaAlSiO4 at pressures up to ~45 GPa using synchrotron-based, single-crystal X-ray diffraction. Using the high-pressure lattice parameters, we also determined a third-order Birch-Murnaghan equation of state, with V0 = 241.6(1) Å3, KT0 = 220(4) GPa, and KT0′ = 2.6(3) for Fe-free CF, and V0 = 244.2(2) Å3, KT0 = 211(6) GPa, and KT0′ = 2.6(3) for Fe-bearing CF. The addition of Fe into CF-NaAlSiO4 resulted in a 10 ± 5% decrease in the stiffest direction of linear compressibility along the c-axis, leading to stronger elastic anisotropy compared with the Fe-free CF phase. The NaO8 polyhedra volume is 2.6 times larger and about 60% more compressible than the octahedral (Al,Si)O6 sites, with K0NaO8 = 127 GPa and K0(Al,Si)O6 ~304 GPa. Raman spectra of the pure CF-type NaAlSiO4 sample shows that the pressure coefficient of the mean vibrational mode, 1.60(7) cm–1/GPa, is slightly higher than 1.36(6) cm−1/GPa obtained for the Fe-bearing CF-NaAlSiO4 sample. The ability of CF-type phases to contain incompatible elements such as Na beyond the stability field of jadeite requires larger and less-compressible NaO8 polyhedra. Detailed high-pressure crystallographic information for the CF phases provides knowledge on how large alkali metals are hosted in alumina framework structures with stability well into the lowermost mantle.

Geochemistry & Geophysics↗
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↗
DOE OSTI2023

Computation of high-frequency magnetoelastic waves in layered materials

Here, the direct calculation of magnetoelastic wave dispersion in layered media is presented using an efficient, accurate computational technique. The governing, coupled equations for elasticity and magnetism, the Navier and Landau-Lifshitz equations, respectively, are linearized to form a quadratic eigenvalue problem that determines a complex web of wave-number–frequency dispersion branches and their corresponding mode profiles. Numerical discretization of the eigenvalue problem via a spectral collocation method (SCM) is employed to determine the complete dispersion maps for both a single, finite-thickness magnetic layer and a finite magnetic-nonmagnetic double-layer. The SCM, previously used to study elastic waves in nonmagnetic media, is fast, accurate, and adaptable to a variety of sample configurations and geometries. Emphasis is placed on the extremely high-frequency regimes being accessed in ultrafast magnetism experiments. The dispersion maps and modes provide insight into how energy propagates through the coupled system, including how energy can be transferred between elastic- and magnetic-dominated waves as well as between different layers. The numerical computations for a single layer are further understood by a simplified analytical calculation in the high-frequency, exchange-dominated regime where the resonance condition required for energy exchange (an anticrossing) between quasi-elastic and quasi-magnetic dispersion branches is determined. Nonresonant interactions are shown to be well approximated by the dispersion of uncoupled elastic and magnetic waves. The methods and results provide fundamental theoretical tools to model and understand current and future magnetic devices powering spintronic innovation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Rheological and Structural Properties of k‐Carrageenan/Xanthan Gum Gummies Architected With a New Natural Additive

ABSTRACT Biohybrids (BHs)based on bentonite and anthocyanins can be produced as alternatives to synthetic colorants. However, no information is available on the interaction of biohybrid with food system. This research studies the rheological properties of gummies containing k‐carrageenan and xanthan gum incorporated with a BH. The incorporation of BH (≤ 2% w/w) increased the flow consistency index and shear thinning effect in melted gummies at 85°C. During cooling from 90°C to 25°C, the presence of BH altered the sol–gel transition temperature typical in carrageenan systems, probably by the formation of carrageenan complexes. At 25°C, all gummies had a solid‐like response with elastic modulus ( G ') higher than viscous modulus ( G "). G ' values increased linearly with BH concentration, and the linear viscoelastic region and breaking strain decreased with BH concentration, indicating that this natural additive acted as a reinforcing material.

Valencia, German Ayala [Department of Chemical and↗
DOE OSTI2024

Wind Turbine Aeroelastic Stability in OpenFAST

Wind turbines are growing in size and increasingly suffer from aeroelastic instabilities. Unfortunately, numerical models often show inconsistent results during verification studies. We address this gap by first introducing novel linearization capabilities within the open-source aero-hydro-servo-elastic framework OpenFAST. Next, a code-to-code benchmark study is presented that compares modal parameters between OpenFAST and HAWCStab2 for a land-based version of the International Energy Agency 15-MW reference wind turbine modeled with quasi-steady aerodynamics. The two solvers are in strong agreement except for discrepancies in the second rotor flapwise modes. The differences are attributed to the torsional flexibility of the tower, which is assumed torsionally stiff in the OpenFAST model. Work is ongoing to close this modeling gap. The aeroelastic stability of a low-specific-power land-based wind turbine is also investigated. The impact of design choices is discussed, high-lighting how narrow the margins are between a stable design and an unstable design.

17 WIND ENERGY↗
DOE OSTI2022

Effects of Model Uncertainties in Underground Chemical Explosions on Far-field Results

We used the CTH shock physics code to simulate the explosion of an 18-t chemical explosive at a depth of 250 m. We used the CTH in the two-dimensional axisymmetric (cylindrical) geometry (2DC) and most simulations included fully tamped explosions in wet tuff. Our study focused on parametric studies of three of the traditional strength models available in CTH, namely, geologic-yield, elastic perfectly-plastic von Mises, and Johnson-Cook strength (flow stress) models. We processed CTH results through a code that generates Reduced Displacement Potential (RDP) histories for each simulation. Since RDP is the solution of the linear wave equation in spherical coordinates, it is mainly valid at far-enough distance from the explosion the elastic radius. Among various parameters examined, we found the yield strength to have the greatest effect on the resulting RDP, where the peak RDP reduces almost linearly in log-log space as the yield strength increases. Moreover, an underground chemical explosion results in a cavity whose final diameter is inversely proportional to the material yield strength, i.e., as the material's yield strength increases the resulting final cavity radius decreases. Additionally, we found the choice of explosive material (COMP-C4 versus COMP-B) has minor effects on the peak RDP, where denser COMP-C4 shows higher peak RDP than the less dense COMP-B by a factor of ~1.1. In addition to wet tuff, we studied explosions in dry tuff, salt, and basalt, for a single strength model and yield strength value. We found wet tuff has the highest peak RDP value, followed by dry tuff, salt, and basalt. 2DC simulations of explosions in 11 m radius spherical, hemispherical, and cylindrical cavities showed the RDP signals have much lower magnitude than tamped explosions, where the cavity explosions mimicked nearly decoupled explosions.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2022

A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

58 GEOSCIENCES↗
DOE OSTI2022

Ab initio investigation of elastic properties of dilute Cu alloys for high-gradient accelerating structures

The elastic properties and solid solution strengthening of several Cu alloys that may be promising for high gradient applications were studied using ab initio calculations based on density functional theory. Specifically, the temperature and solute concentration dependence of elastic properties including bulk and shear moduli, thermal expansion coefficient, and heat capacity of these alloys were calculated. We found that the bulk and shear moduli are linearly dependent on the solute concentration in the dilute regime and change significantly for different Cu alloys. The thermal expansion coefficient and heat capacity are comparatively less sensitive to solute concentrations. Solid solution strengthening is observed to be linearly correlated with the change of lattice parameter of Cu. This work provides an in-depth understanding of the elastic properties of Cu alloys, which is key for predicting their performance at high gradient in normal conducting accelerating structures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Evaluation of dual-weighted residual and machine learning error estimation for projection-based reduced-order models of steady partial differential equations

Projection-based reduced-order models (pROMs) show great promise as a means to accelerate many-query applications such as forward error propagation, solving inverse problems, and design optimization. In order to deploy pROMs in the context of high-consequence decision making, accurate error estimates are required to determine the region(s) of applicability in the parameter space. The following paper considers the dual-weighted residual (DWR) error estimate for pROMs and compares it to another promising pROM error estimate, machine learned error models (MLEM). Here, we show how DWR can be applied to ROMs and then evaluate DWR on two partial differential equations (PDEs): a two-dimensional linear convection–reaction–diffusion equation, and a three-dimensional static hyper-elastic beam. It is shown that DWR is able to estimate errors for pROMs extrapolating outside of their training set while MLEM is best suited for pROMs used to interpolate within the pROM training set.

42 ENGINEERING↗
DOE OSTI2021

Magnetoelastic standing waves induced in UO 2 by microsecond magnetic field pulses

Significance UO 2 is of particular interest due to its use as a nuclear fuel, its unusually poor thermal conductivity, and its uniquely strong magnetoelastic properties characterized by a linear coupling between strain and magnetization. We clearly observe mechanical standing waves and elastic resonances in a single crystal, induced by microsecond magnetic field pulses. The phase and amplitude of the mechanical resonances are modulated by the piezomagnetic switching behavior, offering a sensitive probe for otherwise hard to detect magnetic dynamics in UO 2 . We highlight the importance of magnetic domains in the attenuation of the mechanical resonances and the robustness of the antiferromagnetic state in UO 2 against externally applied magnetic fields to 150 T.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Model-independent determination of the dipole response of 66 Zn using quasimonoenergetic and linearly polarized photon beams

Background: Photon strength functions are an important ingredient in calculations relevant for the nucleosynthesis of heavy elements. The relation to the photoabsorption cross section allows to experimentally constrain photon strength functions by investigating the photoresponse of atomic nuclei. Purpose: Here, we determine the photoresponse of 66 Zn in the energy region of 5.6 MeV to 9.9 MeV and analyze the contribution of the 'elastic' decay channel back to the ground state. In addition, for the elastic channel electric and magnetic dipole transitions were separated. Methods: Nuclear resonance fluorescence experiments were performed using a linearly polarized quasi-monoenergetic photon beam at the High Intensity γ-ray Source. Photon beam energies from 5.6 to 9.9 MeV with an energy spread of about 3% were selected in steps of 200–300 keV. Two high purity germanium detectors were used for the subsequent γ-ray spectroscopy. Results: Full photoabsorption cross sections are extracted from the data making use of the monoenergetic character of the photon beam. For the ground-state decay channel, the average contribution of electric and magnetic dipole strengths is disentangled. The average branching ratio back to the ground state is determined as well. Conclusions: The new results indicate lower cross sections when compared to the values extracted from a former experiment using bremsstrahlung on 66 Zn. In the latter, the average branching ratio to the ground state is estimated from statistical-model calculations in order to analyze the data. Corresponding estimates from statistical-model calculations underestimate this branching ratio compared to the values extracted from the present analysis, which would partly explain the high cross sections determined from the bremsstrahlung data.

59 ≤ A ≤ 89↗
DOE OSTI2021

Report on FY 2021 testing in support of integrating Alloy 800H and 2.25Cr-1Mo for EPP analysis

The use of simplified Elastic-Perfectly Plastic (EPP) design methods avoids the complexities and limitations of the elastic analysis approach in evaluating strain limits and creep-fatigue damage for Class A metallic coolant boundary components in ASME Boiler and Pressure Vessel Code, Section III, Division 5, Subsection HB, Subpart B. No stress classification and linearization are required in the EPP methods. Further, the strain limits and creep-fatigue damage evaluation procedures based on elastic analysis results were developed for the reactor operating conditions where uncoupled, rate-independent plasticity and stationary creep can be used to approximate the deformation in the component adequately. But at very high temperatures, the deformation is rate-dependent, and it can only be approximated adequately by unified viscoplastic model. Thus, for example, the elastic analysis approach cannot be used to evaluate strain limits and creep-fatigue damage for Alloy 617 above 650 °C. The EPP methods remove the temperature restriction. They can be used to evaluate strain limits and creep-fatigue damage for the full range of temperatures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Ultrasound Characterization of Plastic Bonded Explosives With Varied Binder Content and Density

Advanced explosive formulations and methods of production require rapid characterization for both processes and materials. Here, we present results for ultrasound characterization of a plastic bonded explosive (PBX) formulation as a function of binder weight percent and pressing density using through-transmission sound speed measurements. The explosive in this study was composed of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) combined with the fluoropolymer binder Kynar Ultraflex B. Binder content ranged from 2.6 wt.% to 11.5 wt.% with associated explosive formulation densities from 1.69 g/cc to 1.89 g/cc. We find linear and nonlinear changes in longitudinal and shear sound speeds, peak frequencies, and measured elastic properties for the different PBX formulations. These measured changes can differentiate the tested formulations, highlighting the potential for ultrasound as a nondestructive characterization method for newly-formulated PBXs. We also demonstrate through-transmission measurements using different coupling media and present results from an ultrasound technique called resonant ultrasound spectroscopy (RUS). Both transmission and RUS approaches may allow for rapid and even automated characterization of PBX components suitable for high-throughput manufacturing operations.

Engineering↗
DOE OSTI2022

Sensitivity of the U.S. economy to oil prices controlling for domestic production and imports

This paper investigates the sensitivity of U.S. economic performance to oil price changes, accounting for changes in the domestic petroleum supply-demand balance over the last decade. Herein, a non-linear (asymmetric) autoregressive distributed lag (NARDL) model is used to estimate the U.S. GDP elasticity with respect to the oil price, controlling for oil production, consumption and imports, and macroeconomic variables. The positive and negative components of the oil price both have statistically significant long-run impacts on the real U.S. GDP. The parameter estimates imply that a 1% positive and permanent oil price shock, all else the same, would have a long-run impact of –0.045% on the U.S. economy (an elasticity of –0.045), but short-run parameters on the positive oil price terms are not significant at the 10% level. The long-run parameter on the negative oil price component term implies an elasticity of –0.034. Thus, controlling for domestic U.S. oil production and oil trade, the long-run oil price elasticity of the U.S. GDP remains within the range of estimates from previous studies. The results also show that domestic oil production and consumption have short-run impacts on the U.S. GDP. The potential extent of interactions among these variables, and implications for the net economic impacts, under an oil price shock are subjects of future research.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Viscoelastic necking dynamics between attractive microgels

Hypothesis: Microgels can deform and interpenetrate and display colloid/polymer duality. The effective interaction of microgels in the collapsed state is governed by the interplay of polymer–solvent interfacial tension and bulk elasticity. A connecting neck is shown to mediate microgel interaction, but its temporal evolution has not been addressed. Here we hypothesize that the necking dynamics of attractive microgels exhibits liquid-like or solid-like behavior over different time and length scales. Experiments: We simulate the merging and pinching of attractive microgels with different crosslinking densities in explicit solvent using dissipative particle dynamics. The temporal coalescence dynamics of microgels is investigated and compared with simple liquid and polymeric droplets. We model the neck growth on long time scales using Maxwell model of polymer relaxation and compare the theoretical prediction with simulation data. The mechanical strength of the neck is characterized systematically via simulated pinch-off of microgels by steered molecular dynamics. Findings: We evidence a crossover in the coalescence dynamics reflecting the viscoelastic signature of microgels. In contrast to the common knowledge that viscoelastic materials respond elastically on short time scales, the early expansion of the microgel neck exhibits a linear behavior, similar to the viscous coalescence of liquid droplets. However, the late regime with arrested dynamics resembles sintering of solid particles. Through an analytical model relating microgel dynamics to neck growth, we show that the long-term behavior is governed by stress relaxation of the polymers in the neck region and predict an exponential decay in the rate of growth, which agrees favorably with the simulation. Different from coalescence, the thread thinning in microgel breakup primarily highlights its polymeric characteristics.

42 ENGINEERING↗
DOE OSTI2020

Reconsidering astrophysical constraints on macroscopic dark matter

Macroscopic dark matter—“macros”—refers to a broad class of alternative candidates to particle dark matter with still unprobed regions of parameter space. These candidates would transfer energy primarily through elastic scattering with approximately their geometric cross section. For sufficiently large cross sections, the linear energy deposition could produce observable signals if a macro were to pass through compact objects such as white dwarfs or neutron stars in the form of thermonuclear runaway, leading to a type Ia supernova or superburst, respectively. Here, we update the constraints from white dwarfs. These are weaker than previously inferred in important respects because of more careful treatment of the passage of a macro through the white dwarf and greater conservatism regarding the size of the region that must be heated to initiate runaway. On the other hand, we place more stringent constraints on macros at a low cross section, using new data from the Montreal White Dwarf Database. New constraints are inferred from the low mass x-ray binary 4U 1820-30, in which more than a decade passed between successive superbursts. Updated microlensing constraints are also reported.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗