Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix elements”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Analytic techniques for solving the transport equations in electroweak baryogenesis

We develop an efficient method for solving transport equations, particularly in the context of electroweak baryogenesis. It provides fully-analytical results under mild approximations and can also test semi-analytical results, which are applicable in more general cases. Key elements of our method include the reduction of the second-order differential equations to first order, representing the set of coupled equations as a block matrix of the particle densities and their derivatives, identification of zero modes, and block decomposition of the matrix. We apply our method to calculate the baryon asymmetry of the Universe (BAU) in a Standard Model effective field theory framework of complex Yukawa couplings to determine the sensitivity of the resulting BAU to modifications of various model parameters and rates, and to estimate the effect of the commonly-used thin-wall approximation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

A high-performance implementation of atomistic spin dynamics simulations on x86 CPUs

Atomistic spin dynamics simulations provide valuable information about the energy spectrum of magnetic materials in different phases, allowing one to identify instabilities and the nature of their excitations. However, the time cost of evaluating the dynamical correlation function S(q,t) increases quadratically as the number of spins N, leading to significant computational effort, making the simulation of large spin systems very challenging. In this work, here, we propose to use a highly optimized general matrix multiply (GEMM) subroutine to calculate the dynamical spin-spin correlation function that can achieve near-optimal hardware utilization. Furthermore, we fuse the element-wise operations in the calculation of S(q,t) into the in-house GEMM kernel, which results in further performance improvements of 44\% - 71\% on several relatively large lattice sizes when compared to the implementation that uses the GEMM subroutine in OpenBLAS, which is the state-of-the-art open source library for Basic Linear Algebra Subroutine (BLAS).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Investigation of crossed-twin structure formation in magnesium and magnesium alloys

In this work, the effect of alloying addition on the propensity for twin-twin interactions to transform into crossed-twin structures in magnesium alloys is investigated. A full-field elasto-viscoplastic fast Fourier transform (EVP-FFT) framework combined with a discrete twin model and dislocation density-based hardening law for slip strengths is used to calculate the micromechanical fields in the crystalline matrix around the interacting twins. AZ31 and MgLi alloys are selected along with pure Mg to study the influence of plastic anisotropy in connection with alloying elements. These alloys were selected since their plastic anisotropy measure, which is defined as the ratio between the critical resolved shear stress for pyramidal $\langle$c +a$\rangle$ and basal $\langle$a$\rangle$ slip modes, spanned a wide range. To quantify the role of twin thicknesses, we probe a range of impinging twin thicknesses while fixing the recipient twin thickness. The analysis reveals that: (i) the local driving stress for crossed twin structure formation generated from the interaction of the two twins is lower in a low plastically anisotropic material, like a MgLi alloy, than a high plastically anisotropic material like pure Mg and (ii) the critical impinging twin thickness needed to form the crossed twin structure in pure Mg, AZ31 and MgLi alloys is ~0.5, ~0.75 and ~1.5 times the recipient twin thickness. We propose a relationship between the tendency for crossed twin structure formation and the experimentally observed higher ductility in MgLi alloys compared to pure Mg. Finally, one key implication of the findings is that crossed- twin structure formation can be hindered and the ductility of magnesium alloy thereby improved by properly choosing alloying elements that lower the slip strength for pyramidal $\langle$c + a$\rangle$ slip.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Shock compaction heating and collisional processes in the production of type 3 ordinary chondrites: Lessons from the (nearly) unique L3 chondrite melt breccia Northwest Africa 8709

Abstract Northwest Africa (NWA) 8709 is a rare example of a type 3 ordinary chondrite melt breccia and provides critical information for the shock compaction histories of chondrites. An L3 protolith for NWA 8709 is inferred on the basis of oxygen isotope composition, elemental composition, diverse mineral chemistry, and overall texture. NWA 8709 is among the most strongly shocked type 3 chondrites known, and experienced complete melting of the matrix and partial melting of chondrules. Unmelted phases underwent FeO reduction and partial homogenization, with reduction possibly occurring by reaction of olivine and low‐Ca pyroxene with an S‐bearing gas that was produced by vaporization. Chondrules and metal grains became foliated by uniaxial compaction, and during compression, chondrules and fragments became attached to form larger clumps. This process, and possibly also melt incorporation into chondrules to cause “inflation,” may have contributed to anomalously large chondrule sizes in NWA 8709. The melt breccia character is attributed to strong shock affecting a porous precursor. Data‐model comparisons suggest that a precursor with 23% porosity that was impacted by a 3 km/s projectile could have produced the meteorite. The rarity of other type 3 ordinary chondrite melt breccias implies that the immediate precursors to such chondrites were lower in porosity than the NWA 8709 precursor, or experienced weaker shocks. Altogether, the data imply a predominantly “quiet” dynamical environment to form most type 3 ordinary chondrites, with compaction occurring in a series of relatively weak shock events.

Ruzicka, Alex M.↗
DOE OSTI2024

Performance Improvements of the Griffin Solvers in FY24

The Griffin code is a MOOSE-based reactor physics application jointly developed by Idaho National Laboratory and Argonne National Laboratory under the Department of Energy Office of Nuclear Energy Nuclear Energy Advanced Modeling and Simulation Program. This fiscal year, we have made significant efforts to improve the performance of transport solver options and cross-section generation for the efficient use of Griffin in advanced reactor applications. For the HFEM-PN solver, the residual evaluations of HFEM kernels were optimized by utilizing the pre- computed averaged cross sections for individual elements. Numerical integration involving the evaluation of basis functions at quadrature points was bypassed by facilitating precomputed element mass matrices for response matrices. Red-black iterations were improved by introducing a new generalized minimum residual based solver. The memory usage of response matrix storage was significantly reduced by applying basis function rotations on interfaces and calculating volumetric odd-parity moments on the fly. Additionally, the adjoint flux and transient calculation capabilities of the HFEM-PN solver were successfully implemented and verified using the TWIGL benchmark problem. For the DFEM-SN solver, memory footprint and computation time were significantly reduced by not treating angular flux vectors as the MOOSE nonlinear system vectors. Specifically for IQS, scalar adjoint weighting was introduced to further eliminate angular adjoint flux storage in the MOOSE auxiliary system. It was demonstrated through the three-dimensional Advanced Burner Test Reactor core problem that the memory usage for transient calculations with the IQS method was reduced by over 7.5× compared to before the optimizations. For the self-shielding application programming interface, a new double-heterogeneity treatment method, named the Bell Function-Based Analytic Two-Region Slowing Down Method, was developed to efficiently flux-volume homogenize TRISO particles with the matrix. Additionally, optimizations were made to hyper- fine group (HFG) slowing down calculations by pretabulating collision probability coefficients and grouping isotopes, significantly reducing the computational time for calculating scattering sources per HFG. Lastly, the pin power reconstruction module was extended to account for temporal behavior in a microreactor analysis problem, specifically for a control drum transient. Verification tests for each of these improvements demonstrated significant performance enhancements and memory reduction.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Design and Development of Novel Equiatomic Refractory Multi Principal elemental Alloys Based on MoNbTi System for Use in Irradiation Environments

Multi Principal Elemental Alloys (MPEA) have emerged as promising materials for next-generation nuclear reactors due to their exceptional irradiation resistance. Eight equiatomic MPEA based on the MoNbTi system, comprising of elements with low thermal neutron absorption cross-sections were explored using a combined approach employing empirical parameter estimations, and CALPHAD simulations by which the phases and elemental segregation observed in all the alloys in the as-cast state were predicted. Solution heat treatment at 1500°C transformed five alloys into single-phase matrix materials, enhancing homogeneity and reducing hardness. The densities of the alloys ranged between 6.47 to 7.68 g/cm3, hardness between 472 and 656 VHN, Young’s modulus between 142 GPa to 169 GPa, shear modulus between 54 GPa and 62 GPa, bulk modulus between 117 GPa to 194 GPa and Poisson’s ratio between 0.3 to 0.35. The in-situ high temperature Xray diffraction results, differential scanning calorimetry and dilatometry results up to 1000°C suggested the high temperature phase stability of the MPEA. Subsequent ageing heat treatment at 800 and 1000 oC for 96 hours revealed significant secondary phase precipitation in MoNbTiZr, MoNbTiZrV, and MoNbTiCrA. Oxidation studies at 800 oC for 24 hours in air revealed superior oxidation resistance and cubic rate law dependence in Cr containing MPEA especially in MoNbTiCrAl, while severe mass gain resulting in total disintegration and exfoliation in MoNbTiZr and MoNbTiZrV. This comprehensive study underscores the potential of novel MPEA as promising materials for advanced nuclear reactor applications, shedding light on their microstructural control, mechanical properties, thermal stability, and oxidation resistance.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Influence of sizing concentration on strength, stiffness, and porosity in textile grade carbon fiber (TCF)-Epoxy composites: Revealing inverse trends

The effect of fiber sizing (i.e., surface treatment) concentration (0 %, 1.36 %, 1.52 %, 1.94 %, and 2.13 %) on the mechanical properties (tensile, flexural, interlaminar shear strength (ILSS), and low velocity impact) of textile grade carbon fiber (TCF)-epoxy composite is examined. An inverse relationship between the strength and stiffness of the composite is observed with increased sizing concentration. The root mean square (RMS) roughness of the fiber surface increased from 17.8 nm (unsized) to 22.7 nm with 2.13 % sizing concentration. It was found that the tensile strength increased by 131 % from 221.4 ± 18.5 MPa (unsized) to 510.8 ± 28.05 MPa (for 1.36 % sizing) and further by 155 %–563.7 ± 14.95 MPa at 2.13 % sizing. On the contrary, the stiffness is initially increased by 126 % from 33.52 ± 7.80 GPa (unsized) to 75.9 ± 3.21 GPa (for 1.36 % sizing) but reduced with further increase in the sizing concentration. A single fiber pull-out test is simulated using the finite element method to validate the reverse trend in strength and stiffness. The varying sizing concentration is simulated by introducing an interface of varying thickness between fibers and matrix. Simulation results confirm that a thicker interface, corresponding to a higher sizing concentration, decreases interfacial shear stress, enhancing material strength while decreasing stiffness. The reverse trend in strength and stiffness with the sizing concentration aligns with experimental observations. In conclusion, the present study emphasizes the importance of sizing concentration for mechanical properties and provide a design criterion for customized high-strength and high-stiffness applications.

Porosity↗
DOE OSTI2020

Three-colour digital holographic microscopy based on a Mach – Zender interferometer and a solid-state RGB laser

We have developed a facility for recording colour digital holograms involving a Mach – Zender interferometer and a pulsed solid-state RGB laser, which radiates simultaneously at three wavelengths: 451, 532, and 634 nm. These wavelengths are obtained in the conversion of the radiation of a pulsed Nd{sup 3+} : YAG laser (1064 nm), which operates in the regime of intracavity parametric oscillation at a wavelength of 1570 nm, with the subsequent conversion of its frequencies to the red, green, and blue spectral regions in nonlinear KTP-crystal elements. The combination of this three-wavelength pulsed light source with the Mach – Zender interferometer provides the simultaneous recording of three spectral digital holograms with one matrix RGB radiation detector. Pulsed three-colour digital holography shows promise for the analysis of fast processes and rapidly varying colour phase objects. (paper)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Progress on Light-Element Standard Cross Sections at Los Alamos [Slides]

Presenters identified the following statuses of the light-element standards: Extension of N-N scattering analysis to 200 MeV will likely come in time for ENDF/B VIII.1 and new data are to be added even in the energy range up to 100 MeV; for n- 3 He, a more recent 4 He R-matrix analysis from which updated n+ 3 He cross sections could be obtained for ENDF/B VIII.1; for n- 6 Li, more work is required to get a reasonable fit to all data up to E n = 8 MeV. Of particular concern are the t+ 4 He scattering data at E t > 14 MeV. For n- 10 B, new work on the 11 B system has not yet been started and there is much new data to add for n+ 10 B reactions, and maybe also for &#945 + 7 Li reactions. The proliferation of excited-state channels will limit the upper neutron energy of the analysis to a few MeV. For n- 12 C, new elastic and inelastic scattering data from Vanhoy support the results of the present analysis. More work is needed to extend the fit successfully above 6.5 MeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

A weighted Shifted Boundary Method for free surface flow problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods and was recently introduced for the Poisson, linear advection/diffusion, Stokes, Navier-Stokes, acoustics, and shallow-water equations. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we extend the SBM to the simulation of incompressible Navier-Stokes flows with moving free-surfaces, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach prevents spurious pressure oscillations in time, which would otherwise be produced if the total active fluid volume were to change abruptly over a time step. In fact, the proposed weighted SBM method induces small mass (i.e., volume) conservation errors, which converge quadratically in the case of piecewise-linear finite element interpolations, as the grid is refined. Finally, we present an extensive set of two- and three-dimensional tests to demonstrate the robustness and accuracy of the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A Poromechanical Framework for Internal Interactions Induced by Solid Inclusions

The framework of poromechanics is generalized to simulate the multiscale behavior of porous media subjected to internal loadings stemming from the growth of solid inclusions. This generalization is designed to enable the study of anisotropic internal stress generation from solid growth within the pores, while recovering isotropic fluid-induced loading as a particular case. For this purpose, a mathematical strategy to define constitutive tensors in a thermodynamically consistent form is proposed, thus offering new opportunities for determining the poromechanical properties of a porous solid through advanced experimentation or micromechanical models. The framework is specialized by means of established elastic solutions for single pore–matrix interaction, as well as through homogenization schemes considering the interaction among congruent pores. In particular, the second Eshelby solution and the Tanaka–Mori–Benveniste homogenization scheme are used to derive a microporoelastic model. At an elemental scale, the model is tested under mixed control conditions by replicating different scenarios of geomaterial testing. In addition, the model characteristics are outlined with reference to inelastic microscopic loadings replicating chemo-mechanical forcing, such as expansive crystal formation. Through a series of parametric analyses, it is shown that the microstructure of the pores significantly influences the properties of porous media. Most notably, it is shown that the effects of a solid forming within the pores depend in a highly nonlinear fashion on the constitutive characteristics of the inhomogeneities and can therefore not be readily quantified or predicted without models capturing the diverse multiscale interactions among pores, inhomogeneities, and matrix.

chemomechanics↗
DOE OSTI2021

Unraveling Cu Chemical Signature in CdTe by Spectral Fluorescence Mapping

X-ray absorption spectroscopy (XAS) has been shown to be a powerful tool to unravel the chemical environment of a given atom within a matrix. When used in correlative X-ray microscopy approaches, XAS allows one to probe with nanoscale precision regions of particular interest in an absorber. Herein, we use X-ray absorption near edge structure (XANES) to evaluate the chemical environment of Cu atoms within a CdTe solar cell. The reconstruction of XANES spectra from XRF maps have unfolded 2D maps of Cu chemical structures. In this work, we found that most Cu atoms exist in Cu 2 Te and Cu 1.4 Te phase. Moreover, we found traces of CuTe, Cu 2 O, CuO, Cu 2 S, CuS, and metallic Cu phase. Investigating Cu chemical structures at different performing areas, we found no observable correlation between Cu chemical structures and electrical performance. This approach allows tracking of Cu chemical structures along with electrical performance and elemental distribution simultaneously, with high spatial resolution in a statistically practical way.

CdTe↗
DOE OSTI2024

A high-order Shifted Interface Method for Lagrangian shock hydrodynamics

Here, we present a new method for two-material Lagrangian hydrodynamics, which combines the Shifted Interface Method (SIM) with a high-order Finite Element Method. Our approach relies on an exact (or sharp) material interface representation, that is, it uses the precise location of the material interface. The interface is represented by the zero level-set of a continuous high-order finite element function that moves with the material velocity. This strategy allows to evolve curved material interfaces inside curved elements. By reformulating the original interface problem over a surrogate (approximate) interface, located in proximity of the true interface, the SIM avoids cut cells and the associated problematic issues regarding implementation, numerical stability, and matrix conditioning. Accuracy is maintained by modifying the original interface conditions using Taylor expansions. We demonstrate the performance of the proposed algorithms on established numerical benchmarks in one, two and three dimensions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Evaluation of Light-Element Reactions in the Resolved Resonance Region

Light-element reactions at low energies in the resolved resonance region are important for a range of applications in basic and applied sciences including nuclear reactors, nonproliferation, cultural heritage, forensics and environmental control, rare event investigations and nuclear astrophysics. In this paper, we report on an effort to evaluate charged-particle cross sections in the resolved resonance region and produce evaluated nuclear data files for further processing and inclusion in evaluated data libraries. We discuss the open issues in R-matrix calculations as we extend to higher energies, such as dealing with the rapidly growing number of open channels and merging with the regime of smooth cross sections described by the statistical model, and present attempts to address these issues in neutron-induced reactions relevant to nuclear reactor applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

A Physically-Based Model for Thermo-Oxidative and Hydrolytic Aging of Elastomers

A computationally efficient model is proposed to capture the loss of mechanical performance due to chemical aging that are formed as the competition of chain scission and cross-link formation/dissolution, such as thermo-oxidative aging or hydrolytic aging. The model should be considered an extension of our recent models [1, 2, 3] which further simplifies the matrix behavior based on the assumption of independence of environmental and mechanical damage. The model uses this assumption to reduce the necessary material parameters needed to model constitutive and inelastic behavior of elastomers during aging. To this end, the model can provide accurate predictions of the material performance with the significantly fewer number of fitting parameters. The model is relevant for all decay mechanisms formed by the occurrence of two simultaneous micro-mechanisms; (i) formation/reduction of the cross-links, and (ii) chain scission, both of which are present in thermo-oxidation and hydrolytic aging. Assuming the alteration of the chain density along the aging trajectory is identical to the peroxide cross-link density for thermo-oxidation, and the change of the average molecular weight for hydrolysis, the strain energy of polymer matrix can be rewritten as a function of deformation, deformation history, storage time and aging temperature. Next, the modified network alteration model is formulated for implementation into Finite Element (FE) simulations. The model is built on the presumption of homogeneous and consistent oxygen/water absorption and thus is mainly relevant for relatively thin samples exposed to environmental loads for a long time. The proposed model includes only six physically inspired material parameters. Therefore, while it is computationally efficient, it shows good agreement with own experimental data, which performed on various range of accelerated aging temperatures and times. With respect to its computational efficiency, simplicity, accuracy, and interpret-ability, the model is the right choice for advanced implementations in FE programs.

42 ENGINEERING↗
DOE OSTI2020

Formation of multicomponent alloy particles in doped ceria under I 2+ ion irradiation and thermal annealing

Multicomponent alloys resembling the e-phase particles in nuclear fuel have the potential to immobilize 99 Tc in a stable waste form. Formation of such multicomponent alloys under extreme conditions is studied in a non-activated model system: CeO 2 doped with 2 wt% Mo, 1.5 wt% Ru, 0.75 wt% Pd, 0.5 wt% Re and 0.25 wt% Rh. The doped CeO 2 films were irradiated with I 2+ ions (610°C, 10 16 and 5×10 16 I 2+ /cm 2 ). For selected samples post-irradiation heat treatment was conducted (900°C, 1100°C). Analytical TEM revealed the formation of 5-20 nm precipitates containing Pd, Mo, and Re sometimes associated with cavities. Iodine was found to segregate to grain boundaries and cracks. After heat treatment at 1100°C, the CeO 2 matrix had recrystallized. Precipitates containing Pd, Mo, Re, and Ru were observed near the interface with the polycrystalline YSZ substrate. It follows that Pd is the most mobile element, followed by Mo and Re. Finally, while irradiation promotes precipitation, high-temperature effects are more significant.

36 MATERIALS SCIENCE↗
DOE OSTI2022

237 Np Mössbauer Isomer Shifts: A Lesson About the Balance of Static and Dynamic Electron Correlation in Heavy Element Complexes

A large set of neptunium compounds with different oxidation states (III to VII) was assembled to study the Mössbauer isomer shift by wave function calculations and better understand covalency in f-elements complexes. The contact density approach was used to calculate the isomer shift using complete active space self-consistent field (CASSCF) multiconfiguration wave functions, as well as matrix product states [from Density Matrix Renormalization Group (DMRG) algorithms] for large active spaces. Dynamic correlation effects for the isomer shifts were treated via CASPT2 energy derivatives with respect to the nuclear radius. The CASSCF calculations appear to produce different orbital overlocalization errors for low and high Np oxidation states. For compounds with low Np oxidation numbers, the errors can be attributed to the overlocalization of the 5f orbitals. For the compounds with high Np oxidation numbers, the main errors arise from an overlocalization of ligand orbitals and concomitant to weak donation bonding. Attempts to mitigate the overlocalization errors with large active spaces using DMRG were only partially successful, showing that explicit treatment of dynamic correlation is necessary for accurate predictions of Mössbauer isomer shifts. The CASPT2 calculations perform very satisfactorily. For a subset of Np compounds, both static and dynamic correlation effects were substantial. Furthermore, a rational active space selection based on orbital entanglement diagrams proved beneficial for determining the optimal reference wave function.

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

Development of Direct Injection/Ionization Mass Spectrometry Methods for Whole Molecule Characterization

The objective of this work is to adapt ambient ionization mass spectrometry (AMS) techniques for the rapid analysis of intact uranium complexes, stable strontium, cerium, and explosive compounds. The methods used were “soft ionization” techniques, which facilitate the detection of whole molecule complexes. The soft ionization mass spectrometry (MS) techniques that were investigated include paper spray ionization (PSI), matrix-assisted ionization (MAI), electrospray ionization (ESI) and direct analysis in real time (DART). For the first time, PSI-MS methods were successfully developed for whole molecule uranium-containing analytes (uranyl acetate, uranyl nitrate, and uranyl-tributylphosphate complexes). This was also the first demonstration of uranium complex detection and characterization and one of the few examples of inorganic analysis using MAI techniques. Proof of concept experiments also putatively identified matrix-derived ions and ion complexes that have not previously been described in the literature. Additionally, PSI-MS on cotton swipe samples doped with a multi-element standard containing μg levels of U, Bi, Pb, Cd, Fe, and Zn were directly analyzed without purification, representing a major improvement over existing methods. Both PSI and MAI methods demonstrated limits of detection (LODs) in the 10 - 100’s ng for various uranyl species within a range of 10’s ppm - 100’s ppb, dependent on analytical method and analyte species. AMS methods were also developed for other inorganics, including Ce and Sr, and organic explosive residues to address specific challenges in environmental monitoring and forensics. Further refinement and qualification of the AMS techniques developed within this effort would lead to significant cost reduction and timeliness by facilitating the triage and queueing of samples for subsequent more sensitive and time-consuming analyses.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗