Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fuel lattices”

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

Experimental confirmation of first-principles thermal conductivity in Zirconium-doped ThO 2

The degradation of thermal conductivity in advanced nuclear fuels due to the accumulation of fission products and irradiation-induced defects is inevitable, and must be considered as part of safety and efficiency analyses of nuclear reactors. Here, this study examines the thermal conductivity of a zirconium-doped ThO 2 crystal, synthesized via the hydrothermal method using a spatial domain thermoreflectance technique. Zirconium is one of the soluble fission products in oxide fuels that can effectively scatter heat-carrying phonons in the crystalline lattice of fuel. Thus, thermal property measurements of zirconium-doped ThO 2 single crystals provide insights into the effects of substitutional zirconium doping, isolated from extrinsic factors such as grain boundary scattering. The experimental results are compared with first-principles calculations of the lattice thermal conductivity of ThO 2 , employing an iterative solution of the Peierls-Boltzmann transport equation. Additionally, the non-perturbative Green's function methodology is utilized to compute phonon-point defect scattering rates, accounting for local distortions around point defects, including mass difference changes, interatomic force constants, and structural relaxation. The congruence between the predicted results from first-principles calculations and the measured temperature-dependent thermal conductivity validates the computational methodology. Furthermore, the methodologies employed in this study enable systematic investigations of thermal conductivity reduction by fission products, potentially leading to the development of more accurate fuel performance codes.

36 - MATERIALS SCIENCE↗

Reactor physics benchmark of Westinghouse PWR core design suite for high burnup/high enrichment fuel - Part 1: Pin and lattice

Westinghouse has performed a comprehensive set of code-to-code benchmark comparisons to corroborate application of its state-of-the-art NEXUS-based code suite to high burnup and high enrichment (HB/HE) fueled PWR core design and reload analysis. This paper focuses on analysis of the predictions for pin and lattice designs representative of HB/HE core reloads from the Westinghouse advanced lattice code PARAGON2, the Monte-Carlo code SERPENT2 and MPACT, the deterministic neutronic component of VERA. A companion paper presented at this conference analyzes results obtained with the Westinghouse ANC9 core simulator, using lattice data provided by PARAGON2 with the NEXUS cross-section methodology, and MPACT for 2D and 3D core configurations incorporating HB/HE fuel reloads to achieve cycle length extension to 24-month in high-power density PWRs. The results of the pin and lattice benchmark presented in this paper show that the neutronic predictions from PARAGON2 for HB/HE fueled PWRs are in remarkable agreement with predictions from both SERPENT2 and MPACT including for depletion to high fuel burnup values, and employing a set of challenging lattices which feature combined use of burnable absorbers to obtain the required reactivity hold-down to support PWR transition to HB/HE fuel and 24-month operational cycles. As corroborated by this analysis and the analysis contained in the companion paper, the Westinghouse NEXUS-based PWR core analysis package can be confidently applied by industry for design analysis of PWR core reloads adopting HB/HE fuel and properly support the industry commercialization plans for in-reactor deployment. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Nuclear Fusion Reactions in Deuterated Metals

Nuclear fusion reactions of D-D are examined in an environment comprised of high density cold fuel embedded in metal lattices in which a small fuel portion is activated by hot neutrons. Such an environment provides for enhanced screening of the Coulomb barrier due to conduction and shell electrons of the metal lattice, or by plasma induced by ionizing radiation (γ quanta). We show that neutrons are far more efficient than energetic charged particles, such as light particles (e−, e+) or heavy particles (p, d, α) in transferring kinetic energy to fuel nuclei (D) to initiate fusion processes. It is well-known that screening increases the probability of tunneling through the Coulomb barrier. Electron screening also significantly increases the probability of large- versus small-angle Coulomb scattering of the reacting nuclei to enable subsequent nuclear reactions via tunneling. This probability is incorporated into the astrophysical factor S(E). Aspects of screening effects to enable calculation of nuclear reaction rates are also evaluated, including Coulomb scattering and localized heating of the cold fuel, primary D-D reactions, and subsequent reactions with both the fuel and the lattice nuclei. The effect of screening for enhancement of the total nuclear reaction rate is a function of multiple parameters including fuel temperature and the relative scattering probability between the fuel and lattice metal nuclei. Screening also significantly increases the probability of interaction between hot fuel and lattice nuclei increasing the likelilhood of Oppenheimer-Phillips processes opening a potential route to reaction multiplication. We demonstrate that the screened Coulomb potential of the target ion is determined by the nonlinear Vlasov potential and not by the Debye potential. In general, the effect of screening becomes important at low kinetic energy of the projectile. We examine the range of applicability of both the analytical and asymptotic expressions for the well-known electron screening lattice potential energy Ue, which is valid only for E >> Ue (E is the energy in the center of mass reference frame). We demonstrate that for E ≤ Ue, a direct calculation of Gamow factor for screened Coulomb potential is required to avoid unreasonably high values of the enhancement factor f (E) by the analytical—and more so by the asymptotic—formulas.

Vladimir Pines↗

Assessment of Core Physics Characteristics of Extended Enrichment and Higher Burnup LWR Fuels using the Polaris/PARCS Two-Step Approach (Vol. 2: BWR Fuel)

Nuclear fuel vendors and utilities are currently investigating changes to fuel contents and fuel designs for more economical and safer reactor operations. Extending cycle lengths beyond 18-month cycles for pressurized water reactors (PWRs) and 24-month cycles for boiling water reactors (BWRs) requires extending fuel enrichments beyond the current 5 wt % 235U limit. Therefore, low-enriched uranium plus (LEU+) fuel is expected to be used in current light-water reactor fleets in the near term. LEU+ is a subset of high-assay low-enriched uranium (HALEU) and is a term to describe fuel enrichments above 5% up to 10%. A series of studies were conducted at Oak Ridge National Laboratory (ORNL) to compare low-enriched uranium (LEU) with LEU+ fuel with respect to isotopic fuel content, lattice parameters, and core physics to identify any challenges in operation, storage, and transportation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Thermal and hydraulic performance of volumetrically heated triply periodic minimal surface heaters

Advanced heat transfer surfaces, mainly Triply Periodic Minimal Surfaces (TPMSs) have great potential to increase heat transfer performance over traditional heat transfer surfaces due to their tortuous flow path and higher surface area-volume ratio. The emergence of additive manufacturing of nuclear fuel motivates experimental investigation of such advanced geometries in the context of a nuclear reactor core. Here, this study employed a novel method to explore the performance of gyroid and diamond TPMS geometries. Volumetrically heated heaters were additively manufactured from conductive polymer filament, and airflow was used to simulate the convective heat transfer from the core within a nuclear reactor. The mass flow and volumetric generation rates were varied, local Nusselt number correlations were obtained using embedded thermocouples to measure the solid and fluid temperature, and friction factor correlations were developed using differential pressure measurements. The gyroid TPMS presented a friction factor 1.5 times higher than the diamond TPMS and 90 times that of the rod bundle. TPMS geometries showed Nusselt numbers 8–10 times higher than that of the rod bundle, also maintaining colder internal temperatures, which suggested that they could sustain 250–275 % higher power density than the rod bundle for a given centerline temperature. Hence, having the potential to greatly increase the safety margins and power output of a nuclear reactor core with similar volume.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

MOOGLE: A Multi-Objective Optimization tool for three-dimensional nuclear fuel assembly design

MOOGLE is a new genetic algorithm methodology for the three-dimensional design of nuclear fuel assemblies. MOOGLE uses common fuel rod types as the decision variable to develop a suite of three-dimensional fuel assemblies to provide optimized solutions to the design problem. Pressurized Water Reactor (PWR) fuel assemblies were optimized using IFBA and gadolinium (Gd 2 O 3 ) as burnable poisons in order to compare how burnable poison choice affects optimization results. Boiling Water Reactor (BWR) fuel bundles were also optimized using three unique fuel rod palettes to study how the size of the design space affects optimization results. Burnable poison analysis showed that utilizing IFBA and Gd 2 O 3 as burnable poisons produced the best and widest range of optimized solutions. Further, BWR fuel bundle optimization results indicate that the inclusion of additional fuel rod types produced a wider solution space but did not improve optimization results for regions explored using fewer unique fuel rods. These tests demonstrate MOOGLE's ability to analyze the tradeoffs between the inclusion of different fuel elements and their effects on assembly performance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

SCALE Lattice Physics Code Assessments of Accident Tolerant Fuel

This report highlights accident-tolerant fuel (ATF) code assessment activities performed under Project NRCHQ- 60-17-T-0017, Lattice Physics Enhancements and Assessment. ATF covers a broad range of advanced fuel and clad designs for light water reactors (LWRs) to enhance performance under several accident conditions. Several ATF concepts are anticipated to be deployed as lead test rods (LTRs) or lead test assemblies (LTAs) within the next five years. The purpose of this work is to assess the predictive capabilities of NRC neutronics codes that underpin various licensing calculations. ATF designs use different fuel and clad materials compared to standard UO 2 fuel and zirconium-alloy claddings (hereafter UO 2 -Zry). These new materials and geometrical designs need to be assessed to quantify the impact of nuclear data uncertainties on quantities of interest (QOIs) in licensing calculations and the impact of modeling approximations which may be valid for UO 2 -Zry but not for ATF. This report outlines a systematic approach for ATF neutronics code assessment which includes sensitivity and uncertainty (S/U) analysis of nuclear data, identification of experimental benchmark and gaps for code validation, investigation of modeling approximations, and code-to-code comparisons of calculated QOIs against high-fidelity reference continuous energy (CE) Monte Carlo (MC) calculations. This report focuses on the assessment of the SCALE/Polaris lattice physics code for reactor safety analysis. Polaris lattice physics calculations generate few group (FG) cross sections for PARCS full-core calculations. (Full-core analysis, spent fuel analysis, and severe accident analysis will be performed in future work.) The selected ATF concepts for this report include Cr 2 O 3 and Al 2 O 3 -Cr 2 O 3 -doped UO 2 fuel, U 3 Si 2 fuel, FeCrAl cladding, SiC cladding, and Cr-coated cladding.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Scale effects on core design, fuel costs, and spent fuel volume of pressurized water reactors

The desire to improve the economic competitiveness and deployment pace of nuclear energy through modularization, manufacturing, and series production had led to the development of smaller size reactors. As the standard 17x17 fuel technology is mainly maintained in the pressurized water reactors (PWRs) category, this translates into a lower number of fuel assemblies in the core and sometimes a reduced fuel height. To assess the impact of such scale change in core design on fuel cycle cost and spent fuel volume, a scoping analysis tool is developed based on infinite lattice calculations, leakage, fuel management reduced models, and levelized unit cost of electricity (LCOE) estimate. As such, cost dynamics driven by fuel specific power, burnup, core leakage, feed, cycle length, fuel assembly height as well as uranium market data are captured with consistent set of assumptions and analysis methods. A selection of 5 reactor designs representative of leading PWR developers is assessed and compared. Pursuing higher specific powers and optimal burnups are highlighted as the main fuel cost reduction drivers, nevertheless, practical limitations and opportunities must be evaluated to establish the feasibility of such enhanced fuel operation. In consequence, a detailed core design is performed using SIMULATE3 code for 5 PWR variations including natural and forced coolant circulation modes, two reactor scales, power densities of 73, 112, and 123 kW/l and higher discharge burnups. Design and optimization are performed at the lattice level, for the reflector, and at the core loading level. Satisfactory steady-state operation including power distribution, coolant operating limits, and reactivity requirements are analyzed and reported in this paper. The fuel economics of the detailed core designs confirm the scoping analysis findings. Despite the unlocked power uprates in small PWRs, the achievable burnup for a given fuel specific power requires more enrichment and shorter fuel height results in higher fabrication costs per mass of fuel, which makes scaling down core size a more expensive endeavor on the fuel cycle front. Spent fuel volumes are reported for the PWRs designed in this paper. Furthermore, these volumes are driven by the core average discharge burnup regardless of the scale in consideration. Additional cost and core performance aspects related to heavy reflector gains, fuel-reflector substitution, and disposal cost policy in the U.S. are examined.

42 ENGINEERING↗

Evidence of Xe-incorporation in the bubble superlattice in irradiated U-Mo fuel

For years, researchers have reported competing ideas on the formation of fission gas bubble super lattices in metallic fuels, which was postulated to be comprised of elemental xenon (Xe) atoms. However, the chemical and physical arrangement of these elements within this gas bubble superlattice has not been verified. In this contribution, Xe was chemically profiled using atomic resolution scanning transmission electron microscopy (STEM) and atom probe tomography (APT) techniques in irradiated uranium-molybdenum (U-Mo) fuels. These complementary techniques provide conclusive evidence of Xe presence in bubble superlattice up to fission densities of 4.5×10 21 fissions/cm 3 and provide quantitative assessment of Xe content within the bubble superlattice. Based on the results of simultaneous imaging and spectroscopy using STEM, the chemical composition of the Xe bubble superlattice present in fission densities less than 4.5x10 21 fissions/cm 3 in U-Mo fuel was measured and found to be up to 8±1.3 atomic %. APT data complements STEM findings on Xe distribution in irradiated U-Mo fuel sample. APT showed a chemical composition of 0.7±0.1 atomic%. Finally, this study concludes by providing conclusive evidence that Xe is not only present within an atomically arranged bubble superlattice but provides fundamental insight the physical state of Xe in low enriched U-Mo monolithic fuel.

36 MATERIALS SCIENCE↗

Review—Development of Highly Active and Durable Hybrid Compressive Platinum Lattice Catalysts for Polymer Electrolyte Membrane Fuel Cells: Mathematical Modeling and Experimental Work

This review provides a comprehensive overview on the development of highly active and durable platinum catalysts with ultra-low Pt loadings for polymer electrolyte membrane fuel cells (PEMFCs) through a combined mathematical modeling and experimental work. First, simulation techniques were applied to evaluate the validity of the Tafel approximation for the calculation of the mass activity (MA) and specific activity (SA). A one-dimensional agglomeration model was developed and solved to understand the effects of exchange current density, porosity, agglomerate size, Nafion ® film thickness, and Pt loading on the MA and SA. High porosity (> 60%) and agglomerations at high Pt loadings cause the loss of the Tafel approximation and consequently the decrease in MA and SA. A new structure parameter was introduced to estimate the real porous structure using the fractal theory. The volumetric catalyst density was corrected by the fractal dimension (measured by Hg porosimetry), which gave a good agreement with the experimental values. The loading-dependent Tafel equation was then derived, which contains both the utilization and the non-linear scaling factor. Second, activated carbon composite support (ACCS) with optimized surface area, porosity, pore size, and pore size distribution was developed. The hydrophilic/hydrophobic ratio, structural properties (amorphous/crystalline ratio), and the number of active sites were optimized through metal-catalyzed pyrolysis. Stability of ACCS and Pt/ACCS were evaluated using an accelerated stress test (AST). The results indicated that Pt/ACCS showed no significant loss of MA and power density after 5,000 cycles at 1.0–1.5 V, while the commercial Pt/C catalysts showed drastic losses of MA and power density. Finally, monolayers of compressed Pt (core–shell-type Pt 3 Co 1 ) catalysts were structured by diffusing Co atoms (previously embedded in ACCS) into Pt. Compressive Pt lattice (Pt * ) catalysts were synthesized through an annealing procedure developed at the University of South Carolina (USC). The Pt * /ACCS catalyst showed high initial power density (rated) of 0.174 g Pt kW −1 and high stability (24 mV loss) at 0.8 A cm −2 after 30,000 cycles (0.6–1.0 V). The outstanding performance of Pt * /ACCS is due to the synergistic effect of ACCS and compressive Pt * lattice.

25 ENERGY STORAGE↗

Electrochemical loading enhances deuterium fusion rates in a metal target

Nuclear fusion research for energy applications aims to create conditions that release more energy than required to initiate the fusion process1. To generate meaningful amounts of energy, fuels such as deuterium need to be spatially confined to increase the collision probability of particles2, 3–4. We therefore set out to investigate whether electrochemically loading a metal lattice with deuterium fuel could increase the probability of nuclear fusion events. Here we report a benchtop fusion reactor that enabled us to bombard a palladium metal target with deuterium ions. These deuterium ions undergo deuterium–deuterium fusion reactions within the palladium metal. We showed that the in situ electrochemical loading of deuterium into the palladium target resulted in a 15(2)% increase in deuterium–deuterium fusion rates. This experiment shows how the electrochemical loading of a metal target at the electronvolt energy scale can affect nuclear reactions at the megaelectronvolt energy scale.

Chen, Kuo-Yi↗

On-the-fly response function generation method for composite coarse mesh

The hybrid stochastic deterministic transport code COMET, based on the incident response expansion theory, is used to model reactor cores with high fidelity and formidable computational speed. COMET models a reactor core using a library of incident flux response expansion coefficients that are pre computed for all the unique lattice cells (e.g., fuel assemblies, reflector blocks, etc.) in the core. In order to further improve its computational efficiency in pre-calculating the response library a new response function generation method is developed to compute the response functions for the composite coarse meshes made of a smaller set of unique lattices on the fly within the COMET's deterministic transport core sweep. The efficiency is achieved by eliminating a number of unique lattices that can be made up from the reduced set of unique meshes on the fly. The numerical process consists of the following steps. First, the boundary condition on composite coarse mesh boundaries is projected onto the expansion basis to compute the incident flux moments on external surfaces of all the basic (reduced set of unique) coarse meshes. Secondly, the deterministic sweeping solver in COMET is used to converge on the outgoing/incoming flux expansion moments crossing interfaces between the basic coarse meshes. Thirdly, the response functions for the composite coarse meshes are constructed as a superposition on the fly. The new response function generation method was tested on 88 composite coarse meshes consisting of CANDU fuel bundles and moderator blocks. It was found that response functions generated by the new method agree very well with those generated by direct Monte Carlo calculations. The average and maximum relative differences in the surface-to-surface response coefficients computed by the two methods are 0.10% and 0.20%, respectively. Similarly, the average and maximum relative differences in the response fission densities are 0.13% and 0.43%, respectively. These discrepancies are within one standard deviation of the stochastic uncertainties. The new method is five times faster than the original direct Monte Carlo method. The size of the response function library for the new method is five times smaller than that for the original method, leading to significantly less requirement for the computer hard drive space and memory. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Numerical investigation and experimental validation of water condensation in the gas diffusion layer with different properties

Liquid water in porous media hinders the transport of reactant to catalyst layer, where the electrochemical reactions occur, which in turn, reduces performance and lifetime of a proton exchange membrane fuel cell (PEMFC). In this study, two types of gas diffusion layer (GDL) materials are used to study the two-phase saturation distribution in a PEMFC. The correlations between the effective and anisotropic transport properties and porosity of the two GDLs are solved by the pore scale model. In addition, the anisotropic liquid water permeability and the relationship between saturation pressure and capillary pressure is determined by using the Lattice-Boltzmann method. Fuel cell performance as well as liquid water distribution in the GDL using neutron radiography are obtained for validation. Under wet conditions, the cell with Freudenberg GDL performs better than that using Toray GDL, especially at high current densities. The results of the two-phase model simulation show that the peak water saturation for Toray GDL occurs in the catalyst layer and can reach 40%, while the peak water saturation for Freudenberg GDL occurs inside the GDL directly under the ribs and only reaches 25% saturation. The combined results provide key insights to enable high power density operation of a PEMFC through GDL material optimization.

42 ENGINEERING↗

Transition Core Modeling for Extended Enrichment & Accident-Tolerant Fuels Using Polaris/PARCS

Commercial light water reactor (LWR) operators and fuel vendors are currently interested in increasing the low-enriched uranium (LEU) fuel enrichments from the current limit of 5.0 $^w/_o$ $^{235}U$ up to 10 $^ w/_o$ $^{235}U$ (referred to as "LEU+") in their current fleets; they are also interested in using accident-tolerant fuel (ATF) with both LEU and LEU+ fuel. This report aims to identify modeling challenges and accuracy concerns in transition core analysis using the SCALE Polaris lattice physics code and U.S. Nuclear Regulatory Commission core simulator PARCS. At the time this study was started, no publicly available LEU+ core designs existed for boiling water reactor (BWR) or pressurized water reactor (PWR) systems. Therefore, fuel lattices were shuffled within a multi-assembly model to mimic neutronically challenging lattice combinations seen in transition cores, such as a fresh LEU+ lattice next to depleted LEU lattices. In addition to multi-assembly models, whole-core BWR transition core calculations were performed for ATF and LEU+ fuel using an existing Hatch-1 Cycle 3 core model. A whole-core BWR model was chosen due to the more heterogeneous core designs compared to those for a PWR core. Since the original core is an old checkerboard core design and no core or fuel design optimization was performed for the modeled fuel types, these core calculations were intended only to provide: (1) Comparisons of core characteristics of interest, such as the pin power distributions and peaking factors, Doppler temperature coefficients (DTCs), and control blade worths (CBWs) under challenging core designs, (2) Identification of reactor physics challenges in modeling LEU+ and ATF cores, and (3) A stress test for the Polaris/PARCS two-step modeling approach, including characterization of the relative accuracy for predicting characteristics of interest such as pin power distributions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Ab Initio Evaluation of Uranium Carbide S(α,β) and Thermal Neutron Cross Sections

Uranium Carbide (UC) is a nuclear fuel material which offers better neutron economy and lower fuel-cycle costs compared to conventional mixed-oxide fuels. UC’s lattice binding and dynamical properties impact thermal neutron scattering and low temperature epithermal resonance absorption. The Thermal Scattering Law (TSL) describes the scattering system available energy and momentum transfer states. There is no TSL evaluation for UC in the ENDF/B-VIII.0 database; herein, ab-initio lattice dynamics (AILD) techniques are invoked to calculate the phonon spectrum for UC using spin-orbit-coupling density functional theory (DFT). The TSLs, inelastic and elastic thermal scattering cross sections for Uranium and Carbon in UC, respectively, are calculated in FLASSH for use in higher fidelity reactor design calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Investigating Liquid Water Transport in Different Pore Structure of Gas Diffusion Layers for PEMFC Using Lattice Boltzmann Method

Proton exchange membrane fuel cells (PEMFC) require a gas diffusion layer (GDL) to aid in the transport of liquid fuel to the catalyst layer. In this work, direct modeling using the Lattice Boltzmann Method (LBM) was applied to X-ray CT scans of four different carbon gas diffusion layers to understand the mass transport properties through the samples. Three injection orientations were used to study local saturation levels, water evolution through the sample, and mass transport behavior at breakthrough conditions. The LBM, combined with computational fluid dynamic modeling techniques, can accurately predict liquid saturation at the macro and micro scale, which provides more insight into the mass transport phenomena through the GDL. The change of pore structure and orientation in both the in-plane and through-plane determines the path that liquid water must take, which could aid or impact PEMFC performance. The outcomes from this work will also benefit any research that needs knowledge of internal mass transport qualities of gas diffusion media.

Sepe, M.↗

Numerical implementation of the MOOSE subchannel module (SCM) algorithm

The subchannel module of moose (SCM), previously referred to as Pronghorn-SC, is a subchannel code designed to resolve single-phase flow fields and calculate the relevant flow variables, in nuclear reactor fuel-pin assemblies. The assemblies it models are: water-cooled, with bare fuel-pins in a square ducted, quadrilateral lattice and liquid-metal/water-cooled, with bare/wire-wrapped fuel pins in a hexagonal ducted, triangular lattice. Previous publications have presented the development, validation and verification of SCM. This work presents an comparative overview of the different solvers implemented within SCM. This includes the solver algorithm description, and a performance comparison. The two test cases chosen to demonstrate the solver performance were taken from the PSBT enthalpy mixing benchmark and the ORNL-19 pin benchmark.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

SIMULATeQCD: A simple multi-GPU lattice code for QCD calculations

The rise of exascale supercomputers has fueled competition among GPU vendors, driving lattice QCD developers to write code that supports multiple APIs. Moreover, new developments in algorithms and physics research require frequent updates to existing software. These challenges have to be balanced against constantly changing personnel. At the same time, there is a wide range of applications for HISQ fermions in QCD studies. This situation encourages the development of software featuring a HISQ action that is flexible, high-performing, open source, easy to use, and easy to adapt. In this technical paper, we explain the design strategy, provide implementation details, list available algorithms and modules, and show key performance indicators for SIMULATeQCD, a simple multi-GPU lattice code for large-scale QCD calculations, mainly developed and used by the HotQCD collaboration. The code is publicly available on GitHub.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗