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 109 records · Page 6

The Cosmic Neutrino Background is within Reach of Future Neutrino Telescopes

The cosmic neutrino background (C$ν$B) can be boosted to high energies due to scatterings with energetic cosmic rays (CRs) across cosmological scales. Previous calculations focused on neutral current incoherent and coherent elastic scatterings of cosmic-ray protons off relic neutrinos. However, charged current interactions and deep inelastic scatterings are also expected to occur, which enhances the boosted relic neutrino fluxes on Earth. Here, we compute the \textit{total} diffuse boosted cosmic neutrino background (DBC$ν$B) arising from CRs at all redshifts in the Universe, accounting for neutral current and charged current elastic and deep inelastic scatterings. We find that IceCube already places an upper limit on the cosmic neutrino background overdensity in cosmological scales of ~$\mathcal{O}(100-1000)$ at $E_ν=10^{10}$ GeV, for a lightest neutrino mass of $m_ν \gtrsim 0.1$ eV. We further show that IceCube-Gen2 could test $\mathcal{O}(1-10)$ C$ν$B overdensities, and the combination of $10$ future neutrino telescopes with similar sensitivity would allow us to test the $Λ$CDM expected C$ν$B density for a lightest neutrino mass compatible with the KATRIN bound.

Herrera, Gonzalo [MIT, MKI; Harvard U.; Virginia T↗

Comparison of the Computational Efficiency of the Original Versus Reformulated High-Fidelity Generalized Method of Cells

The High-Fidelity Generalized Method of Cells (HFGMC) micromechanics model has recently been reformulated by Bansal and Pindera (in the context of elastic phases with perfect bonding) to maximize its computational efficiency. This reformulated version of HFGMC has now been extended to include both inelastic phases and imperfect fiber-matrix bonding. The present paper presents an overview of the HFGMC theory in both its original and reformulated forms and a comparison of the results of the two implementations. The objective is to establish the correlation between the two HFGMC formulations and document the improved efficiency offered by the reformulation. The results compare the macro and micro scale predictions of the continuous reinforcement (doubly-periodic) and discontinuous reinforcement (triply-periodic) versions of both formulations into the inelastic regime, and, in the case of the discontinuous reinforcement version, with both perfect and weak interfacial bonding. The results demonstrate that identical predictions are obtained using either the original or reformulated implementations of HFGMC aside from small numerical differences in the inelastic regime due to the different implementation schemes used for the inelastic terms present in the two formulations. Finally, a direct comparison of execution times is presented for the original formulation and reformulation code implementations. It is shown that as the discretization employed in representing the composite repeating unit cell becomes increasingly refined (requiring a larger number of sub-volumes), the reformulated implementation becomes significantly (approximately an order of magnitude at best) more computationally efficient in both the continuous reinforcement (doubly-periodic) and discontinuous reinforcement (triply-periodic) cases.

Arnold, Steven M↗

Iterative solution of the multistream electron transport equation. I - Comparison with laboratory beam injection experiments

The Neumann iteration method presently used for solving the electron transport equation in which energy, attitude, and pitch angle are independent variables is fast, and can compute numerical point-response-function solutions of the electron transport equation. Because both the inelastic cross sections and angular elastic cross sections of the model are empirically based, the solutions obtained represent a test of compatibility between various sets of cross sections and energy deposition measurements. The use of a numerical quadrature based on analytic phase function forms yields accurate phase function integrals at low computational cost.

Porter, H. S.↗

String breaking in the heavy quark limit with scalable circuits

Quantum simulations of non-Abelian gauge theories require efficient mappings onto quantum computers and practical state preparation and measurement procedures. A truncation of the Hilbert space of non-Abelian lattice gauge theories with matter in the heavy quark limit is developed. This truncation is applied to SU(2) lattice gauge theory in 1+1D to map the theory efficiently onto a quantum computer. Scalable variational circuits are found to prepare the vacuum and single meson states. It is also shown how these state preparation circuits can be used to perform measurements of the number of mesons produced during the system’s time evolution. A state with a single qq¯ pair is prepared on quantum hardware and the inelastic production of qq¯ pairs is observed using 104 qubits on IBM’s Heron quantum computer ibm_torino.

Ciavarella, Anthony N↗

Perturbative QCD contribution to transverse single spin asymmetries in the Drell-Yan process and SIDIS

In a previous publication [Benić , ], we have computed the perturbative QCD contribution to transverse single spin asymmetries (SSAs) in semi-inclusive deep inelastic scattering (SIDIS) involving the g T ( x ) distribution. In this paper, we first present a more efficient derivation of the asymmetries which is applicable to both transverse and longitudinal SSAs and correct some inconsistencies in our previous calculation. We then adapt the method to compute transverse SSAs in the Drell-Yan process proportional to g T ( x ) and its gluonic counterpart and discuss the crossing symmetry between the results for SIDIS and the Drell-Yan process. Finally, we present numerical results for various asymmetries measurable at the EIC, RHIC, COMPASS, and Fermilab (SpinQuest), including also part of the genuine twist-three corrections to g T ( x ) from a recent global analysis. We find that the asymmetries can reach percent-level magnitude, if the kinematics predominantly probes the large- x (valence) region of the polarized proton, but remain at subpercent levels otherwise. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nonlinear structural and life analyses of a combustor liner

Three-dimensional, nonlinear, finite element structural analyses were performed for a simulated aircraft combustor liner specimen in order to assess the capability of nonlinear analyses using classical inelastic material models to represent the thermoplastic-creep response of the component. In addition, the computed stress-strain history at the critical location was input into life prediction methods in order to evaluate the ability of these procedures to predict crack initiation life. It is concluded that: (1) elastic analysis is adequate for obtaining strain range and critical location; (2) inelastic analyses did not accurately represent cyclic behavior of materials; and (3) none of the crack initiation life prediction methods were satisfactory.

Kaufman, A.↗

Multiple Concentric Cylinder Model (MCCM) user's guide

A user's guide for the computer program mccm.f is presented. The program is based on a recently developed solution methodology for the inelastic response of an arbitrarily layered, concentric cylinder assemblage under thermomechanical loading which is used to model the axisymmetric behavior of unidirectional metal matrix composites in the presence of various microstructural details. These details include the layered morphology of certain types of ceramic fibers, as well as multiple fiber/matrix interfacial layers recently proposed as a means of reducing fabrication-induced, and in-service, residual stress. The computer code allows efficient characterization and evaluation of new fibers and/or new coating systems on existing fibers with a minimum of effort, taking into account inelastic and temperature-dependent properties and different morphologies of the fiber and the interfacial region. It also facilitates efficient design of engineered interfaces for unidirectional metal matrix composites.

Williams, Todd O.↗

High-fidelity micromechanical modeling of the effects of defects on damage and creep behavior in single tow ceramic matrix composite

Despite the superior properties of ceramic matrix composites (CMCs), their fabrication process generates inevitable defects with high density that significantly influence material integrity and residual useful life. Recent efforts have focused on CMC property prediction and investigation of their different inelastic mechanisms, but very limited work exists on understanding the influence of defects on inelastic responses of CMCs. Here, we introduce a three-dimensional micromechanics computational framework that includes experimentally informed material microstructure and architectural variabilities to investigate CMC response in the presence of manufacturing-induced defects. A developed microstructure generation algorithm is used to generate the representative volume elements in the micromechanics models from complex material morphology information obtained from extensive characterization studies of C/SiNC and SiC/SiNC CMCs. A fracture mechanics-informed matrix damage model is reformulated, in which the model considers the growth of porosity and microcracks in the as-received material. A progressive fiber damage model as well as Orthotropic viscoplasticity creep formulation are also utilized for the study. This methodology is implemented within the micromechanics framework to investigate the complex temperature- and time-dependent load transfer and damage mechanisms of CMCs under operation loading conditions. The developed framework explores the influence of as-received defects on the damage and creep behavior in service conditions. It also provides new insights into the effect of size, shape, distribution, and location of these defects on material response. The framework is then calibrated with unidirectional CMC minicomposite results available in the literature.

36 MATERIALS SCIENCE↗

The importance of phonon spectrum measurements for accurate prediction of thermal scattering laws [Slides]

Although approximate phonon spectra (e.g. Debye model) are sufficient to reproduce integral quantities (e.g. total INS cross section) in reasonable agreement with the measured data, more accurate phonon spectra are needed to achieve agreement with the differential data. SANS component is essential for reproducing the total scattering cross section; this is the case for nanodiamonds and other materials where SANS components are present (such as porous materials like nuclear graphite). A format proposal for SANS data in ENDF or GNDS is needed. Therefore, differential (e.g. phonon spectrum) and integral (e.g. total cross section) measurements are necessary in the process of creating an ENDF evaluation. We have successfully developed a method by which ENDF evaluations can be directly compared to the relevant experimentally measured quantities. INS experiments at SNS, ILL, ISIS, and, in the future, ESS, as well as total cross section measurements at the VESUVIO instrument at ISIS, are accessible through user programs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Three-dimensional inelastic analysis for hot section components, BEST 3D code

The goal is the development of an alternative stress analysis tool, distinct from the finite element method, applicable to the engineering analysis of gas turbine engine structures. The boundary element method was selected for this development effort on the basis of its already demonstrated applicability to a variety of geometries and problem types characteristic of gas turbine engine components. Major features of the BEST3D computer program are described, and some of the significant developments carried out as part of the Inelastic Methods Contract are outlined.

Wilson, Raymond B.↗

Potential bill impacts of dynamic electricity pricing on California utility customers

The rapid growth of renewable generation is creating challenges for the California grid in the form of the “duck curve,” with increasingly steep ramping required for conventional generation resources in the morning and evening, and growing curtailment of solar resources in midday periods. Time-varying electricity tariffs have received considerable attention as a tool to address these challenges, with a renewed recent focus on the potential for dynamic tariffs that vary to reflect conditions on the grid in near-real time. Consideration of dynamic tariffs may raise concerns about the financial impact on utility customers, especially for those who have limited flexibility to modify their electricity consumption in response. Specific areas of concern include electricity bills, bill volatility, and equity implications related to cost shifting among customer groups. In this paper we leverage smart meter data for more than 400,000 California utility customers, spanning residential, commercial, industrial, and agricultural customers, to assess potential customer bill impacts arising from a multi-component dynamic tariff . Specifically, we compute impacts on customer bills and bill volatility under the assumption of fully inelastic demand, i.e., where customers do not change their consumption patterns in response to the tariff. We also assess various approaches designing subscription load shapes that customers can pre-purchase as a hedge that may provide a measure of protection against large negative impacts, while still incentivizing the modification of loads on the margin. We compare and contrast the relative impacts on different customer classes and discuss benefits and pitfalls of different dynamic tariff structures and subscription load shapes.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Advances in 3-D Inelastic Analysis Methods for hot section components

3-D Inelastic Analysis Methods are described. These methods consist of a series of new computer codes embodying a progression of mathematical models (mechanics of materials, specialty finite element, boundary element) for streamlined analysis of: (1) combustor liners, (2) turbine blades, and (3) turbine vanes. These models address the effects of high temperatures and thermal/mechanical loadings on the local (stress/strain) and global (dynamics, buckling) structural behavior of the three selected components. Three computer codes, referred to as MOMM (Mechanics of Materials Model), MHOST (MARC-Hot Section Technology), and BEST (Boundary Element Stress Technology), have been developed and are briefly described in this paper.

Chamis, Christos C.↗

User's manual to the ACTION computer code

The form and interpretation of input and output data are defined. The Analysis of Crash Transients in Inelastic and Geometrically Nonlinear structures program (ACTION) performs nonlinear transient response analysis of structures subjected to time varying loads, allowing for nonlinear, time independent material properties and large geometry changes.

Kamat, M. P.↗

Inelastic, nonlinear analysis of stiffened shells of revolution by numerical integration

This paper describes the latest addition to the STARS system of computer programs, STARS-2P, for the plastic, large deflection analysis of axisymmetrically loaded shells of revolution. The STARS system uses a numerical integration scheme to solve the governing differential equations. Several unique features for shell of revolution programs that are included in the STARS-2P program are described. These include orthotropic nonlinear kinematic hardening theory, a variety of shell wall cross sections and discrete ring stiffeners, cyclic and nonproportional mechanical and thermal loading capability, the coupled axisymmetric large deflection elasto-plastic torsion problem, an extensive restart option, arbitrary branching capability, and the provision for the inelastic treatment of smeared stiffeners, isogrid, and waffle wall constructions. To affirm the validity of the results, comparisons with available theoretical and experimental data are presented.

Levine, H. S.↗

Scaling the SciDAC QuantOm Workflow

As part of the Scientific Discovery through Advanced Computing (SciDAC) program, the Quantum Chromodynamics Nuclear Tomography (QuantOM) project aims to analyze data from Deep Inelastic Scattering (DIS) experiments conducted at Jefferson Lab and the upcoming Electron Ion Collider. The DIS data analysis is performed on an event-level by combining the input from theoretical and experimental nuclear physics into a single, composable workflow. The optimization itself (I.e. fitting the experimental data with theoretical predictions) is carried out by a machine / deep learning algorithm. The size of the acquired DIS data as well as the complexity of the workflow itself require that the analysis is performed across multiple GPUs on high performance computing systems, such as Polaris at Argonne National Laboratory. This presentation discusses the novelties and challenges that came along with parallelizing this workflow. Recent results are compared to common distributed training techniques.

Lersch, Daniel↗

Data-Driven Constitutive Model for the Inelastic Response of Metals: Application to 316H Steel

Here, predictions of the mechanical response of structural elements are conditioned by the accuracy of constitutive models used at the engineering length-scale. In this regard, a prospect of mechanistic crystal-plasticity-based constitutive models is that they could be used for extrapolation beyond regimes in which they are calibrated. However, their use for assessing the performance of a component is computationally onerous. To address this limitation, a new approach is proposed whereby a surrogate constitutive model (SM) of the inelastic response of 316H steel is derived from a mechanistic crystal plasticity-based polycrystal model tracking the evolution of dislocation densities on all slip systems. The latter is used to generate a database of the expected plastic response and dislocation content evolution associated with several instances of creep loading. From the database, a SM is developed. It relies on the use of orthogonal polynomial regression to describe the evolution of the dislocation content. The SM is then validated against predictions of the dead load creep response given by the polycrystal model across a range of temperatures and stresses. When the SM is used to predict the response of 316H during complex non monotonic loading, extrapolating to new loading conditions, it is found that predictions compare particularly well against those from the physics-based polycrystal model.

36 MATERIALS SCIENCE↗

Direct quarkonium-plus-gluon production in DIS in the color glass condensate

We compute the differential cross-section for direct quarkonium production accompanied by a gluon in high-energy deep inelastic scattering (DIS) at small- x . We employ the Non-Relativistic QCD factorization framework, focusing on the S -wave contribution to the formation of the quarkonium, and including both color singlet and octet contributions. Our short distance coefficients for the production of the heavy quark pair are obtained within the Color Glass Condensate effective field theory. Our results pave the way towards the next-to-leading order computation of direct quarkonium in DIS, as well as the study of azimuthal correlations of direct quarkonium and jet.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Inclusive prompt photon-jet correlations as a probe of gluon saturation in electron-nucleus scattering at small x

We compute the differential cross-section for inclusive prompt photon+quark production in deeply inelastic scattering of electrons off nuclei at small x (e + A DIS) in the framework of the Color Glass Condensate effective field theory. The result is expressed as a convolution of the leading order (in the strong coupling α s ) impact factor for the process and universal dipole matrix elements, in the limit of hard photon transverse momentum relative to the nuclear saturation scale Q s,A (x). We perform a numerical study of this process for the kinematics of the Electron-Ion Collider (EIC), exploring in particular the azimuthal angle correlations between the final state photon and quark. We observe a systematic suppression and broadening pattern of the back-to-back peak in the relative azimuthal angle distribution, as the saturation scale is increased by replacing proton targets with gold nuclei. Our results suggest that photon+jet final states in inclusive e + A DIS at high energies are in general a promising channel for exploring gluon saturation that is complementary to inclusive and diffractive dijet production. They also provide a sensitive empirical test of the universality of dipole matrix elements when compared to identical measurements in proton-nucleus collisions. However because photon+jet correlations at small x in EIC kinematics require jet reconstruction at small k ⊥ , it will be important to study their feasibility relative to photon-hadron correlations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗