Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ELASTIC BODY”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Assessing Critical Conditions for Scour Near Obstructions using Bed Shear, Particle Onset of Motion Balances, and CFD-DEM Modeling of Granular Beds

Computational Fluid Dynamics combined with a Discrete Element Method is one of the computational methods that can be used to model multiphase flows. In this method various phases, gas and liquid or solid, are present in the same computational domain. The local averaged Navier–Stokes equations determine the flow of the continuous phase fluid and are solved using the traditional CFD finite volume approach. DEM is based on a Lagrangian formulation, which solves the equations of motion, expressed in ordinary differential equations, for representative particles as they move in space and time. The interactions between the continuous fluid phase and discrete solid phase are modeled with the use of Newton’s laws of motion via drag force. The particles interact with each other and with the boundaries of the fluid continuum, and the resulting contact forces are included in the equations of motion. The properties of solid particles and boundaries are treated as elastic bodies, with specified density, elastic modulus, and Poisson’s ratio. Particle shapes may vary from single spherical particles to more complex-shaped composite particles. The particles may be introduced into the domain by random or structured injection at a point, surface, or volume, depending on the application. More details on the formulation can be found in the Simcenter STAR-CCM+ User’s Manual and OpenFOAM website.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Diffuse energy transport and coda-wave interferometry for resonant transmission between reverberant structures

Approximate analytic and numerical theories are developed with which to model, and compare with laboratory measurements, the diffuse ultrasonic field transmitted from one elastic body to another through a resonant coupling. Particular focus is on the sensitivity of the transmitted field to perturbations in the natural frequency of the coupling, and on the manifestation of such perturbations as time delays in the second body. Immediate application is to measurements of slow dynamic elastic nonlinearity. It is found that the developed approximate theories do a good job of predicting the time dependence of mean square transmitted signal and signal spectra. They also predict the time delays, both their erratic character and their means. Analytic dependence of these time shifts on the fractional frequency perturbation df/f is derived.

42 ENGINEERING↗
DOE OSTI2020

A Generalized Stress Inversion Approach With Application to Residual Stress Estimation

In this work develop a generalized stress inversion technique (or the generalized inversion method) capable of recovering stresses in linear elastic bodies subjected to arbitrary cuts. Specifically, given a set of displacement measurements found experimentally from digital image correlation (DIC), we formulate a stress estimation inverse problem as a partial differential equation-constrained optimization problem. We use gradient-based optimization methods, and we accordingly derive the necessary gradient and Hessian information in a matrix-free form to allow for parallel, large-scale operations. By using a combination of finite elements, DIC, and a matrix-free optimization framework, the generalized inversion method can be used on any arbitrary geometry, provided that the DIC camera can view a sufficient part of the surface. We present numerical simulations and experiments, and we demonstrate that the generalized inversion method can be applied to estimate residual stress.

42 ENGINEERING↗
DOE OSTI2021

Impact of three-body forces on elastic nucleon-nucleus scattering observables

In a previous series of papers we investigated the domain of applicability of chiral potentials to the construction of a microscopic optical potential (OP) for elastic nucleon-nucleus scattering. The OP was derived at the first order of the spectator expansion of the Watson multiple scattering theory and its final expression was a folding integral between the nucleon-nucleon ( N N ) t matrix and the nuclear density of the target. In the calculations N N and three-nucleon ( 3 N ) chiral interactions were used for the target density and only the N N interaction for the N N t matrix. The purpose of this study is to achieve another step towards the calculation of a more consistent OP introducing the 3 N force also in the dynamic part of the OP. The full treatment of the 3 N interaction is beyond our present capabilities. Thus, here it is approximated with a density dependent N N interaction obtained after the averaging over the Fermi sphere. In practice, in our model the 3 N force acts as a medium correction of the bare N N interaction used to calculate the t matrix. Even if the 3 N force is treated in an approximate way, this method naturally extends our previous model of the OP and allows a direct comparison of our present and previous results. We consider as case studies the elastic scattering of nucleons off C 12 and O 16 . We present results for the differential cross section and the spin observables for different values of the projectile energy. From the comparison with the experimental data and with the results of our previous model we assess the importance of the 3 N interaction in the dynamic part of the OP. Conclusions: Our analysis indicates that the contribution of the 3 N force in the t matrix is small for the differential cross section and it is sizable for the spin observables, in particular, for the analyzing power. We find that the two-pion exchange term is the major contributor to the 3 N force. A chiral expansion order-by-order analysis of the scattering observables confirms the convergence of our results at the next-to-next-to-next-to-leading-order, as already established in our previous work.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Kohn anomaly and elastic softening in body-centered cubic molybdenum at high pressure

Transition metals in body-centered cubic (bcc) structures under compression can display several novel physical properties because of their complex electronic structures and electron-phonon interactions. Here, we used inelastic x-ray scattering experiments in a diamond-anvil cell up to ∼45 GPa and density-functional theory calculations up to 210 GPa to investigate the phonon dispersions, and electronic and elastic properties of single-crystal molybdenum (Mo). Our results show a pressure-induced Kohn anomaly at 𝑞∼0.5 along the [ξ00] direction in the longitudinal acoustic mode at ∼45 GPa; this anomaly is triggered by the pressure-enhanced Fermi-surface nesting effect. Theoretical calculations show that electron redistributions in the 𝑠-to-𝑑 orbitals of bcc-Mo contribute to the shear modulus anomaly at ∼50 GPa. In contrast, the Young's modulus anomaly in bcc-Mo at ∼210 GPa results from a Lifshitz-type electronic topological transition. In conclusion, our results shed light on the complex electronic behaviors that are associated with macroscopic elastic properties in typical bcc 𝑑-block transition metals under compression.

3-dimensional systems↗
DOE OSTI2021

Deuteron-Deuteron Elastic and Three and Four-Body Breakup Scattering Using the Faddeev-Yakubovskii Equations

The deuteron-deuteron elastic and three and four-body breakup scattering cross section have been calculated using the Faddeev-Yakubovskii (FY) chain-of-partition momentum-space equations. In this calculation the initial two-cluster potential is split into separable and non-separable components, and the effective potential is reduced to elastic two-body and three and four-body breakup open-channels and a closed-channel many-body contribution. The closed-channel contribution is determined by minimizing a variational bound. The Coulomb interaction was included by expanding the initial and final Coulomb states in a Coulomb-Sturmian basis. The three sets of chain-of-partition integral equations were solved for the elastic and three and four-body breakup scattering amplitudes. The calculations were performed for the S = 2 spin/L = 0 angular-momentum state. The elastic and double-breakup calculations were performed for energies up to E = 5.48 MeV, while the single-breakup calculations were performed for energies up to E = 4.17 MeV. In the case of elastic scattering, the calculated scattering length of 5 a dd = 7.8 ± 0.3 fm is in good agreement with a FY cluster reduction calculation. The calculated phase shift is smaller than that predicted by the resonating group model and this difference is believed to be due to the differences in the potential and calculation methods. The breakup cross sections were calculated as a function of initial deuteron momentum and fragmented-deuteron momentum. Here, the d+d → d+n+p cross sections were compared with neutron yield measurements and, while the measurements also included the L > 0 components, the general features were consistent. Estimates of the calculational uncertainties/bias are provided.

42 ENGINEERING↗
DOE OSTI2020

Excitons in two-dimensional metal-halide perovskites

Georgia Tech School of Chemistry and Biochemistry Special Summer Seminar Series, Prof Carlos Silva; While polarons --- charges bound to a lattice deformation induced by electron-phonon coupling --- are primary photoexcitations at room temperature in bulk metal-halide hybrid organic-inorganic perovskites (HOIP), excitons --- Coulomb-bound electron-hole pairs --- are the stable quasi-particles in their two-dimensional (2D) analogues. Here we address the fundamental question: are polaronic effects consequential for excitons in 2D-HIOPs? Based on our recent work, we argue that polaronic effects are manifested intrinsically in the exciton spectral structure, which is comprised of multiple non-degenerate resonances with constant inter-peak energy spacing. We highlight measurements of population and dephasing dynamics that point to the apparently deterministic role of polaronic effects in excitonic properties. We contend that an interplay of long-range and short-range exciton-lattice couplings give rise to exciton polarons, a character that fundamentally establishes their effective mass and radius, and consequently, their quantum dynamics. Given this complexity, a fundamentally far-reaching issue is how Coulomb-mediated many-body interactions---elastic scattering such as excitation-induced dephasing, inelastic exciton bimolecular scattering, and multi-exciton binding---depend upon the specific exciton-lattice coupling within the structured excitation lineshape. We measure the intrinsic and density-dependent exciton dephasing rates of the multiple excitons and their dependence on temperature by means of two-dimensional coherent excitation spectroscopy. We find that diverse excitons display distinct intrinsic dephasing rates mediated by phonon scattering involving different effective phonons, and contrasting rates of exciton-exciton elastic scattering. These findings establish specifically the consequence of distinct lattice dressing on exciton many-body quantum dynamics, which critically define fundamental optical properties that underpin photonics and quantum optoelectronics.

solar cells, metal-halide perovskites, excitons, p↗
DOE OSTI2021

Crack Stability in Breached Fuel

This report describes the fracture mechanics formalism to evaluate the stability of cracks in a fuel cladding. Recognized consensus-body linear elastic fracture mechanics (LEFM) was applied to identify the crack instability length, or the length at which unstable mechanical crack extension would occur, as a function of pellet swelling loading (radial strain and fraction conversion from UO 2 to U 3 O 8 ) at the local cracked cladding region for two postulated fracture toughness (K IC ) values of the cladding. The crack opening displacement (COD) and the crack opening area (COA) were also identified. This analysis informs evaluations for crack extension and the potential for pellet debris loss from the fuel rod for cases of pellet oxidation in dry storage canisters where inadvertent residual water may undergo radiolysis causing oxidizing conditions to pellets exposed to the canister environment through breached cladding.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Crack Stability in Breached Fuel

This report describes the fracture mechanics formalism to evaluate the stability of cracks in a fuel cladding. Recognized consensus-body linear elastic fracture mechanics (LEFM) was applied to identify the crack instability length, or the length at which unstable mechanical crack extension would occur, as a function of pellet swelling loading (radial strain and fraction conversion from UO 2 to U 3 O 8 ) at the local cracked cladding region for two postulated fracture toughness (K IC ) values of the cladding. The crack opening displacement (COD) and the crack opening area (COA) were also identified. This analysis informs evaluations for crack extension and the potential for pellet debris loss from the fuel rod for cases of pellet oxidation in dry storage canisters where inadvertent residual water may undergo radiolysis causing oxidizing conditions to pellets exposed to the canister environment through breached cladding.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

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↗
DOE OSTI2019

A charge optimized many-body potential for iron/iron-fluoride systems

Here, a classical interatomic potential for iron/iron-fluoride systems is developed in the framework of the charge optimized many-body (COMB) potential. This interatomic potential takes into consideration the effects of charge transfer and many-body interactions depending on the chemical environment. The potential is fitted to a training set composed of both experimental and ab initio results of the cohesive energies of several Fe and FeF 2 crystal phases, the two fluorine molecules F 2 and the F 2 –1 dissociation energy curve, the Fe and FeF 2 lattice parameters of the ground state crystalline phase, and the elastic constants of the body centered cubic Fe structure. The potential is tested in an NVT ensemble for different initial structural configurations as the crystal ground state phases, F 2 molecules, iron clusters, and iron nanospheres. In particular, we model the FeF 2 /Fe bilayer and multilayer interfaces, as well as a system of square FeF 2 nanowires immersed in an iron solid. It has been shown that there exists a reordering of the atomic positions for F and Fe atoms at the interface zone; this rearrangement leads to an increase in the charge transfer among the atoms that make the interface and put forward a possible mechanism of the exchange bias origin based on asymmetric electric charge transfer in the different spin channels.

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

A 3D Full Stress Tensor Model for Oklahoma

Abstract The stress tensor is an important property for upper crustal studies such as those that involve pore fluids and earthquake hazards. At tectonic plate scale, plate boundary forces and mantle convection are the primary drivers of the stress field. In many local settings (10–100 s of km and <10 km depth) in tectonic plate interiors, we can simplify by assuming a constant background stress field that is perturbed by local heterogeneity in density and elasticity. Local stress orientation and sometimes magnitude can be estimated from earthquake and borehole‐based observations when available. Modeling of the local stress field often involves interpolating sparse observations. We present a new method to estimate the 3D stress field in the upper crust and demonstrate it for Oklahoma. We created a 3D material model by inverting multiple types of geophysical observations simultaneously. Integrating surface‐wave dispersion, local travel times and gravity observations produces a model of P‐wave velocity, S‐wave velocity, and density. The stress field can then be modeled using finite element simulations. The simulations are performed using our simplified view of the local stress field as the sum of a constant background stress field that is perturbed by local density and elasticity heterogeneity and gravitational body forces. An orientation of N82°E, for the maximum compressive tectonic force, best agrees with previously observed stress orientations and faulting types in Oklahoma. The gravitational contribution of the horizontal stress field has a magnitude comparable to the tectonic contribution for the upper 5 km of the subsurface.

58 GEOSCIENCES↗
DOE OSTI2021

Shock melting and the hcp-bcc phase boundary of Mg under dynamic loading

The high-temperature and high-pressure response of magnesium has been investigated through shock-release experiments performed up to shock melting. The longitudinal and bulk sound speeds of Mg are reported along the Hugoniot from 25 to 56 GPa and used to determine the elastic properties and Grüneisen parameter. The hexagonal close-packed (hcp)–body-centered cubic (bcc) phase transition is marked by a reduction in the determined shear wave speed. Thermal softening is observed to begin between 40 and 44 GPa, with incipient melt at 55.5 GPa, in close agreement with previous diffraction measurements under laser shock loading. Examination of the release profiles showed that two different responses were observed, depending on the peak stress. When plotted in a stress-energy phase diagram, the two responses are observed to form separate lines that intersect the static hcp-bcc phase transition and incipient melt. The results indicate that the hcp-bcc transition occurs on the Hugoniot at 28.4 GPa and place the hcp-bcc-liquid triple point ~20 GPa. In conclusion, a multiphase equation of state is developed which places the melt boundary below the previously reported static measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Universal behavior in fragmenting brittle, isotropic solids across material properties

A bonded particle model is used to explore how variations in the material properties of brittle, isotropic solids affect critical behavior in fragmentation. To control material properties, a model is proposed which includes breakable two- and three-body particle interactions to calibrate elastic moduli and mode I and mode II fracture toughnesses. In the quasistatic limit, fragmentation leads to a power-law distribution of grain sizes which is truncated at a maximum grain mass that grows as a nontrivial power of system size. In the high-rate limit, truncation occurs at a mass that decreases as a power of increasing rate. A scaling description is used to characterize this behavior by collapsing the mean-square grain mass across rates and system sizes. Consistent scaling persists across all material properties studied, although there are differences in the evolution of grain size distributions with strain as the initial number of grains at fracture and their subsequent rate of production depend on Poisson's ratio. Importantly, this evolving granular structure is found to induce a unique rheology where the ratio of the shear stress to pressure, an internal friction coefficient, decays approximately as the logarithm of increasing strain rate. The stress ratio also decreases at all rates with increasing strain as fragmentation progresses and depends on elastic properties of the solid.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Oil-Pressure Based Apparatus for In-Situ High-Energy Synchrotron X-Ray Diffraction Studies During Biaxial Deformation

Background: Understanding biaxial loading response at the microstructural level is crucial in helping better design sheet manufacturing processes and calibrate/validate material deformation models. Objective: The objective of this work was to develop a low-cost testing apparatus to probe, with sufficient spatial resolution, the micro-mechanical response of a sheet material in-situ under biaxial loading conditions. Methods: The testing apparatus fabricated as a part of this study operates in a similar fashion to a standard bulge test and uses oil pressure to generate biaxial loading conditions. This biaxial testing apparatus was operated within a synchrotron beamline to characterize the mechanical response of a flash-processed steel sheet using in-situ high-energy X-ray diffraction (XRD) measurements. Further, the GSAS-II package was utilized to develop a workflow for the analysis of the large volume of diffraction data acquired. The workflow was then used to extract the peak position, width, and integrated intensity of the XRD peaks corresponding to the major body-centered cubic phase. Results: The equi-biaxial nature of the loading in the measured area was independently corroborated using experimental (XRD) and simulation (finite element analysis) methods. Furthermore, we discuss the evolution of elastic strain in the major body-centered cubic phase as a function of applied oil pressure and location on the steel sheet. Conclusions: A key advantage of the biaxial apparatus fabricated in this synchrotron study is demonstrated using the results obtained for the flash-processed steel sheet – i.e., mapping the lattice plane-dependent response to biaxial loading for a relatively large sample area in a spatially resolved manner.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Regional-Scale 3D Ground-Motion Simulations of M w 7 Earthquakes on the Hayward Fault, Northern California Resolving Frequencies 0–10 Hz and Including Site-Response Corrections

Large earthquake ground-motion simulations in 3D Earth models provide constraints on site-specific shaking intensities but have suffered from limited frequency resolution and ignored site response in soft soils. In this work we report new regional-scale 3D simulations for moment magnitude 7.0 scenario earthquakes on the Hayward Fault, northern California with SW4. Simulations resolved significantly broader band frequencies (0–10 Hz) than previous studies and represent the highest resolution simulations for any such earthquake to date. Seismic waves were excited by a kinematic rupture following Graves and Pitarka (2016) and obeyed wave propagation in a 3D Earth model with topography from the U.S. Geological Survey (USGS) assuming a minimum shear wavespeed, VSmin, of 500 m/s. We corrected motions for linear and nonlinear site response for the shear wavespeed, VS, from the USGS 3D model, using a recently developed ground-motion model (GMM) for Fourier amplitude spectra . At soft soil locations subjected to strong shaking, the site-corrected intensities reflect the competing effects of linear amplification by low V S material, reduction of stiffness during nonlinear deformation, and damping of high frequencies. Sites with near-surface V S of 500 m/s or greater require no linear site correction but can experience amplitude reduction due to nonlinear response. Averaged over all sites, we obtained reasonable agreement with empirical ergodic median GMMs currently used for seismic hazard and design ground motions (epsilon less than 1), with marked improvement at soft sedimentary sites. At specific locations, the simulated shaking intensities show systematic differences from the GMMs that reveal path and site effects not captured in these ergodic models. Results suggest how next generation regional-scale earthquake simulations can provide higher spatial and frequency resolution while including effects of soft soils that are commonly ignored in scenario earthquake ground-motion simulations.

58 GEOSCIENCES↗
DOE OSTI2023

Experimental and Phenomenological Investigations of the MiniBooNE Anomaly

This thesis covers a range of experimental and theoretical efforts to elucidate the origin of the $4.8\sigma$ MiniBooNE low energy excess (LEE). We begin with the follow-up MicroBooNE experiment, which took data along the BNB from 2016 to 2021. This thesis specifically presents MicroBooNE's search for $\nu_e$ charged-current quasi-elastic (CCQE) interactions consistent with two-body scattering. The two-body CCQE analysis uses a novel reconstruction process, including a number of deep-learning-based algorithms, to isolate a sample of $\nu_e$ CCQE interaction candidates with $75\%$ purity. The analysis rules out an entirely $\nu_e$-based explanation of the MiniBooNE excess at the $2.4\sigma$ confidence level. We next perform a combined fit of MicroBooNE and MiniBooNE data to the popular $3+1$ model; even after the MicroBooNE results, allowed regions in $\Delta m^2$-$\sin^2 2_{\theta_{\mu e}}$ parameter space exist at the $3\sigma$ confidence level. This thesis also demonstrates that the MicroBooNE data are consistent with a $\overline{\nu}_e$-based explanation of the MiniBooNE LEE at the $<2\sigma$ confidence level. Next, we investigate a phenomenological explanation of the MiniBooNE excess combining the $3+1$ model with a dipole-coupled heavy neutral lepton (HNL). It is shown that a 500 MeV HNL can accommodate the energy and angular distributions of the LEE at the $2\sigma$ confidence level while avoiding stringent constraints derived from MINER$\nu$A elastic scattering data. Finally, we discuss the Coherent CAPTAIN-Mills experiment--a 10-ton light-based liquid argon detector at Los Alamos National Laboratory. The background rejection achieved from a novel Cherenkov-based reconstruction algorithm will enable world-leading sensitivity to a number of beyond-the-Standard Model physics scenarios, including dipole-coupled HNLs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗