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.

Efficient generation of one-group cross sections for coupled Monte Carlo depletion calculations

Coupled Monte Carlo depletion systems provide a versatile and an accurate tool for analyzing advanced thermal and fast reactor designs for a variety of fuel compositions and geometries. The main drawback of Monte Carlo-based systems is a long calculation time imposing significant restrictions on the complexity and amount of design-oriented calculations. This paper presents an alternative approach to interfacing the Monte Carlo and depletion modules aimed at addressing this problem. The main idea is to calculate the one-group cross sections for all relevant isotopes required by the depletion module in a separate module external to Monte Carlo calculations. Thus, the Monte Carlo module will produce the criticality and neutron spectrum only, without tallying of the individual isotope reaction rates. The onegroup cross section for all isotopes will be generated in a separate module by collapsing a universal multigroup (MG) cross-section library using the Monte Carlo calculated flux. Here, the term "universal" means that a single MG cross-section set will be applicable for all reactor systems and is independent of reactor characteristics such as a neutron spectrum; fuel composition; and fuel cell, assembly, and core geometries. This approach was originally proposed by Haeck et al. and implemented in the ALEPH code. Implementation of the proposed approach to Monte Carlo burnup interfacing was carried out through the BGCORE system. One-group cross sections generated by the BGCORE system were compared with those tallied directly by the MCNP code. Analysis of this comparison was carried out and led to the conclusion that in order to achieve the accuracy required for a reliable core and fuel cycle analysis, accounting for the background cross section (σ0) in the unresolved resonance energy region is essential. An extension of the one-group cross-section generation model was implemented and tested by tabulating and interpolating by a simplified σ0 model. A significant improvement of the one-group cross-section accuracy was demonstrated.

The stochastic implicit Euler method - A stable coupling scheme for Monte Carlo burnup calculations

Existing Monte Carlo burnup codes use various schemes to solve the coupled criticality and burnup equations. Previous studies have shown that the coupling schemes of the existing Monte Carlo burnup codes can be numerically unstable. Here we develop the Stochastic Implicit Euler method - a stable and efficient new coupling scheme. The implicit solution is obtained by the stochastic approximation at each time step. Our test calculations demonstrate that the Stochastic Implicit Euler method can provide an accurate solution to problems where the methods in the existing Monte Carlo burnup codes fail. © 2013 Elsevier Ltd. All rights reserved.

A global assessment of manufacturing: Economic development, energy use, carbon emissions, and the potential for energy efficiency and materials recycling

We present in two parts an assessment of global manufacturing. In the first part, we review economic development, pollution, and carbon emissions from a country perspective, tracking the rise of China and other developing countries. The results show not only a rise in the economic fortunes of the newly industrializing nations, but also a significant rise in global pollution, particularly air pollution and CO2 emissions largely from coal use, which alter and even reverse previous global trends. In the second part, we change perspective and quantitatively evaluate two important technical strategies to reduce pollution and carbon emissions: energy efficiency and materials recycling. We subdivide the manufacturing sector on the basis of the five major subsectors that dominate energy use and carbon emissions: (a) iron and steel, (b) cement, (c) plastics, (d) paper, and (e) aluminum. The analysis identifies technical constraints on these strategies, but by combined and aggressive action, industry should be able to balance increases in demand with these technical improvements. The result would be high but relatively flat energy use and carbon emissions. The review closes by demonstrating the consequences of extrapolating trends in production and carbon emissions and suggesting two options for further environmental improvements, materials efficiency, and demand reduction. © 2013 by Annual Reviews. All rights reserved.

The distribution and economic losses of alien species invasion to China

Invasive alien species have become one of the most serious environmental issues in the world. Data of taxon, origin, pathway, and environmental impacts of invasive alien microorganisms, invertebrates, amphibians and reptiles, fish, birds, mammals, herbs, trees, and, marine organisms in terrestrial, aquatic, and marine ecosystems of China were analyzed during 2001 and 2003, based on literature retrieval and field survey. There were 283 invasive alien species in China, and the number of species of invasive alien microorganisms, aquatic plants, terrestrial plants, aquatic invertebrates, terrestrial invertebrates, amphibians and reptiles, fish, and mammals were 19, 18, 170, 25, 33, 3, 10, and 5, respectively. The proportion of invasive alien species originated from America, Europe, Asia, Africa, and Oceania were 55.1, 21.7, 9.9, 8.1, and 0.6%, respectively. Methods for estimation of direct economic losses to agriculture, forestry, stockbreeding, fishery, road and water transportation, storage, water conservancy, environment and public facilities, and human health were established. Methods for estimation of indirect economic losses caused by invasive alien species to service functions of forest ecosystems, agricultural ecosystems, grassland ecosystems, and wetland ecosystems were also established. The total economic losses caused by invasive alien species to China were to the time of USD 14.45 billion, with direct and indirect economic losses accounting for 16.59% and 83.41% of total economic losses, respectively.

Regulation of gamete release in the economic brown seaweed Hizikia fusiforme (Phaeophyta)

Gamete release is an essential event in artificial seeding of the economic brown seaweed, Hizikia fusiforme. Mass egg release occurred in the dark, with few eggs being discharged in the light. Release of eggs was elicited with eight practical salinity units (one PSU = 1 g sea salts l(-1)) and was inhibited by salinity levels > 32 PSU. Egg release was optimal at 23 degrees C, and was decreased by 72% in agitated seawater compared to unstirred seawater. Inhibitors of photosynthesis and ions channels suppressed egg release, indicating that this process was physiologically associated with photosynthetic activity and ion transport.

Photosynthetic characteristics of the economic brown seaweed Hizikia fusiforme (Sargassaceae, Phaeophyta), with special reference to its "leaf" and receptacle

Photosynthetic responses to irradiance and temperature of "leaves" and receptacles were compared in February ( vegetative stage) and May ( reproductive stage) in the seaweed, Hizikia fusiforme ( Harvey) Okamura (Sargassaceae, Phaeophyta) from Nanao Island, Shantou, China. Irradiance-saturated photosynthesis (P-max) was significantly higher in receptacles than in "leaves" on a fresh weight basis, and that of "leaves" was greater in May than in February at ambient seawater temperatures. The optimum temperature for P-max was 30 degrees C for both "leaves" and receptacles, being 5 - 10 degrees C higher than the ambient seawater temperature. The apparent photosynthetic efficiencies were greater in receptacles than in "leaves" within the tested temperature range of 10 - 40 degrees C. The irradiance for saturating photosynthesis for both "leaves" and receptacles was temperature-dependent, with the highest values ( about 200 mu mol photons m(-2) s(-1)) at 30 degrees C.

Photosynthetic utilization of inorganic carbon in the economic brown alga, Hizikia Fusiforme (Sargassaceae) from the South China Sea

The mechanism of inorganic carbon (C-i) acquisition by the economic brown macroalga, Hizikia fusiforme (Harv.) Okamura (Sargassaceae), was investigated to characterize its photosynthetic physiology. Both intracellular and extracellular carbonic anhydrase (CA) were detected, with the external CA activity accounting for about 5% of the total. Hizikia fusiforme showed higher rates of photosynthetic oxygen evolution at alkaline pH than those theoretically derived from the rates of uncatalyzed CO2 production from bicarbonate and exhibited a high pH compensation point (pH 9.66). The external CA inhibitor, acetazolamide, significantly depressed the photosynthetic oxygen evolution, whereas the anion-exchanger inhibitor 4,4'-diisothiocyano-stilbene-2,2'-disulfonate had no inhibitory effect on it, implying the alga was capable of using HCO3- as a source of C-i for its photosynthesis via the mediation of the external CA. CO2 concentrations in the culture media affected its photosynthetic properties. A high level of CO2 (10,000 ppmv) resulted in a decrease in the external CA activity; however, a low CO2 level (20 ppmv) led to no changes in the external CA activity but raised the intracellular CA activity. Parallel to the reduction in the external CA activity at the high CO2 was a reduction in the photosynthetic CO2 affinity. Decreased activity of the external CA in the high CO2 grown samples led to reduced sensitiveness of photosynthesis to the addition of acetazolamide at alkaline pH. It was clearly indicated that H. fusiforme, which showed CO2-limited photosynthesis with the half-saturating concentration of C-i exceeding that of seawater, did not operate active HCO3- uptake but used it via the extracellular CA for its photosynthetic carbon fixation.

Defects in gallium nitride nanowires: First principles calculations

Atomic configurations and formation energies of native defects in an unsaturated GaN nanowire grown along the [001] direction and with (100) lateral facets are studied using large-scale ab initio calculation. Cation and anion vacancies, antisites, and interstitials in the neutral charge state are all considered. The configurations of these defects in the core region and outermost surface region of the nanowire are different. The atomic configurations of the defects in the core region are same as those in the bulk GaN, and the formation energy is large. The defects at the surface show different atomic configurations with low formation energy. Starting from a Ga vacancy at the edge of the side plane of the nanowire, a N-N split interstitial is formed after relaxation. As a N site is replaced by a Ga atom in the suboutermost layer, the Ga atom will be expelled out of the outermost layers and leaves a vacancy at the original N site. The Ga interstitial at the outmost surface will diffuse out by interstitialcy mechanism. For all the tested cases N-N split interstitials are easily formed with low formation energy in the nanowires, indicating N-2 molecular will appear in the GaN nanowire, which agrees well with experimental findings.

Magnetic coupling properties of mn-doped ZnO nanowires: First-principles calculations

Based on the density functional theory, we study the magnetic coupling properties of Mn-doped ZnO nanowires. For the nanowires with passivated surfaces, the antiferromagnetic state is found and the Mn atoms have a clustering tendency. When the distance between two Mn atoms is large, the system energetically favors the paramagnetic or spin-glass state. For the nanowires with unpassivated surfaces, the ferromagnetic (FM) coupling states appear between the two nearest Mn atoms, and the zinc vacancies can further stabilize the FM states between them. The electrons with enough concentration possibly mediate the FM coupling due to the negative exchange splitting of conduction band minimum induced by the s-d coupling, which could be useful in nanomaterial design for spintronics. (C) 2008 American Institute of Physics.