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.

Award Winner in the Young Investigator Category, 2014 Society for Biomaterials Annual Meeting and Exposition, Denver, Colorado, April 16-19, 2014: Periodically perforated core-shell collagen biomaterials balance cell infiltration, bioactivity, and mechanical properties.

Orthopedic tissue engineering requires biomaterials with robust mechanics as well as adequate porosity and permeability to support cell motility, proliferation, and new extracellular matrix (ECM) synthesis. While collagen-glycosaminoglycan (CG) scaffolds have been developed for a range of tissue engineering applications, they exhibit poor mechanical properties. Building on previous work in our lab that described composite CG biomaterials containing a porous scaffold core and nonporous CG membrane shell inspired by mechanically efficient core-shell composites in nature, this study explores an approach to improve cellular infiltration and metabolic health within these core-shell composites. We use indentation analyses to demonstrate that CG membranes, while less permeable than porous CG scaffolds, show similar permeability to dense materials such as small intestine submucosa (SIS). We also describe a simple method to fabricate CG membranes with organized arrays of microscale perforations. We demonstrate that perforated membranes support improved tenocyte migration into CG scaffolds, and that migration is enhanced by platelet-derived growth factor BB-mediated chemotaxis. CG core-shell composites fabricated with perforated membranes display scaffold-membrane integration with significantly improved tensile properties compared to scaffolds without membrane shells. Finally, we show that perforated membrane-scaffold composites support sustained tenocyte metabolic activity as well as improved cell infiltration and reduced expression of hypoxia-inducible factor 1α compared to composites with nonperforated membranes. These results will guide the design of improved biomaterials for tendon repair that are mechanically competent while also supporting infiltration of exogenous cells and other extrinsic mediators of wound healing.

Superposition codes for mismatched decoding

An achievable rate is given for discrete memoryless channels with a given (possibly suboptimal) decoding rule. The result is obtained using a refinement of the superposition coding ensemble. The rate is tight with respect to the ensemble average, and can be weakened to the LM rate of Hui and Csiszár-Körner, and to Lapidoth's rate based on parallel codebooks. © 2013 IEEE.

Second-order rate region of constant-composition codes for the multiple-access channel

This paper presents an achievable second-order rate region for the discrete memoryless multiple-access channel. The result is obtained using a random-coding ensemble in which each user's codebook contains codewords of a fixed composition. It is shown that this ensemble performs at least as well as i.i.d. random coding in terms of second-order asymptotics, and an example is given where a strict improvement is observed. © 2013 IEEE.

Monitoring and preventing numerical oscillations in 3D simulations with coupled Monte Carlo codes

Previous studies have reported that different schemes for coupling Monte Carlo (MC) neutron transport with burnup and thermal hydraulic feedbacks may potentially be numerically unstable. This issue can be resolved by application of implicit methods, such as the stochastic implicit mid-point (SIMP) methods. In order to assure numerical stability, the new methods do require additional computational effort. The instability issue however, is problem-dependent and does not necessarily occur in all cases. Therefore, blind application of the unconditionally stable coupling schemes, and thus incurring extra computational costs, may not always be necessary. In this paper, we attempt to develop an intelligent diagnostic mechanism, which will monitor numerical stability of the calculations and, if necessary, switch from simple and fast coupling scheme to more computationally expensive but unconditionally stable one. To illustrate this diagnostic mechanism, we performed a coupled burnup and TH analysis of a single BWR fuel assembly. The results indicate that the developed algorithm can be easily implemented in any MC based code for monitoring of numerical instabilities. The proposed monitoring method has negligible impact on the calculation time even for realistic 3D multi-region full core calculations. © 2014 Elsevier Ltd. All rights reserved.