HTGR fuel element depletion benchmark: Stage three results

Recently, a new numerical benchmark exercise for High Temperature Gas Cooled Reactor (HTGR) fuel depletion was defined. The purpose of this benchmark is to provide a comparison basis for different codes and methods applied to the burnup analysis of HTGRs. The benchmark specifications include three different models: (1) an infinite lattice of tristructural isotropic (TRISO) fuel particles, (2) an infinite lattice of fuel pebbles, and (3) a prismatic fuel including fuel and coolant channels. In this paper, we present the results of the third stage of the benchmark obtained with MCNP based depletion code BGCore and deterministic lattice code HELIOS 1.9. The depletion calculations were performed for three-dimensional model of prismatic fuel with explicitly described TRISO particles as well as for two-dimensional model, in which double heterogeneity of the TRISO particles was eliminated using reactivity equivalent physical transformation (RPT). Generally, good agreement in the results of the calculations obtained using different methods and codes was observed.

HTGR fuel element depletion benchmark: Stage one results

In this paper, we reported the results of the first stage of HTGR fuel element depletion benchmark obtained with BGCore and HELIOS depletion codes. The results of the k-inf are generally in good agreement. However, significant deviation in concentrations of several nuclides between MCNP based and HELIOS codes was observed.

The impact of sand slugs against beams and plates: Coupled discrete particle/finite element simulations

The impact of a slug of dry sand particles against a metallic sandwich beam or circular sandwich plate is analysed in order to aid the design of sandwich panels for shock mitigation. The sand particles interact via a combined linear-spring-and-dashpot law whereas the face sheets and compressible core of the sandwich beam and plate are treated as rate-sensitive, elastic-plastic solids. The majority of the calculations are performed in two dimensions and entail the transverse impact of end-clamped monolithic and sandwich beams, with plane strain conditions imposed. The sand slug is of rectangular shape and comprises a random loose packing of identical, circular cylindrical particles. These calculations reveal that loading due to the sand is primarily inertial in nature with negligible fluid-structure interaction: the momentum transmitted to the beam is approximately equal to that of the incoming sand slug. For a slug of given incoming momentum, the dynamic deflection of the beam increases with decreasing duration of sand-loading until the impulsive limit is attained. Sandwich beams with thick, strong cores significantly outperform monolithic beams of equal areal mass. This performance enhancement is traced to the "sandwich effect" whereby the sandwich beams have a higher bending strength than that of the monolithic beams. Three-dimensional (3D) calculations are also performed such that the sand slug has the shape of a circular cylindrical column of finite height, and contains spherical sand particles. The 3D slug impacts a circular monolithic plate or sandwich plate and we show that sandwich plates with thick strong cores again outperform monolithic plates of equal areal mass. Finally, we demonstrate that impact by sand particles is equivalent to impact by a crushable foam projectile. The calculations on the equivalent projectile are significantly less intensive computationally, yet give predictions to within 5% of the full discrete particle calculations for the monolithic and sandwich beams and plates. These foam projectile calculations suggest that metallic foam projectiles can be used to simulate the loading by sand particles within a laboratory setting. © 2013 Elsevier Ltd.

Study on small-strain behaviours of methane hydrate sandy sediments using discrete element method

Methane hydrate bearing soil has attracted increasing interest as a potential energy resource where methane gas can be extracted from dissociating hydrate-bearing sediments. Seismic testing techniques have been applied extensively and in various ways, to detect the presence of hydrates, due to the fact that hydrates increase the stiffness of hydrate-bearing sediments. With the recognition of the limitations of laboratory and field tests, wave propagation modelling using Discrete Element Method (DEM) was conducted in this study in order to provide some particle-scale insights on the hydrate-bearing sandy sediment models with pore-filling and cementation hydrate distributions. The relationship between shear wave velocity and hydrate saturation was established by both DEM simulations and analytical solutions. Obvious differences were observed in the dependence of wave velocity on hydrate saturation for these two cases. From the shear wave velocity measurement and particle-scale analysis, it was found that the small-strain mechanical properties of hydrate-bearing sandy sediments are governed by both the hydrate distribution patterns and hydrate saturation. © 2013 AIP Publishing LLC.

Lateral response of monopiles using centrifuge testing and finite element analysis

Offshore wind capacity is expected to grow exponentially over the next decade resulting in the production of a considerable amount of renewable energy. Monopiles are currently the most popular type of foundation for supporting offshore wind turbines in shallow to medium depth waters. In this paper, the load-deformation response of a 3.8 m diameter monopile installed in soft clays when subjected to axial and lateral loading is investigated using centrifuge testing and soil pore-fluid coupled three-dimensional finite element analysis. Monopile deformation is principally assessed in terms of its lateral displacements and bending moments. Its behaviour as a short rigid pile is discussed using concepts such as its rotation at mudline and the pile depth at which pivoting occurs. © 2014 Taylor & Francis Group.

Adapting data processing to compare model and experiment accurately: A discrete element model and magnetic resonance measurements of a 3d cylindrical fluidized bed

Discrete element modeling is being used increasingly to simulate flow in fluidized beds. These models require complex measurement techniques to provide validation for the approximations inherent in the model. This paper introduces the idea of modeling the experiment to ensure that the validation is accurate. Specifically, a 3D, cylindrical gas-fluidized bed was simulated using a discrete element model (DEM) for particle motion coupled with computational fluid dynamics (CFD) to describe the flow of gas. The results for time-averaged, axial velocity during bubbling fluidization were compared with those from magnetic resonance (MR) experiments made on the bed. The DEM-CFD data were postprocessed with various methods to produce time-averaged velocity maps for comparison with the MR results, including a method which closely matched the pulse sequence and data processing procedure used in the MR experiments. The DEM-CFD results processed with the MR-type time-averaging closely matched experimental MR results, validating the DEM-CFD model. Analysis of different averaging procedures confirmed that MR time-averages of dynamic systems correspond to particle-weighted averaging, rather than frame-weighted averaging, and also demonstrated that the use of Gaussian slices in MR imaging of dynamic systems is valid. © 2013 American Chemical Society.

Structural tunability and switchable exciton emission in inorganic-organic hybrids with mixed halides

Room-temperature tunable excitonic photoluminescence is demonstrated in alloy-tuned layered Inorganic-Organic (IO) hybrids, (C12H 25NH3)2PbI4(1-y)Br4y (y = 0 to 1). These perovskite IO hybrids adopt structures with alternating stacks of low-dimensional inorganic and organic layers, considered to be naturally self-assembled multiple quantum wells. These systems resemble stacked monolayer 2D semiconductors since no interlayer coupling exists. Thin films of IO hybrids exhibit sharp and strong photoluminescence (PL) at room-temperature due to stable excitons formed within the low-dimensional inorganic layers. Systematic variation in the observed exciton PL from 510 nm to 350 nm as the alloy composition is changed, is attributed to the structural readjustment of crystal packing upon increase of the Br content in the Pb-I inorganic network. The energy separation between exciton absorption and PL is attributed to the modified exciton density of states and diffusion of excitons from relatively higher energy states corresponding to bromine rich sites towards the lower energy iodine sites. Apart from compositional fluctuations, these excitons show remarkable reversible flips at temperature-induced phase transitions. All the results are successfully correlated with thermal and structural studies. Such structural engineering flexibility in these hybrids allows selective tuning of desirable exciton properties within suitable operating temperature ranges. Such wide-range PL tunability and reversible exciton switching in these novel IO hybrids paves the way to potential applications in new generation of optoelectronic devices. © 2013 AIP Publishing LLC.

Finite Element studies on indentation size effect using a higher order strain gradient theory

In this work, a Finite Element implementation of a higher order strain gradient theory (due to Fleck and Hutchinson, 2001) has been used within the framework of large deformation elasto-viscoplasticity to study the indentation of metals with indenters of various geometries. Of particular interest is the indentation size effect (ISE) commonly observed in experiments where the hardness of a range of materials is found to be significantly higher at small depths of indentation but reduce to a lower, constant value at larger depths. That the ISE can be explained by strain gradient plasticity is well known but this work aims to qualitatively compare a gamut of experimental observations on this effect with predictions from a higher order strain gradient theory. Results indicate that many of the experimental observations are qualitatively borne out by our simulations. However, areas exist where conflicting experimental results make assessment of numerical predictions difficult. © 2012 Elsevier Ltd. All rights reserved.

Finite element analysis of the ring rolling process with integrated closed-loop control

In FEA of ring rolling processes the tools' motions usually are defined prior to simulation. This procedure neglects the closed-loop control, which is used in industrial processes to control up to eight degrees of freedom (rotations, feed rates, guide rolls) in real time, taking into account the machine's performance limits as well as the process evolution. In order to close this gap in the new simulation approach all motions of the tools are controlled according to sensor values which are calculated within the FE simulation. This procedure leads to more realistic simulation results in comparison to the machine behaviour. © 2012 CIRP.

Efficient parametric uncertainty analysis within the hybrid Finite Element/Statistical Energy Analysis method

This paper is concerned with the development of efficient algorithms for propagating parametric uncertainty within the context of the hybrid Finite Element/Statistical Energy Analysis (FE/SEA) approach to the analysis of complex vibro-acoustic systems. This approach models the system as a combination of SEA subsystems and FE components; it is assumed that the FE components have fully deterministic properties, while the SEA subsystems have a high degree of randomness. The method has been recently generalised by allowing the FE components to possess parametric uncertainty, leading to two ensembles of uncertainty: a non-parametric one (SEA subsystems) and a parametric one (FE components). The SEA subsystems ensemble is dealt with analytically, while the effect of the additional FE components ensemble can be dealt with by Monte Carlo Simulations. However, this approach can be computationally intensive when applied to complex engineering systems having many uncertain parameters. Two different strategies are proposed: (i) the combination of the hybrid FE/SEA method with the First Order Reliability Method which allows the probability of the non-parametric ensemble average of a response variable exceeding a barrier to be calculated and (ii) the combination of the hybrid FE/SEA method with Laplace's method which allows the evaluation of the probability of a response variable exceeding a limit value. The proposed approaches are illustrated using two built-up plate systems with uncertain properties and the results are validated against direct integration, Monte Carlo simulations of the FE and of the hybrid FE/SEA models. © 2013 Elsevier Ltd.

A fixed-grid b-spline finite element technique for fluid-structure interaction

We present a fixed-grid finite element technique for fluid-structure interaction problems involving incompressible viscous flows and thin structures. The flow equations are discretised with isoparametric b-spline basis functions defined on a logically Cartesian grid. In addition, the previously proposed subdivision-stabilisation technique is used to ensure inf-sup stability. The beam equations are discretised with b-splines and the shell equations with subdivision basis functions, both leading to a rotation-free formulation. The interface conditions between the fluid and the structure are enforced with the Nitsche technique. The resulting coupled system of equations is solved with a Dirichlet-Robin partitioning scheme, and the fluid equations are solved with a pressure-correction method. Auxiliary techniques employed for improving numerical robustness include the level-set based implicit representation of the structure interface on the fluid grid, a cut-cell integration algorithm based on marching tetrahedra and the conservative data transfer between the fluid and structure discretisations. A number of verification and validation examples, primarily motivated by animal locomotion in air or water, demonstrate the robustness and efficiency of our approach. © 2013 John Wiley & Sons, Ltd.

How polluted is the Yangtze River? Water quality downstream from the Three Gorges Dam

The concentrations of major anions and cations, nitrogen and phosphorus, dissolved and particulate trace elements, and organic pollutants were determined for the middle and lower reaches of the Yangtze River (Changjiang) from below the Three Gorges Dam (TGD) to the mouth at Shanghai in November 2006. The concentration of dissolved inorganic phosphate (DIP) was constant at a low level of 6-8 mu gP/L, but the concentration of nitrate (NO3-) approximately doubled downstream and was closely correlated with K+. This translated to a daily load of well over 1000 It of dissolved inorganic nitrogen (DIN) at Datong. The average concentrations of dissolved Pb (0.078 +/- 0.023 mu g/L), Cd (0.024 +/- 0.009 mu g/L), Cr(0.57 +/- 0.09 mu g/L), Cu (1.9 +/- 0.7 mu g/L), and Ni (0.50 +/- 0.49 mu g/L) were comparable with those in other major world rivers, while As (3.3 +/- 1.3 mu g/L) and Zn (1.5 +/- 0.6 mu g/L) were higher by factors of 5.5 and 2.5, respectively. The trace element contents of suspended particles of As (31 +/- 28 mu g/g), Pb (83 +/- 34 mu g/g), and Ni (52 +/- 16 mu g/g) were close to maximum concentrations recommended for rivers by the European Community (EC). The average concentrations of Cd (2.6 +/- 1.6 mu g/g), Cr (185 +/- 102 mu g/g), Cu (115 +/- 106 mu g/g), and Zn (500 +/- 300 mu g/g) exceeded the EC standards by a factor of two, and Hg (4.4 +/- 4.7 mu g/g) by a factor of 4 to 5. Locally occurring peak concentrations exceed these values up to fourfold, among them the notorious elements As, Hg, and Tl. All dissolved and particulate trace element concentrations were higher than estimates made twenty years ago [Zhang, J., Geochemistry of trace metals from Chinese river/estuary systems: an overview. Estuar Coast Shelf Sci 1995; 41: 631-658.]. The enormous loads of anthropogenic pollutants disposed to the river were diluted by the large water discharge of the Yangtze even during the lowest flow resulting in the relatively low concentration levels of trace elements and organic pollutants observed. We estimated loads of e.g. As, Pb and Ni to the East China Sea to be about 4600 kg As d(-1), 3000 kg Pb d(-1), and 2000 kg Ni d(-1). About 6000 t d(-1) of dissolved organic carbon (DOC) was delivered into the sea at the time of our cruise. We tested for 236 organic pollutants, and only the most infamous were found to be barely above detection limits. We estimated that the load of chlorinated compounds, aromatic hydrocarbons, phenols, and PAHs were between 500 and 3500 kg d(-1). We also detected eight herbicides entering the estuary with loads of 5-350 kg d(-1). The pollutant load, even when at low concentrations, are considerable and pose an increasing threat to the health of the East China Sea ecosystem. (c) 2008 Elsevier B.V. All rights reserved.

An important pitfall of pseudo-static finite element analysis

Finite Element (FE) pseudo-static analysis can provide a good compromise between simplified methods of dynamic analysis and time domain analysis. The pseudo-static FE approach can accurately model the in situ, stresses prior to seismic loading (when it follows a static analysis simulating the construction sequence) is relatively simple and not as computationally expensive as the time domain approach. However this method should be used with caution as the results can be sensitive to the choice of the mesh dimensions. In this paper two simple examples of pseudo-static finite element analysis are examined parametrically, a homogeneous slope and a cantilever retaining wall, exploring the sensitivity of the pseudo-static analysis results on the adopted mesh size. The mesh dependence was found to be more pronounced for problems with high critical seismic coefficients values (e.g. gentle slopes or small walls), as in these cases a generalised layer failure mechanism is developed simultaneously with the slope or wall mechanism. In general the mesh width was found not to affect notably the predicted value of critical seismic coefficient but to have a major impact on the predicted movements. © 2012 Elsevier Ltd.

Novel fluid grid and voidage calculation techniques for a discrete element model of a 3D cylindrical fluidized bed

A discrete element model (DEM) combined with computational fluid dynamics (CFD) was developed to model particle and fluid behaviour in 3D cylindrical fluidized beds. Novel techniques were developed to (1) keep fluid cells, defined in cylindrical coordinates, at a constant volume in order to ensure the conditions for validity of the volume-averaged fluid equations were satisfied and (2) smoothly and accurately measure voidage in arbitrarily shaped fluid cells. The new technique for calculating voidage was more stable than traditional techniques, also examined in the paper, whilst remaining computationally-effective. The model was validated by quantitative comparison with experimental results from the magnetic resonance imaging of a fluidised bed analysed to give time-averaged particle velocities. Comparisons were also made between theoretical determinations of slug rise velocity in a tall bed. It was concluded that the DEM-CFD model is able to investigate aspects of the underlying physics of fluidisation not readily investigated by experiment. © 2014 The Authors.