SEARCH · Engineering Papers

Results for “Matrix”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,441 records · Page 80

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Predictive Tools for Customizing Heat Treatment of Additively Manufactured Aerospace Components

Laser-bed powder fusion (LBPF) additive manufacturing is increasingly being used to produce components of complex geometries using the Ni-base superalloy Inconel 718. The composition and the microstructure of the alloy are currently well optimized for wrought components made using conventional manufacturing processes such as rolling, forging, extrusion, etc. The attractive mechanical properties of the alloy result from the underlying austenitic matrix with fine equiaxed grains, and a high density and uniform distribution of the precipitation hardening phase, γ". Heat treatment steps such as homogenization, solutioning and aging are well documented for the wrought alloy. However, when the same wrought alloy compositions are used for the additive manufacturing (AM) processes, the asprocessed microstructure is significantly different, because of the different thermal history associated with LBPF, including rapid solidification and multiple temperature excursions that lead to multiple re-melting and reheating in the solid state. Rapid solidification introduces potential non-equilibrium effects at the moving solid-liquid interfaces that impact the extent of solute segregation, as well as the morphology of the dendritic grains that form. In order to recover the target mechanical properties, AM components have to undergo post-process heat treatments. However, such heat treatments have to be custom designed for the AM process and the component geometry because of the expected vast differences in the microstructure at various locations of a component with complex geometry. The homogenization and precipitation steps should be optimized for the component so that target mechanical properties can be obtained throughout the part. The objective of this research is to utilize High Performance Computing in phase field simulations of microstructure evolution during post-processing of AM components. The physics-based modeling will be beneficial in reducing the experimental effort required for heat treatment process selection, optimization, and certification, thus leading to a significant reduction in energy consumption for AM and post-processing heat treatment. The optimization study will help identify heat treatments steps that are critical for development of a final desired microstructure with the minimum energy input. This combined with shortening of the production cycle (time-to-market) by reducing the number of failed parts (property targets), and reduction in the number of iterations for process optimization, will enable 30-40% savings in the energy costs. Phase field simulations of the degree of homogenization and the effect of local matrix composition on the nucleation and growth of competing precipitating phases were performed using the Microstructure Evolution Using Massively Parallel Phase Field Simulations code developed in-house at the Oak Ridge National Laboratory. The simulations were able to successfully capture the kinetics of nucleation and growth, and morphologies of various precipitating phases as a function of local matrix compositions and composition gradients characteristic of local microstructures arising from location-dependent variations in the thermal conditions. Future work will involve extending the simulations to a length scale consisting of multiple dendrites, so that the effect of homogenization on the coarsening of the dendrites can be simulated and used as an additional input to the optimization of the heat treatment process.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Composing preconditioners for multiphysics PDE systems with applications to Generalized MHD

New patch smoothers or relaxation techniques are developed for solving linear matrix equations coming from systems of discretized partial differential equations (PDEs). One key linear solver challenge for many PDE systems arises when the resulting discretization matrix has a near null space that has a large dimension, which can occur in generalized magnetohydrodynamic (GMHD) systems. Patch-based relaxation is highly effective for problems when the null space can be spanned by a basis of locally supported vectors. The patch-based relaxation methods that we develop can be used either within an algebraic multigrid (AMG) hierarchy or as stand-alone preconditioners. These patch-based relaxation techniques are a form of well-known overlapping Schwarz methods where the computational domain is covered with a series of overlapping sub-domains (or patches). Patch relaxation then corresponds to solving a set of independent linear systems associated with each patch. In the context of GMHD, we also reformulate the underlying discrete representation used to generate a suitable set of matrix equations. In general, deriving a discretization that accurately approximates the curl operator and the Hall term while also producing linear systems with physically meaningful near null space properties can be challenging. Unfortunately, many natural discretization choices lead to a near null space that includes non-physical oscillatory modes and where it is not possible to span the near null space with a minimal set of locally supported basis vectors. Further discretization research is needed to understand the resulting trade-offs between accuracy, stability, and ease in solving the associated linear systems.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Radial Deconsolidation and Leach-Burn-Leach of AGR-3/4 Compacts 1-3, 4-3, 10-1, and 10-2

The Advanced Gas Reactor (AGR) Fuel Development and Qualification Program’s third and fourth irradiation experiments (AGR-3/4), originally planned as separate tests, were combined into one test train for irradiation in the Advanced Test Reactor at Idaho National Laboratory (INL). The irradiation test began on December 14, 2011 and ended on April 12, 2014 (Collin 2016). The originally planned AGR-3 and AGR-4 irradiation experiments were both focused on obtaining data on fission product transport to support modeling improvement. The AGR-3 experimental plan was focused on gaseous and metallic fission product release from the kernels and diffusion in the coatings during irradiation and post-irradiation safety testing. The AGR-4 experimental plan was focused on diffusivities and sorptivities in the compact matrix and reactor graphite (Petti et al. 2005). These two goals were combined in the AGR-3/4 irradiation, which consisted of twelve independently monitored capsules that each contained four AGR-3/4 compacts in a single stack surrounded by an inner ring of matrix or graphite and an outer ring of graphite. There were two capsule types: a standard capsule and a fuel body in which the outer graphite ring included a floor and cap that fully enclosed the fuel (Stempien et al. 2018a). The fuel body design supported post-irradiation safety testing of the intact fuel and ring assembly to provide data on fission product transport and release from matrix and graphite at accident temperatures (Demkowicz 2017)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Solid Phase Gradient Alloying Method via ShAPE

Alloying has been used to confer desirable properties to various metallic systems since the bronze age. For example, Ni, Mn, and Cr are alloyed into an Fe matrix to make steel, which has improved strength and corrosion resistance (as displayed in Fig. 1). In addition to the type of alloying elements, the amount of alloying elements added into a matrix is also very critical for a material’s mechanical property, physical property, machinability, and cost. Conventionally, new alloys are discovered by combining alloy precursors and elemental constituents into a mixture and melting them together. The solidified product or ingot represents a usually non-homogeneous combination of the chemistry with a microstructure dictated by the physics of solidification. It is time and energy intensive to optimize the alloying amounts via the casting process and hence, economically not conducive to rapid alloy discovery. In addition, the cast microstructure seen in the as-cast ingot is often not the ideal, nor even desired structure for optimized performance. Cast materials often need further processing such as homogenization, annealing, and deformation work put in through rolling, forging, extruding, etc. to have favorable microstructures and properties. A faster method to discover new alloy combinations and evaluate them in their worked or processed form is needed. Recently, Laser Engineered Net Shaping was applied to fabricate gradient compositional material for alloy designing. However, the mechanical properties of this gradient alloy were poor due to the existence of impurities, oxidation and cracks that are inherent in the melt-solidification process. Because of these limitations, we propose to use a solid-phase process (ShAPE) to create bulk materials (extrudates) that vary in composition from one end of the extruded solid to the other. Various manufacturing methods for alloy design including conventional casting method, laser based combinatorial method and current solid phase gradient alloying method are displayed and compared in Fig. 2. ShAPE machine and a schematic of ShAPE process are displayed in Fig. 3. Starting material is processed by a rotating and plunging die to form an extrudate. We are proposing through this methodology to invent a new, bulk-scale combinatorial technique. With this novel technique, alloying element can be dissolved into the matrix with a continued gradient without bulk melting. Formation of inter-metallics and defects originating from melt processing can be avoided. In addition, the ShAPE process creates the mixed alloy chemistry, and meanwhile subjects the material to severe plastic strain, which induces the favorable “worked” microstructure. Products from the ShAPE process can be tested in hardness directly providing a path to rapid evaluation of mechanical properties. The ability to create graded structures using friction stir processing (FSP) has been demonstrated, however using the ShAPE process we believe will lead to much faster and higher fidelity results. When combined with high throughput screening methods of physical, mechanical and microstructural property characterization, efficiency and accuracy of alloy design can be significantly improved.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Determination of Average Burnup in AGR-3/4 Compacts 1-4 and 7-4

The Advanced Gas Reactor (AGR) Fuel Development and Qualification Program third and fourth irradiation experiments (AGR-3/4), originally planned as separate tests, were combined in one test train for irradiation in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL). The irradiation test began on December 14, 2011, and ended on April 12, 2014 (Collin 2016). The originally planned AGR-3 and AGR-4 irradiation experiments were both focused on obtaining data on fission product transport to support the improvement of modeling. The AGR-3 experimental plan was focused on gaseous and metallic fission product release from the kernels and diffusion in the coatings during irradiation and post-irradiation safety testing. The AGR-4 experimental plan was focused on diffusivities and sorptivities in the compact matrix and reactor graphite. These two goals were combined in the AGR-3/4 irradiation, which consisted of 12 independently monitored capsules that each contained four AGR-3/4 compacts in a single stack surrounded by an inner ring of matrix or graphite and an outer ring of graphite. Two capsule types were used: a standard capsule and a so-called fuel body, in which the outer graphite ring included a floor and cap that fully enclosed the fuel. The fuel body design supported post-irradiation safety testing of the intact fuel and ring assembly to provide data on fission product transport and release from the matrix and graphite at accident temperatures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

From Fossil Fuels to Green Mines: Understanding Hydrogen-Shale Interactions for Secure Underground Storage

Hydrogen-powered mining trucks are rapidly gaining traction as a key driver of decarbonization in the mining industry. This surge necessitates a decisive shift towards renewable energy sources, primarily for hydrogen production. However, the inherent intermittency of renewables poses a challenge for continuous operations, demanding efficient storage solutions for hydrogen fuel. Hydrogen, derived from renewables, requires large-scale storage, with underground hydrogen storage (UHS) in depleted oil and gas reservoirs projected to provide this capacity. An impermeable seal (caprock) overlying depleted reservoirs is crucial for preventing leakage and ensuring secure storage. Shale, a potential caprock, has been recognized for its effectiveness in containing natural gas for geological time-periods. This study investigates the impact of hydrogen exposure on the matrix permeability of Marcellus shale, a key factor influencing its sealing capacity. A novel time-resolved approach was employed, measuring matrix permeability after each exposure to understand the dynamic interactions between hydrogen and the shale matrix. Samples were exposed to hydrogen in pressure reactions at realistic reservoir conditions at varying time durations. Preliminary results reveal an initial increase in permeability, while subsequent exposures show a decrease, indicating potential microstructural changes.

08 HYDROGEN↗

DOE OSTI2024

Sister Rod Destructive Examinations (FY23): Appendix J: Leaching of High Burnup Used Nuclear Fuel in Deionized Water

The leaching experiment aims to understand the trends in the radiolysis-enhanced dissolution of HBU SNF when exposed to water (e.g., in-reactor or in-pool cladding failures). Specimens from a baseline M5-clad rod and a heat-treated M5 rod were cut from the fractured CIRFT specimens and placed in 100 mL deionized water for a period of 128 days. Both radial and axial sections were cut to provide different surface areas of fuel in contact with the leachate. During the four-month exposure period, aliquot samples of the leachate were analyzed using gamma spectroscopy and inductively coupled plasma - mass spectrometry (ICP-MS). The analysis quantified the amount of fuel leached into the solutions and provided individual isotopic release fractions (of 30+ isotopes) which were compared as a function of time and surface area of the fuel exposed. Consistent with existing literature, the leaching followed a trend in which isotopes of certain elements such as Cs and Mo were among the first species in the matrix to dissolve, and with the highest release rates. This was followed by a gradual matrix dissolution consisting of uranium and other actinides and a slower-than-matrix release from some isotopes, including Ru and Rh. It was also observed that the circumferential samples having less exposed fuel surface area (to begin with) leached more than the axial samples for a majority of the isotopes during the timespan of the study. Previous experiments were conducted by decladding or exposing the pellet-clad surface completely, but the samples used in this study retained the cladding. One possible explanation for the higher leaching rate of the circumferential samples is that the pellet-clad interface, which has a greater density of grain boundaries and defects, may be the most vulnerable area to leaching of fuel in the presence of water.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2024

Efficient Scalable Contact Network Generation from Population Data

Modeling the contacts among a population is critical to understanding the dynamics of a disease outbreak. Contact networks, where nodes are individuals and edges are contacts among them, are used to represent these complex individual-level interactions. In this work, we are given the daily activity schedules of an urban population that represent the activity location and time of individuals in a population during a single twenty four hour period over multiple days. Using collocation to determine contact between individuals, our goal is to extract hourly contact networks from large-scale activity data. We improve upon the existing adjacency matrix-based method by implementing our custom sparse matrix multiplication algorithm. Starting with a Python implementation, we achieve a 1600x speed up in the computation with a fast custom designed sparse matrix multiplier algorithm implemented in the C++ language. This work is central to future parallel designs of the problem.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Upcycling of CFRP Waste: Viable Eco-friendly Chemical Recycling and Manufacturing of Novel Repairable and Recyclable Composites

The rapid growth of the carbon fiber-reinforced polymer (CFRP) composite market has driven researchers to find value-added applications for outdated prepregs, manufacturing scraps, and end-of-life components. Currently, most CFRP waste is incinerated or landfilled, which squanders its residual value and burdens the environment. Various mechanical, thermal, and chemical approaches have been attempted to recover the carbon fiber, polymer matrix, or both. However, these current practices are often costly, energy-intensive, and generate secondary waste and pollution. To address these shortcomings, this project aims to develop a viable and sustainable chemical recycling technology for CFRP waste that can efficiently and cost-effectively decompose the polymer matrix and manufacture new recyclable composites from both the recovered carbon fiber (rCF) and decomposed matrix polymer (DMP).

36 MATERIALS SCIENCE↗

DOE OSTI2025

Poisson Equation for a (General) Homogeneous d-Dimensional Ellipsoid with Applications to Beam Envelope Tracking

This note describes the solution of the free-space Poisson equation in the interior of a $d$-dimensional homogeneous ellipsoid, and the associated space charge fields. An explicit formula (\ref{Sformula}) is provided that relates the $d\times d$ matrix describing the space charge (quadratic) potential to the $d\times d$ covariance matrix of the ellipsoid. For the cases $d=2$ and $d=3$, this result is used to determine the linear map corresponding to a space charge kick, that may be used to push the beam $6\times 6$ covariance matrix during envelope tracking. The treatment of upright ellipsoids for $d=2$ and $d=3$ is well-represented in the literature. However, the approach taken here emphasizes a general ellipsoid with arbitrary correlations in any dimension. The Appendix provides a general solution of the free-space Poisson equation in dimension $d$ for a source distribution with ellipsoidal symmetry.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

In-Silico Design of Next Generation Cellulose-Derived Packaging Materials (CRADA Final Report)

Developing sustainable solutions for single-use packaging is an important objective to combat the environmental crisis of plastics pollution. Most embodiments of cellulose-based packaging materials, including CellophaneTM, are completely biodegradable in both terrestrial and marine environments. However, petroleum-derived alternatives offer some performance advantages for metrics such as moisture barriers and mechanical properties. This project leverages molecular dynamics simulation to investigate how molecular modifications to cellulose-based polymer assemblies impact their material properties. An important performance criterion for the modified materials was to retain biodegradability; thus, modifications by naturally occurring, biodegradable additives were the focus of this study. Specifically, we developed models with xylan and lignin of varying monomeric compositions into the cellulose matrix. The mechanical properties were investigated by performing stress-strain simulations, and the water barrier and hydrophobicity were investigated by simulating the water contact angle. Our findings indicate that the incorporation of xylan into the cellulose matrix tends to increase the mechanical properties with an optimal loading of ~27 wt%. We also predict that orienting the nanoscale directionality of the xylan chains such that they are perpendicular to the cellulose fibrils will dramatically increase mechanical strength. In contrast, the incorporation of lignin tends to weaken the composite at all loadings investigated. Simulations of water contact angle predicted that coating polymers on the surface of the cellulose assembly creates a more hydrophobic surface than incorporating them throughout the matrix. Of the coatings investigated, lignin resulted in the most hydrophobic surface, followed by pectin and keratin, which both imparted modest increases in hydrophobicity. Future experimental work done by Futamura will focus on designing material prototypes to capitalize on the predictions of performance enhancement obtained from molecular modeling. While substantial progress was made by the simulations performed in this project, there still exists a vast parameter space that we were unable to investigate, including branching, functional group decoration, and degree of polymerization of polymer additives. However, the methods developed in this initial investigation will facilitate more rapid evaluation of the impact of molecular characteristics on the performance of biopolymer composite materials and thereby accelerate future materials discovery efforts in this area.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Develop a new integrated macro→micro←nano (MMN) multiscale modeling framework to optimize high strength aluminum alloys and processes for vehicle light-weighting

Bending tests provide a means to study plane strain fracture performance of 6000 series aluminum alloys. Metrics from bending tests have been correlated with self-pierce riveting (SPR) performance of a high strength AA6111 automotive aluminum alloy in previous works. Using the ORNL HPC resources, this project developed an innovative macro→micro←nano (MMN) multiscale microstructure-based finite element (FE) code to further understanding of the relationship between microstructure and fracture properties of high-strength 6000-series alloys. This work started with microstructural characterization in both mesoscale and nanoscale and bending performance characterization of AA6111 HS2-T6 alloy at Ford, the MMN framework was applied to this alloy to simulate 3-point VDA bending. From the results of the macro-modeling of 3-point VDA bending, the critical region of fracture was identified, and the region geometry was used to construct the micro-model. The fracture criterion of micron-scale precipitates and aluminum matrix which contains submicron and nano particles (AL-SMP-NP) in the micro-model was calibrated and validated by comparing simulated and measured bending results. With the AL-SMP-NP fracture strain obtained, the fracture strain of Al-matrix containing nano particles (ALNP) will similarly be determined by a submicron scale-model using an edge-constrained FE modeling approach developed by Hu et al. With the ALNP fracture strain obtained, the fracture strain of Al-matrix containing no particles will similarly be determined by a nano scale-model using an edge-constrained FE modeling approach. After the MMN framework is built and fracture criterion calibrated, nano-model FE simulations with virtual microstructures was performed to obtain a reduced order model (ROM) of the fracture criterion of the ALNP as a function of volume fraction, size, and distribution of the nanoparticle. This nano→submicron→micro modeling part allows exploration of the influence of different material nanostructures from different process conditions on the bending properties within the multiscale bending simulation framework and the ROM of material bendability as a function of nanoparticle size and shape was established. This obtained reduced order model (ROM) could help guide the design and selection of materials to improve existing SPR process models that could replace trial-and-error rivet/die selection and help to design new rivet/die combinations capable of robustly joining new higher strength 6000 5 series alloys in automotive body structures. This would enable lightweighting of Ford vehicles leading to greater fuel efficiency and reduce manufacturing time and energy.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Pion and Kaon decays from lattice QCD and |Vud|, |Vus|

We review lattice calculations of the form factors for Kaon semileptonic decays and the leptonic decay constants of Kaons and Pions. We discuss the determinations of the CKM matrix elements $|V_{ud}|$ and $|V_{us}|$ using the form factors and experimental results. We use the CKM matrix elements to test first-row unitarity of the CKM matrix.

Simone, J. N. [Fermilab] (ORCID:0000000185153337)↗

DOE OSTI2025

FY24 Task 5: Leachate Disposition

Directly feeding sludge solids to the high-level waste (HLW) Waste Treatment Plant represents an alternative flowsheet seeking to initiate sludge processing as soon as possible. Key processing functions previously captured during baseline pretreatment operations include leaching and washing prior to solids concentration. These operations should be considered in the potential direct feed flowsheets to maximize waste feed loading, minimize HLW volume, and mitigate corrosion challenges associated with vitrification of high phosphate and fluoride concentrations. Additionally, single-shell tank (SST) retrievals and waste transfers to double-shell tanks (DSTs) in a direct feed flowsheet would likely also benefit from some level of leaching, washing, and solids concentration in order to reduce DST space and mission duration. These operations could occur in a new facility or potentially in available DSTs. If washing and leaching are utilized, an effective disposition pathway for the wash water and leachate solutions are needed. Three target species that benefit significantly from leaching and washing are phosphate, fluoride and aluminum. Phosphate (PO 4 3- ) and fluoride (F - ) can contribute substantially to the amount of carrier fluid needed for dissolution, and the resulting volume of liquid generated. Disposition of this retrieval solution should be evaluated in order to prevent crystallization of these anions throughout system processing. Since there is a high probability that any retrieval solutions will be at or near their PO 4 3- and F - solubility limits, evaporation or blending with a high Na supernate (>3.5 M) is not recommended for the wash water streams without a method to remove precipitants prior to solution disposal. Additionally, aluminum present in the southeast quadrant of the Hanford site represents roughly 60% of the waste solids in the initial processing tanks. These aluminum solids are in the form of gibbsite (Al(OH) 3 ) and can pose significant challenges for processing due to the fast-settling times and high solids loading associated with these materials. Easily remediated by caustic addition to the solids, these wash solutions could be processed through crystalline silicotitanate (CST) ion exchange columns to prepare the supernate solutions for disposition. The current target for feed conditions to the Low Activity Waste (LAW) melter are waste streams that contain nominally 5-6 M Na. Fractions within the tanks contain upwards of 0.2 M phosphate and fluoride in solution at 3.5 M Na. Concentrating these solutions above 5 M Na would result in an exceedance of the solubility limits, and potential for uncontrolled precipitation of the phosphate and fluoride crystal material. The resulting crystalline salt material is typically sodium fluoride phosphate, also referred to as natrophosphate (Na 7 FPO 4 ·19H 2 O). Salt phases are of importance in tank waste due to their chemical reactivity, which can result in precipitation, dissolution, or transformation, impacting any downstream processes (Bolling et al. 2020, Russell, Snow, and Peterson 2010). Salt generation and precipitation could pose challenges by causing system plugging and melter corrosion if left in the supernate stream, or limit sodium molarity of the supernate that would be accepted without incident in waste operations. To understand the impact of this salt generation, the crystallization of natrophosphate in multiple simulant feed matrices was studied to understand the implications of various tank waste supernate chemistries. Three matrices were examined: high PO 4 3- /low F - , low PO 4 3- /high F - , and an average matrix. Subsequent testing was performed with the average matrix with the inclusion of CsNO 3 , and a final run with the average matrix including CsNO 3 and a 137 Cs spike for tracer purposes.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

Printable Fiber Reinforced Cement Composites – Feasibility Study

Additive manufacturing is enabling the manufacturability of structures with previously unattainable complexity or functionality, and there is growing interest in additive manufacturing of “printed” concrete structures. The focus of this Phase 1 Technical Collaboration (TC) project was to evaluate feasibility of printing hybrid cement composite structures reinforced with textile carbon fibers (tCF). This project leverages other (non-IACMI) projects on cement formulation and printing process development, as well as on the production process for tCF. This project’s primary focus was to explore cement composite mix design with textile carbon fibers to be manufactured by MonteFibre (TC partner) and evaluate suitable fiber-matrix interface or sizing for cement composites working with Michelman (TC Partner). This project supports IACMI’s goal of reducing the cost and embodied energy of carbon fiber composites. Cost is one of the fundamental challenges to carbon fiber reinforced cement composites. Cement is an extremely inexpensive material (approximately $\$$0.05/lb). Adding 1 wt% of conventional carbon fiber to cement quadruples its cost. Therefore, the need to use low-cost carbon fiber and ensure that the additional cost of the carbon fiber has a greater cost benefit to the final product. This was the first preliminary evaluation to integrate tCF reinforcement in cement composites, and such potential tCF utilization should significantly reduce materials cost. Cement composite production is energy and emissions intensive, thus by strengthening it less material will be required. Hence, the embodied energy and production time of the resulting structures will be reduced. Additionally, integrating these new materials into additive processes can enable selective use of the material in high load or stress areas. It is noteworthy that past work in this field of fiber reinforced cement composites did not consider the optimization of fiber-matrix interface using suitable sizing. Carbon fiber reinforcement offers potential added benefits of thermal conductivity (which affects cure rate) and flow behavior that could provide opportunities for site specific utilization of carbon fiber on hybrid cement structures (e.g. use the fiber reinforcement on outer surfaces to enhance strength and modulus and then infiltrating the internal structures with conventional concrete). MonteFibre was the industry lead and planned on supplying the tCF for this project. However, during the short Phase-1 duration of this project, MonteFibre was unable to produce tCF for this project due to manufacturing plant being off-line throughout the course of the project. The project team decided to pursue an alternate option which involved demonstrating printable concrete with steel fibers by the ORNL lead, Dr. Brian Post. The University of Tennessee collaboration team focused on evaluating the suitable chemical sizing for carbon fibers working with Michelman and also developed methods for material characterization of cement-based composites to evaluate the material response for compression, shear, flexure, and tension. The two milestones for University of Tennessee, Knoxville were realized related to identification of one sizing suitable for carbon fiber reinforced cement composite and developing data associated with mechanical behavior of unreinforced (neat) and carbon fiber reinforced cement composites. ORNL could not complete the task of carbon fiber reinforced printed cement composites due to the reasons mentioned earlier, but was able to replace tCF with steel fibers to demonstrate the feasibility of printing with fiber reinforced cement composites. The Project team reviewed possible sizing chemistry available in collaboration with Michelman for use on carbon fiber reinforcement in cement composites and concrete applications. Our initial goal was to identify a sizing most promising for formulation with textile carbon fibers (tCF) to deliver excellent mechanical properties in composite material state. Since tCF was not available for this project as originally envisioned, the team continued this task to identify a suitable sizing for carbon fiber applications by applying such sizing to lower cost carbon fibers currently available commercially from Zoltek called Panex fibers. At a future time this can be optimized for textile carbon fibers from Montefibre. The bulk of previous work on carbon fiber reinforced cement has neglected the importance of fiber-matrix adhesion on mechanical properties of the cement composite and identifying this missing link was an important accomplishment for future research. Tensile behavior of fiber reinforced concrete is important to evaluate in order to realize the dream of concrete products that do not need reinforcing steel. Important sample preparation and testing procedures were addressed in this study and it was concluded that substantial improvements in tensile behavior, without compromising compressive strength, and improved ductility can result from the use of carbon fiber reinforcement.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Developments in Performance and Portability for MadGraph5_aMC@NLO

Event generators simulate particle interactions using Monte Carlo techniques, providing the primary connection between experiment and theory in experimental high energy physics. These software packages, which are the first step in the simulation worflow of collider experiments, represent approximately 5 to 20% of the annual WLCG usage for the ATLAS and CMS experiments. With computing architectures becoming more heterogeneous, it is important to ensure that these key software frameworks can be run on future systems, large and small. In this contribution, recent progress on porting and speeding up the Madgraph5_aMC@NLO event generator on hybrid architectures, i.e. CPU with GPU accelerators, is discussed. The main focus of this work has been in the calculation of scattering amplitudes and "matrix elements", which is the computational bottleneck of an event generation application. For physics processes limited to QCD leading order, the code generation toolkit has been expanded to produce matrix element calculations using C++ vector instructions on CPUs and using CUDA for NVidia GPUs, as well as using Alpaka, Kokkos and SYCL for multiple CPU and GPU architectures. Performance is reported in terms of matrix element calculations per time on NVidia, Intel, and AMD devices. The status and outlook for the integration of this work into a production release usable by the LHC experiments, with the same functionalities and very similar user interfaces as the current Fortran version, is also described.

Valassi, Andrea↗

DOE OSTI2023

Electroweak box diagrams on the lattice for pion and neutron decay

CKM matrix is unitary by construction in the standard model(SM). The recent analyses on the first row of CKM matrix show ≈3σ tension with unitarity. Nonperturbative calculations of the radiative corrections can reduce the theory uncertainty in CKM matrix elements. Here we compute the electroweak box contribution to the pion and kaon β decays using seven Nf = 2 + 1 + 1 HISQ-Clover lattice with various pion mass and lattice spacing. The continuum and chiral limit is taken using the leading dependence on M π and a, where M π extrapolation is taken to the physical pion mass and SU(3) symmetric mass for pion and kaon box contribution, respectively. Our results are $\square^{VA}_{γW}|_{π}$ = 2.820(28) × 10 –3 and $\square^{VA}_{γW}|_{K}$ = 2.384(17) × 10 –3 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Improved quantum algorithms for linear and nonlinear differential equations

We present substantially generalized and improved quantum algorithms over prior work for inhomogeneous linear and nonlinear ordinary differential equations (ODE). Specifically, we show how the norm of the matrix exponential characterizes the run time of quantum algorithms for linear ODEs opening the door to an application to a wider class of linear and nonlinear ODEs. In [1], a quantum algorithm for a certain class of linear ODEs is given, where the matrix involved needs to be diagonalizable. The quantum algorithm for linear ODEs presented here extends to many classes of non-diagonalizable matrices including singular matrices. The algorithm here is also exponentially faster than the bounds derived in [1] for certain classes of diagonalizable matrices. Our linear ODE algorithm is then applied to nonlinear differential equations using Carleman linearization (an approach taken recently by us in [2]). The improvement over that result is two-fold. First, we obtain an exponentially better dependence on error. This kind of logarithmic dependence on error has also been achieved by [3], but only for homogeneous nonlinear equations. Second, the present algorithm can handle any sparse matrix (that models dissipation) if it has a negative log-norm (including non-diagonalizable matrices), whereas [2] and [3] additionally require normality.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗