Some examples of multiple solutions in limiting flaw calculations using PD6493:1991 fracture mechanics analyses

The use of PC-based PD6493:1991 fracture assessment procedures has revealed that, under certain circumstances, flaws of different dimensions may be found as being limiting or critical for identical applied conditions. The main causes for multiple solutions are a steep applied stress gradient, residual stress relaxation and flaw re-characterisation. This work uses several case studies to illustrate some of the circumstances under which multiple solutions occurs.

CALCULATIONS OF INLET DISTORTION INDUCED COMPRESSOR FLOWFIELD INSTABILITY.

Calculations are presented predicting the onset of flow instability for a multistage low speed axial compressor operating in circumferentially distorted inlet flow. The most important feature of the model used is that it attempts to properly account for the fluid dynamic interaction between the spoiled and unspoiled sectors of the compressor. The calculations show that there is an approximate stability criterion, the annulus averaged slope of the compressor pressure rise characteristic equal to zero, that is valid whenever the dynamics of the compressor distorted flowfield can be considered independent of the compressor environment. This approximate criterion is used to investigate the relationship between the present model and the 'parallel compressor' model. Further calculations are performed to investigate cases of interest when the dynamics of the compressor flowfield are coupled to the environment. Resonant cases and cases when the distortion is unsteady are studied.

Discrete element calculations of the impact of a sand column against rigid structures

Discrete particle simulations of column of an aggregate of identical particles impacting a rigid, fixed target and a rigid, movable target are presented with the aim to understand the interaction of an aggregate of particles upon a structure. In most cases the column of particles is constrained against lateral expansion. The pressure exerted by the particles upon the fixed target (and the momentum transferred) is independent of the co-efficient of restitution and friction co-efficient between the particles but are strongly dependent upon the relative density of the particles in the column. There is a mild dependence on the contact stiffness between the particles which controls the elastic deformation of the densified aggregate of particles. In contrast, the momentum transfer to a movable target is strongly sensitive to the mass ratio of column to target. The impact event can be viewed as an inelastic collision between the sand column and the target with an effective co-efficient of restitution between 0 and 0.35 depending upon the relative density of the column. We present a foam analogy where impact of the aggregate of particles can be modelled by the impact of an equivalent foam projectile. The calculations on the equivalent projectile are significantly less intensive computationally and yet give predictions to within 5% of the full discrete particle calculations. They also suggest that "model" materials can be used to simulate the loading by an aggregate of particles within a laboratory setting. © 2012 Elsevier Ltd. All rights reserved.

Measurements and calculations of transport AC loss in second generation high temperature superconducting pancake coils

Theoretical and experimental AC loss data on a superconducting pancake coil wound using second generation (2 G) conductors are presented. An anisotropic critical state model is used to calculate critical current and the AC losses of a superconducting pancake coil. In the coil there are two regions, the critical state region and the subcritical region. The model assumes that in the subcritical region the flux lines are parallel to the tape wide face. AC losses of the superconducting pancake coil are calculated using this model. Both calorimetric and electrical techniques were used to measure AC losses in the coil. The calorimetric method is based on measuring the boil-off rate of liquid nitrogen. The electric method used a compensation circuit to eliminate the inductive component to measure the loss voltage of the coil. The experimental results are consistent with the theoretical calculations thus validating the anisotropic critical state model for loss estimations in the superconducting pancake coil. © 2011 American Institute of Physics.

Numerical stability of the predictor-corrector method in Monte Carlo burnup calculations of critical reactors

Monte Carlo burnup codes use various schemes to solve the coupled criticality and burnup equations. Previous studies have shown that the simplest methods, such as the beginning-of-step and middle-of-step constant flux approximations, are numerically unstable in fuel cycle calculations of critical reactors. Here we show that even the predictor-corrector methods that are implemented in established Monte Carlo burnup codes can be numerically unstable in cycle calculations of large systems. © 2013 Elsevier Ltd. All rights reserved.

Comparison among MCNP-based depletion codes applied to burnup calculations of pebble-bed HTR lattices

The double-heterogeneity characterising pebble-bed high temperature reactors (HTRs) makes Monte Carlo based calculation tools the most suitable for detailed core analyses. These codes can be successfully used to predict the isotopic evolution during irradiation of the fuel of this kind of cores. At the moment, there are many computational systems based on MCNP that are available for performing depletion calculation. All these systems use MCNP to supply problem dependent fluxes and/or microscopic cross sections to the depletion module. This latter then calculates the isotopic evolution of the fuel resolving Bateman's equations. In this paper, a comparative analysis of three different MCNP-based depletion codes is performed: Montburns2.0, MCNPX2.6.0 and BGCore. Monteburns code can be considered as the reference code for HTR calculations, since it has been already verified during HTR-N and HTR-N1 EU project. All calculations have been performed on a reference model representing an infinite lattice of thorium-plutonium fuelled pebbles. The evolution of k-inf as a function of burnup has been compared, as well as the inventory of the important actinides. The k-inf comparison among the codes shows a good agreement during the entire burnup history with the maximum difference lower than 1%. The actinide inventory prediction agrees well. However significant discrepancy in Am and Cm concentrations calculated by MCNPX as compared to those of Monteburns and BGCore has been observed. This is mainly due to different Am-241 (n,γ) branching ratio utilized by the codes. The important advantage of BGCore is its significantly lower execution time required to perform considered depletion calculations. While providing reasonably accurate results BGCore runs depletion problem about two times faster than Monteburns and two to five times faster than MCNPX. © 2009 Elsevier B.V. All rights reserved.