Coupled neutronic thermo-hydraulic analysis of full PWR core with Monte-Carlo based BGCore system

BGCore reactor analysis system was recently developed at Ben-Gurion University for calculating in-core fuel composition and spent fuel emissions following discharge. It couples the Monte Carlo transport code MCNP with an independently developed burnup and decay module SARAF. Most of the existing MCNP based depletion codes (e.g. MOCUP, Monteburns, MCODE) tally directly the one-group fluxes and reaction rates in order to prepare one-group cross sections necessary for the fuel depletion analysis. BGCore, on the other hand, uses a multi-group (MG) approach for generation of one group cross-sections. This coupling approach significantly reduces the code execution time without compromising the accuracy of the results. Substantial reduction in the BGCore code execution time allows consideration of problems with much higher degree of complexity, such as introduction of thermal hydraulic (TH) feedback into the calculation scheme. Recently, a simplified TH feedback module, THERMO, was developed and integrated into the BGCore system. To demonstrate the capabilities of the upgraded BGCore system, a coupled neutronic TH analysis of a full PWR core was performed. The BGCore results were compared with those of the state of the art 3D deterministic nodal diffusion code DYN3D (Grundmann et al.; 2000). Very good agreement in major core operational parameters including k-eff eigenvalue, axial and radial power profiles, and temperature distributions between the BGCore and DYN3D results was observed. This agreement confirms the consistency of the implementation of the TH feedback module. Although the upgraded BGCore system is capable of performing both, depletion and TH analyses, the calculations in this study were performed for the beginning of cycle state with pre-generated fuel compositions. © 2011 Published by Elsevier B.V.

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.

Use of axially graded burnable boron for hot-spot temperature reduction in a pressurized water reactor core

Shortly after the loading of a pressurized water reactor (PWR) core, the axial power distribution in fresh fuel has already attained the characteristic, almost flat shape, typical of burned fuel. At beginning of cycle (BOC), however, the axial distribution is centrally peaked. In assemblies hosting uniform burnable boron rods, this BOC peaking is even more pronounced. A reduction in the axial peaking is today often achieved by shortening the burnable boron rods by some 30 cm at each edge. It is shown that a two-zone grading of the boron rod leads, in a representative PWR cycle, to a reduction of the hot-spot temperature of approximately 70 °C, compared with the nongraded case. However, with a proper three-zone grading of the boron rod, an additional 20 °C may be cut off the hot-spot temperature. Further, with a slightly skewed application of this three-zone grading, an additional 50 °C may be cut off. The representative PWR cycle studied was cycle 11 of the Indian Point 2 station, with a simplification in the number of fuel types and in the burnup distribution. The analysis was based on a complete three-dimensional burnup calculation. The code system was ELCOS, with BOXER as an assembly code for the generation of burnup-dependent cross sections and SILWER as a three-dimensional core code with thermal-hydraulic feedback.

Application of axial discontinuity factors for high conversion water reactors

The paper shows that generating cross sections using three-dimensional geometry and application of axial discontinuity factors are essential requirements for obtaining accurate prediction of criticality and zone average reaction rates in highly heterogeneous RBWR-type systems using computer codes based on diffusion theory approximation. The same methodology as presented here will be used to generate discontinuity factors for each axial interface between fuel assembly zones to ensure preservation of reaction rates in each zone and global multiplication factor. The use of discontinuity factors and three-dimensional cross sections may allow for a coarser energy group structure which is desirable to simplify and speed up transient calculations.

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.

Generation of SFR few-group constants using the Monte Carlo code serpent

In this study, the Serpent Monte Carlo code was used as a tool for preparation of homogenized few-group cross sections for the nodal diffusion analysis of Sodium cooled Fast Reactor (SFR) cores. Few-group constants for two reference SFR cores were generated by Serpent and then employed by nodal diffusion code DYN3D in 2D full core calculations. The DYN3D results were verified against the references full core Serpent Monte Carlo solutions. A good agreement between the reference Monte Carlo and nodal diffusion results was observed demonstrating the feasibility of using Serpent for generation of few-group constants for the deterministic SFR analysis.

Average corotation of line segments near a point and vortex identification

An easy-to-interpret kinematic quantity measuring the average corotation of material line segments near a point is introduced and applied to vortex identification. At a given point, the vector of average corotation of line segments is defined as the average of the instantaneous local rigid-body rotation over "all planar cross sections" passing through the examined point. The vortex-identification method based on average corotation is a one-parameter, region-type local method sensitive to the axial stretching rate as well as to the inner configuration of the velocity gradient tensor. The method is derived from a well-defined interpretation of the local flow kinematics to determine the "plane of swirling" and is also applicable to compressible and variable-density flows. Practical application to direct numerical simulation datasets includes a hairpin vortex of boundary-layer transition, the reconnection process of two Burgers vortices, a flow around an inclined flat plate, and a flow around a revolving insect wing. The results agree well with some popular local methods and perform better in regions of strong shearing. Copyright © 2013 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Multi-scale permeability of murine bone measured by nanoindentation

The determination of lacunar-canalicular permeability is essential to understand the mechano-transduction mechanism of bone. Murine models are widely used to investigate skeletal growth and regulation, but the value of lacunar-canalicular permeability is still unclear. To address this question, a poroelastic analysis based on nanoindentation data was used to calculate the lacunar-canalicular permeability of wild type C57BL/6 mice of 12 months. Cross-sections of three tibiae were indented using spherical fluid cell indenter tips of two sizes. Results suggest that the value of lacunar-canalicular intrinsic permeability of B6 female murine tibia is in the order of 10 -24 m2. The distribution of the values of intrinsic permeability suggests that with larger contact sizes, nanoindentation alone is capable of capturing the multi-scale permeability of bone. Multi-scale permeability of bone measured by nanoindentation will lead to a better understanding of the role of fluid flow in mechano-transduction. © 2013 American Society of Civil Engineers.

Hybrid core carbon fiber composite sandwich panels: Fabrication and mechanical response

Carbon fiber reinforced polymer (CFRP) composite sandwich panels with hybrid foam filled CFRP pyramidal lattice cores have been assembled from a carbon fiber braided net, 3D woven face sheets and various polymeric foams, and infused with an epoxy resin using a vacuum assisted resin transfer process. Sandwich panels with a fixed CFRP truss mass have been fabricated using a variety of closed cell polymer and syntactic foams, resulting in core densities ranging from 44-482kgm-3. The through thickness and in-plane shear modulus and strength of the cores increased with increasing foam density. The use of low compressive strength foams within the core was found to result in a significant reduction in the compressive strength contributed by the CFRP trusses. X-ray tomography led to the discovery that the trusses develop an elliptical cross-section shape during pressure assisted resin transfer. The ellipticity of the truss cross-sections increased, and the lattice contribution to the core strength decreased as the foam density was reduced. Micromechanical modeling was used to investigate the relationships between the mechanical properties and volume fractions of the core materials and truss topology of the hybrid core. The specific strength and moduli of the hybrid cores lay between those of the CFRP lattices and foams used to fabricate them. However, their volumetric and gravimetric energy absorptions significantly exceeded those of the materials from which they were fabricated. They compare favorably with other lightweight energy absorbing materials and structures. © 2013.

Increasing allowable deformation criteria through application of level II LTSM approach

In the central part of the Delft railway tunnel project, an underground railway station is being built at very close distance to the existing station building, which is still in operation. Although elaborate sensitivity analyses were made, some unforeseen deformations were encountered during the first phases of the execution process. Especially the installation of temporary sheet pile walls as well as the installation of a huge amount of grout anchor piles resulted in deformations exceeding the predicted final deformations as well as the boundary values defined by a level I limiting tensile strain method (LTSM) approach. In order to ensure the execution process, supplementary analyses were made to predict future deformations, and this for multiple cross sections. These deformations were implemented into a finite element model of the masonry of the building in order to define probable crack formation. This Level II LTSM approach made it possible to increase the initially foreseen deformation criteria and the continuation of the works. Design steps, design models and monitoring results will be explained within this paper.

Characteristics of oxide-free InGaAlAs layers grown by narrow stripe selective MOVPE

Narrow stripe selective MOVPE has been used to grow high quality oxide-free InGaAlAs layers on an InP substrate patterned with SiO2 masks at optimized growth conditions. Mirror-like surface morphologies and abrupt cross sections are obtained in all samples without spike growth at the mask edge. For the narrow stripe selectively grown InGaAlAs layers with a mesa width of about 1.2 mu m, a bandgap wavelength shift of 70 nm, a photoluminescence (PL) intensity of more than 80% and a PL full width at half maximum (FWHM) of less than 60 meV are obtained simultaneously with a small mask width variation from 0 to 40 mu m. The characteristics of the thickness enhancement ratio and the PL spectrum dependence on the mask width are presented and explained by considering both the migration effect from a masked region and the lateral vapour diffusion effect.

Fabrication of InGaAlAs MQW buried heterostructure lasers by narrow stripe selective MOVPE

Fabrication of InGaAlAs MQW buried heterostructure (BH) lasers by narrow stripe selective MOVPE is demonstrated in this paper. High quality InGaAlAs MQWs were first grown by narrow stripe selective MOVPE without any etching process and assessed by analysing the cross sections and PL spectrums of the InGaAlAs MQWs. Furthermore, BHs were fabricated for the InGaAlAs MQW lasers by a developed unselective regrowth method, instead of conventional selective regrowth. The InGaAlAs MQW BH lasers exhibit good device characteristics, with a high internal differential quantum efficiency of 85% and a low internal loss of 6.7 cm(-1). Meanwhile, narrow divergence angles of the far field pattern are obtained for the fabricated lasers.

Deep centers investigations of P-HEMT functional materials of ultra-high-speed microstructures grown by MBE

The deep centers of high electron mobility transistor (HEMT) and pseudomorphic-HEMT (P-HEMT) functional materials of ultra-high-speed microstructures grown by MBE are investigated using deep level transient spectroscopy (DLTS) technique. DLTS spectra demonstrate that midgap states, having larger concentrations and capture cross sections, are measured in n-AlGaAs layers of HEMT and P-HEMT structures. These states may correlate strongly with oxygen content of n-AlGaAs layer. At the same time, one can observe that the movement of DX center is related to silicon impurity that is induced by the strain in AlGaAs layer of the mismatched AlGaAs/InGaAs/GaAs system of P-HEMT structure. The experimental results also show that DLTS technique may be a tool of optimization design of the practical devices.

Spectroscopic properties of anisotropic absorption in a neodymium-doped YAlO3 laser crystal

Absorption spectra of YAlO3:Nd for the three crystallographic axes are investigated at room temperature, The spectral strengths indicate that the absorption coefficient of YAlO3:Nd is anisotropic. The anisotropy of the local electric field acting on the rare-earth ion in a laser crystal is considered, An extended Judd-Ofelt theory is applied to calculate the absorption cross sections and oscillator strengths of the electric-dipole transitions in the different principal directions. Three groups of the phenomenological parameters are derived from a least-squares-fitting procedure. We also analyze theoretically the anisotropy of the optical absorption of YAlO3:Nd crystal in detail. (C) 1997 American Institute of Physics.

Deep centers in AlGaAs/GaAs GRIN-SCH SQW laser structures grown by MBE and MOCVD

The deep centers in AlGaAs/GaAs graded index-separate confinement heterostructure single quantum well (GRIN-SCHSQW) laser structures grown by MBE and MOCVD have been investigated using deep level transient spectroscopy (DLTS) technique, The majority and minority carrier DLTS spectra show that the deep (hole and electron) traps (Hi and E3), having large capture cross sections and concentrations, are observed in the graded n-AlxGa1-xAs layer of laser structures in addition to the well-known DX centers. For laser structures grown by MBE, the deep hole trap H1 and the deep electron trap E3 may be spatially localized in the interface regions of discontinuous variation Al mole fraction of the n-AlxGa1-xAs layer with x = 0.20-0.43. For laser structures grown by MOCVD, the deep electron trap E3 may be spatially localized in the n-AlxGa1-xAs layer with x = 0.18-0.30, and the DX center may be spatially localized in the interface regions of discontinuous variation Al mole fraction of the AlxGa1-xAs layer with x = 0.22-0.30.