Search NASA⌕ Search

Engineering topics

Freels, James D.

Publications and source records attributed to Freels, James D..

Deflection predictions of involute-shaped fuel plates using a fully-coupled numerical approach

This paper describes the modeling and simulation of fluid structure interactions (FSI) of involute-shaped fuel plates used in nuclear research reactors. We believe this to be the first time that this type of application is described in the literature using a fully-coupled, and monolithic, finite element approach. The simulations are validated against plate deflection data for the conceptual design of the Advanced Neutron Source Reactor (ANSR), which was envisioned to be the world’s most powerful nuclear research reactor for neutron scattering and other applications, but was ultimately never completed. The high performance of the ANSR creates a bounding envelope for involute-shaped research reactors such as that used in the High Flux Isotope Reactor (HFIR) at the Oak Ridge National Laboratory (ORNL) which is undergoing research for the conversion from highly-enriched uranium (HEU) to low-enriched uranium (LEU) fuel. As such, the findings from the present FSI analyses carried out herein for the ANSR plates provide good guidelines and inform designers what should be expected for the next generation of plates in the HFIR. Furthermore it is shown herein that the current approach can accurately capture the leading-edge deflections of the involute-shaped plates and simulations can predict the ‘S-shaped’ deflection of the first mode instilling confidence in the methodology.

42 ENGINEERING↗

Predicting Large Deflections of Multiplate Fuel Elements Using a Monolithic FSI Approach

As part of the Global Threat Reduction Initiative, the Oak Ridge National Laboratory is evaluating conversion of fuel for the High Flux Isotope Reactor (HFIR) from high-enriched uranium to low-enriched uranium. Currently, multiphysics simulations that model fluid-structure interaction phenomena are being performed to ensure the safety of the reactor with the new fuel type. A monolithic solver that fully couples fluid and structural dynamics is used to model deflections in the new design. A classical experiment is chosen to validate the capabilities of the current solver and the method. Here, a single-plate simulation with various boundary conditions as well as a five-plate simulation are presented. Finally, use of the monolithic solver provides stable solutions for the large deflections and the tight coupling of the fluid and structure and the maximum deflections are captured accurately.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Thermal Safety Analyses for the Production of Plutonium-238 at the High Flux Isotope Reactor

There has been a considerable effort over the previous few years to demonstrate and optimize the production of plutonium-238 ( 238 Pu) at the High Flux Isotope Reactor (HFIR). This effort has involved resources from multiple divisions and facilities at the Oak Ridge National Laboratory (ORNL) to demonstrate the fabrication, irradiation, and chemical processing of targets containing neptunium-237 ( 237 Np) dioxide (NpO 2 )/aluminum (Al) cermet pellets. A critical preliminary step to irradiation at the HFIR is to demonstrate the safety of the target under irradiation via documented experiment safety analyses. The steady-state thermal safety analyses of the target are simulated in a finite element model with the COMSOL Multiphysics code that determines, among other crucial parameters, the limiting maximum temperature in the target. Safety analysis efforts for this model discussed in the present report include: (1) initial modeling of single and reduced-length pellet capsules in order to generate an experimental knowledge base that incorporate initial non-linear contact heat transfer and fission gas equations, (2) modeling efforts for prototypical designs of partially loaded and fully loaded targets using limited available knowledge of fabrication and irradiation characteristics, and (3) the most recent and comprehensive modeling effort of a fully coupled thermo-mechanical approach over the entire fully loaded target domain incorporating burn-up dependent irradiation behavior and measured target and pellet properties, hereafter referred to as the production model. These models are used to conservatively determine several important steady-state parameters including target stresses and temperatures, the limiting condition of which is the maximum temperature with respect to the melting point. The single pellet model results provide a basis for the safety of the irradiations, followed by parametric analyses in the initial prototypical designs that were necessary due to the limiting fabrication and irradiation data available. The calculated parameters in the final production target model are the most accurate and comprehensive, while still conservative. Over 210 permutations in irradiation time and position were evaluated, and are supported by the most recent inputs and highest fidelity methodology. The results of these analyses show that the models presented in this report provide a robust and reliable basis for previous, current and future experiment safety analyses. In addition, they reveal an evolving knowledge of the steady-state behavior of the NpO 2 /Al pellets under irradiation for a variety of target encapsulations and potential conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Automating Nuclear-Safety-Related SQA Procedures with Custom Applications

Nuclear safety-related procedures are rigorous for good reason. Small design mistakes can quickly turn into unwanted failures. Researchers at Oak Ridge National Laboratory worked with COMSOL to define a simulation app that automates the software quality assurance (SQA) verification process and provides results in less than 24 hours.

97 MATHEMATICS AND COMPUTING↗

Research and Development of Multiphysics Models in Support of the Conversion of the High Flux Isotope Reactor to Low Enriched Uranium Fuel

The findings presented in this report are results of a five year effort led by the RRD Division of the ORNL, which is focused on research and development toward the conversion of the High Flux Isotope Reactor (HFIR) fuel from high-enriched uranium (HEU) to low-enriched uranium (LEU). This report focuses on the tasks accomplished by the University of Tennessee Knoxville (UTK) team from the Department of Mechanical, Aerospace, and Biomedical Engineering (MABE) that provided expert support in multiphysics modeling of complex problems associated with the LEU conversion of the HFIR reactor. The COMSOL software was used as the main computational modeling tool, whereas Solidworks was also used in support of computer-aided-design (CAD) modeling of the proposed LEU fuel design. The UTK research has been governed by a statement of work (SOW), which was updated annually to clearly define the specific tasks reported herein. Ph.D. student Isaac T. Bodey has focused on heat transfer and fluid flow modeling issues and has been aided by his major professor Dr. Rao V. Arimilli. Ph.D. student Franklin G. Curtis has been focusing on modeling the fluid-structure interaction (FSI) phenomena caused by the mechanical forces acting on the fuel plates, which in turn affect the fluid flow in between the fuel plates, and ultimately the heat transfer, is also affected by the FSI changes. Franklin Curtis has been aided by his major professor Dr. Kivanc Ekici. M.Sc. student Adam R. Travis has focused two major areas of research: (1) on accurate CAD modeling of the proposed LEU plate design, and (2) reduction of the model complexity and dimensionality through interdimensional coupling of the fluid flow and heat transfer for the HFIR plate geometry. Adam Travis is also aided by his major professor, Dr. Kivanc Ekici. We must note that the UTK team, and particularly the graduate students, have been in very close collaboration with Dr. James D. Freels (ORNL technical monitor and mentor) and have benefited greatly from his leadership and expertise in COMSOL modeling of complex physical phenomena. Both UTK and ORNL teams have used COMSOL releases 3.4 through 5.0 inclusive (with particular emphasis on 3.5a, 4.3a, 4.3b, and 4.4) for most of the work described in this report, except where stated otherwise. Just as in the performance of the research, each of the respective sections has been originally authored by respective authors. Therefore, the reader will observe a contrast in writing style throughout this document.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

3D COMSOL Simulations for Thermal Deflection of HFIR Fuel Plate in the "Cheverton-Kelley" Experiments

Three dimensional simulation capabilities are currently being developed at Oak Ridge National Laboratory using COMSOL Multiphysics, a finite element modeling software, to investigate thermal expansion of High Flux Isotope Reactor (HFIR) s low enriched uranium fuel plates. To validate simulations, 3D models have also been developed for the experimental setup used by Cheverton and Kelley in 1968 to investigate the buckling and thermal deflections of HFIR s highly enriched uranium fuel plates. Results for several simulations are presented in this report, and comparisons with the experimental data are provided when data are available. A close agreement between the simulation results and experimental findings demonstrates that the COMSOL simulations are able to capture the thermal expansion physics accurately and that COMSOL could be deployed as a predictive tool for more advanced computations at realistic HFIR conditions to study temperature-induced fuel plate deflection behavior.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Low-Enriched Uranium Fuel Conversion Activities for the High Flux Isotope Reactor, Annual Report for FY 2011

This report describes progress made during FY11 in ORNL activities to support converting the High Flux Isotope Reactor (HFIR) from high-enriched uranium (HEU) fuel to low-enriched uranium (LEU) fuel. Conversion from HEU to LEU will require a change in fuel form from uranium oxide to a uranium-molybdenum (UMo) alloy. With both radial and axial contouring of the fuel foil and an increase in reactor power to 100 MW, calculations indicate that the HFIR can be operated with LEU fuel with no degradation in performance to users from the current levels achieved with HEU fuel. Studies are continuing to demonstrate that the fuel thermal safety margins can be preserved following conversion. Studies are also continuing to update other aspects of the reactor steady state operation and accident response for the effects of fuel conversion. Technical input has been provided to Oregon State University in support of their hydraulic testing program. The HFIR conversion schedule was revised and provided to the GTRI program. In addition to HFIR conversion activities, technical support was provided directly to the Fuel Fabrication Capability program manager.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Complex Geometry Creation and Turbulent Conjugate Heat Transfer Modeling

The multiphysics capabilities of COMSOL provide the necessary tools to simulate the turbulent thermal-fluid aspects of the High Flux Isotope Reactor (HFIR). Version 4.1, and later, of COMSOL provides three different turbulence models: the standard k-{var_epsilon} closure model, the low Reynolds number (LRN) k-{var_epsilon} model, and the Spalart-Allmaras model. The LRN meets the needs of the nominal HFIR thermal-hydraulic requirements for 2D and 3D simulations. COMSOL also has the capability to create complex geometries. The circular involute fuel plates used in the HFIR require the use of algebraic equations to generate an accurate geometrical representation in the simulation environment. The best-estimate simulation results show that the maximum fuel plate clad surface temperatures are lower than those predicted by the legacy thermal safety code used at HFIR by approximately 17 K. The best-estimate temperature distribution determined by COMSOL was then used to determine the necessary increase in the magnitude of the power density profile (PDP) to produce a similar clad surface temperature as compared to the legacy thermal safety code. It was determined and verified that a 19% power increase was sufficient to bring the two temperature profiles to relatively good agreement.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Preliminary Multiphysics Analyses of HFIR LEU Fuel Conversion using COMSOL

The research documented herein was performed by several individuals across multiple organizations. We have previously acknowledged our funding for the project, but another common thread among the authors of this document, and hence the research performed, is the analysis tool COMSOL. The research has been divided into categories to allow the COMSOL analysis to be performed independently to the extent possible. As will be seen herein, the research has progressed to the point where it is expected that next year (2011) a large fraction of the research will require collaboration of our efforts as we progress almost exclusively into three-dimensional (3D) analysis. To the extent possible, we have tried to segregate the development effort into two-dimensional (2D) analysis in order to arrive at techniques and methodology that can be extended to 3D models in a timely manner. The Research Reactors Division (RRD) of ORNL has contracted with the University of Tennessee, Knoxville (UTK) Mechanical, Aerospace and Biomedical Engineering Department (MABE) to perform a significant fraction of this research. This group has been chosen due to their expertise and long-term commitment in using COMSOL and also because the participating students are able to work onsite on a part-time basis due to the close proximity of UTK with the ORNL campus. The UTK research has been governed by a statement of work (SOW) which clearly defines the specific tasks reported herein on the perspective areas of research. Ph.D. student Isaac T. Bodey has focused on heat transfer, fluid flow, modeling, and meshing issues and has been aided by his major professor Dr. Rao V. Arimilli and is the primary contributor to Section 2 of this report. Ph.D student Franklin G. Curtis has been focusing exclusively on fluid-structure interaction (FSI) due to the mechanical forces acting on the plate caused by the flow and has also been aided by his major professor Dr. Kivanc Ekici and is the primary contributor to Section 4 of this report. The HFIR LEU conversion project has also obtained the services of Dr. Prashant K. Jain of the Reactor & Nuclear Systems Division (RNSD) of ORNL. Prashant has quickly adapted to the COMSOL tools and has been focusing on thermal-structure interaction (TSI) issues and development of alternative 3D model approaches that could yield faster-running solutions. Prashant is the primary contributor to Section 5 of the report. And finally, while incorporating findings from all members of the "COMSOL team" (i.e., the team) and contributing as the senior COMSOL leader and advocate, Dr. James D. Freels has focused on the 3D model development, cluster deployment, and has contributed primarily to Section 3 and overall integration of this report. The team has migrated to the current release of COMSOL at version 4.1 for all the work described in this report, except where stated otherwise. Just as in the performance of the research, each of the respective sections has been originally authored by the respective authors. Therefore, the reader will observe a contrast in writing style throughout this document.

COMSOL↗

Exploiting New Features of COMSOL Version 4 on Conjugate Heat Transfer Problems

Users of COMSOL Multiphysics at version 3.5a and earlier have enjoyedmany features that have provided not only a good user experience atthe GUI interface, but also the capability to solve many classes ofproblems in a consistent manner with the physics being simulated. Withthe new release version 4.0 and later (4+) of COMSOL, the user isprovided a dramatic new interface from which to interact, and many newfeatures ``under the hood'' for solving problems more efficiently and with evengreater accuracy and consistency than before. This paper will exploreseveral of these new version 4+ features for the conjugate heattransfer class of problems. Our environment is challenging in that wedemand high-quality solutions for nuclear-reactor systems and themodels tend to become large and difficult to solve. Areasinvestigated include turbulence modeling, distributed parallelprocessing, solver scaling, and opengl graphics issues in a Linuxcomputing environment.

99 GENERAL AND MISCELLANEOUS↗

Upgrading the HFIR Thermal-Hydraulic Legacy Code Using COMSOL

Modernization of the High Flux Isotope Reactor (HFIR) thermal-hydraulic (TH) design and safety analysis capability is an important step in preparation for the conversion of the HFIR core from a high enriched uranium (HEU) fuel to a low enriched uranium (LEU) fuel. Currently, an important part of the HFIR TH analysis is based on the legacy Steady State Heat Transfer Code (SSHTC), which adds much conservatism to the safety analysis. The multi-dimensional multi-physics capabilities of the COMSOL environment allow the analyst to relax the number and magnitude of conservatisms, imposed by the SSHTC, to present a more physical model of the TH aspect of the HFIR.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

2D Thermal Hydraulic Analysis and Benchmark in Support of HFIR LEU Conversion using COMSOL

The research documented herein was funded by a research contract between the Research Reactors Division (RRD) of Oak Ridge National Laboratory (ORNL) and the University of Tennessee, Knoxville (UTK) Mechanical, Aerospace and Biomedical Engineering Department (MABE). The research was governed by a statement of work (SOW) which clearly defines nine specific tasks. This report is outlined to follow and document the results of each of these nine specific tasks. The primary goal of this phase of the research is to demonstrate, through verification and validation methods, that COMSOL is a viable simulation tool for thermal-hydraulic modeling of the High Flux Isotope Reactor (HFIR) core. A secondary goal of this two-dimensional phase of the research is to establish methodology and data base libraries that are also needed in the full three-dimensional COMSOL simulation to follow. COMSOL version 3.5a was used for all of the models presented throughout this report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗