Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational inelasticity”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

An inelastic analysis of a welded aluminum joint

Butt-weld joints are most commonly designed into pressure vessels which then become as reliable as the weakest increment in the weld chain. In practice, weld material properties are determined from tensile test specimen and provided to the stress analyst in the form of a stress versus strain diagram. Variations in properties through the thickness of the weld and along the width of the weld have been suspect but not explored because of inaccessibility and cost. The purpose of this study is to investigate analytical and computational methods used for analysis of welds. The weld specimens are analyzed using classical elastic and plastic theory to provide a basis for modeling the inelastic properties in a finite-element solution. The results of the analysis are compared to experimental data to determine the weld behavior and the accuracy of prediction methods. The weld considered in this study is a multiple-pass aluminum 2219-T87 butt weld with thickness of 1.40 in. The weld specimen is modeled using the finite-element code ABAQUS. The finite-element model is used to produce the stress-strain behavior in the elastic and plastic regimes and to determine Poisson's ratio in the plastic region. The value of Poisson's ratio in the plastic regime is then compared to experimental data. The results of the comparisons are used to explain multipass weld behavior and to make recommendations concerning the analysis and testing of welds.

Vaughan, R. E.↗
DOE OSTI2023

Effects of anisotropy and disorder on the superconducting properties of niobium

We report results for the superconducting transition temperature and anisotropic energy gap for pure niobium based on Eliashberg’s equations and electron and phonon band structures computed from density functional theory. The electronic band structure is used to construct the Fermi surface and calculate the Fermi velocity at each point on the Fermi surface. The phonon bands are in excellent agreement with inelastic neutron scattering data. The corresponding phonon density of states and electron–phonon coupling define the electron–phonon spectral function, α 2 F(p, p'; ω), and the corresponding electron–phonon pairing interaction, which is the basis for computing the superconducting properties. The electron–phonon spectral function is in good agreement with existing tunneling spectroscopy data except for the spectral weight of the longitudinal phonon peak at ℏω LO = 23 meV. We obtain an electron–phonon coupling constant of λ = 1.057, renormalized Coulomb interaction μ* = 0.218, and transition temperature T c = 9.33 K. The corresponding strong-coupling gap at T = 0 is modestly enhanced, Δ 0 = 1.55 meV, compared to the weak-coupling BCS value $Δ^{\text {wc}}_{0}$ = 1.78 k B T c = 1.43 meV. The superconducting gap function exhibits substantial anisotropy on the Fermi surface. We analyze the distribution of gap anisotropy and compute the suppression of the superconducting transition temperature using a self-consistent T-matrix theory for quasiparticle-impurity scattering to describe niobium doped with non-magnetic impurities. We compare these results with experimental results on niobium SRF cavities doped with nitrogen impurities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2020

Solution of the Nonlinear High-Fidelity Generalized Method of Cells Micromechanics Relations via Order-Reduction Techniques

The High-Fidelity Generalized Method of Cells (HFGMC) is one technique, distinct from traditional finite element approaches, for accurately simulating nonlinear composite material behavior. In this work, the HFGMC global system of equations for doubly-periodic repeating unit cells with nonlinear constituents has been reduced in size through the novel application of a Petrov-Galerkin Proper Orthogonal Decomposition order-reduction scheme in order to improve its computational efficiency. Order-reduced models of an E-glass/Nylon 12 composite led to a 4.8 to 6.3x speedup in the equation assembly/solution runtime while maintaining model accuracy. This corresponded to a 21 to 38 percent reduction in total runtime. The significant difference in assembly/solution and total runtimes was attributed to the evaluation of integration point inelastic field quantities; this step was identical between the unreduced and order-reduced models. Nonetheless, order-reduced techniques offer the potential to significantly improve the computational efficiency of multiscale calculations.

Ricks, Trenton M.↗
NASA NTRS2001

Micromechanics-Based Inelastic Finite Element Analysis Accomplished Via Seamless Integration of MAC/GMC

A critical issue in the micromechanics-based analysis of composite structures becomes the availability of a computationally efficient homogenization technique: one that is 1) Capable of handling the sophisticated, physically based, viscoelastoplastic constitutive and life models for each constituent; 2) Able to generate accurate displacement and stress fields at both the macro and the micro levels; 3) Compatible with the finite element method. The Generalized Method of Cells (GMC) developed by Paley and Aboudi is one such micromechanical model that has been shown to predict accurately the overall macro behavior of various types of composites given the required constituent properties. Specifically, the method provides "closed-form" expressions for the macroscopic composite response in terms of the properties, size, shape, distribution, and response of the individual constituents or phases that make up the material. Furthermore, expressions relating the internal stress and strain fields in the individual constituents in terms of the macroscopically applied stresses and strains are available through strain or stress concentration matrices. These expressions make possible the investigation of failure processes at the microscopic level at each step of an applied load history.

Arnold, Steven M.↗
NASA NTRS1980

Three dimensional inelastic finite element analysis of laminated composites

Formulations of the inelastic response of laminated composites to thermal and mechanical loading are used as the basis for development of the computer NALCOM (Nonlinear Analysis of Laminated Composites) computer program which uses a fully three dimensional isoparametric finite element with 24 nodes and 72 degrees of freedom. An incremental solution is performed with nonlinearities introduced as pseudoloads computed for initial strains. Equilibrium iteration may be performed at every step. Elastic and elastic-plastic response of boron/epoxy and graphite/epoxy graphite/epoxy and problems of curing 0/90 sub s Gr/Ep laminates with and without circular holes are analyzed. Mechanical loading of + or - 45sub s Gr/Ep laminates is modeled and symmetry conditions which exist in angle-ply laminates are discussed. Results are compared to experiments and other analytical models when possible. All models are seen to agree reasonably well with experimetnal results for off-axis tensile coupons. The laminate analyses show the three dimensional effects which are present near holes and free corners.

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

Consistent Pl Analysis of Aqueous Uranium-235 Critical Assemblies

The lethargy-dependent equations of the consistent Pl approximation to the Boltzmann transport equation for slowing down neutrons have been used as the basis of an IBM 704 computer program. Some of the effects included are (1) linearly anisotropic center of mass elastic scattering, (2) heavy element inelastic scattering based on the evaporation model of the nucleus, and (3) optional variation of the buckling with lethargy. The microscopic cross-section data developed for this program covered 473 lethargy points from lethargy u = 0 (10 Mev) to u = 19.8 (0.025 ev). The value of the fission neutron age in water calculated here is 26.5 square centimeters; this value is to be compared with the recent experimental value given as 27.86 square centimeters. The Fourier transform of the slowing-down kernel for water to indium resonance energy calculated here compared well with the Fourier transform of the kernel for water as measured by Hill, Roberts, and Fitch. This method of calculation has been applied to uranyl fluoride - water solution critical assemblies. Theoretical results established for both unreflected and fully reflected critical assemblies have been compared with available experimental data. The theoretical buckling curve derived as a function of the hydrogen to uranium-235 atom concentration for an energy-independent extrapolation distance was successful in predicting the critical heights of various unreflected cylindrical assemblies. The critical dimensions of fully water-reflected cylindrical assemblies were reasonably well predicted using the theoretical buckling curve and reflector savings for equivalent spherical assemblies.

Fieno, Daniel↗
DOE OSTI2020

Predictive Model of Charge Mobilities in Organic Semiconductor Small Molecules with Force-Matched Potentials

Charge mobility of crystalline organic semiconductors (OSC) is limited by local dynamic disorder. Recently, the charge mobility for several high mobility OSCs, including TIPS-pentacene, were accurately predicted from a density functional theory (DFT) simulation constrained by the crystal structure and the inelastic neutron scattering spectrum, which provide direct measures of the structure and the dynamic disorder in the length scale and energy range of interest. However, the computational expense required for calculating all of the atomic and molecular forces is prohibitive. In this work, we demonstrate the use of density functional tight binding (DFTB), a semiempirical quantum mechanical method that is 2 to 3 orders of magnitude more efficient than DFT. We show that force matching a many-body interaction potential to DFT derived forces yields highly accurate DFTB models capable of reproducing the low-frequency intricacies of experimental inelastic neutron scattering (INS) spectra and accurately predicting charge mobility. We subsequently predicted charge mobilities from our DFTB model of a number of previously unstudied structural analogues to TIPS-pentacene using dynamic disorder from DFTB and transient localization theory. The approach we establish here could provide a truly rapid simulation pathway for accurate materials properties prediction, in our vision applied to new OSCs with tailored properties.

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

Elucidating correlated defects in metal organic frameworks using theory-guided inelastic neutron scattering spectroscopy

Metal organic frameworks (MOFs) that incorporate metal oxide cluster nodes, exemplified by UiO-66, have been widely studied, especially in terms of their deviations from the ideal, defect-free crystalline structures. Although defects such as missing linkers, missing nodes, and the presence of adventitious synthesis-derived node ligands (such as acetates and formates) have been proposed, their exact structures remain unknown. Previously, it was demonstrated that defects are correlated and span multiple unit cells. The highly specialized techniques used in these studies are not easily applicable to other MOFs. Thus, there is a need to develop new experimental and computational approaches to understand the structure and properties of defects in a wider variety of MOFs. Here, we show how low-frequency phonon modes measured by inelastic neutron scattering (INS) spectroscopy can be combined with density functional theory (DFT) simulations to provide unprecedented insights into the defect structure of UiO-66. We are able to identify and assign peaks in the fingerprint region (<100 cm –1 ) which correspond to phonon modes only present in certain defective topologies. Specifically, this analysis suggests that our sample of UiO-66 consists of predominantly defect-free fcu regions with smaller domains corresponding to a defective bcu topology with 4 and 2 acetate ligands bound to the Zr 6 O 8 nodes. Importantly, the INS/DFT approach provides detailed structural insights (e.g., relative positions and numbers of acetate ligands) that are not accessible with microscopy-based techniques. Furthermore, the quantitative agreement between DFT simulations and the experimental INS spectrum combined with the relative simplicity of sample preparation, suggests that this methodology may become part of the standard and preferred protocol for the characterization of MOFs, and, in particular, for elucidating the structure defects in these materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Disentangling the role of bond lengths and orbital symmetries in controlling $T_c$ of optimally doped YBa 2 Cu 3 O 7

Optimally doped YBa 2 Cu 3 O 7 (YBCO) has a high critical temperature, at 92 K. It is largely believed that Cooper pairs form in YBCO and other cuprates because of spin fluctuations, but the issue and the detailed mechanism are far from settled. In the present work, we employ a state-of-the-art first-principles ability to compute both the low- and high-energy spin fluctuations in optimally doped YBCO. We benchmark our results against recent inelastic neutron scattering and resonant inelastic x-ray scattering measurements. Further, we use strain as an external parameter to modulate the spin fluctuations and superconductivity. We disentangle the roles of barium-apical oxygen hybridization, interlayer coupling, and orbital symmetries by applying an idealized strain, and also a strain with a fully relaxed structure. We show that shortening the distance between Cu layers is conducive to enhanced Fermi surface nesting, which increases spin fluctuations and drives up $T_c$. However, when the structure is fully relaxed, electrons flow to the $d_{z^2}$ orbital as a consequence of a shortened Ba-O bond, which is detrimental for superconductivity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

braggedgemodeling

braggedgemodeling (bem) is an open-source Python package for modeling neutron Bragg-edge imaging. It computes the wavelength-dependent total neutron cross-section of a material (coherent and incoherent elastic, coherent and incoherent inelastic scattering, and absorption) from its crystal structure, and implements the March-Dollase texture model and the Jorgensen peak profile, supporting quantitative analysis of energy-resolved neutron imaging data (phase, stress/strain, and texture). Published in the Journal of Open Source Software (2018).

Lin, Jiao [Oak Ridge National Laboratory (ORNL), O↗
DOE OSTI2024

Recombination of 𝐵 𝑐 mesons in ultrarelativistic heavy-ion collisions

High-energy heavy-ion collisions have been suggested as a favorable environment for the production of 𝐵 𝑐 mesons, due to a much larger abundance of charm and bottom quarks compared to elementary reactions. Motivated by recent CMS data for 𝐵$^{+}_{𝑐}$ production in Pb-Pb(5.02 TeV) collisions at the Large Hadron Collider (LHC), we deploy a previously developed transport approach for charmonia and bottomonia to evaluate the kinetics of 𝐵 𝑐 mesons throughout the fireball formed in these reactions. The main inputs to our approach are two transport parameters: the 𝐵 𝑐 's reaction rate and equilibrium limit. Both quantities are determined by previous calculations via a combination of charm and bottom sectors. In-medium binding energies of 𝐵 𝑐 mesons are calculated from a thermodynamic 𝑇 matrix with a lattice-QCD constrained potential, and figure in their inelastic reaction rates. Temperature-dependent equilibrium limits include charm- and bottom-quark fugacities based on their initial production. We compute the centrality dependence of inclusive 𝐵 𝑐 production and transverse-momentum (𝑝 𝑇 ) spectra using two different recombination models: instantaneous coalescence and resonance recombination. The main uncertainty in the resulting nuclear modification factors, 𝑅 AA , is currently associated with the 𝐵 𝑐 cross section in elementary 𝑝⁢𝑝 collisions, caused by the uncertainty in the branching ratio for the 𝐵$^{−}_{𝑐}$ → 𝐽/𝜓⁢𝜇 − ⁢$\overline{𝜈}$ decay. Our results indicate a large enhancement of the 𝑅 AA at low 𝑝 𝑇 , with significant regeneration contributions up to 𝑝 𝑇 ≃ 20GeV. Comparisons to CMS data are carried out but firm conclusions will require a more accurate value of the branching ratio, or alternative channels to measure the 𝐵 𝑐 production in 𝑝⁢ collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2002

Nuclear model calculations and their role in space radiation research

Proper assessments of spacecraft shielding requirements and concomitant estimates of risk to spacecraft crews from energetic space radiation requires accurate, quantitative methods of characterizing the compositional changes in these radiation fields as they pass through thick absorbers. These quantitative methods are also needed for characterizing accelerator beams used in space radiobiology studies. Because of the impracticality/impossibility of measuring these altered radiation fields inside critical internal body organs of biological test specimens and humans, computational methods rather than direct measurements must be used. Since composition changes in the fields arise from nuclear interaction processes (elastic, inelastic and breakup), knowledge of the appropriate cross sections and spectra must be available. Experiments alone cannot provide the necessary cross section and secondary particle (neutron and charged particle) spectral data because of the large number of nuclear species and wide range of energies involved in space radiation research. Hence, nuclear models are needed. In this paper current methods of predicting total and absorption cross sections and secondary particle (neutrons and ions) yields and spectra for space radiation protection analyses are reviewed. Model shortcomings are discussed and future needs presented. c2002 COSPAR. Published by Elsevier Science Ltd. All right reserved.

NASA Center JSC↗
NASA NTRS1972

Production cross sections for Lee-Wick massive electromagnetic bosons and for spin-zero and spin-one W bosons at high energies.

Production cross sections for three types of hypothetical particles are calculated in the presented paper. Several (Z, Z') cases were studied corresponding to elastic scattering off protons and neutrons (either free or embedded within a Fermi sea), coherent scattering off a nucleus, and inelastic scattering off a proton (in which case Z' denotes a nucleon resonance or hadronic system in the continuum). Detailed structure-function data are used to improve the accuracy of the inelastic scattering calculation. Results of calculations are given for beam energies between 50 and 10,000 GeV, and masses between 5 and 40 GeV for the massive Lee-Wick spin-1 boson. Cross sections were computed for resonant and semiweak processes. The production cross section of spin-zero weak intermediate bosons was found to be at least one order of magnitude smaller than for spin-1 weak bosons in nearly all regions of interest. The production cross section of spin-zero weak intermediate bosons for inelastic scattering off protons compares with that for elastic scattering in the regions of interest. In the case of massive spin-1 bosons and spin-1 weak intermediates, the main contribution to total production cross section off protons is elastic.

Linsker, R.↗
DOE OSTI2025

Open-source library for performance-portable neutrino reaction rates: Application to neutron star mergers

A realistic and detailed description of neutrinos in binary neutron star (BNS) mergers is essential to build reliable models of such systems. To this end, we present bns_nurates, a novel open-source numerical library designed for the efficient on-the-fly computation of neutrino interactions, with particular focus on regimes relevant to BNS mergers. bns_nurates targets a higher level of accuracy and realism in the implementation of commonly employed reactions by accounting for relevant microphysics effects on the interactions, such as weak magnetism and mean field effects. It also includes the contributions of inelastic neutrino scattering off electrons and positrons and (inverse) nucleon decays. Finally, it offers a way to reconstruct the neutrino distribution function in the framework of moment-based transport schemes. As a first application, we compute both energy-dependent and energy-integrated neutrino emissivities and opacities for conditions extracted from a BNS merger simulation with m1 transport scheme. We find some qualitative differences in the results when considering the impact of the additional relevant reactions and of microphysics effects. For example, neutrino-electron/positron scattering reactions are important for the energy exchange of heavy-type neutrinos as they do not undergo semileptonic charged-current processes, when μ± are not accounted for. Moreover, weak magnetism and mean field effects can significantly modify the contribution of β processes for electron-type (anti)neutrinos, increasing at the same time the importance of (inverse) neutron decays. Here, the improved treatment for the reaction rates also modifies the conditions at which neutrinos decouple from matter in the system, potentially affecting their emission spectra.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Accelerating resonant spectroscopy simulations using multishifted biconjugate gradient

Resonant spectroscopies, which involve intermediate states with finite lifetimes, provide important insights into collective excitations in quantum materials that are otherwise inaccessible. However, theoretical understanding in this area is often limited by the numerical challenges of solving Kramers-Heisenberg-type response functions for large-scale systems. To address this, we introduce a multishifted biconjugate gradient algorithm that exploits the shared structure of Krylov subspaces across spectra with varying incident energies, effectively reducing the computational complexity to that of linear spectroscopies. Both mathematical proofs and numerical benchmarks confirm that this algorithm substantially accelerates spectral simulations, achieving constant complexity independent of the number of incident energies, while ensuring accuracy and stability. This development provides a scalable, versatile framework for simulating advanced spectroscopies in quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Development of a MOOSE-based crystal plasticity model with irradiation defect evolution for irradiation creep in 316

Irradiation creep and irradiation swelling are two of the lifetime limiting factors for structural materials in nuclear reactors. These mechanical effects are driven by irradiation defect evolution and the interaction of those defects with dislocations in the microstructure. We present here a coupled cluster dynamics – crystal plasticity approach to model irradiation swelling and creep behavior in 316 SS. The time-dependent evolution of irradiation defects is calculated with a cluster dynamics approach and passed to the crystal plasticity model to compute the dislocation evolution. We show the impact of the irradiation defect evolution on the stress state in the material, which drives inelastic deformation through dislocation motion. The inelastic deformation in the 316 SS is dependent on the dose rate, where the inelastic deformation driven by the early-stage irradiation defect evolution determines the mechanical behavior of the 316 SS.

316 Stainless Steels↗
NASA NTRS1983

Advances and trends in structural and solid mechanics; Proceedings of the Symposium, Washington, DC, October 4-7, 1982

The mechanics of materials and material characterization are considered, taking into account micromechanics, the behavior of steel structures at elevated temperatures, and an anisotropic plasticity model for inelastic multiaxial cyclic deformation. Other topics explored are related to advances and trends in finite element technology, classical analytical techniques and their computer implementation, interactive computing and computational strategies for nonlinear problems, advances and trends in numerical analysis, database management systems and CAD/CAM, space structures and vehicle crashworthiness, beams, plates and fibrous composite structures, design-oriented analysis, artificial intelligence and optimization, contact problems, random waves, and lifetime prediction. Earthquake-resistant structures and other advanced structural applications are also discussed, giving attention to cumulative damage in steel structures subjected to earthquake ground motions, and a mixed domain analysis of nuclear containment structures using impulse functions.

Noor, A. K.↗