Search NASA⌕ Search

SEARCH · Search NASA

Results for “Scattering 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.

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

At least 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Impact of moderate pump–Stokes chirp on femtosecond coherent anti-Stokes Raman scattering spectra

The effects of inducing moderate chirp in the pump and Stokes beams in chirped-probe-pulse femtosecond coherent anti-Stokes Raman scattering (CARS) spectroscopy are discussed. This is an important issue for measurements in high-pressure systems where the CARS beams must be transmitted through glass of significant thickness to reach the probe volume of interest. The effects were investigated both experimentally, by inserting disks of SF-11 glass into the pump and Stokes beam paths, and theoretically, by incorporating pulse chirp into our time-dependent density matrix (TDDM) simulations. Experimentally, we found that inducing moderate pump and Stokes chirp produced significant spectral narrowing of the nonresonant four-wave mixing contribution to the CARS signal, and this allowed us to control the resonant excitation bandwidth to eliminate interferences from neighboring species. The temperature measurement accuracy and precision were essentially unchanged with respect to the unchirped pump and Stokes case. The effect of moderate pump and Stokes chirp on the narrowing of the Raman excitation efficiency envelope was investigated by solving the full set of TDDM through direct numerical integration. The numerical efficiency of the solutions was enhanced dramatically by developing a parallel version of the TDDM code. The calculated Raman excitation efficiency profile was incorporated in our phenomenological spectral fitting codes.

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

Symbol alphabets from plabic graphs II: rational letters

Symbol alphabets of n-particle amplitudes in N = 4 super-Yang-Mills theory are known to contain certain cluster variables of G(4, n) as well as certain algebraic functions of cluster variables. The first paper arXiv:2007.00646 in this series focused on n = 8 algebraic letters. In this paper we show that it is possible to obtain all rational symbol letters (in fact all cluster variables) by solving matrix equations of the form C Z = 0 if one allows C to be an arbitrary cluster parameterization of the top cell of G + (n-4, n).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Asymptotic behavior of meson transition form factors

One of the open issues in evaluations of the contribution from hadronic light-by-light scattering to the anomalous magnetic moment of the muon (g - 2)µ concerns the role of heavier scalar, axial-vector, and tensormeson intermediate states. The coupling of axial vectors to virtual photons is suppressed for small virtualities by the Landau–Yang theorem, but otherwise there are few rigorous constraints on the corresponding form factors. In this paper, we first derive the Lorentz decomposition of the two-photon matrix elements into scalar functions following the general recipe by Bardeen, Tung, and Tarrach. Based on this decomposition, we then calculate the asymptotic behavior of the meson transition form factors from a light-cone expansion in analogy to the asymptotic limits for the pseudoscalar transition form factor derived by Brodsky and Lepage. Finally, we compare our results to existing data as well as previous models employed in the literature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Single-ion magnetism in the extended solid-state: insights from X-ray absorption and emission spectroscopy

Large single-ion magnetic anisotropy is observed in lithium nitride doped with iron. The iron sites are two-coordinate, putting iron doped lithium nitride amongst a growing number of two coordinate transition metal single-ion magnets (SIMs). Uniquely, the relaxation times to magnetisation reversal are over two orders of magnitude longer in iron doped lithium nitride than other 3d-metal SIMs, and comparable with high-performance lanthanide-based SIMs. To understand the origin of these enhanced magnetic properties a detailed characterisation of electronic structure is presented. Access to dopant electronic structure calls for atomic specific techniques, hence a combination of detailed single-crystal X-ray absorption and emission spectroscopies are applied. Together K-edge, L 2,3 -edge and Kβ X-ray spectroscopies probe local geometry and electronic structure, identifying iron doped lithium nitride to be a prototype, solid-state SIM, clean of stoichiometric vacancies where Fe lattice sites are geometrically equivalent. Extended X-ray absorption fine structure and angular dependent single-crystal X-ray absorption near edge spectroscopy measurements determine Fe I dopant ions to be linearly coordinated, occupying a D 6h symmetry pocket. The dopant engages in strong 3dπ-bonding, resulting in an exceptionally short Fe–N bond length (1.873(7) Å) and rigorous linearity. It is proposed that this structure protects dopant sites from Renner–Teller vibronic coupling and pseudo Jahn–Teller distortions, enhancing magnetic properties with respect to molecular-based linear complexes. The Fe ligand field is quantified by L 2,3 -edge XAS from which the energy reduction of 3d z 2 due to strong 4s mixing is deduced. Quantification of magnetic anisotropy barriers in low concentration dopant sites is inhibited by many established methods, including far-infrared and neutron scattering. We deduce variable temperature L 3 -edge XAS can be applied to quantify the J = 7/2 magnetic anisotropy barrier, 34.80 meV (~280 cm –1 ), that corresponds with Orbach relaxation via the first excited, M J = ±5/2 doublet. The results demonstrate that dopant sites within solid-state host lattices could offer a viable alternative to rare-earth bulk magnets and high-performance SIMs, where the host matrix can be tailored to impose high symmetry and control lattice induced relaxation effects.

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

A proposed method for addressing large unphysical uncertainties in mubar

Mubar uncertainties in Section MF 34 MT 2 of ENDF/B-VIII.0 data are sometimes too large. Physical limitations of the bounds of mubar limit the maximum uncertainty to be < 1.0 and usually $\ll$ 1.0. Two ad hoc methods are proposed to artificially constrain the randomly sampled values of mubar to allowable values. The application of NJOY’s mubar covariance matrix to the “sandwich rule” is also discussed and a consistent way to apply sensitivity coefficients is described. Mubar uncertainties for the Jezebel k eff were found to be 156 pcm – a value comparable to other evaluated cross section uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Two-magnon bound states in the Kitaev model in a [111] field

It is now well established that the Kitaev honeycomb model in a magnetic field along the [111] direction harbors an intermediate gapless quantum spin liquid (QSL) phase sandwiched between a gapped non-Abelian QSL at low fields $H < H_{c1}$ and a partially polarized phase at high fields $H > H_{c2}$. Here, we analyze the low-field and high-field phases and phase transitions in terms of single- and two-magnon excitations using exact diagonalization and density matrix renormalization group methods. In this work, we find that the energy to create a bound state of two-magnon $Δ_p$ becomes lower than the energy to create a single spin flip $Δ_s$ near $H_{c2}$. In the entire Kitaev spin liquid, $Δ_p < Δ_s$ and both gaps vanish at $H_{c2}$. We make testable predictions for magnon pairing that could be observable in Raman scattering measurements on Kitaev QSL candidate materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Signatures of the attractive interaction in spin spectra of one-dimensional cuprate chains

Identifying the minimal model for cuprates is crucial for explaining the high- T c pairing mechanism. Recent photoemission experiments have suggested a significant near-neighbor attractive interaction V in cuprate chains, favoring pairing instability. To determine its strength, we systematically investigate the dynamical spin structure factors S ( q , ω ) using the density matrix renormalization group. Our analysis quantitatively reveals a notable softening in the two-spinon continuum, particularly evident in the intense spectrum at large momentum. This softening is primarily driven by the renormalization of the superexchange interaction, as determined by a comparison with the slave-boson theory. We also demonstrate the feasibility of detecting this spectral shift in thin-film samples using resonant inelastic x-ray scattering. Therefore, this provides a distinctive fingerprint for the attractive interaction, motivating future experiments to unveil essential ingredients in cuprates. Published by the American Physical Society 2024

Shen, Zecheng↗
DOE OSTI2025

Structure and Sulfur: Tuning the Viscoelastic and Surface Properties of Natural Keratin Fibers

Natural keratin fibers, such as wool, possess a complex hierarchical structure that governs their mechanical properties and surface energy. However, the extent to which these characteristics are influenced by combined contributions of structural variations (e.g., fiber diameter, intermediate filament (IF) packing) and chemical composition (e.g., disulfide bond density) remains poorly understood. In this study, we investigate wool fibers from five sheep breeds (Merino, Polwarth, Cheviot, Eider, and Devon) to elucidate how these factors influence viscoelasticity and surface interactions. Using a multimodal approach integrating interfacial and bulk characterization methods, including inverse gas chromatography (IGC), atomic force microscopy-infrared spectroscopy (AFM-IR), X-ray photoelectron spectroscopy (XPS), uniaxial tensile testing, and synchrotron small-angle X-ray scattering (SAXS), we show that the nanometer-thick 18-methyleicosanoic acid (18-MEA) layer is consistently present across all wool types and plays a key role in governing hydrophobicity and surface heterogeneity. A controlled isothermal treatment at 200 °C, designed to cleave disulfide bonds, results in a nearly 40% reduction in specific surface area across all fiber types, accompanied by a significant decrease in tensile strength and 80% reduction in elongation at break for Merino and Devon wool, but limited influence on the mechanical properties of Eider fibers. Furthermore, rate-dependent tensile testing within the elastic regime reveals distinct viscoelastic responses among the fiber types, suggesting that the sulfur-rich protein matrix surrounding IFs and its structure contribute actively to stress partitioning. Altogether, when combined with conclusions from SAXS measurements of IF spacing, our work offers compelling insights into the role of the keratin-associated protein (KAP) matrix in shaping wool fiber mechanics. Differences in mechanical behavior among wool types, despite similar IF spacing or sulfur content, highlight the importance of matrix composition and cross-linking density, suggesting that the molecular architecture of the KAP network may be a dominant factor in determining fiber performance.

X-ray scattering↗
DOE OSTI2020

Characterization of the Sulfur Saturated Versions of the High-Chromium Matrix Glasses

In this report, the Savannah River National Laboratory provides chemical analysis of a series of sulfur saturated melt versions of simulated nuclear waste glasses, and chemical analysis of the wash solutions resulting from the preparation of these glasses. The glasses were selected and fabricated by the Pacific Northwest National Laboratory as part of a broader study of the influence of glass composition on chemical durability, sulfur retention, and other properties. The resulting data will be used in the development of improved property/composition models to support operation of the Hanford Tank Waste Treatment and Immobilization Plant. Some degree of scatter among the Al₂O₃, B₂O₃, CaO, K₂O, Na₂O, and SiO₂ measurements was noted for the study glasses. The measured concentrations of B₂O₃, chlorine, and fluorine were below the targeted values for most of the study glasses, likely because of volatility during the multiple melting steps. The measured concentrations of Cr₂O₃ and K₂O were generally low relative to the targeted values. These components were identified in the wash solutions, indicating partitioning to the excess sulfate phase that was washed from the glasses prior to analyses. As expected, the measured concentrations of SO₃ in the glasses were higher than targeted due to the use of the sulfur saturation method in fabricating these glasses. Minor scatter among the triplicate measurements of some of the analytes of the wash solutions was noted. The wash solutions contained mainly sulfur and sodium, with moderate concentrations of chromium and potassium. Revision 1 of this report corrects the omission of Al₂O₃ values from Table A-4 in Appendix A.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Doping dependence of the electron-phonon coupling in two families of bilayer superconducting cuprates

While electron-phonon coupling (EPC) is crucial for Cooper pairing in conventional superconductors, its role in high-Tc superconducting cuprates is debated. Here, using resonant inelastic x-ray scattering at the oxygen K edge, we study the EPC in Bi 2 Sr 2 CaCu 2 O 8+δ (Bi2212) and Nd 1+x Ba 2-x Cu 3 O 7-δ (NBCO) at different doping levels ranging from heavily underdoped (p=0.07) to overdoped (p=0.21). We analyze the data with a localized Lang-Firsov model that allows for the coherent excitations of two phonon modes. While electronic band dispersion effects are non-negligible, we are able to perform a study of the relative values of EPC matrix elements in these cuprate families. In the case of NBCO, the choice of the excitation energy allows us to disentangle modes related to the CuO chains and the CuO 2 planes. Combining the results from the two families, we find the EPC strength decreases with doping at q∥=(-0.25, 0) r.l.u., but has a nonmonotonic trend as a function of doping at smaller momenta. This behavior is attributed to the screening effect of charge carriers. We also find that the phonon intensity is enhanced in the vicinity of the charge-density-wave excitations while the extracted EPC strength appears to be less sensitive to their proximity. By performing a comparative study of two cuprate families, we are able to identify general trends in the EPC for the cuprates and provide experimental input to theories invoking a synergistic role for this interaction in d -wave pairing.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Probing Carbon Mineralization Mechanisms in Pore and Bulk Fluids by Harnessing Architected Calcium Silicates

The ability to synthesize materials with well-controlled pore structures gives us unprecedented control over probing fluid interactions with reactive interfaces and advancing calibrated insights into coupled chemo-morphological interactions. One of the primary challenges in developing crystalline silicate materials lies in achieving ordered pore structures. Existing approaches of producing amorphous mesoporous metal silicates via sol–gel methods and heat-treatment of these materials to produce crystalline phases cause the pore structures in the amorphous phases to collapse. To overcome this challenge, carbon coating of amorphous mesoporous calcium silicate particles is carried out to retain the pore structure, while the material is heated to produce crystalline calcium silicate with calcium sulfate inclusions. The pore diameter in these materials is about 3.9 nm, with a surface area and a pore volume of 28.75 m 2 /g and 0.092 cm 2 /g, respectively. The mechanisms of carbon mineralization are investigated by reacting architected calcium silicate with 1 M Na 2 CO 3 and monitoring the evolution in the structural phases using operando wide-angle X-ray scattering (WAXS) measurements. Formation of stable calcium carbonate polymorph or calcite and metastable calcium carbonate polymorph or vaterite in pore and bulk fluids, respectively, resulting from the reaction between Na 2 CO 3 and CaSiO 3 , are noted. The mechanisms associated with carbon mineralization are delineated using ReaxFF molecular dynamics (MD) simulations. The surface dissolution reaction is initiated by 2H + ions that replace a Ca 2+ ion in the Ca–silicate matrix. Ca 2+ ions in the solution initially react with water to form calcium hydroxide and eventually form calcium (bi)carbonate. A slow and gradual increase in the formation of sodium silicate in the solution resulting from the reactions of silicic acid or the silicon dioxide reaction with sodium hydroxide is noted. When carbon mineralization occurs in environments bearing interfacial fluids, calcite is the dominant calcium carbonate polymorph, as determined using experiments with pore fluids and molecular-scale simulations. In conclusion, these studies provide fundamental insights into the mechanisms underlying the carbon mineralization of calcium silicate informed by experiments and molecular-scale simulations.

Calcium↗
DOE OSTI2022

Local conformations and heterogeneities in structures and dynamics of isotactic polypropylene adsorbed onto carbon fiber

Carbon fiber (CF) reinforced polymers (CFRPs) have experienced widespread use in various industries. One of the most important parameters that controls the macroscopic property of CFRPs is the interface between a polymer matrix and CF. There is growing evidence to suggest the formation of a bound polymer layer (BPL), i.e., polymer chains that physically adsorb on a filler surface. However, this interface is always in contact with the thicker part of a polymer matrix, rendering its understanding a difficult task. Therein, we use CF-reinforced isotactic polypropylene (iPP) as a rational CFRP. To characterize the BPL on the CF surface, we extracted it from the CFRP using solvent-rinsing with p-xylene. The physical and thermal properties of the BPL were characterized by differential scanning calorimetry and thermogravimetric analysis, while its microscopic structures and dynamics were probed by small-angle neutron scattering and quasi-elastic neutron scattering (QENS) techniques. Subsequently, we employed atomistic molecular dynamics (MD) simulations to complement the QENS results above the bulk melting temperature and reveal details that were experimentally inaccessible. In this work, we observed that the degree of crystallinity of the BPL was quite lower than the bulk, while the melting temperature of the BPL remained the same as the bulks. Within the given length and time scales probed by QENS and MD, we also observed that most of the bound chains were mobile, with the formation of a high-density region (less than 1 nm in thickness) near the CF surface. The segmental dynamics of the bound chains probed by both QENS and MD were also much faster than those of the free chains, possibly due to the presence of a free surface region at the topmost surface of the BPL. Furthermore, the MD results demonstrated that the backbone chains and side groups lie nearly flat on the CF surface, which is the driving force for the flattening process of the iPP BPL to overcome the conformational entropy loss in the total free energy.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Electronic Structure, Spectroscopy and Correlation Effects in Novel Materials

This project was targeted on electronic structure, spectroscopic studies, and correlation effects in a variety of novel materials of great current interest. Spectroscopies resolved highly in momentum, energy or spatial dimensions are playing a key role in unraveling the nature of the ground state and excitation properties in wide classes of novel materials. The seminal insights thus obtained are of critical importance not only for answering some of the fundamental questions facing condensed matter physics and materials science today, but also for understanding and thus helping to design and develop new materials with desirable properties, which will continue to be the key to the survival of mankind and its energy needs as a technological society long into the future. However, spectroscopies do not provide a direct map of electronic states, but act as a very complex ‘filter’ or ‘mapping’ of the underlying spectrum. This connection between the electronic states and measured spectra—called the ‘matrix element effect’—is in general an extremely complex function of the phase space of the experiment (e.g. energy/polarization of photons in photoemission), presenting both a challenge and an opportunity. So motivated, this project pursued techniques for realistic treatment of electronic spectra of a wide variety of materials, which served as a prelude to formulating and implementing methodologies for making direct connection with various spectroscopies such as ARPES, STS/STM, and inelastic light scattering. Specific systems are topological materials, including 2D ultrathin films beyond graphene, novel superconductors, and nanocrystals, among others. Our goal was to exploit the strengths of various spectroscopies to piece together the most complete picture of electronic states in systems of current interest, enabling direct and sharpened confrontations with theoretical models, and also to help advance the reach of various spectroscopies.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A Measurement of the Pion-Energy Dependence of Muon Neutrino Charged-Current Scattering to Final States With One Charged Pion in NOvA

The study of neutrino oscillations is a main priority for particle physics as the most immediately tractable lever on physics beyond the Standard Model. In particular, more insight into violation of the combined symmetry of charge-conjugation plus parity could yield clues to the origin of matter-antimatter asymmetry, and in some theoretical frameworks the neutrino mass could give insights into dark matter. NOvA is a long-baseline accelerator neutrino experiment with both a near and far detector that seeks to measure several of the parameters of the neutrino mixing matrix, as well as carry out a broad program of additional physics. NOvA has innovated and developed a variety of techniques in the space of neutrino physics, including expanding the use of machine learning techniques in reconstruction. NOvA has also set the stage for the next generation US-based long baseline experiment, DUNE. This dissertation details the creation of a new neutrino interaction vertex reconstruction package for NOvA, which offers enormous improvements in accuracy above the previous vertexer it replaces. This is accomplished using a Convolutional Visual Network trained on large datasets of simulated events. The vertexer is validated thoroughly against data. This dissertation also presents a cross-section measurement for $\nu_\mu + N \rightarrow \mu^- + 1\pi^\pm + X$ (where X does not include additional charged pions) binned in pion kinetic energy, a challenging measurement in general and for NOvA in particular. Measurements like this provide crucial inputs for neutrino interaction models utilized by neutrino oscillation experiments. This measurement in particular is important to the energy ranges relevant for DUNE in a band not well-covered by any other experiment.

Ewart, Erin [Indiana U., Bloomington (main)]↗
DOE OSTI2023

Compositional Modifications to Alter and Suppress Laves Phases in Al x CrMoTa y Ti Alloys

The development of refractory complex concentrated alloys in the Al–Cr–Mo–Ta–Ti alloy system is reported. Alloys with modified Al and Ta concentrations are designed using CALPHAD tools and produced via arc melting and characterized in both as-cast and annealed forms. Properties of the alloys, nature of the microstructures, and phase transformation behavior are described via X-ray diffraction, microstructural characterization, microhardness, and differential scanning calorimetry. Two alloys, namely, Al 0.25 CrMoTa 0.8 Ti and Al 0.75 CrMoTa 0.8 Ti, are represented by a body-centered-cubic matrix phase after annealing, along with a secondary Cr–Ta Laves phase of the C15 and C14 polytypes, respectively. In as-cast and annealed forms, the Al 0.75 CrMoTa 0.45 Ti alloy comprises a single-bcc phase. Microhardness of the Laves phase containing alloys demonstrates susceptibility to cracking, whereas the Al 0.75 CrMoTa 0.45 Ti alloy displays high specific hardness, signs of ductility as evidenced by slip traces near indentations, and minimal scatter of hardness values.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Nanoscale Ion Transport Enhances Conductivity in Solid Polymer-Ceramic Lithium Electrolytes

The predictive design of flexible and solvent-free polymer electrolytes for solid-state batteries requires an understanding of the fundamental principles governing the ion transport. In this work, we establish a correlation among the composite structures, polymer segmental dynamics, and lithium ion (Li + ) transport in a ceramic-polymer composite. Elucidating this structure–property relationship will allow tailoring of the Li + conductivity by optimizing the macroscopic electrochemical stability of the electrolyte. The ion dissociation from the slow polymer segmental dynamics was found to be enhanced by controlling the morphology and functionality of the polymer/ceramic interface. The chemical structure of the Li + salt in the composite electrolyte was correlated with the size of the ionic cluster domains, the conductivity mechanism, and the electrochemical stability of the electrolyte. Polyethylene oxide (PEO) filled with lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) or lithium bis(fluorosulfonyl) imide (LiFSI) salts was used as a matrix. A garnet electrolyte, aluminum substituted lithium lanthanum zirconium oxide (Al-LLZO) with a planar geometry, was used for the ceramic nanoparticle moieties. Further, the dynamics of the strongly bound and highly mobile Li + were investigated using dielectric relaxation spectroscopy. The incorporation of the Al-LLZO platelets increased the number density of more mobile Li + . The structure of the nanoscale ion-agglomeration was investigated by small-angle X-ray scattering, while molecular dynamics (MD) simulation studies were conducted to obtain the fundamental mechanism of the decorrelation of the Li + in the LiTFSI and LiFSI salts from the long PEO chain.

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

Enabling attractive-repulsive potentials in binary-collision-approximation monte-carlo codes for ion-surface interactions

Abstract Binary Collision Approximation (BCA) codes for ion-material interactions, such as SRIM, Tridyn, F-TRIDYN, and SDtrimSP, have historically been limited to screened Coulomb potentials even at low energies due to the difficulty in numerically solving the Distance of Closest Approach (DOCA) problem for attractive-repulsive potentials. Techniques such as direct n-body simulation or modifications to Newton’s method are either prohibitively costly or not guaranteed to work for all potentials. Advanced rootfinding techniques, such as companion matrix solvers, offer a solution. For many attractive-repulsive potentials, however, a companion matrix cannot be used directly, because there is no way to put the associated functions into a monomial basis form. A complementary technique is proxy rootfinding—by finding the best-fit polynomial approximant of a function, the zeros of the approximant can be guaranteed to be close to the zeros of the function. Using the Chebyshev basis and grid offers additional guarantees with regards to the quality of the approximation, the speed of convergence, and the avoidance of Runge’s phenomenon. By finding Chebyshev interpolants and using the Chebyshev-Frobenius companion matrix, the zeros of any real function on a bounded domain can be found. Here we show that using an Adaptive Chebyshev Proxy Rootfinder with Automatic Subdivision (ACPRAS) with appropriate scaling functions, numerical issues presented by attractive-repulsive potentials, including those of scale, can be handled. Using these techniques, we show that it is possible to include any physically reasonable interatomic potential in a BCA code, and to guarantee correctness of the resulting scattering angle calculations.

Materials Science↗
DOE OSTI2025

Code Coverage Status of the ARC Code PERSENT

The Argonne Reactor Code (ARC) software system supports users in their fast reactor design goals by providing neutronic, thermal-hydraulic, and structural analysis capabilities. PERSENT fulfills the role of generating reactivity coefficients for a given time point of a REBUS calculation usable in a point kinetics based safety analysis capability. PERSENT also provides a sensitivity coefficient capability on eigenvalue, reactivity worth, and several other key coefficients that are used in the follow-on safety analysis. Given a co-variance matrix, PERSENT can carry out the uncertainty quantification to indicate the amount of error in the reactivity coefficients derived from the errors in the cross section measurements. With continued improvement of computational resources, many of the geometry modeling capabilities in DIF3D that were primarily used in low order schemes are not really needed anymore. Today, the diffusion and transport capabilities of DIF3D-VARIANT are primarily used in the reactor design process with some scattered usage of DIF3D-FD and DIF3D-Nodal. PERSENT is part of the ARC code system and is built around DIF3D-VARIANT and the flux solution it provides. The purpose of the present work is to identify a set of test problems for PERSENT and assess the code coverage of PERSENT for those test problems. PERSENT treats the DIF3D executable as an external executable and thus the code coverage considerations only need to focus on the PERSENT source code and only a fraction of the connected modules in the existing ARC software library. The goal is to document what parts of the existing PERSENT code are touched by the set of test problems and which are not. Because the verification work done on PERSENT was focused on the most common uses of PERSENT for fast reactor analysis, the code coverage assessment of those capabilities is the highest priority. This will ensure that nothing is being missed by the existing verification test problems that users of PERSENT rely upon. The code coverage analysis of PERSENT was performed with the Code Coverage Tool of the Intel Fortran compiler which requires modifications to the compilation of PERSENT. The detailed coverage tables are given for each submodule of PERSENT. Most of the uncovered parts/files could be easily ignored because they are either for error message and debugging output or not needed by PERSENT today. Only a few uncovered parts of PERSENT deserve extending the verification test suite.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗