Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Serpent”

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 37 records · Page 2

Modeling Enhancements, Cross-Section Generation Updates, and Benchmarking with Shift

This technical report documents the modeling enhancements, cross-section generation updates, and bench marking with the Shift Monte Carlo code performed under the US Department of Energy Nuclear Energy Advanced Modeling and Simulation Program in FY 2024. The work performed included several modeling enhancements, such as integration of cross-section generation in Titan and the ability to produce microscopic multigroup cross sections with Shift. Benchmarking of the cross sections produced by Shift and the two-step workflow with Griffin was performed for three problems: the Advanced Breeder Test Reactor, a generic pebble bed reactor, and a TRISO heat pipe microreactor. Comparisons of results from these benchmark problems were done with Serpent, OpenMC, and Griffin. These enhancements provide a robust foundation for applying Shift for both reference and two-step neutronics analysis for advanced reactor simulation.

97 MATHEMATICS AND COMPUTING

Validation of Numerical Tools for Calculating Reactivity Feedback in Sodium Fast Reactors Using SEFOR Experimental Data

The Southwest Experimental Fast Oxide Reactor (SEFOR) was an experimental sodium-cooled fast breeder reactor operated from 1969 to 1972 with experiments designed to measure Doppler reactivity feedback in a wide temperature range from around 350 °F to temperatures approaching the melting point of mixed oxide fuel of around 5000 °F, providing valuable data for code validations. Co-supported by the Department of Energy (DOE) Fast Reactor Program (FRP) and the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, the SEFOR benchmark project focused on using the experimental data to validate numerical tools that are used in industry and academia to design and license sodium-cooled fast reactors (SFRs). By the end of FY-25, substantial progress was achieved in the SEFOR benchmark study. A variety of numerical tools commonly used for modeling SFRs were applied to develop models for SEFOR core configurations I-D, I-E, I-I, and I-J. These included Monte Carlo codes such as MCNP, Serpent, and Shift; deterministic codes such as the legacy Argonne Reactor Computation (ARC) suite and the high-fidelity NEAMS code Griffin; and the system analysis code SAS4A/SASSYS-1 (SAS). Using these models, both SEFOR zero-power experiments and power-ascending tests were successfully simulated. Comparisons were performed against experimental measurements of core criticalities, reflector worth, kinetics parameters (Λ/βeff), isothermal reactivity feedback (from 350 °F to 760 °F at zero power), and power-ascending reactivity feedback (as power increased from 0.4 MW to 17 MW). In general, these comparisons demonstrated very good agreement between numerical results and experimental data. In Fiscal Year 26 (FY-26), the SEFOR benchmark project will continue to address the modeling issues identified in FY-25. Effort will focus on the simulation of reactivity insertion transients in SEFOR core II using the ARC/SAS model. Future work will also focus on incorporating BISON into the SEFOR core modeling process to enable the first Multiphysics simulations of the isothermal tests based on the MOOSE framework.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS

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.

HTGR Multiphysics Application Drivers FY26 Updates

This report summarizes FY26 progress under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program's high-temperature gas-cooled reactor (HTGR) application driver work, covering a wide range of activities such as code validation and multi-physics code assessment. 1) A detailed SAM model of the High-Temperature Engineering Test Reactor (HTTR) was developed using a unique-block grouping approach, with an extended parallel thermal network method to capture block-to-block conduction and radiation heat transfer, and applied to steady-state simulations of the HTTR 30~MW and 9~MW cases. 2) In another activity, SAM's newly implemented multi-component gas flow model was validated against the Natural convection Shutdown heat removal Test Facility (NSTF) argon ingress experiment, correctly capturing the density-driven suppression and thermal recovery of natural circulation observed when argon is introduced into the air-cooled Reactor Cavity Cooling System (RCCS) loop. 3) For the OECD/NEA High Temperature Test Facility (HTTF) benchmark, we co-led the international benchmark activities as well as the OECD/NEA final benchmark report to be released at the end of this year. 4) Finally, the coupled Griffin-SAM modeling capability for pebble-bed HTGRs was advanced by verifying the Griffin neutronics solution against Serpent Monte Carlo for a realistic non-uniform temperature distribution, resolving several deficiencies in the SAM-to-Griffin temperature transfer scheme, and enabling distinct fuel kernel, moderator, and coolant temperatures for cross section feedback. These new features were demonstrated in a PBR load-following transient.

Lee, Alvin

On the Statistical Uncertainty of Monte Carlo-Calculated Scattering Sensitivities

Sensitivity coefficients calculated with Monte Carlo codes are widely used for nuclear data uncertainty quantification in the modeling and simulation of complex 3D reactor systems. This study systematically compares sensitivity coefficients and associated statistical uncertainties for the multiplication factor and fuel temperature reactivity across multiple Monte Carlo codes (SCALE/KENO, SCALE/Shift, MCNP, and Serpent) using simple models representing light-water reactors and advanced reactor concepts. For multiplication factor sensitivities, statistical uncertainties are generally acceptable, although scattering sensitivities show significantly larger statistical uncertainties than, for example, fission and capture reactions. Fuel temperature reactivity sensitivities show significantly larger statistical uncertainties across all reactions. Elastic scattering sensitivities are the most problematic: all Monte Carlo codes fail to resolve energy-dependent coefficients, and they produce dramatically different energy-collapsed values. Critically, the use of these sensitivity coefficients in nuclear data uncertainty propagation leads to reduced statistical uncertainties in individual uncertainty contributions. This can lead to the masking of unusable sensitivity coefficients and producing misleading uncertainty results. The findings of this study show that new or enhanced methods are needed to improve Monte Carlo elastic scattering sensitivity calculations. Additionally, this study shows the relevance of verifying sensitivity coefficients through direct perturbation calculations for individual nuclide reactions, instead of only for total cross sections as commonly done.

Bostelmann, Rike [ORNL] (ORCID:0000000165968088)

Geobotanical techniques for discriminating serpentine rock types in Western United States

The effectiveness of airborne multispectral imagery for identifying vegetation indicative of serpentine soils, which usually contain deposits of Ni, Cr, and Co, was examined with ground truth data. Analyses were carried out in SW Oregon, the Jasper Ridge of the Santa Cruz Mountains, and the Sierra Nevada foothills. Airborne studies were performed with 10-channel scanner, a IR scanner, and with false color composites, while ground studies classified the vegetation and mineral contents of the soil. The airborne studies were successfully used to discern serpentine soils from nonultramafic rock types, although difficulties were encountered in separating serpentine from nonserpentinized peridotite. A lack of deciduous trees on serpentized land was noted.

Mouat, D.

Spectroscopy of Moses Rock Kimberlite Diatreme

Three types of remote sensing data (Airborne Imaging Spectroscopy (AIS), NS001, Zeiss IR-photographs) were obtained for the Moses Rock kimberlite dike in southern Utah. The goal is to identify and characterize the mantle derived mafic component in such volcanic features. The Zeiss and NS001 images provide information on the regional setting and allow units of the dike to be distinguished from surrounding material. A potential unmapped satellite dike was identified. The AIS data provide characterizing information of the surface composition of the dike. Serpentized olivine-bearing soils are (tentatively) identified from the AIS spectra for a few areas within the dike.

Pieters, C. M.

Spectroscopy of Moses Rock dike using remote sensing

Zeiss IR-photographs, NS0001 (TM simulator) and airborne imaging spectrometer (AIS) data were obtained for the Moses Rock kimberlite dike in southern Utah to identify and characterize the distinctive mafic mineralogy of the dike as well as the surrounding sedimentary rocks. The Zeiss and NS001 images provide information on the regional setting and allow units of the dike to be distinguished from the sediments. The AIS data are narrow images obtained in 128 near-infrared channels and provide characterizing information on the surface composition through. Three distinct spectroscopic units were found which have been tentatively identified as serpentized olivine-bearing soils found in the dike and two types of gypsum bearing soils found in the surrounding sedimentary soils.

Mustard, J. F.

Alteration of Al-rich inclusions inside amoeboid olivine aggregates in the Allende meteorite

The primary phases of Al-rich inclusions in amoeboid olivine aggregates have undergone alteration reactions with the solar nebular gas. The simplest interpretation of the present observations is that melilite was the first primary phase to disappear with falling temperature, and was replaced by grossular + anorthite + feldspathoids, followed by fassaite; spinel was the last phase to be altered. Thermodynamic calculations suggest that Na-rich phlogopite could have formed at about 470 K and chlorite at about 328 K at a water fugacity of 0.000001, which is that of a gas of solar composition in this temperature range. The olivine around Al-rich inclusions is not serpentized, indicating the cessation of gas-solid equilibrium above 274 K.

Hashimoto, Akihiko

Abundance and distribution of ultramafic microbreccia in Moses Rock dike - Quantitative application of mapping spectroscopy

Data from the Airborne Imaging Spectrometer were used to map the distribution and abundance of the serpentized ultramafic microbreccia (SUM) component in the Moses Rock dike, which is a Tertiary diatreme located on the Colorado Plateau in Utah. The geologic setting and composition of Moses Rock dike are discussed together with its texture and the relationship to the bedrock of surface materials. These observations along with laboratory spectroscopic data are used to interpret surface mineralogy of the dike and the surrounding regions from the imaging spectometer data. The spatial distribution and the abundance of the primary surface components were calculated using a nonlinear model for the mixing of spectra from multicomponent surfaces. The derived SUM distribution and abundance data support McGetchin's (1968) model for the emplacement of Moses Rock dike as a fluidized system.

Mustard, John F.

Software engineering activities at SEI (Software Engineering Institute)

Prototyping was shown to ease system specification and implementation, especially in the area of user interfaces. Other prototyping approaches do not allow for the evolution of the prototype into a production system or support maintenance after the system is fielded. A set of goals is presented for a modern user interface environment and Serpent, a prototype implementation that achieves these goals, is described.

Chittister, Clyde

Carbon Isotope Systematics in Mineral-Catalyzed Hydrothermal Organic Synthesis Processes at High Temperature and Pressures

Observation of methane in the Martian atmosphere has been reported by different detection techniques. Reduction of CO2 and/or CO during serpentization by mineral surface catalyzed Fischer-Tropsch Type (FTT) synthesis may be one possible process responsible for methane generation on Mars. With the evidence a recent study has discovered for serpentinization in deeply buried carbon rich sediments, and more showing extensive water-rock interaction in Martian history, it seems likely that abiotic methane generation via serpentinization reactions may have been common on Mars. Experiments involving mineral-catalyzed hydrothermal organic synthesis processes were conducted at 750 C and 5.5 Kbars. Alkanes, alcohols and carboxylic acids were identified as organic compounds. No "isotopic reversal" of delta C-13 values was observed for alkanes or carboxylic acids, suggesting a different reaction pathway than polymerization. Alcohols were proposed as intermediaries formed on mineral surfaces at experimental conditions. Carbon isotope data were used in this study to unravel the reaction pathways of abiotic formation of organic compounds in hydrothermal systems at high temperatures and pressures. They are instrumental in constraining the origin and evolution history of organic compounds on Mars and other planets.

Fu, Qi

Carbon Isotope Measurements of Experimentally-Derived Hydrothermal Mineral-Catalyzed Organic Products by Pyrolysis-Isotope Ratio Mass Spectrometry

We report results of experiments to measure the C isotope composition of mineral catalyzed organic compounds derived from high temperature and high pressure synthesis. These experiments make use of an innovative pyrolysis technique designed to extract and measure C isotopes. To date, our experiments have focused on the pyrolysis and C isotope ratio measurements of low-molecular weight intermediary hydrocarbons (organic acids and alcohols) and serve as a proof of concept for making C and H isotope measurements on more complicated mixtures of solid-phase hydrocarbons and intermediary products produced during high temperature and high pressure synthesis on mineral-catalyzed surfaces. The impetus for this work stems from recently reported observations of methane detected within the Martian atmosphere [1-4], coupled with evidence showing extensive water-rock interaction during Martian history [5-7]. Methane production on Mars could be the result of synthesis by mineral surface-catalyzed reduction of CO2 and/or CO by Fischer-Tropsch Type (FTT) reactions during serpentization reactions [8,9]. Others have conducted experimental studies to show that FTT reactions are plausible mechanisms for low-molecular weight hydrocarbon formation in hydrothermal systems at mid-ocean ridges [10-12]. Further, recent experiments by Fu et al. [13] focus on examining detailed C isotope measurements of hydrocarbons produced by surface-catalyzed mineral reactions. Work described in this paper details the experimental techniques used to measure intermediary organic reaction products (alcohols and organic acids).

Socki, Richard A.

Carbon Isotopes of Alkanes in Hydrothermal Abiotic Organic Synthesis Processes at High Temperatures and Pressures: An Experimental Study

Observation of methane in the Martian atmosphere has been reported by different detection techniques [1-4]. With more evidence showing extensive water-rock interaction in Martian history [5-7], abiotic formation by Fischer-Tropsch Type (FTT) synthesis during serpentization reactions may be one possible process responsible for methane generation on Mars [8, 9]. While the experimental studies performed to date leave little doubt that chemical reactions exist for the abiotic synthesis of organic compounds by mineral surface-catalyzed reactions [10-12], little is known about the reaction pathways by which CO2 and/or CO are reduced under hydrothermal conditions. Carbon and hydrogen isotope measurements of alkanes have been used as an effective tool to constrain the origin and reaction pathways of hydrocarbon formation. Alkanes generated by thermal breakdown of high molecular weight organic compounds have carbon and hydrogen isotopic signatures completely distinct from those formed abiotically [13-15]. Recent experimental studies, however, showed that different abiogenic hydrocarbon formation processes (e.g., polymerization vs. depolymerization) may have different carbon and hydrogen isotopic patterns [16]. Results from previous experiments studying decomposition of higher molecular weight organic compounds (lignite) also suggested that pressure could be a crucial factor affecting fractionation of carbon isotopes [17]. Under high pressure conditions, no experimental data are available describing fractionation of carbon isotope during mineral catalyzed FTT synthesis. Thus, hydrothermal experiments present an excellent opportunity to provide the requisite carbon isotope data. Such data can also be used to identify reaction pathways of abiotic organic synthesis under experimental conditions.

Fu, Qi

Hydrogen Isotope Measurements of Organic Acids and Alcohols by Pyrolysis-GC-MS-TC-IRMS: Application to Analysis of Experimentally Derived Hydrothermal Mineral-Catalyzed Organic Products

We report results of experiments to measure the H isotope composition of organic acids and alcohols. These experiments make use of a pyroprobe interfaced with a GC and high temperature extraction furnace to make quantitative H isotope measurements. This work compliments our previous work that focused on the extraction and analysis of C isotopes from the same compounds [1]. Together with our carbon isotope analyses our experiments serve as a "proof of concept" for making C and H isotope measurements on more complex mixtures of organic compounds on mineral surfaces in abiotic hydrocarbon formation processes at elevated temperatures and pressures. Our motivation for undertaking this work stems from observations of methane detected within the Martian atmosphere [2-5], coupled with evidence showing extensive water-rock interaction during Mars history [6-8]. Methane production on Mars could be the result of synthesis by mineral surface-catalyzed reduction of CO2 and/or CO by Fischer-Tropsch Type (FTT) reactions during serpentization [9,10]. Others have conducted experimental studies to show that FTT reactions are plausible mechanisms for low-molecular weight hydrocarbon formation in hydrothermal systems at mid-ocean ridges [11-13]. Our H isotope measurements utilize an analytical technique combining Pyrolysis-Gas Chromatograph-Mass Spectrometry-High Temperature Conversion-Isotope Ratio Mass Spectrometry (Py-GC-MS-TC-IRMS). This technique is designed to carry a split of the pyrolyzed GC-separated product to a Thermo DSQII quadrupole mass spectrometer as a means of making qualitative and semi-quantitative compositional measurements of separated organic compounds, therefore both chemical and isotopic measurements can be carried out simultaneously on the same sample.

Socki, Richard A.

Summary Report Of The FY25 Reactor Physics Verification And Validation Exercises In The Advanced Reactor Technologies - Gas-cooled Reactor Program

Valdiation and verification of numerical tools is critical for ensuring reasonable predictions for design scoping, licensing, and safety analsyis. In this report, two reactor physics verification and validation exercises are presented. The first of these exercises focuses on burnup analysis with data from the Advanced Gas Reactor (AGR) program. Simulations are performed with Monte Carlo N-Particle (MCNP) and are compared with the experimental measurements for the AGR 1 and 2 experiments that utilize both UCO and UO2 fuel. The second exercises utilizes data from the HTR-Proteus experiments to perform reactor physics validation. Specifications of the experimental facility are provdied, along with a demonstration of initial modeling efforts in Serpent for one of the determistic packing experiments. Both cases are part of the Generation-IV international forum (GIF) Very High-Temperature Reactor (VHTR) Computational Methods, Validation, and Benchmarking (CMVB) program, an international collaborative organization dedicated to the verification and validation of High-Temperature Gas-Cooled Reactor (HTGR) analysis. Participation in the CMVB allows the US Department of Energy (DOE) to leverage these existing validation activities to provide extra value through benchmarking activities with other CMVB members.

and Benchmarking (CMVB) program

The HTR-Proteus Benchmark: Analysis and Use as a Verification and Validation Case

This presentation outlines the evaluation and application of the HTR-Proteus benchmark as a verification and validation (V&V) case for advanced reactor modeling tools. The work supports the U.S. Department of Energy’s HALEU Availability Program (HAP) and the joint DOE/NRC DNCSH project, which aims to reduce criticality safety uncertainties in commercial-scale HALEU fuel cycle and transportation systems. The HTR-Proteus experiments, conducted at the Paul Scherrer Institute, provide high-fidelity data for TRISO-fueled, graphite-moderated pebble bed reactors with high neutron leakage—conditions relevant to HALEU transport scenarios. This study focuses on Cores 4.2 and 4.3 of the HTR-Proteus benchmark, analyzing key sources of uncertainty including pebble packing, TRISO particle positioning, and core height. Using Project Chrono for realistic pebble geometries and Serpent, SHIFT, and MCNP for neutronics simulations, the study quantifies the impact of these uncertainties on the effective multiplication factor (keff). Results show that a sample size of 110 pebble configurations is sufficient to converge keff, with ±30 pcm uncertainty due to packing randomness. TRISO positioning and core height variations also significantly influence keff, highlighting the importance of detailed modeling in V&V efforts. The benchmark serves as a valuable test case for validating the Griffin reactor physics code and improving confidence in HALEU system simulations.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS