Core design options for high conversion BWRs operating in Th- 233U fuel cycle

Several options of fuel assembly design are investigated for a BWR core operating in a closed self-sustainable Th-233U fuel cycle. The designs rely on an axially heterogeneous fuel assembly structure consisting of a single axial fissile zone "sandwiched" between two fertile blanket zones, in order to improve fertile to fissile conversion ratio. The main objective of the study was to identify the most promising assembly design parameters, dimensions of fissile and fertile zones, for achieving net breeding of 233U. The design challenge, in this respect, is that the fuel breeding potential is at odds with axial power peaking and the core minimum critical power ratio (CPR), hence limiting the maximum achievable core power rating. Calculations were performed with the BGCore system, which consists of the MCNP code coupled with fuel depletion and thermo-hydraulic feedback modules. A single 3-dimensional fuel assembly having reflective radial boundaries was modeled applying simplified restrictions on the maximum centerline fuel temperature and the CPR. It was found that axially heterogeneous fuel assembly design with a single fissile zone can potentially achieve net breeding, while matching conventional BWR core power rating under certain restrictions to the core loading pattern design. © 2013 Elsevier B.V. All rights reserved.

Neutronic optimization in high conversion Th-233U fuel assembly with simulated annealing

This paper reports on fuel design optimization of a PWR operating in a self sustainable Th-233U fuel cycle. Monte Carlo simulated annealing method was used in order to identify the fuel assembly configuration with the most attractive breeding performance. In previous studies, it was shown that breeding may be achieved by employing heterogeneous Seed-Blanket fuel geometry. The arrangement of seed and blanket pins within the assemblies may be determined by varying the designed parameters based on basic reactor physics phenomena which affect breeding. However, the amount of free parameters may still prove to be prohibitively large in order to systematically explore the design space for optimal solution. Therefore, the Monte Carlo annealing algorithm for neutronic optimization is applied in order to identify the most favorable design. The objective of simulated annealing optimization is to find a set of design parameters, which maximizes some given performance function (such as relative period of net breeding) under specified constraints (such as fuel cycle length). The first objective of the study was to demonstrate that the simulated annealing optimization algorithm will lead to the same fuel pins arrangement as was obtained in the previous studies which used only basic physics phenomena as guidance for optimization. In the second part of this work, the simulated annealing method was used to optimize fuel pins arrangement in much larger fuel assembly, where the basic physics intuition does not yield clearly optimal configuration. The simulated annealing method was found to be very efficient in selecting the optimal design in both cases. In the future, this method will be used for optimization of fuel assembly design with larger number of free parameters in order to determine the most favorable trade-off between the breeding performance and core average power density.

Comparison of square and hexagonal fuel lattices for high conversion PWRs

This paper reports on an investigation into fuel design choices of a pressurized water reactor operating in a self-sustainable Th- 233U fuel cycle. In order to evaluate feasibility of this concept, two types of fuel assembly lattices were considered: square and hexagonal. The hexagonal lattice may offer some advantages over the square one. For example, the fertile blanket fuel can be packed more tightly reducing the blanket volume fraction in the core and potentially allowing to achieve higher core average power density. The calculations were carried out with Monte-Carlo based BGCore code system and the results were compared to those obtained with Serpent Monte-Carlo code and deterministic transport code BOXER. One of the major design challenges associated with the SB concept is high power peaking due to the high concentration of fissile material in the seed region. The second objective of this work is to estimate the maximum achievable core power density by evaluation of limiting thermal hydraulic parameters. The analysis showed that both fuel assembly designs have a potential of achieving net breeding. Although hexagonal lattice was found to be somewhat more favorable because it allows achieving higher power density, while having breeding performance comparable to the square lattice case. © Carl Hanser Verlag München.

Investigation of fuel assembly design options for high conversion thorium fuel cycle in PWRs

This study explores the basic possibility of achieving a self-sustainable Th-U233 fuel cycle that can be adopted in the current generation of Pressurized Water Reactors. This study outlines some fuel design strategies to achieve (or to approach as closely as possible) a sustainable fuel cycle. Major design tradeoffs in the core design are discussed. Preliminary neutronic analysis performed on the fuel assembly level with BOXER computer code suggests that net breeding of U233 is feasible in principle within a typical PWR operating envelope. However, some reduction in the core power density and/or shorter than typical fuel cycle length would most likely be required in order to achieve such performance.

Axial discontinuity factors for the nodal diffusion analysis of high conversion BWR cores

High conversion LWRs concepts typically rely on a heterogeneous core configuration, where fissile zones are interspersed with fertile blanket zones in order to achieve a high conversion ratio. Modeling such a heterogeneous structure of these cores represents a significant challenge to the conventional reactor analysis methods. It was recently suggested to overcome such difficulties, in particular, for the case of axially heterogeneous reduced moderation BWRs, by introducing an additional set of discontinuity factors in axial direction at the interfaces between fissile and fertile fuel assembly zones. However, none of the existing nodal diffusion core simulators have the capability of accounting for discontinuity of homogeneous nodal fluxes in axial direction since the fuel composition of conventional LWRs is much more axially uniform. In this work, we modified the nodal diffusion code DYN3D by introducing such a capability. The new version of the code was tested on a series of reduced moderation BWR cases with Th-U233 and U-Pu-MA fuel. The library of few-group homogenized cross sections and the data required for the calculation of discontinuity factors were generated using the Monte Carlo transport code Serpent. The results obtained with the modified version of DYN3D were compared with the reference Monte Carlo solutions and were found to be in good agreement. The current analysis demonstrates that high conversion LWRs can in principle be modeled using existing nodal diffusion core simulators. © 2013 Elsevier Ltd. All rights reserved.

Low-order modelling of the response of ducted flames in annular geometries

The adoption of lean premixed prevaporised combustion systems can reduce NOx emissions from gas turbines, but unfortunately also increases their susceptibility to thermoacoustic instabilities. Initially, acoustic waves can produce heat release fluctuations by a variety of mechanisms, often by perturbing the equivalence ratio. If correctly phased, heat release fluctuations can subsequently generate more acoustic waves, which at high amplitude can result in significant structural damage to the combustor. The prediction of this phenomenon is of great industrial interest. In previous work, we have coupled a physics based, kinematic model of the flame with a network model to provide the planar acoustic response necessary to close the feedback loop and predict the onset and amplitude of thermoacoustic instabilities in a lab-scale, axisymmetric single burner combustor. The advantage of a time domain approach is that the modal interaction, the influence of harmonics, and flame saturation can be investigated. This paper extends this approach to more realistic, annular geometries, where both planar and circumferential modes must be considered. In lean premixed prevaporised combustors, fluctuations in equivalence ratio have been shown to be a dominant cause of unsteady combustion. These can occur, for example, due to velocity perturbations in the premix ducts, which can lead to equivalence ratio fluctuations at the fuel injectors, which are subsequently convected downstream to the flame surfaces. Here, they can perturb the heat release by locally altering the flame speed, enthalpy of combustion, and, indirectly, the flame surface area. In many gas turbine designs, particularly aeroengines, the geometries are composed of a ring of premix ducts linking a plenum and an annular combustor. The most unstable modes are often circumferential modes. The network model is used to characterise the flow response of the geometry to heat fluctuations at an appropriate location, such as the fuel injectors. The heat release at each flame holder is determined in the time domain using the kinematic flame model derived, as a function of the flow perturbations in the premix duct. This approach is demonstrated for an annular ring of burners on a in a simple geometry. The approach is then extended to an industrial type gas turbine combustor, and used to predict the limit cycle amplitudes. Copyright © 2012 by ASME.

On the effects of turbulence modelling in design optimisation for high-lift devices

Aerodynamic shape optimisation is being increasingly utilised as a design tool in the aerospace industry. In order to provide accurate results, design optimisation methods rely on the accuracy of the underlying CFD methods applied to obtain aerodynamic forces for a given configuration. Previous studies of the authors have highlighted that the variation of the order of accuracy of the CFD solver with a fixed turbulence model affects the resulting optimised airfoil shape for a single element airfoil. The accuracy of the underlying CFD model is even more relevant in the context of high-lift configurations where an accurate prediction of flow is challenging due to the complex flow physics involving transition and flow separation phenomena. This paper explores the effect of the fidelity of CFD results for a range of turbulence models within the context of the computational design of aircraft configurations. The NLR7301 multi-element airfoil (main wing and flap) is selected as the baseline configuration, because of the wealth of experimental an computational results available for this configuration. An initial validation study is conducted in order to establish optimal mesh parameters. A bi-objective shape optimisation problem is then formulated, by trying to reveal the trade-off between lift and drag coefficients at high angles of attack. Optimisation of the airfoil shape is performed with Spalart-Allmaras, k - ω SST and k - o realisable models. The results indicate that there is consistent and complementary impact to the optimum level achieved from all the three different turbulence models considered in the presented case study. Without identifying particular superiority of any of the turbu- lence models, we can say though that each of them expressed favourable influence towards different optimality routes. These observations lead to the exploration of new avenues for future research. © 2012 AIAA.

On the effects of turbulence modelling in design optimisation for high-lift devices

Aerodynamic shape optimisation is being increasingly utilised as a design tool in the aerospace industry. In order to provide accurate results, design optimisation methods rely on the accuracy of the underlying CFD methods applied to obtain aerodynamic forces for a given configuration. Previous studies of the authors have highlighted that the variation of the order of accuracy of the CFD solver with a fixed turbulence model affects the resulting optimised airfoil shape for a single element airfoil. The accuracy of the underlying CFD model is even more relevant in the context of high-lift configurations where an accurate prediction of flow is challenging due to the complex flow physics involving transition and flow separation phenomena. This paper explores the effect of the fidelity of CFD results for a range of turbulence models within the context of the computational design of aircraft configurations. The NLR7301 multi-element airfoil (main wing and flap) is selected as the baseline configuration, because of the wealth of experimental an computational results available for this configuration. An initial validation study is conducted in order to establish optimal mesh parameters. A bi-objective shape optimisation problem is then formulated, by trying to reveal the trade-off between lift and drag coefficients at high angles of attack. Optimisation of the airfoil shape is performed with Spalart-Allmaras, k - ω SST and k - ε realisable models. The results indicate that there is consistent and complementary impact to the optimum level achieved from all the three different turbulence models considered in the presented case study. Without identifying particular superiority of any of the turbu- lence models, we can say though that each of them expressed favourable influence towards different optimality routes. These observations lead to the exploration of new avenues for future research. © 2012 by the authors.

Collapse of a composite beam made from ultra high molecular-weight polyethylene fibres

Hot-pressed laminates with a [0/90]48 lay-up, consisting of 83% by volume of ultra high molecular-weight polyethylene (UHMWPE) fibres, and 17% by volume of polyurethane (PU) matrix, were cut into cantilever beams and subjected to transverse end-loading. The collapse mechanisms were observed both visually and by X-ray scans. Short beams deform elastically and collapse plastically in longitudinal shear, with a shear strength comparable to that observed in double notch, interlaminar shear tests. In contrast, long cantilever beams deform in bending and collapse via a plastic hinge at the built-in end of the beam. The plastic hinge is formed by two wedge-shaped microbuckle zones that grow in size and in intensity with increasing hinge rotation. This new mode of microbuckling under macroscopic bending involves both elastic bending and shearing of the plies, and plastic shear of the interface between each ply. The double-wedge pattern contrasts with the more usual parallel-sided plastic microbuckle that occurs in uniaxial compression. Finite element simulations and analytical models give additional insight into the dominant material and geometric parameters that dictate the collapse response of the UHMWPE composite beam in bending. Detailed comparisons between the observed and predicted collapse responses are used in order to construct a constitutive model for laminated UHMWPE composites. © 2013 Elsevier Ltd.

High-level scene structure using visibility and occlusion

We demonstrate a new method for extracting high-level scene information from the type of data available from simultaneous localisation and mapping systems. We model the scene with a collection of primitives (such as bounded planes), and make explicit use of both visible and occluded points in order to refine the model. Since our formulation allows for different kinds of primitives and an arbitrary number of each, we use Bayesian model evidence to compare very different models on an even footing. Additionally, by making use of Bayesian techniques we can also avoid explicitly finding the optimal assignment of map landmarks to primitives. The results show that explicit reasoning about occlusion improves model accuracy and yields models which are suitable for aiding data association. © 2011. The copyright of this document resides with its authors.