SEARCH · Engineering Papers
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.
This fact sheet describes what Automatic Receptacle Controls (ARCs) are, how they work, ARC benefits and savings, and last, how they are being incorporated into energy codes. This fact sheet is meant to be a resource for energy codes adopters and influencers to help encourage them to adopt ARC provisions in their energy code.
Certain embodiments of the present invention are configured to facilitate analyzing computer code more efficiently. For example, by conducting a first level abstraction (e.g., symbolic interpretation and algebraic simplification) and a second level abstraction (e.g., generalization) of the computer code, the analysis may more accurately account for variations in the code that may occur as a result of register renaming, instruction reordering, choice of instructions, etc. while minimizing the cost of computations required to perform the analysis.
The present work describes the latest advances and progress in the development of the OpenMC Monte Carlo code, an open-source code originating from the Massachusetts Institute of Technology. First, an overview of the development workflow of OpenMC is given. Various enhancements to the code such as real-time XML input validation, state points, plotting, OpenMP threading, and coarse mesh finite difference acceleration are described. (authors)
In this paper, the work that has been done to implement variance reduction techniques in a three dimensional, multi group Monte Carlo code - Tortilla, that works within the frame work of the commercial deterministic code - Attila, is presented. This project is aimed to develop an integrated Hybrid code that seamlessly takes advantage of the deterministic and Monte Carlo methods for deep shielding radiation detection problems. Tortilla takes advantage of Attila's features for generating the geometric mesh, cross section library and source definitions. Tortilla can also read importance functions (like adjoint scalar flux) generated from deterministic calculations performed in Attila and use them to employ variance reduction schemes in the Monte Carlo simulation. The variance reduction techniques that are implemented in Tortilla are based on the CADIS (Consistent Adjoint Driven Importance Sampling) method and the LIFT (Local Importance Function Transform) method. These methods make use of the results from an adjoint deterministic calculation to bias the particle transport using techniques like source biasing, survival biasing, transport biasing and weight windows. The results obtained so far and the challenges faced in implementing the variance reduction techniques are reported here. (authors)
LATC is an essential part of COSINE code package, which stands for Core and System Integrated Engine for design and analysis. LATC performs 2D multi-group assembly transport calculation and generates few group constants and the required cross-section data for CORE, the core simulator code. LATC is designed to have the capability of modeling the API 000 series assemblies. The development is a continuously improved process. Currently, LATC uses well-proven technology to achieve the key functions. In the next stage, more advanced methods and modules will be implemented. At present, WIMS and WIMS improved format library could be read in LATC code. For resonance calculation, equivalent relation with rational approximations is utilized. For transport calculation, two options are available. One choice is collision probability method in cell homogenization while discrete coordinate method in assembly homogenization, the other is method of characteristics in assembly homogenization directly. For depletion calculation, an improved linear rate 'constant power' depletion method has been developed. (authors)
DYNSUB is a high-fidelity coupled code system consisting of the reactor simulator DYN3D and the sub-channel code SUBCHANFLOW. It describes nuclear reactor core behavior with pin-by-pin resolution for both steady-state and transient scenarios. In the course of the coupled code system's active development, super-homogenization (SPH) and generalized equivalence theory (GET) discontinuity factors may be computed with and employed in DYNSUB to compensate pin level homogenization errors. Because of the largely increased numerical problem size for pin-by-pin simulations, DYNSUB has bene fitted from HPC techniques to improve its numerical performance. DYNSUB's coupling scheme has been structurally revised. Computational bottlenecks have been identified and parallelized for shared memory systems using OpenMP. Comparing the elapsed time for simulating a PWR core with one-eighth symmetry under hot zero power conditions applying the original and the optimized DYNSUB using 8 cores, overall speed up factors greater than 10 have been observed. The corresponding reduction in execution time enables a routine application of DYNSUB to study pin level safety parameters for engineering sized cases in a scientific environment. (authors)
Monte Carlo codes are becoming more widely used for reactor analysis. Some of these applications involve the generation of diffusion theory parameters including macroscopic cross sections and diffusion coefficients. Two approximations used to generate diffusion coefficients are assessed using the Monte Carlo code MC21. The first is the method of homogenization; whether to weight either fine-group transport cross sections or fine-group diffusion coefficients when collapsing to few-group diffusion coefficients. The second is a fundamental approximation made to the energy-dependent P1 equations to derive the energy-dependent diffusion equations. Standard Monte Carlo codes usually generate a flux-weighted transport cross section with no correction to the diffusion approximation. Results indicate that this causes noticeable tilting in reconstructed pin powers in simple test lattices with L2 norm error of 3.6%. This error is reduced significantly to 0.27% when weighting fine-group diffusion coefficients by the flux and applying a correction to the diffusion approximation. Noticeable tilting in reconstructed fluxes and pin powers was reduced when applying these corrections. (authors)
The nuclear fuel cycle may be interpreted as a network or graph, thus allowing methods from formal graph theory to be used. Nodes are often idealized as nuclear fuel cycle facilities (reactors, enrichment cascades, deep geologic repositories). With the advent of modern object-oriented programming languages - and fuel cycle simulators implemented in these languages - it is natural to define a class hierarchy of facility types. Bright is a quasi-static simulator, meaning that the number of material passes through a facility is tracked rather than natural time. Bright is implemented as a C++ library that models many canonical components such as reactors, storage facilities, and more. Cyclus is a discrete time simulator, meaning that natural time is tracked through out the simulation. Therefore a robust, dependent type system was developed to enable inter-operability between Bright and Cyclus. This system is capable of representing any fuel cycle facility. Types declared in this system can then be used to automatically generate code which binds a facility implementation to a simulator front end. Facility model wrappers may be used either internally to a fuel cycle simulator or as a mechanism for inter-operating multiple simulators. While such a tool has many potential use cases it has two main purposes: enabling easy performance of code-to-code comparisons and the verification and the validation of user input.
Accurate cross sections are of great importance to Monte Carlo particle transport calculation codes which use pointwise cross sections. It is well known that cross sections have much dependence on temperature on account of Doppler effect which is particularly vital for resonance absorption cross sections. In general, pre-prepared continuous-energy libraries are used in Monte Carlo codes, which is feasible in many applications. However, in multi-physics calculations where neutronics and thermal hydraulics are coupled, cross sections at a fine interval of temperatures are desired, which brings a great challenge to computer memory. Consequently, different methods about on the fly Doppler broadening for cross sections are proposed, which can generate cross sections at the specific temperature points during neutron simulations. These methods can be roughly classified into two categories: one is to calculate the required cross sections by different algorithms, such as interpolation using a high order function expansion, multipole representation, and different integration algorithms based on Doppler broadening equation. The other, for instance, target motion sampling (TMS), does not calculate the Doppler broadening cross sections but takes into account Doppler effect by rejection technique. In this paper, several on the fly Doppler broadening algorithms based on Doppler broadening equation are implemented in the Reactor Monte Carlo code RMC, and then the accuracy and efficiency are discussed and compared. In this paper, several algorithms for on the fly Doppler broadening based on Doppler broadening equation have been implemented in RMC and the focus falls on the comparison of the Gauss Legendre quadrature and Gauss Hermite quadrature. Gauss Legendre quadrature is more accurate than Gauss Hermite quadrature but also more time consuming. Base temperature has a significant impact on Gauss Hermite quadrature while increasing base temperature to room temperature (300 K) can give satisfactory results. In addition, the proportion of overhead of Doppler broadening treatment in total computational time tends to decrease with an increase number of nuclides in fuel. For depletion calculation, both algorithms are in good agreement with reference and the relative errors of K{sub inf} are within 5 times of standard deviation. As for total computational time, the ratio between Gauss Legendre quadrature or Gauss Hermite quadrature and reference decreases and reaches 9.6 or 3.4 at last, respectively. Four important isotopes are chosen to compare the atom density and all relative errors are within 0.25%. After comparisons of different algorithms and different base temperatures, it can be found that Gauss Hermite quadrature broadening from 300 K can perform on the fly cross section Doppler broadening with least total computational time and acceptable precision. (authors)
The in-pile behavior of light-water-reactor (LWR) fuels is governed by the physics from various time and spatial scales. Accordingly, multi-scale modeling and simulation are needed to accurately predict the fuel performance at the engineering scale. To improve the traditional fuel performance modeling, the Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program targets at developing multi-scale, multiphysics modeling tools for in-pile fuel behavior prediction under various operating or accident scenarios. In traditional fuel performance modeling, the history of fuel is represented by a single parameter, fuel burn-up. Burn-up measures the energy that has been extracted from fuels. However, it is not an adequate measure of fuel history as it ignores the flux effect, which is critical for the irradiation-induced microstructural evolution in fuels. Moreover, the traditional fuel performance models are usually fitted to experimental data. Their capabilities of being extrapolated to new operating conditions and being extended to new types of fuels are strongly limited by the empirical or semi-empirical nature. To improve the burn-up based models, a microstructure based fuel approach, the MOOSE-BISON-MARMOT (MBM) tool-kit, is being developed under the fuel product line (FPL) of the NEAMS program for predictive fuel performance modeling. In the MBM approach, the evolution of fuel microstructure and changes in fuel properties are described by the MARMOT code. The transient fuel properties obtained from MARMOT are then fed into the BISON code for engineering scale fuel performance modeling. Both BISON and MARMOT are developed within the MOOSE framework, which was developed at Idaho National Laboratory. In this summary, the development and the capabilities of the MARMOT code will be introduced. In the past few years, several materials models have been developed including: (1) a grain growth model that describes the grain size evolution as a function of time, temperature and fuel porosity, considering the effects of impurity and porosity dragging and anisotropic grain boundary energy, (2) a fission gas release model that describes the fission gas diffusion in grain interior, fission gas bubble formation in grain interior and at grain boundaries, and fission gas release to the fuel-cladding gap, (3) a fracture model that describes the fracture initiation and propagation in grain interior and along grain boundaries, and (4) a thermal conductivity model that describes the change in thermal conductivity due to grain growth, gas bubble formation, fracture and radiation damage. Several other models are also being developed in MARMOT for hydride formation in Zr-based cladding, radiation damage and grain size subdivision. In most cases, these materials models are coupled with each other. For instance, to accurate predict the transient thermal conductivity, information on grain size, fission gas density, porosity and fracture is required. Further development is on going to improve all these models. These models are developed primarily for UO{sub 2} fuels used in LWRs. They can be extended to other fuels with possibly new development. All models in MARMOT will be rigorously verified and validated before being applied for fuel performance modeling. A version of MARMOT was recently released in September 2015 with an assessment report. (authors)
In this study, a preliminary structural assessment on a graphite component of a VHTR was done in accordance with the semi-probabilistic assessment procedures of the latest 2015 edition of the ASME Boiler and Pressure Vessel Code, Section III, Div. 5, Subpart HHA-Graphite Materials. Establishment of the Weibull distribution for the tensile strength is essential to set up the stress limits based on the probability of failure of the graphite materials. All parameters required in the ASME code were identified and evaluated for a nuclear graphite material, IG-110. The assessment of design level was done on one of fuel elements in a Korean VHTR design by two methodologies in the ASME code: (1) Simplified assessment and (2) Full assessment. The structural integrity of the fuel element was assured by both assessment methods. Also, the effect of the number of integration volume groups on the probability of failure in the full assessment was investigated and a recommendation for the conservative results was suggested. (authors)
A new neutron transport lattice physics code, cosLATC, which is an essential part in COSINE (Core and System Integrated Engine) code package, has been developed by SPICRI. It is mainly used for 2-D transport calculations on Pressurized Water Reactor (PWR) assemblies or pin cells. cosLATC generates few group constants for the core physics code cosCORE in COSINE package. With regards to cross-section data, WIMS-D 69 and 281 energy groups library is employed currently. cosLATC also has capability of resonance treatment, depletion and fuel temperature calculation. cosLATC uses modular tracking and parallelization technique in MOC transport calculation, which can save lots of memory and running time. In this paper, MOC module in cosLATC are introduced, the accuracy and efficiency of the module are verified by C5G7 benchmark and a PWR typical assembly. (authors)
Predictive modeling and simulation of the performance of a nuclear reactor and its fuel is a challenging task because of the large number of coupled physical phenomena that must be addressed. In addition to this intrinsic complexity, model uncertainty must be accounted for in any analysis if the model will be used to facilitate design or operational decisions which may impact safety and performance. Rigorous, structured uncertainty analyses are performed by first characterizing the model's input uncertainties and, by propagating the uncertainty through the model in order to estimate the output's uncertainty. This paper presents the state of advancement of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) project to characterize the uncertainty of the CFD code Nek5000 using the Dakota package. Nek5000 is a high order spectral element computational fluid dynamics (CFD) code developed at Argonne National Laboratory for high resolution spectral-filtered large eddy simulations (LESs) and unsteady Reynolds averaged Navier-Stokes (RANS) simulations. The Dakota package developed at Sandia National Laboratory can be integrated with many scientific and engineering codes to facilitate efficient, effective uncertainty quantification (UQ) and sensitivity analyses. Dakota is designed to interface with external packages/codes such as Nek5000 through a script-driven interface called fork. Such interface provides maximum flexibility to the users for pre- and post-processing tasks preceding and following each function evaluation. (authors)
System codes such as TRACE and RELAP5 have the option of adding boron to be dissolved in the water. In principle, the density and the viscosity of the water are not affected when the amounts of dissolved boron are very small. However, if the boron concentration is significant, water density and viscosity will be affected. Version 5.890 of the TRACE code has a new boron model with the option of modifying the water density and viscosity to account for the presence of dissolved boron. The input variable ISOLUT is used for this purpose. If ISOLUT=0, no boron is dissolved in the water. If ISOLUT=1, boron is dissolved in the water, and the liquid water density and viscosity are not modified by the dissolved boron. If ISOLUT=2, boron is dissolved in the water, and the liquid water density and viscosity are modified by the dissolved boron. This paper evaluates some features of this new boron model in the TRACE code.
The MSRs is a liquid-fuel reactor in which uranium, fission products and actinides are dissolved in the molten salt. The molten salt reactors are featured by circulating fuel, low pressure and high temperature. The MSRs meets many of the future goals of nuclear technology, in particular for what concerns as improved sustainability and unique characteristics in terms of actinide burning. The fluid nature of fuel in MSRs gives extra flexibility in reactor design, fuel fabrication and recycling. But because of the application of fluid fuel, MSRs have several special characteristics compared with the traditional solid-fuel reactors. Firstly, the transport of the delayed neutron precursors by the fluid flow reduces the contribution of delayed neutrons to the chain reaction.. Secondly, the fact that the fuel is dissolved in the coolant rather than separated from the coolant by the claddings results in a much stronger coupling phenomenon between the neutronics and thermal-hydraulics. Thirdly, the shape of the fuel is determined by the container because of the characteristics of fluid, the unstructured meshes are required to model the complex geometry accurately. However most of the conventional reactor codes were developed for solid-fuel reactor and therefore not capable of taking liquid fuel flow into account. To address all these particular modeling issues of MSRs, a multi-physics code named TANSY has been developed in our laboratory. In this paper, the benchmark problem based on the MSFR design was calculated to verify the code.
The R2Smesh ('Rigorous 2-step') code system, developed by KIT for the determination of shut-down dose rate distribution on mesh grid, can be utilized to provide photon radiation dose field around activated components which are moved from the irradiation site in fusion reactor to another external location. R2Smesh couples MCNP transport calculations and FISPACT inventory calculations through suitable interfaces utilizing MCNP mesh tally capability. Neutron flux and the decay gamma source are provided with high resolution on the mesh grid superimposed to the real geometry to take into account the spatial variations of the flux and the decay gamma source distributions. This feature enables the decay gamma source distribution to be transferred from the irradiation site in the reactor to any other external location for the determination of photon flux and dose rate distributions. In this work, we applied the feature of the portable decay gamma source in R2Smesh code system to a dose rate analysis for a tunnel connecting the ITER tokamak complex building (TCB) and the hot cell facility building (HCFB). This paper presents the use of the portability feature of the decay gamma source in R2Smesh interface code system developed by KIT for a dose rate analysis of the ITER tunnel connecting TCB and HCFB. This application shows that the decay gamma source of the activated components could be transferred successfully from one coordinate system to another coordinate system and superimposed to the geometry of the components. This is very useful when the activated components need to be transported from the reactor to somewhere else for maintenance, testing or storage and the dose rate analysis is requested. Another portable feature, which could convert the 3D decay gamma source distribution generated by R2Smesh system to the MCNP 'SDEF' card source description and avoid the use of source sampling subroutine, is under consideration. With this feature, one decay gamma source result for the specific reactor components could be shared by different institutions who need to use the same decay gamma source for different analysis. (authors)
This paper presents a joint work conducted at CEA Saclay and JAEA Tokai aimed at comparing the Monte Carlo codes TRIPOLI-4{sup R} and MVP on a selection of ICSBEP benchmarks. Our goal is to establish a common set of Monte Carlo input decks, as a basis for rigorous inter-code comparison in criticality-safety. As a reference, we will use the MCNP Criticality Validation Suite: other Monte Carlo developers might easily join this effort in the future. For the purpose of inter-code comparison, the TRIPOLI-4{sup R} and MVP input decks have been exactly translated from those of MCNP, without any further assumptions. Both TRIPOLI-4{sup R} and MVP have been run with the same ENDF/B-VII.0 evaluated nuclear data and, as far as possible, the same simulation options as in the original LANL work (same initial source, same number of active and discarded cycles and neutrons per cycle). (authors)
For application in criticality safety, numerical code systems require profound validation against experimental data. In terms of burn-up credit, this refers to the inventory determination as well as the subsequent criticality calculation. Within this framework, a dedicated, conservative estimation of uncertainties, bias and bias uncertainties for the resulting calculated multiplication factor k{sub eff} is indispensable. However, measures to reduce excess conservatism may become desirable, depending on the application case. For this reason GRS implemented a methodology to validate a calculation code system for criticality safety analyses for systems with irradiated nuclear fuel. The work at hand demonstrates this methodology by means of the GRS home-built inventory determination system KENOREST and ORNL's well-known SCALE 6 CSASS criticality safety calculation sequence. A number of 40 radiochemical analysis samples taken from SFCOMPO and elsewhere have been calculated to derive isotopic correction factors in different ways, depending on the required level of conservatism and the calculational efforts to be carried out. This includes bounding methods as well as different Monte Carlo sampling approaches. Furthermore, from the ICSBEP a number of 515 critical benchmark experiment configurations comprising water moderated low enriched UO{sub 2} fuel rods and water moderated mixed oxide (MOX) fuel rods have been selected and evaluated. Important absorber and structure materials are also covered. Various trending and sensitivity and uncertainty analysis methods have been applied, again to derive conservative estimates on the various contributions to k{sub eff} uncertainties. The overall methodology is applicable to low enriched, water moderated PWR uranium dioxide fuels in a burn-up range of about 10 to 60 GWd/tHM, and validates the reactivity contributions of the major uranium and plutonium isotopes. The methodology itself has no limitation which requires the exclusion of further isotopes as e.g. fission products. However, its application is constrained by the experimental database being publicly available. After the potential provision of an extended database in the future, the range of applicability can easily be extended to more nuclides. In addition, it is not limited to the code systems being used here but can also be applied using different numerical tools. For a generic PWR spent fuel pool rack and different enrichment and burn-up combinations, the impact of the current level of validation on k{sub eff} is demonstrated and quantified. (authors)