Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer Codes Computer Programs”

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 1,801 records · Page 100

Gyrokinetic understanding of the edge pedestal transport driven by resonant magnetic perturbations in a realistic divertor geometry

Self-consistent simulations of neoclassical and electrostatic turbulent transport in a DIII-D H-mode edge plasma under resonant magnetic perturbations (RMPs) have been performed using the global total-f gyrokinetic particle-in-cell code x-point gyrokinetic code (XGC), in order to study density pump-out and electron heat confinement. The RMP field is imported from the extended magneto-hydrodynamics code M3D-C1, taking into account the linear two-fluid plasma response. With both neoclassical and turbulence physics considered together, the XGC simulation reproduces two key features of experimentally observed edge transport under RMPs: increased radial particle transport in the pedestal region that is sufficient to account for the experimental pump-out rate and suppression of the electron heat flux in the steepest part of the edge pedestal. In the simulation, the density fluctuation amplitude of modes moving in the electron diamagnetic direction increases due to interaction with RMPs in the pedestal shoulder and outward, while the electron temperature fluctuation amplitude decreases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

MCNP6 Results for the Phase III Sensitivity Benchmark of the OCED/NEA Expert Group on Uncertainty Analysis for Criticality Safety Assessment

Within the last decade, there has been increasing interest in the calculation of cross section sensitivity coefficients of k eff for integral experiment design and uncertainty analysis. The OECD/NEA has an Expert Group devoted to Sensitivity and Uncertainty Analysis within the Working Party for Nuclear Criticality Safety. This expert group has developed benchmarks to assess code capabilities and performance for doing sensitivity and uncertainty analysis. Phase III of a set of sensitivity benchmarks evaluates capabilities for computing sensitivity coefficients. MCNP6 has the capability to compute cross section sensitivities for k eff using continuous-energy physics. To help verify this capability, results for the Phase III benchmark cases are generated and submitted to the Expert Group for comparison. The Phase III benchmark has three cases: III.1, an array of MOX fuel pins, III.2, a series of infinite lattices of MOX fuel pins with varying pitches, and III.3 two spheres with homogeneous mixtures of UF 4 and polyethylene with different enrichments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Terry Turbopump Expanded Operating Band Modeling and Simulation Efforts in Fiscal Year 2021 Extended Period of Performance (Final Report)

This report documents the progress made under the Terry Turbine Expanded Operating Band (TTEXOB) program's modeling and simulation (MODSIM) initiative at Sandia National Laboratories (SNL ). It describes the US Federal Fiscal Year 2021 (FY21) extended period-of-performance MODSIM work completed since the closure of FY20 with due reference to the Texas A&M University (TAMU) hybrid milestone 5/6 experimental program. This work, which falls under Milestone 7 of the program, provides a counterpart to the various experiments. The overall TTEXOB program and its milestone-based approach are described in the program's Summary Plan. Details of the individual milestone test plans can be found in the corresponding detailed test plan, e.g. the Milestone 3 and 4 Detailed Test Plan. SNL MODISM is conducted alongside experiments performed at TAMU, and SNL technical staff regularly consults with TAMU on the experimental program. In FY21, MELCOR code models and capabilities were exercised in two different contexts: experimental comparisons to the TAMU ZS-1 and GS-2, and stand-alone analyses of a station black-out (SBO) scenario in a generic boiling water reactor (BWR). Code to experiment comparisons met with fair success when turbine losses were well characterized as for the ZS-1 turbine. Both deterministic and Bayesian calibration processes were used to find a recommended turbine torque multiplier for ZS-1 type turbines. This process could be repeated for GS-2 type turbines if GS-2 losses were better understood. Stand-alone generic BWR SBO calculations revealed that three different modes of self-regulating turbopump behavior may be observed depending on certain modeling parameters and choices having to do with turbine nozzles. Aspects of this predicted behavior may have been observed in TAMU GS-2 experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

eqtools: Modular, extensible, open-source, cross-machine Python tools for working with magnetic equilibria

As plasma physics research for fusion energy transitions to an increasing emphasis on cross-machine collaboration and numerical simulation, it becomes increasingly important that portable tools be developed to enable data from diverse sources to be analyzed in a consistent manner. This paper presents eqtools, a modular, extensible, open-source toolkit implemented in the Python programming language for handling magnetic equilibria and associated data from tokamaks. eqtools provides a single interface for working with magnetic equilibrium data, both for handling derived quantities and mapping between coordinate systems, extensible to function with data from different experiments, data formats, and magnetic reconstruction codes, replacing the diverse, non-portable solutions currently in use. Moreover, while the open-source Python programming language offers a number of advantages as a scripting language for research purposes, the lack of basic tokamak-specific functionality has impeded the adoption of the language for regular use. Implementing equilibrium-mapping tools in Python removes a substantial barrier to new development in and porting legacy code into Python. In this paper, we introduce the design of the eqtools package and detail the workflow for usage and expansion to additional devices. Additionally, the implementation of a novel three-dimensional spline solution (in two spatial dimensions and in time) is also detailed. Finally, verification and benchmarking for accuracy and speed against existing tools are detailed. Wider deployment of these tools will enable efficient sharing of data and software between institutions and machines as well as self-consistent analysis of the shared data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

MCNP® Code Version 6.3.1: Verification & Validation Testing (Rev.1)

This report describes the verification and validation testing performed on MCNP® code version 6.3.1. The purpose of this report is to act as a compendium of test suite descriptions and results. Accordingly the document is divided into two parts. Part I (this part) describes each test suite in the following sections. Part II provides the results of testing each suite and comparisons to experimental and/or alternative computational results, as appropriate.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Improvements to PHISICS/RELAP5-3D{sup C} Capabilities for Simulating HTGRs

Accurate, high-resolution full core simulations of prismatic high temperature reactors for burnup and transient analysis remain a challenge. The cores are geometrically large, but some features are quite small (tri-structural isotropic (TRISO) particles, burnable poison pellets). The treatment of both neutron scattering in graphite and resonance capture are complex and not adequately captured using the methods traditionally used in High Temperature Gas-cooled Reactors (HTGR) and commonly used in Light Water Reactors (LWR). For transient analysis, temperature feedback is a function of the TRISO fuel form, but full core models cannot resolve phenomena at this scale without averaging over space and energy. For burnup calculations, a second level of heterogeneity must be resolved to capture the local effects of fuel and burnable poisons, while accurately propagating their effects through and between blocks. The HTGR Methods and Simulation group has been utilizing the INL-developed PHISICS/RELAP5-3D{sup C} (P/R), code suite for participation in the Organization for Economic Cooperation and Development (OECD) Modular High Temperature Gas Cooled Reactor (MHTGR) Benchmark, funded by the US Department of Energy (DOE) under the Advanced Reactor Technologies (ART) program. The MHTGR-350 benchmark consists of several lattice, steady-state and transient exercises, and the P/R code has been successfully applied to most of the cases defined. In addition to this code-to-code verification activity, it is planned that P/R will be validated using data from several Loss of Forced Cooling (LOFC) tests performed at the High Temperature Test Reactor (HTTR) in a joint effort with Japan Atomic Energy Agency (JAEA) under the Civil Nuclear Energy Working Group (CNWG) project. These verification and validation simulations will identify code and model discrepancies, sensitivities, and establish the current state-of-the-art in HTGR simulation. The Department of Astronautical, Electrical and Energy Engineering (DIAEE), 'Sapienza' University has been involved in the development of P/R since 2014 through various graduate and PhD student projects. In the framework of this collaboration the P/R optimization for application on the OECD MHTGR-350 and HTTR models will be carried out.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Accelerated application development: The ORNL Titan experience

The use of computational accelerators such as NVIDIA GPUs and Intel Xeon Phi processors is now widespread in the high performance computing community, with many applications delivering impressive performance gains. However, programming these systems for high performance, performance portability and software maintainability has been a challenge. In this paper we discuss experiences porting applications to the Titan system. Titan, which began planning in 2009 and was deployed for general use in 2013, was the first multi-petaflop system based on accelerator hardware. To ready applications for accelerated computing, a preparedness effort was undertaken prior to delivery of Titan. In this paper we report experiences and lessons learned from this process and describe how users are currently making use of computational accelerators on Titan.

97 MATHEMATICS AND COMPUTING↗

The Jefferson Lab High Power Broadband THz Facility

Jefferson Lab operates a broadband terahertz (THz) facility, whose light is generated by 300 fs FWHM 135 pC bunches of electrons at an energy of 125 MeV as they traverse a magnetic compression chicane prior to entering the optical cavity of the Jefferson Lab Free Electron Laser (FEL)[1-3]. The repetition rate of the electron bunches, and hence that of the light pulses, can be as high as 75 MHz continuous. Approximately 1 microJoule of energy is emitted for each pulse, which implies an average power of up to 75W, a peak power of several MW and electric fields of MV/m. Niche applications of this unique source are pump-probe dynamics, both in and out of equilibrium, and large-scale imaging. The beamline itself is designed to deliver an f/8 beam 2mm vertical × 3mm horizontal in size and 60% vertically polarized from a subtended angle of ~200 milliradians horizontal by ~135 milliradians vertical in the accelerator. The first mirror is copper-braid cooled. The philosophy of the optical design is to use a relay optics configuration to transport the beam via reflection off metal mirrors through a series of focal points, with the final one being at a diamond window. The optical scheme was developed using the Synchrotron Radiation Workshop (SRW)[4]. This code performs a full calculation of the electric field from a relativistic electron. It does not handle multiparticle coherent enhancement[5], but it is significantly different from all the other synchrotron radiation calculations because it retains a term called the Coulomb term in which the electric field falls off as 1/R2. In addition to calculating the electric field and the intensity, the code also allows one to propagate the field through the optical focusing elements and apertures. The beam exiting the final diamond window in the THz user lab passes into a series of vacuum chambers with considerable flexibility for attenuation, collimation, and/or beam splitting for imaging, spectroscopy and 2 photon pump-probe measurements. In this talk we will discuss the details of the philosophy, optical design and operating experience of this facility and also present ideas for discussion. We deeply appreciate the excellent support of our colleagues without which these facilities would not have been built and operated. This work was supported by the U.S. Dept. of Energy under contract DE-AC05-06-OR23177, and the Commonwealth of Virginia. The FEL facility is supported by the Office of Naval Research. 1. S. V. Benson et al, ?The VUV/IR/THz Free Electron Laser Program at Jefferson Lab?, Nucl. Instr. and Meth. A649 9 (2011). 2. G.L. Carr, M.C. Martin, W.R. McKinney, K. Jordan, G.R. Neil and G.P. Williams ?High Power Terahertz Radiation from Relativistic Electrons?, Nature 420 153-156 (2002). 3. J.M. Klopf et al. ?The Jefferson Lab high power THz user facility?, Nucl. Instrum. Methods A582, 114-116 (2007). 4. O. Chubar, P. Elleaume, "Accurate And Efficient Computation Of Synchrotron Radiation In The Near Field Region", proc. of the EPAC98 Conference, 22-26 June 1998, p.1177-1179. 5. C. J. Hirschmugl, M. Sagurton and G. P. Williams, ?Multiparticle coherence calculations for synchrotron radiation emission? Physical Review A44, 1316 (1991).

Klopf, John↗

Assessment of the Griffin Reactor Multiphysics Application Using the Empire Micro Reactor Design Concept

In late 2019, INL and ANL agreed to jointly develop the reactor physics code named Griffin based on the integration of the two code suites, MAMMOTH/Rattlesnake (INL) and MC2 - 3/PROTEUS (ANL). Griffin is being developed based on the MOOSE framework and MOOSE quality assurance procedures. This decision was made to be able to allow DOE-NE to efficiently invest funding to this area and to provide effective and timely support for existing and potential users; the latter includes industry and government organizations who are developing various types of advanced reactors in the near and long term. Since MAMMOTH/Rattlesnake has been developed based on the MOOSE framework, the INL/ANL Griffin development team agreed to build Griffin beginning with a merger of MAMMOTH and Rattlesnake into a single code and moving forward by implementing capabilities from the PROTEUS suite into Griffin. Moving forward, both ANL and INL efforts are equally invested in the Griffin project, with management support, to provide an advanced reactor multiphysics tool to assist in reactor design, optimization, and safety analysis. Much work remains in moving Griffin forward to migrate PROTEUS capabilities and to optimize performance to meet user needs. The main objective of this work is to assess the current status of Griffin capabilities in terms of performance and accuracy, to determine priorities for PROTEUS migration, and to identify capabilities and features to improve for supporting the code integration effort. For this assessment, the Empire micro reactor problem that was developed in the ARPA-E MEITNER program was selected as an advance reactor concept of interest to the technical community. The Empire reactor problem was expanded from its original incomplete specification to be a small heat-pipe-cooled micro reactor core with ~113 cm radius and 70 cm in height, composed of 18 fuel assemblies, 12 control drums, and beryllium radial and axial reflectors. In the current model, using 5 cm axial reflectors specified in the original Empire assembly model, more than 10% of neutrons leak axially and through the empty center safety hole, as well as through heat pipe channels in fuel assembly elements that extend through the top reflector region. Several calculation models of the core were defined for systematic assessment, including 2-D and 3-D fuel assemblies and whole cores with cylindrical boundaries. Cross sections were generated using Serpent 2, and meshes were produced using the Argonne mesh tool or the INL neutronics meshing tools combined with CUBIT. Cross sections and meshes were converted to the ISOXML and Exodus formats, respectively, so that Griffin and PROTEUS could use consistent data for solving the reactor problems. With the prepared cross sections and meshes, PROTEUS was run first to ensure that all input data were correctly generated and input options in terms of angle, mesh, and energy group were accurately determined. Comparisons against Serpent 2 solutions were made in terms of eigenvalue and pin power. The same calculations and comparisons were then conducted using Griffin. For the fuel assembly and whole core problems, the PROTEUS eigenvalues agreed well with reference Serpent 2 solutions within 100 and 30 pcm, respectively, and pin power differences relative to Serpent 2 were overall less than 2.2% and RMS 0.8% for the whole core models. This indicated that all input data were properly prepared. Using the same data, Griffin was run selecting the SAAF-CFEM SN solver with Legendre-Gaussian quadrature and NDA and DSA for acceleration. It was found that the SAAF-CFEM solver of Griffin required finer meshes to achieve eigenvalue and pin power solutions in good agreement with Serpent 2, consequently requiring more memory requirement and longer computation time. On the other hand, the SPH-Diffusion 2-D core calculations performed using Griffin were able to recover the exact eigenvalue from the reference Serpent 2 solutions, resulting in a pin-power distribution with an RMS of 0.6% and maximum absolute difference of less than 1.4%. The runtimes for SPH-Diffusion for the 2-D core were less than 3 minutes on 40 cores. During this evolution of this evaluation, many updates were made in Griffin by the Griffin development team of INL (focusing on software updates) and ANL (reviewing and supporting software updates) to complete this assessment. Observations from the code assessment are presented in the conclusion section of this report, followed by a discussion of recommendations for future work.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

SU-F-T-468: Efficient Scanning Data Analysis for TomoTherapy

Purpose: The standard TomoTherapy 2D scanning system requires manual analysis using spreadsheets or similar software. Manual analysis is time consuming and does not allow efficient use of the medical physicist or timely results that may indicate issues during the acquisition of data. The purpose of this study is to determine if an in-house solution can provide timely analysis of TomoTherapy scan data. Methods: A custom MATLAB program utilizing a streamlined graphical interface was developed to interpret the native scanning data acquired by the TomoTherapy Electrometer Measurement System (TEMS). The gold beam data provided by the manufacturer was used as a baseline for the acquired data. Two separate sets of scan data were used including longitudinal, transverse, and percent depth dose (PDD). The scan data were compared visually in an overlay and analytically for each beam width using the full width half max (FWHM) for longitudinal profiles and difference inside the field core for the transverse profiles as recommended by TG-148. A plot of the PDD percent error at each point along the curve was generated to show agreement. The time to analysis was measured on a laptop running an Intel i5 processor with 8GB of RAM. Results: The software program required on average less than 5 seconds to load and analyze the results. No conversion or reformatting of the data was required. This allowed real-time visualization of the comparison to the gold beam data as well as the parameters recommended by TG-148. The average percent difference was 0.3% for the longitudinal profiles and 0.1% for the transverse profiles. All PDD data was less than 1% different beyond dmax. Conclusion: We have developed a tool to analyze TomoTherapy scan data in real-time. This allows an efficient method to perform beam commissioning, annual QA, and adjust the beam parameters if necessary.

61 RADIATION PROTECTION AND DOSIMETRY↗

Exploiting variability for energy optimization of parallel programs

Here in this paper we present optimizations that use DVFS mechanisms to reduce the total energy usage in scientific applications. Our main insight is that noise is intrinsic to large scale parallel executions and it appears whenever shared resources are contended. The presence of noise allows us to identify and manipulate any program regions amenable to DVFS. When compared to previous energy optimizations that make per core decisions using predictions of the running time, our scheme uses a qualitative approach to recognize the signature of executions amenable to DVFS. By recognizing the "shape of variability" we can optimize codes with highly dynamic behavior, which pose challenges to all existing DVFS techniques. We validate our approach using offline and online analyses for one-sided and two-sided communication paradigms. We have applied our methods to NWChem, and we show best case improvements in energy use of 12% at no loss in performance when using online optimizations running on 720 Haswell cores with one-sided communication. With NWChem on MPI two-sided and offline analysis, capturing the initialization, we find energy savings of up to 20%, with less than 1% performance cost.

97 MATHEMATICS AND COMPUTING↗

Modeling and Simulation for Safeguards [Slides]

The purpose of this talk is to give an overview of the role of modeling and simulation in Safeguards R&D and introduce you to (some of) the tools used. Some definitions are: (1) Modeling - the representation, often mathematical, of a process, concept, or operation of a system, often implemented by a computer program; (2) Simulation - the representation of the behavior or characteristics of one system through the use of another system, especially a computer program designed for the purpose; and (3) Safeguards - the timely detection of diversion of significant quantities of nuclear material. The role of modeling and simulation are: (1) Calculate amounts of material (plant modeling); (2) Calculate signatures of nuclear material etc. (source terms); and (3) Detector performance (radiation transport and detection). Plant modeling software (e.g. FACSIM) gives the flows and amount of material stored at all parts of the process. In safeguards this allow us to calculate the expected uncertainty of the mass and evaluate the expected MUF. We can determine the measurement accuracy required to achieve a certain performance.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Exploring MPI Communication Models for Graph Applications Using Graph Matching as a Case Study

Traditional implementations of parallel graph operations on distributed memory platforms are written using Message Passing Interface (MPI) point to point communication primitives such as Send-Recv (blocking and nonblocking). Apart from this classical model, the MPI community has over the years added other communication models; however, their suitability for handling the irregular traffic workloads typical of graph operations remain comparatively less explored. Our aim in this paper is to study these relatively underutilized communication models of MPI for graph applications. More specifically, we evaluate MPI’s one-sided programming, or Remote Memory Access (RMA), and nearest neighborhood collectives using a process graph topology. There are features in these newer models that are intended to better map to irregular communication patterns, as exemplified in graph algorithms. As a concrete application for our case study, we use distributed memory implementations of a weighted graph matching algorithm to investigate performances of MPI-3 RMA and neighborhood collective operations compared to nonblocking Send- Recv. A matching in a graph is a subset of edges such that no two matched edges are incident on the same vertex. A maximum weight matching is a matching of maximum weight computed as the sum of the weights of matched edges. Algorithms for optimal matching have quadratic complexity and are difficult to parallelize. More importantly, execution of graph matching is dominated by high volume of irregular memory accesses, making it an ideal candidate for studying the effects of various MPI communication models on graph applications at scale. Our neighborhood collectives and RMA implementations yields up to 6× speedup over traditional optimized nonblocking Send-Recv implementations, on thousands of cores of the NERSC Cori supercomputer. We believe the lessons learned from this case study can be adopted to benefit a wider range of graph application codes.

Ghosh, Sayan↗

TUMME: Tsinghua University Minnesota Master Equation program

We report that TUMME is a program for assembling and solving master equations for gas-phase chemical kinetics based on chemically significant eigenmodes. TUMME has interfaces to the Gaussian, Polyrate, and/or MSTor output files that allow the master equation code to obtain the microcanonical flux coefficients needed for the coefficient matrix of the master equation. The flux coefficients for reactions with barriers can be calculated by multi-structural variational transition state theory with small-curvature tunneling (MS-VTST/SCT) or by simpler approximations to this such as conventional transition state theory without tunneling (also called RRKM theory). The flux coefficients for barrierless reactions are provided by a hard-sphere model. TUMME is written in double precision with Python 3; quadruple and octuple precision are also available for some subtasks in C++. The Python code can run in serial or parallel (MP or MPI), and the C++ code can run on a single processor or on multiple processors with OpenMP.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Flattening the Radial Temperature Profile across the Transformational Challenge Reactor Core

The Transformational Challenge Reactor (TCR) program is demonstrating an agile development approach to advanced nuclear reactor design, which has traditionally utilized a linear design process. In leveraging artificial intelligence, additive manufacturing, advanced materials, and cutting-edge modeling and simulation, the TCR program aims to minimize the high cost and lengthy deployment timelines now standard in the nuclear industry. Within a relatively short period of time, a robust and mature advanced gas-cooled reactor was iteratively designed under the TCR program, using these cutting-edge technologies. The TCR is a 3 MWt gas-cooled microreactor fueled with uranium nitride (UN) tristructural isotropic (TRISO) fuel particles. Though manufactured via traditional means, these UN TRISO particles are loaded into additively manufactured silicon carbide (SiC) cans [4]. Once loaded with TRISO particles, the SiC cans are densified using a chemical vapor infiltration process. The additively manufactured SiC enables significantly more freedom in the design of the fuel form than could ever be achieved using traditionally manufactured SiC. The helium coolant, pressurized to 5 MPa, enters the core at 300°C and nominally exits it at 500°C. Typically, the most thermally limiting components in any reactor design are the fuel assemblies in the core center. To provide a wide thermal margin in these central fuel assemblies, the flow may be biased toward the center of the core to more effectively cool these fuel assemblies with more power deposition and flatten the core’s radial temperature distribution. An analytical fluid model of the TCR core was developed to explore methods for biasing the flow away from the cooler outer fuel assemblies and towards the hotter inner ones. Higher-fidelity models developed in STAR-CCM+ 2020.3.1, a computational fluid dynamics code, were then utilized to verify the analytical model’s findings.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Chlorite Dissolution Kinetics at Variable pH and Temperatures up to 275 °C

The objective of this suite of experiments was to develop a useful kinetic dissolution expression for chlorite that would be applicable over a wide range of solution pH and temperature conditions representative of subsurface conditions in natural and/or engineering geothermal reservoirs. The resulting rate equation is dependent on both pH and temperature and utilizes two specific dissolution mechanisms (an “acid” and a “neutral” mechanism). Rate parameters were derived from data collected at LLNL and available published data from Lowson et al. (2007). The form of this rate equation (Equation 1) should be easy to incorporate into most existing reactive transport codes for use in prediction of rock-water interactions in EGS systems. Specifically, we find that the dissolution of chlorite, a sheet-silicate, is relatively slow at elevated temperatures (100-275 °C), compared to other framework silicate minerals for which higher temperature kinetic data are available (e.g., quartz and feldspars; see Palandri and Kharaka, 2004). This finding is in conflict with previously reported high activation energies for chlorite based on extrapolation of low-temperature experimental data to higher temperatures. Additionally, we note that the dissolution of chlorite does not increase under alkaline conditions, as has been noted for other minerals, but rather remains at the same level as that noted for neutral conditions, with only a weak dependence on temperature. Future work includes similar rate measurements and the derivation of useful rate equations for illite, smectite, and biotite from pH 3 to 10 and 100 to 275 °C. Geochemical alteration, changing stress fields, mass transport and heat transfer incorporated into computational models are needed to optimize geothermal energy production for EGS systems. The resulting mineral rate equations from this work can be directly incorporated in modeling efforts to fully assess the impacts of geochemical alteration on long-term fracture permeability for EGS systems across the Geothermal Program. Chlorite rate equations and data have been submitted to the Geothermal Data Repository.

15 GEOTHERMAL ENERGY↗

Improved Cross Section Generation Capability of Griffin in FY22

The Griffin code is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. The code includes a variety of deterministic steady-state transport solvers for fixed source, k-eigenvalue, adjoint, and subcritical multiplication as well as transient solvers for spatial dynamics with the improved quasi-static method. Griffin uses cross section data in the ISOXML format generated from external deterministic or Monte Carlo cross section generation codes. In recent years, the MC 2 -3 modules have been added to Griffin for fast reactor cross section generation, and the self-shielding application programming interface (SSAPI) was implemented in the ISOXML module for thermal reactor cross section generation. The on-the-fly slowing down method and double-heterogeneity treatment have been implemented to SSAPI and verified against particulate fuel-bearing graphite-moderated thermal reactor problems with high accuracy. This year, work has been focused on improving the cross section generation capability of ISOXML and streamlining the cross section generation procedures. In addition, the form function data were added to ISOXML in order to support the pin power reconstruction capability that was newly implemented in Griffin in this fiscal year. To facilitate the cross section generation using MC 2 -3 and SSAPI in Griffin, the cross section generation workflows have been set up for both fast and thermal spectrum reactors. The MOOSE action system tool was devised for fast spectrum problems, and the MOOSE stochastic tool was adopted to the branch calculation procedure for thermal spectrum problems. Meanwhile, to ensure the accuracy of group-constants, the thermal up-scattering kernel calculator accounting for resonance scattering was implemented in ISOXML, demonstrating the accurate computation of a Doppler-broadened scattering kernel of any Legendre order within a reasonable timescale. Other aspects of ISOXML, such as deletion solver and data, documentation, ISOXML file management, and the interface for Mixture, were improved as well.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

CTF Parallel Performance Improvements

CTF is a modernized version of COBRA-TF (Coolant Boiling in Rod Arrays-Two Fluid), which is a thermal hydraulic (T/H) sub-channel code based on a two-phase, two-fluid model that formulates the conservation equations of mass, energy, and momentum for three fields of vapor, continuous liquid, and entrained liquid droplets. Under the Consortium for Advanced Simulation of Light Water Reactors (CASL) program sponsored by the U.S. Department of Energy, CTF is being further developed and improved for modeling and simulation of light water reactors. The code improvements for predicting T/H responses of pressurized water reactor (PWR) cores are focused on its applications to specific challenge problems deemed significant to the industry, such as departure from nucleate boiling (DNB) and crud-induced power shift. The code is also being used for T/H feedback in normal operating condition depletions in PWR and boiling water reactor designs in the reactor core simulator, Virtual Environment for Reactor Applications Core Simulator (VERA-CS), being developed by CASL. CTF modifications under the CASL program include software optimization, new closure models, and parallelization for modeling full reactor core T/H responses. The purpose of this work was to improve the parallel performance of CTF. CTF was originally parallelized via a domain-decomposition strategy using the Message Passing Interface (MPI), in which the basis of a solution domain was a physical assembly in the model. In other words, a user who wanted to perform a parallel solution with CTF had to provide the code one processor per assembly in the model (e.g., 56 cores for a 56 assembly quarter-core model, or 193 cores for a 193 assembly full-core model). This design was rigid, demanding a very specific number of processors for a model; furthermore, when CTF was run in VERA-CS-which typically uses thousands of processors- many processors sat idle during the CTF solve because CTF was unable to use any more processors than there are physical assemblies. This was not necessarily a problem if the CTF execution speed was very fast, but the average time spent in a CTF solve for prototypical PWR operating conditions was significant (i.e., over 20 % of the simulation wall time was being spent in CTF for a standard core depletion simulation). Larger, more computationally intensive simulations performed in support of the CASL DNB challenge problem demonstrated a need for improvements in this original parallelization approach. The specific goal of this work was to allow CTF to use more processors than there are physical assemblies in the model, which would not only speed up CTF and improve the overall efficiency of VERA-CS but also give users more-flexible options for executing their simulations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Refine your search to explore more results.