Fluid-Structure Interaction Modeling of High-Aspect Ratio Nuclear Fuel Plates Using COMSOL.


Book Description

The High Flux Isotope Reactor at the Oak Ridge National Lab is in the research stage of converting its fuel from high-enriched uranium to low-enriched uranium. Due to different physical properties of the new fuel and changes to the internal fuel plate design, the current safety basis must be re-evaluated through rigorous computational analyses. One of the areas being explored is the fluid-structure interaction phenomenon due to the interaction of thin fuel plates (50 mils thickness) and the cooling fluid (water). Detailed computational fluid dynamics and fluid-structure interaction simulations have only recently become feasible due to improved numerical algorithms and advancements in computing technology. For many reasons including the already built-in fluid-structure interaction module, COMSOL has been chosen for this complex problem. COMSOL's ability to solve multiphysics problems using a fully-coupled and implicit solution algorithm is crucial in obtaining a stable and accurate solution. Our initial findings show that COMSOL can accurately model such problems due to its ability to closely couple the fluid dynamics and the structural dynamics problems.




A Parametric Study of a Curved Nuclear Fuel Plate in a Narrow Channel Using Numeric Fluid Structure Interaction Modeling


Book Description

In 2004 the Global Threat Reduction Initiative (GTRI) was established by the National Nuclear Security Administration (NNSA) to quickly identify, secure, remove and/or closely monitor nuclear and radiological materials that pose a high-risk threat to the United States and the international community. Part of GTRI's mission is to convert high performance research reactors and isotope production facilities from their current High Enriched Uranium (HEU) fuel to Low Enriched Uranium (LEU) foil based fuel. In compliance with the conversion portion of the GTRI's mission, the University of Missouri Research Reactor (MURR) is currently trying to convert its reactor fuel. The proposed fuel uses a monolithic U-10Mo foil meat with a zirconium barrier between the aluminum cladding. This is different than the current HEU fuel meat which is comprised of Uranium dispersed in an aluminum matrix in an aluminum cladding. In addition to a change in the physical structure of the fuel, the fuel plate thickness has been significantly decreased. The fuel plates in the MURR reactor are subject to high velocity coolant (water) flow. A decreased thickness in the fuel plate suggests that the rigidity of the fuel plate will decrease as well. With concerns about the hydro-mechanical stability of the newly designed fuel plate being exposed to the high velocity flows in the reactor, there is a need to characterize the structural response of a very thin plate in presence of a velocity flow. Fluid structure interaction (FSI) simulations have been developed to analyze all of the characteristics of a thin fuel plate as the velocity of the water increases across the fuel plate. These models are developed by coupling CFD software, STAR CCM+, with Finite Element Analysis (FEA) software, ABAQUS, to determine the magnitude, location and direction of the deflection of the fuel plate along with other useful metrics to characterize attributes of the flow and movement of the plate. The FSI simulations were designed to replicate the plate thickness, fluid channel geometry and velocities of the experimental set up. The experimental set up consists of a flow loop and curved test section have been constructed for studying plate deflection and channel pressure drop under a variety of fluid flow velocities. The test section consists of two concentric steel cylinders bolted together with spacers between the two to form a fluid channel. The fuel plate inside the test section is clamped at the edges to maintain the axial location of the plate in the fluid channel. The aluminum fuel plate in the test section is 25.5 inches long, 16 mils (0.016 inches) thick and has an arc of 45 degrees. The outer cylinder of the test section has five plexi-glass windows that allow deflection data to be collected at various locations by the laser measurement system. The laser measurement system was fixed at the leading edge of the fuel plate because that is where the maximum deflection has been shown to occur. Since the focus of this study is on the maximum deflection of the fuel plate, the laser measurement system was fixed at that point. The FSI simulations were based on seven different geometries. The first consideration was an ideal geometry that the test section was designed to reflect. This consisted of a 130 mils thick inner fluid channel and a 78 mils thick outer fluid channel. Because the test section was assembled at Argonne National Laboratory (ANL) and then shipped to the University of Missouri, the dimensions of the test section fluid channel geometry was altered during transit. This change in the shape of the fluid channels lead to the need to characterize the shape of the fuel plate and fluid channels again. Since it is physically impossible to characterize the geometry of the entire fuel plate and the azimuthal geometry of the fluid channels of the test section without disassembling it, several geometric options were considered. Three geometries for the shape of the fuel plate were considered. The three fuel plate shapes were based on the distance of the plate based on the distance between the inner radius of the outer cylinder and the location of the fuel plate in the fluid channel. These measurements were obtained by measuring the outer channel gap thickness with a depth micrometer at the pressure tap locations that coincide with the location of the plate. Although this was helpful, it only helped to characterize the shape of the plate between the leading and trailing edges of the plate. This left about eight inches of the fuel plate (4 inches at each end) uncharacterized, leading to the three geometric considerations. The three geometric considerations were combined with two geometric considerations for the fluid channels. An azimuthally constant and varied geometry was considered. This was prompted by the inability to measure the distance of the fluid channels in the azimuthal direction due to the design of the test section. This lead to developing a geometry for the fluid channels that was constant along the azimuthal direction of the fuel plate and another where the fluid channel decreased in thickness as the azimuthal location progressed from the middle of the test section to the edges where the spacers are located. The FSI simulations and the experiments produced deflection results for the 16 mils (0.016 inches) thick fuel plate at velocities ranging from approximately 2 m/s to 4 m/s. Because the fuel plate is so thin, it has been decided that the velocity in the FSI simulations and experiments should be kept relatively low to avoid permanently damaging the aluminum plate in the test section. The As-Built numeric models have been shown to compare well with the experimental results. The As-Built numeric model results also show that the constant azimuthal fluid channel geometry produces a smaller magnitude of deflection for all three fuel plate shapes considered when compared to the models with a varying azimuthal geometry. the FSI simulations "crashed" for all models ran with a velocity beyond 4 m/s. The experimental results showed that maximum deflection results can be obtained beyond 4 m/s. This provides the opportunity to investigate the differences between the experimental and modeling geometries and recalibrate the model to obtain maximum deflection results through the FSI simulations.




An Approach to Modeling Fuel Plate Deformations in Fluid-Structure Interactions


Book Description

As computational power increases, so does the desire to use fluid-structure interaction (FSI) software to design complex systems and components such as fuel plates. Presently, an effort is underway to support the design and qualification of a new nuclear fuel material which is intended for applications in select, plate-type, research and test reactors. Because of the high flow velocities which are experienced within these reactors, a safety parameter of interest is the onset of plate deformation and potential of plate-to-plate collapse. Using FSI software for the entire design process is tempting; however, utilizing FSI software for large, complex systems requires large quantities of computational resources. Currently, computing FSI solutions using off-the-shelf software of components as small as individual fuel plates can take weeks on a desktop computer, thus requiring the use of multiple servers or a cluster to enable a pragmatic solving time. Since computational resources are valuable, the pertinent question to ask is whether or not the resources being used are a necessity. Fluid structure interaction simulations provide a wealth of information regarding flow patterns, but many practical engineering problems do not require such a detailed solution form. If acceptable solutions could be obtained without solving the entire flow field, the required computational resources may be reduced by orders of magnitude. The study detailed herein presents an alternative approach to solving FSI problems using a one dimensional, semi-analytic model derived from first principles. The resulting approach is much less computationally expensive than the alternative of leveraging off-the-shelf FSI software. The FSI model developed herein is specific to solving plate deformations (single plates and arrays of plates). A qualitative benchmark of this study's model is made against an analytic solution ABAQUS. Results are then compared against experimental data collected for a single plate acquired by the University of Missouri, as well as experimental data collected for an array of six plates acquired by at Oregon State University. The outcome of this work has resulted in a new modeling approach for FSI problems applied toward flat plate geometry. Results from this new model approach (while limited) compare reasonably well to commercially available software (ABAQUS) and experimental data over a wide range of flow conditions.




Fluid Structure Interaction of Involute Fuel Plates in the High Flux Isotope Reactor Using a Fully-coupled Numerical Approach


Book Description

This dissertation describes a fully-coupled (FC), finite-element (FE) based, algorithm for modeling and simulation of the fluid-structure interaction (FSI) of involuteshaped fuel plates used in research reactors; specifically the High Flux Isotope Reactor (HFIR) at the Oak Ridge National Laboratory (ORNL). Following a graded approach to code and model validation, a cylinder in cross-flow benchmark is used to establish flow physics as well as properly coupling the FSI phenomena with increasing complexity. As an interim step toward HFIR LEU fuel plate simulations, three experiments are used for validation. The first, performed by Smissaert, is used to envelope large plate deflections and understand the validity of various fluid boundary conditions for single plate comparisons. Continuing with Smissaert's data, a 5-plate simulation is presented showing the first-ever multi-plate simulation using this FC and FE approach. Second, a vibrating plate, presented by Liu et al., is simulated showing the same technique to encompass self-excited, periodic plate deflections. Lastly, an experiment for the conceptual Advanced Neutron Source Reactor (ANSR) using involute plates is utilized to validate the ability of this FC and FE algorithm to predict the deflections of the involute-shaped plates used in the HFIR. The method shown herein accurately captures the established `S-shaped' deflection of the first mode of the involute plate providing guidance that researchers and designers can utilize in the forthcoming design of the next generation of low-enriched uranium (LEU) fuel plates for the HFIR. A `Lessons Learned' section which describes external routine coupling, geometry and meshing guidance, and solver settings used in the computational platform used to perform these FSI simulations is also provided.




Development and Characterization of a Plate Fuel Hotspot Model in COMSOL


Book Description

Aluminum clad plate fuel is common to many high performance water cooled research reactors, including the High Flux Isotope Reactor (HFIR) at ORNL. Fuel manufacturing defects associated with fuel segregation and incomplete bonding of the cladding to the fuel material currently limit the performance of HFIR. A high resolution multi-physics (HRMP) simulation of concurrent fuel segregation and incomplete bonding of fuel cladding is developed in this dissertation. The simulation development begins with a review of legacy modeling of the fuel segregation and cladding non-bond, and then proceeds to identify improvements possible in the HRMP framework. A contact conductance model is selected for the incomplete bonding of cladding to fuel, advancing previous models. A verification of the COMSOL simulation platform used to construct the evaluation model is performed using the method of manufactured solutions, including assessments of the solid conduction modeling domain, the fluid coolant channel domain, and coupled fuel to coolant channel domains. Solution verification is performed with the least squares, grid convergence index approach, and indirect validation of the evaluation model is performed using data generated to establish thermal performance limits for the Advanced Neutron Source reactor. The verification and validation efforts are also extensions to previous work using COMSOL for HRMP modeling of HFIR, and establish numerical and modeling uncertainties. The evaluation model is then employed in a sensitivity assessment of 18 parameters in the evaluation model using Latin Hypercube sampling methods to establish a ranking of parameter importance in predicting four quantities of interest.




Development and Experimental Benchmarking of Numeric Fluid Structure Interaction Models for Research Reactor Fuel Analysis


Book Description

As part of the Global Threat Reduction Initiative (GTRI) reactor conversion program, five U.S. High Performance Research Reactors (HPRRs) are currently studying a novel Low Enriched Uranium (LEU) foil based fuel to replace their current High Enriched Uranium (HEU) dispersion fuel. The proposed changes in the plate design have led to a need to characterize the structural response of the plates in presence of high velocity flow. To generate experiment data, a flow loop and test section have been constructed for studying plate deflection and channel pressure drop under a variety of fluid flow conditions. The work presented here differs from earlier efforts by intentionally offsetting the plate and creating fluid channels of different thickness. This offset effectively simulates manufacturing tolerances of a real fuel assembly. A method for generating ‘As-Built’ numeric models of the experiment geometry is presented. These As-Built numeric models have been shown to dramatically improve matching between experiment and numeric solutions. At higher flow rates, the experiment exhibited a dynamic ‘snap’ behavior that could not be replicated numerically. Additional interrogation of the boundary conditions revealed a possible explanation for this snap, however numeric methods do not yet exist for recreating this behavior. In the numeric and experimental work presented here, plate deflection behavior at low to mid-range flow rates is qualitatively consistent with theoretical expectations.




Computational Mechanics of Fluid-Structure Interaction


Book Description

This book is intended to provide a compilation of the state-of-the-art numerical methods for nonlinear fluid-structure interaction using the moving boundary Lagrangian-Eulerian formulation. Single and two-phase viscous incompressible fluid flows are considered with the increasing complexity of structures ranging from rigid-body, linear elastic and nonlinear large deformation to fully-coupled flexible multibody system. This book is unique with regard to computational modeling of such complex fluid-structure interaction problems at high Reynolds numbers, whereby various coupling techniques are introduced and systematically discussed. The techniques are demonstrated for large-scale practical problems in aerospace and marine/offshore engineering. This book also provides a comprehensive understanding of underlying unsteady physics and coupled mechanical aspects of the fluid-structure interaction from a computational point of view. Using the body-fitted and moving mesh formulations, the physical insights associated with structure-to-fluid mass ratios (i.e., added mass effects), Reynolds number, large structural deformation, free surface, and other interacting physical fields are covered. The book includes the basic tools necessary to build the concepts required for modeling such coupled fluid-structure interaction problems, thus exposing the reader to advanced topics of multiphysics and multiscale phenomena.







Verification and Validation in Scientific Computing


Book Description

Advances in scientific computing have made modelling and simulation an important part of the decision-making process in engineering, science, and public policy. This book provides a comprehensive and systematic development of the basic concepts, principles, and procedures for verification and validation of models and simulations. The emphasis is placed on models that are described by partial differential and integral equations and the simulations that result from their numerical solution. The methods described can be applied to a wide range of technical fields, from the physical sciences, engineering and technology and industry, through to environmental regulations and safety, product and plant safety, financial investing, and governmental regulations. This book will be genuinely welcomed by researchers, practitioners, and decision makers in a broad range of fields, who seek to improve the credibility and reliability of simulation results. It will also be appropriate either for university courses or for independent study.




Modelling Organs, Tissues, Cells and Devices


Book Description

This book presents a theoretical and practical overview of computational modeling in bioengineering, focusing on a range of applications including electrical stimulation of neural and cardiac tissue, implantable drug delivery, cancer therapy, biomechanics, cardiovascular dynamics, as well as fluid-structure interaction for modelling of organs, tissues, cells and devices. It covers the basic principles of modeling and simulation with ordinary and partial differential equations using MATLAB and COMSOL Multiphysics numerical software. The target audience primarily comprises postgraduate students and researchers, but the book may also be beneficial for practitioners in the medical device industry.