Search NASA⌕ Search

Engineering topics

Kumar, S.

Publications and source records attributed to Kumar, S..

At least 19 records

Integration of Online Cross-Section Generation Capability with Depletion and Transient Solvers in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE)-based reactor multiphysics analysis application jointly developed by Argonne and Idaho National Laboratories under the DOENE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In FY25, an online crosssection generation capability based on the Self-Shielding Application Programming Interface (SSAPI) was demonstrated for TRISO-fueled reactor problems under steady-state conditions. This fiscal year, that capability was extended to support depletion and transient multiphysics calculations, enabling high-fidelity analyses that generate self-shielded cross sections on the fly from the actual evolving composition and temperature states rather than from pre-tabulated libraries. For depletion, a two-way coupling was established in which SSAPI computes compact-averaged self-shielded cross sections that the depletion solver then uses to advance the Bateman equations, with the updated compositions returned to SSAPI at each step; the depletion module was refactored to support both library-based and SSAPI-based cross sections, and additional logic was added to track daughter isotopes and to exclude minor isotopes for efficiency. For transient analysis, the SSAPI multigroup library was extended with the kinetics data required for time-dependent calculations, the Improved Quasi-Static (IQS) scheme was coupled with SSAPI, and several supporting capabilities were implemented, including a self-shielding treatment that lets control rods and drums move within a self-shielded model, which had previously been impossible and had ruled out rod- and drum-movement transients with on-the-fly cross sections altogether, a new mixing scheme for delayed-neutron precursor decay constants, a checkpoint-based restart workflow, and performance improvements such as pointwise cross-section interpolation and the bypassing of unnecessary Dancoff factor calculations. The implemented capabilities were verified against Serpent Monte Carlo solutions. For depletion, a prismatic pin-cell problem based on a Next Generation Nuclear Plant (NGNP) Very High Temperature Reactor benchmark showed excellent agreement, with eigenvalue differences within 200 pcm over the entire burnup range (up to 140 MWD/kgU) and fission-product and actinide inventories agreeing to within 0.8% and 2.5%, respectively; a heat-pipe microreactor assembly problem with a much higher fuel loading confirmed the same behavior and quantified the bias introduced when the multigroup equivalence effect is neglected. For transient analysis, a pin-cell problem with a step reactivity insertion and temperature feedback reproduced the analytically expected asymptotic power and showed close agreement between the direct and IQS solutions, and a two-dimensional microreactor core problem with control-drum rotation exercised the new moving-drum self-shielding treatment and demonstrated successful coupling of the online crosssection generation with both the direct and IQS transient methods. The capability was further exercised on a full-core pebble-bed problem, in which Griffin was coupled with the System Analysis Module (SAM) to simulate load-following operation of the gPBR with the Doppler feedback resolved at the TRISO fuel kernel temperature. These developments in Griffin provide a convenient, high-fidelity approach to cross-section generation for advanced thermal reactors with geometrically complex and highly heterogeneous configurations, including TRISO-fueled prismatic and pebble-bed systems, and support steady-state, depletion, and transient multiphysics calculations. They also enable self-shielded cross sections to be evaluated directly at the actual coupled state of the system, thereby establishing a foundation for high-fidelity, fully coupled multiphysics analysis of advanced reactors

Park, H.↗

The Surface-Topography Challenge: A Multi-Laboratory Benchmark Study to Advance the Characterization of Topography

Surface performance is critically influenced by topography in virtually all real-world applications. The current standard practice is to describe topography using one of a few industry-standard parameters. The most commonly reported number is Ra, the average absolute deviation of the height from the mean line (at some, not necessarily known or specified, lateral length scale). However, other parameters, particularly those that are scale-dependent, influence surface and interfacial properties; for example the local surface slope is critical for visual appearance, friction, and wear. The present Surface-Topography Challenge was launched to raise awareness for the need of a multi-scale description, but also to assess the reliability of different metrology techniques. In the resulting international collaborative effort, 153 scientists and engineers from 64 research groups and companies across 20 countries characterized statistically equivalent samples from two different surfaces: a “rough” and a “smooth” surface. The results of the 2088 measurements constitute the most comprehensive surface description ever compiled. We find wide disagreement across measurements and techniques when the lateral scale of the measurement is ignored. Consensus is established through scale-dependent parameters while removing data that violates an established resolution criterion and deviates from the majority measurements at each length scale. Our findings suggest best practices for characterizing and specifying topography. The public release of the accumulated data and presented analyses enables global reuse for further scientific investigation and benchmarking.

42 ENGINEERING↗

The effect of radiation damage on the light yield and uniformity of candidate plastic scintillator tiles for the CMS hadron calorimeter upgrade

A study has been performed to understand the effects of radiation damage on various plastic scintillator tiles considered for a possible upgrade of the hadron calorimeter of the CMS detector. Measurements were made with unirradiated tiles and with tiles that had been irradiated in the CMS collision hall to a dose of 44 kGy. Results are presented for the tiles of different shapes in terms of the energy spectrum, efficiency as a function of the position at which each tile was hit, as well as light yield. All the tiles showed a light reduction of up to about 50%. The tiles with the shape currently used in the CMS detector did not see increased non-uniformity of light collection, while a significant disuniformity was observed for the tiles considered as alternatives.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The Compelling Need for a Mid-Scale Stellarator Facility

In the pursuit of the goal of commercial fusion as an abundant and safe source of energy, the stellarator is a leading concept with compelling attractiveness and demonstrated performance. In this white paper we, as a community of US researchers from Universities, National Laboratories, and Private Industry, involved in studying the stellarator concept, lay out the programmatic and technical motivation for a new mid-size stellarator research facility in the US. This contribution is complementary to several other whitepapers authored by members of our community which address different mid-scale stellarator aspects. A community based technical facility proposal has been prepared by F. Parra, et al: Flexible Stellarator Physics Facility. Two private stellarator companies have submitted proposals supporting the development of a mid-scale stellarator: Thea Energy (C.P.S. Swanson, et al.), Type One Energy (W. Guttenfelder, et al.).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Compelling Need for a Mid-Scale Stellarator Facility

In the pursuit of the goal of commercial fusion as an abundant and safe source of energy, the stellarator is a leading concept with compelling attractiveness and demonstrated performance. A new mid-size stellarator is needed to retire risks and innovate towards a high performance, economically attractive, stellarator Fusion Pilot Plant. In this presentation we, as a community of US researchers from Universities, National Laboratories, and Private Industry, involved in studying the stellarator concept, lay out the programmatic and technical motivation for a new and modern mid-size stellarator research facility. A new mid-scale stellarator is needed to realize the potential predicted by a solid body of theory and simulation along with advances in computational tools for optimization and non-linear turbulence modeling. Notably, it is possible to combine the advantages of the stellarator (steady state, no current drive, no disruptions) with the good confinement regularly achieved in tokamaks. The top priorities for experimental work, and the motivation for a mid-scale stellarator experiment are: turbulence control, non-resonant divertor, MHD stability at large beta, confinement of fast particles, and coil simplification. A new mid-size quasi-symmetric stellarator, built as a user facility, would complement existing research at Wendelstein 7-X and LHD and strongly augment private industry. It would provide a program of innovative research, concept validation, theoretical advancement, and workforce development. Growing support and interest for stellarators by the fusion community and private industry affirms this rationale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Structure of Low-spin States in \(^{45}\)Sc Studied via Coulomb Excitation

The electromagnetic structure of 45 Sc at low excitation energy was investigated via low-energy Coulomb excitation at the Heavy Ion Laboratory (HIL) of the University of Warsaw and at the Inter-University Accelerator Centre (IUAC) in New Delhi. A set of reduced E2, E3, and M1 matrix elements was extracted from the collected data using the GOSIA code. The reduced transition probability B(E2; 11/2 – → 7/2 – ) has been determined, allowing us to deduce the lifetime of the 11/2 – state at 1237 keV. In addition, the upper limit on the reduced transition probability B(E3; 7/2 – → 5/2 + ) has been determined for the first time. Finally, new large-scale shellmodel and beyond-mean-field calculations were performed to interpret the structure of this nucleus.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Magnetotail Variability During Magnetospheric Substorms

Abstract In this work, we present a statistical study of substorms covering a five‐year period 2016–2020. Substorm phases were identified from time series of the SuperMAG AL (SML) index using a list of 5,077 previously identified substorm onsets, the SML peak value marking transition from expansion to recovery phase, and the recovery identified as return to activity less than −100 nT in the SML index. Magnetic field observations from THEMIS, RBSP, and MMS missions were used to study the magnetotail characteristics during the substorm evolution. A superposed epoch analysis indicates that the substorm onset occurs almost simultaneously with a few minutes of uncertainty throughout the magnetotail, ranging from geostationary orbit to 20 RE. The onset in the transition region precedes the ground onset by a few minutes. The peak SML time coincides with the peak of the outer transition region Δ B Z , which suggests that the field‐aligned currents driving the SML activity arise from the outer transition region. Analysis of 2D maps of the tail magnetic field shows that the magnetotail current changes are limited to the center of the tail within | Y | < 10 R E . The substorm recovery is fastest in the inner transition region and lasts longer when moving further out. We did not find major asymmetries in the substorm signatures associated with IMF B Y or B Z .

Kumar, S.↗

A Multiwavelength Investigation of PSR J2229+6114 and its Pulsar Wind Nebula in the Radio, X-Ray, and Gamma-Ray Bands

Abstract G106.3+2.7, commonly considered to be a composite supernova remnant (SNR), is characterized by a boomerang-shaped pulsar wind nebula (PWN) and two distinct (“head” and “tail”) regions in the radio band. A discovery of very-high-energy gamma-ray emission (E γ > 100 GeV) followed by the recent detection of ultrahigh-energy gamma-ray emission (E γ > 100 TeV) from the tail region suggests that G106.3+2.7 is a PeVatron candidate. We present a comprehensive multiwavelength study of the Boomerang PWN (100″ around PSR J2229+6114) using archival radio and Chandra data obtained two decades ago, a new NuSTAR X-ray observation from 2020, and upper limits on gamma-ray fluxes obtained by Fermi-LAT and VERITAS observatories. The NuSTAR observation allowed us to detect a 51.67 ms spin period from the pulsar PSR J2229+6114 and the PWN emission characterized by a power-law model with Γ = 1.52 ± 0.06 up to 20 keV. Contrary to the previous radio study by Kothes et al., we prefer a much lower PWNB-field (B∼ 3μG) and larger distance (d∼ 8 kpc) based on (1) the nonvarying X-ray flux over the last two decades, (2) the energy-dependent X-ray size of the PWN resulting from synchrotron burn-off, and (3) the multiwavelength spectral energy distribution (SED) data. Our SED model suggests that the PWN is currently re-expanding after being compressed by the SNR reverse shock ∼1000 yr ago. In this case, the head region should be formed by GeV–TeV electrons injected earlier by the pulsar propagating into the low-density environment.

Astronomy & Astrophysics↗

Status of PyGriffin Development for Integration into the NEAMS Workbench

The integration of Griffin into Workbench, by way of the PyGriffin code package, was initiated by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY-2022. PyGriffin was developed as a Python wrapper to Griffin to streamline the complex workflow involving mesh, cross section, and Griffin input generation, code execution, and results postprocessing, by improving user experience with high fidelity neutronics analysis through the Workbench GUI interface. PyGriffin can be used as a standalone application or through the PyARC code suite. FY-2023 saw the continued development of PyGriffin to expand its capabilities including advanced post-processing of simulation results and initial development of the Monte Carlo (MC) cross section generation workflow with the Shift MC code. PyGriffin was also approved for open-source software status in FY-2023, initiating plans to move the code repository into an open location. In addition to PyGriffin development, there were many improvements and developments made in PyARC leading to several software releases (latest version is v2.3.0), but the focus of this report is PyGriffin.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Energetic Electron Precipitation Driven by Electromagnetic Ion Cyclotron Waves from ELFIN’s Low Altitude Perspective

We review comprehensive observations of electromagnetic ion cyclotron (EMIC) wave-driven energetic electron precipitation using data collected by the energetic electron detector on the Electron Losses and Fields InvestigatioN (ELFIN) mission, two polar-orbiting low-altitude spinning CubeSats, measuring 50-5000 keV electrons with good pitch-angle and energy resolution. EMIC wave-driven precipitation exhibits a distinct signature in energy-spectrograms of the precipitating-to-trapped flux ratio: peaks at >0.5 MeV which are abrupt (bursty) (lasting ~17 s, or ΔL ~ 0.56) with significant substructure (occasionally down to sub-second timescale). We attribute the bursty nature of the precipitation to the spatial extent and structuredness of the wave field at the equator. Multiple ELFIN passes over the same MLT sector allow us to study the spatial and temporal evolution of the EMIC wave - electron interaction region. Case studies employing conjugate ground-based or equatorial observations of the EMIC waves reveal that the energy of moderate and strong precipitation at ELFIN approximately agrees with theoretical expectations for cyclotron resonant interactions in a cold plasma. Using multiple years of ELFIN data uniformly distributed in local time, we assemble a statistical database of ~50 events of strong EMIC wave-driven precipitation. Most reside at L ~ 5 – 7 at dusk, while a smaller subset exists at L ~ 8 – 12 at post-midnight. The energies of the peak-precipitation ratio and of the half-peak precipitation ratio (our proxy for the minimum resonance energy) exhibit an L-shell dependence in good agreement with theoretical estimates based on prior statistical observations of EMIC wave power spectra. The precipitation ratio’s spectral shape for the most intense events has an exponential falloff away from the peak (i.e., on either side of ~ 1.45 MeV). It too agrees well with quasi-linear diffusion theory based on prior statistics of wave spectra. It should be noted though that this diffusive treatment likely includes effects from nonlinear resonant interactions (especially at high energies) and nonresonant effects from sharp wave packet edges (at low energies). Sub-MeV electron precipitation observed concurrently with strong EMIC wave-driven >1 MeV precipitation has a spectral shape that is consistent with efficient pitch-angle scattering down to ~ 200-300 keV by much less intense higher frequency EMIC waves at dusk (where such waves are most frequent). At ~100 keV, whistler-mode chorus may be implicated in concurrent precipitation. These results confirm the critical role of EMIC waves in driving relativistic electron losses. Nonlinear effects may abound and require further investigation.

79 ASTRONOMY AND ASTROPHYSICS↗

Status of the NEAMS and ARC fast reactor tools integration to the NEAMS Workbench

The Workbench initiative was launched in FY-2017 within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to facilitate the transition from conventional tools to high-fidelity tools. The NEAMS Workbench provides a common user interface for model creation, real-time validation, execution, output processing, and visualization for integrated codes. The integration of the Argonne Reactor Computation (ARC) suite of codes into the NEAMS Workbench through the PyARC module was initiated in FY-2017. The ARC codes, which focus on fast reactor multiphysics analysis, contain both legacy codes like DIF3D and REBUS-3 that were developed with over 30 years of experience, and newer NEAMS additions like MC 2 -3, PERSENT, and PROTEUS. Recent work expended this integration to other tools to support the U.S. fast reactor community such as DASSH for sub-channel thermal-hydraulics, Griffin for high-fidelity deterministic transport calculations, and OpenMC for Monte Carlo simulations (with Shift integration planned for FY-2023). The integration of the “extended ARC” suite of codes into the NEAMS Workbench interface relies on the PyARC and PyGriffin modules to handles the pre- and post-processing of these codes input, and the runtime environment. The PyARC module together with the NEAMS Workbench interface are both released under Open Source Software licenses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

MOOSE Reactor Module Meshing Enhancements to Support Reactor Physics Analysis

The U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) program develops an integrated suite of advanced reactor physics tools built upon the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Each code generally requires an input finite element mesh on which the physics solution is calculated, reported, and transferred to other physics codes. The meshing process is often burdensome for the complex geometries present in reactors due to lack of easy-to-use, open-source meshing tools. To address the bottleneck associated with meshing complex geometries found in nuclear reactors, NEAMS initiated the development of the MOOSE Reactor Module starting in FY21. The Reactor Module builds off the existing MOOSE Mesh System to include targeted meshing capabilities such as the ability to generate hexagonal pin cells, assemblies with ducts, rotating control drums, cores, peripheral zones around a core, as well as the automatic labeling (“reporting IDs”) of pin, assembly, and planar zones to simplify post-processing of results. As a Physics Module in MOOSE, the Reactor Module is open-source, available with any MOOSE installation, directly compatible with MOOSE-based tools, and can be invoked from MOOSEbased applications to generate meshes. Functionality from the Reactor Module has been applied to several advanced reactor concepts to demonstrate user workflow improvements and accuracy. The primary objective of the Reactor Module is to improve useability of MOOSEbased tools by streamlining mesh generation and output inspection processes. During FY22, the functionality of the Reactor Module (and accompanying Mesh System) has been expanded based on user needs. First, the Reactor Geometry Mesh Builder capability developed primarily in FY21 has been refactored and merged to the public MOOSE repository. This capability wraps underlying Reactor Module mesh generators into a “Pin – Assembly – Core” workflow appropriate for conventional Cartesian and hexagonal geometries, and notably assigns material IDs during mesh generation stage and generates only the minimal number of blocks needed in order to reduce computational burden. Biasing and boundary layer options have been added to the base mesh generators as required by thermal hydraulics solvers. The reporting ID functionality has been expanded to differentiate ring-wise and azimuthal sectors within a pin for use with depletion algorithms, and VectorPostProcessor and Reporter objects are now available to integrate solution variables across zones based on ID combinations. Functionality to trim hexagonal meshes along the center or periphery has been developed so users may leverage symmetry and reflective boundary conditions to reduce the mesh size. A flexible and powerful tool to fill the space between two sidesets has been introduced to the framework and can be used for transition layers such as stitching two assemblies together with different numbers of pins, or for complex geometries which do not follow conventional Cartesian/hexagonal patterns. Finally, additional verification problems were performed with NEAMS physics tools in complement with existing NEAMS work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Fast Reactor Physics Model Verification Studies using ARC and PyARC Workflows

PyARC was recently developed at Argonne National Laboratory to automate many of the tasks required in the ARC (Argonne Reactor Computation) fast reactor simulation workflow, from input file generation, code execution, data transfer between ARC codes, and output postprocessing. PyARC will likely be the path forward to train new users of the ARC codes with the goal of wide adoption by the national laboratories, academia, and industry. In particular, for the ANL-JAEA collaboration under the Civil Nuclear Working Group (CNWG) project agreement NE-01, PyARC will be used to model the Joyo and EBR-II reactors for comparisons with measured data and calculated results from JAEA (Task 3: Fast Reactor Fuel and Core). As an additional avenue for verification and validation, this report investigates the use of PyARC towards a variety of existing ARC-based reactor models, in order to understand its efficacy in replicating the behavior of base ARC codes and better understand any limitations within modeling realistic fast reactor problems. To this end, PyARC was used to model the Joyo MKI, RBEC Benchmark-M, PRISM Mod-B, and EBR-II Run 138B cores, and its results were compared to those from existing ARC-based models. It was found that for hexagonal-based geometries PyARC was able to replicate the behavior of ARC codes to within 10 pcm for small reactor cores, and ~150pcm difference in eigenvalue for larger cores. These discrepancies are attributed primarily to differences in local mesh refinement options between ARC and PyARC, which currently cannot be resolved with PyARC’s latest version (1.6.0). In some of these cases, PyARC was used to model steady-state problems with initial core compositions originating from a prior REBUS depletion calculation. While PyARC was not designed to support such steady-state calculations, workarounds were applied to replicate the behavior of ARC-based calculations as closely as possible. Thus, these results demonstrate the wide extent to which they can be applied to fast reactor problems while still providing immense benefit to the user in terms of automating and standardizing common routines within the fast reactor analysis workflow. This study concluded that PyARC will be suitable for modeling the steady-state conditions of the EBR-II and Joyo fast reactors as part of the CNWG project agreement.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Physics demonstration and verification of MOOSE framework reactor module meshing capabilities

The recently developed Reactor module in the open-source MOOSE framework includes finite element meshing capabilities for analysis of common reactor geometries. Capabilities in the Reactor module have been employed to generate meshes for physics applications including a sodium-cooled fast reactor core analysis using Griffin, a fast reactor assembly thermal deformation analysis using MOOSE Tensor Mechanics, and a heat-pipe cooled microreactor coupled analysis using Griffin, Bison, and Sockeye. The process to build these meshes using MOOSE's meshing capabilities is described. Physics simulation results using MOOSE-based meshes have been verified to match results which leverage external meshing software such as Cubit or Argonne's Mesh Tools. MOOSE's Reactor module provides significant advantages compared to the use of external meshing tools when analyzing Cartesian and hexagonal reactor lattices using MOOSE-based applications: accessibility to the end user, low barrier to entry for new users, speed of mesh generation, volume preservation of meshed fuel pins, and simplification of analysis workflow when using MOOSE applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Multiwavelength Observations of the Blazar VER J0521+211 during an Elevated TeV Gamma-Ray State

We report on a long-lasting, elevated gamma-ray flux state from VER J0521+211 observed by VERITAS, MAGIC, and Fermi-LAT in 2013 and 2014. The peak integral flux above 200 GeV measured with the nightly binned light curve is (8.8 ± 0.4) × 10 –7 photons m –2 s –1 , or ~37% of the Crab Nebula flux. Multiwavelength observations from X-ray, UV, and optical instruments are also presented. A moderate correlation between the X-ray and TeV gamma-ray fluxes was observed, and the X-ray spectrum appeared harder when the flux was higher. Using the gamma-ray spectrum and four models of the extragalactic background light (EBL), a conservative 95% confidence upper limit on the redshift of the source was found to be z ≤ 0.31. Unlike the gamma-ray and X-ray bands, the optical flux did not increase significantly during the studied period compared to the archival low-state flux. The spectral variability from optical to X-ray bands suggests that the synchrotron peak of the spectral energy distribution (SED) may become broader during flaring states, which can be adequately described with a one-zone synchrotron self-Compton model varying the high-energy end of the underlying particle spectrum. The synchrotron peak frequency of the SED and the radio morphology of the jet from the MOJAVE program are consistent with the source being an intermediate-frequency-peaked BL Lac object.

79 ASTRONOMY AND ASTROPHYSICS↗

Progress in ITER ECE diagnostic design and integration

The ITER electron cyclotron emission (ECE) diagnostic system has primary roles in providing measurements of the core electron temperature profile and the electron temperature fluctuation associated with the neoclassical tearing modes. The ITER ECE system includes a radial and oblique line-of-sight. Four 43-meter long low-loss transmission lines (TLs) are designed to transmit millimeter wave power in the frequency range of 70–1000 GHz in both X- and O-mode polarization from the port plug to the ECE instrumentation room in the diagnostic building. The measurement instrumentation includes two Fourier transform spectrometer (FTS) systems and two radiometer systems. The Indian Domestic Agency (IN-DA) and United States Domestic Agency share the responsibility. The IN-DA scope excluding instrumentation and control has passed its preliminary design review and is progressing towards the final design review (FDR). In parallel, the diagnostic integration in different areas is ongoing. Furthermore, several captive components for the TLs have passed FDR and will be manufactured for installation in the tokamak building soon. A peer review meeting has been held on the prototype hot calibration source, and its integration and new thermal analysis in the diagnostic shield module are continuing. A prototype TL is being tested. A prototype polarizing Martin-Puplett type FTS, operating in the frequency range 70–1000 GHz, features an in-vacuo fast scanning mechanism and a cryo-cooled dual-channel THz detector system. Its performance has been assessed in detail against ITER requirements.

47 OTHER INSTRUMENTATION↗