Simulating early adoption of alternative fuel vehicles for sustainability

We quantify the conditions that might trigger wide spread adoption of alternative fuel vehicles (AFVs) to support energy policy. Empirical review shows that early adopters are heterogeneous motivated by financial benefits, environmental appeal, new technology, and vehicle reliability. A probabilistic Monte Carlo simulation model is used to assess consumer heterogeneity for early and mass market adopters. For early adopters full battery electric vehicles (BEVs) are competitive but unable to surpass diesels or hybrids due to purchase price premium and lack of charging availability. For mass adoption, simulations indicate that if the purchase price premium of a BEV closes to within 20% of an in-class internal combustion engine (ICE) vehicle, combined with a 60% increase in refuelling availability relative to the incumbent system, BEVs become competitive. But this depends on a mass market that values the fuel economy and CO2 reduction benefits associated with BEVs. We also find that the largest influence on early adoption is financial benefit rather than pro-environmental behaviour suggesting that AFVs should be marketed by appealing to economic benefits combined with pro-environmental behaviour to motivate adoption. Monte Carlo simulations combined with scenarios can give insight into diffusion dynamics for other energy demand-side technologies. © 2012 Elsevier Inc.

Effect of operating conditions and fuel type on crevice hc emissions: Model results and comparison with experiments

A one-dimensional model for crevice HC post-flame oxidation is used to calculate and understand the effect of operating parameters and fuel type (propane and isooctane) on the extent of crevice hydrocarbon and the product distribution in the post flame environment. The calculations show that the main parameters controlling oxidation are: bulk burned gas temperatures, wall temperatures, turbulent diffusivity, and fuel oxidation rates. Calculated extents of oxidation agree well with experimental values, and the sensitivities to operating conditions (wall temperatures, equivalence ratio, fuel type) are reasonably well captured. Whereas the bulk gas temperatures largely determine the extent of oxidation, the hydrocarbon product distribution is not very much affected by the burned gas temperatures, but mostly by diffusion rates. Uncertainties in both turbulent diffusion rates as well as in mechanisms are an important factor limiting the predictive capabilities of the model. However, it seems well suited to sensitivity calculations about a baseline. Copyright © 1999 Society of Automotive Engineers, Inc.

Liquid fuel visualization using laser-induced fluoresence during cold start

The presence of liquid fuel inside the engine cylinder is believed to be a strong contributor to the high levels of hydrocarbon emissions from spark ignition (SI) engines during the warm-up period. Quantifying and determining the fate of the liquid fuel that enters the cylinder is the first step in understanding the process of emissions formation. This work uses planar laser induced fluorescence (PLIF) to visualize the liquid fuel present in the cylinder. The fluorescing compounds in indolene, and mixtures of iso-octane with dopants of different boiling points (acetone and 3-pentanone) were used to trace the behavior of different volatility components. Images were taken of three different planes through the engine intersecting the intake valve region. A closed valve fuel injection strategy was used, as this is the strategy most commonly used in practice. Background subtraction and masking were both performed to reduce the effect of any spurious fluorescence. The images were analyzed on both a time and crank angle (CA) basis, showing the time of maximum liquid fuel present in the cylinder and the effect of engine events on the inflow of liquid fuel. The results show details of the liquid fuel distribution as it enters the engine as a function of crankangle degree, volatility and location in the cylinder. A. semi-quantitative analysis based on the integration of the image intensities provides additional information on the temporal distribution of the liquid fuel flow. © 1998 Society of Automotive Engineers, Inc.

Novel experiment on in-cylinder desorption of fuel from the oil layer

A technique has been developed to measure the desorption and subsequent oxidation of fuel in the oil layer by spiking the oil with liquid fuel and firing the engine on gaseous fuel or motoring with air. Experiments suggest that fuel desorption is not diffusion limited above 50°C and indicated that approximately two to four percent of the cylinder oil layer is fresh oil from the sump. The increase in hydrocarbon emissions is of the order of 100 ppmC1 per 1% liquid fuel introduced into the fresh oil in a methane fired engine at mid-speed and light load conditions. Calculations indicate that fuel desorbing from oil is much more likely to produce hydrocarbon emissions than fuel emerging from crevices. © Copyright 1994 Society of Automotive Engineers, Inc.

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.

BWR fuel design options for self-sustainable Th-233U fuel cycle

In this work, we investigate a number of fuel assembly design options for a BWR core operating in a closed self-sustainable Th-233U fuel cycle. The designs rely on axially heterogeneous fuel assembly structure in order to improve fertile to fissile conversion ratio. One of the main assumptions of the current study was to restrict the fuel assembly geometry to a single axial fissile zone "sandwiched" between two fertile blanket zones. The main objective was to study the effect of the most important design parameters, such as dimensions of fissile and fertile zones and average void fraction, on the net breeding of 233U. The main design challenge in this respect is that the fuel breeding potential is at odds with axial power peaking and therefore limits the maximum achievable core power rating. The calculations were performed with BGCore system, which consists of MCNP code coupled with fuel depletion and thermo-hydraulic feedback modules. A single 3-dimensional fuel assembly with reflective radial boundaries was modeled applying simplified restrictions on maximum central line fuel temperature and Critical Power Ratio. It was found that axially heterogeneous fuel assembly design with single fissile zone can potentially achieve net breeding. In this case however, the achievable core power density is roughly one third of the reference BWR core.

Neutronic optimization in high conversion Th-233U fuel assembly with simulated annealing

This paper reports on fuel design optimization of a PWR operating in a self sustainable Th-233U fuel cycle. Monte Carlo simulated annealing method was used in order to identify the fuel assembly configuration with the most attractive breeding performance. In previous studies, it was shown that breeding may be achieved by employing heterogeneous Seed-Blanket fuel geometry. The arrangement of seed and blanket pins within the assemblies may be determined by varying the designed parameters based on basic reactor physics phenomena which affect breeding. However, the amount of free parameters may still prove to be prohibitively large in order to systematically explore the design space for optimal solution. Therefore, the Monte Carlo annealing algorithm for neutronic optimization is applied in order to identify the most favorable design. The objective of simulated annealing optimization is to find a set of design parameters, which maximizes some given performance function (such as relative period of net breeding) under specified constraints (such as fuel cycle length). The first objective of the study was to demonstrate that the simulated annealing optimization algorithm will lead to the same fuel pins arrangement as was obtained in the previous studies which used only basic physics phenomena as guidance for optimization. In the second part of this work, the simulated annealing method was used to optimize fuel pins arrangement in much larger fuel assembly, where the basic physics intuition does not yield clearly optimal configuration. The simulated annealing method was found to be very efficient in selecting the optimal design in both cases. In the future, this method will be used for optimization of fuel assembly design with larger number of free parameters in order to determine the most favorable trade-off between the breeding performance and core average power density.

Comparison of square and hexagonal fuel lattices for high conversion PWRs

This paper reports on an investigation into fuel design choices of a pressurized water reactor operating in a self-sustainable Th- 233U fuel cycle. In order to evaluate feasibility of this concept, two types of fuel assembly lattices were considered: square and hexagonal. The hexagonal lattice may offer some advantages over the square one. For example, the fertile blanket fuel can be packed more tightly reducing the blanket volume fraction in the core and potentially allowing to achieve higher core average power density. The calculations were carried out with Monte-Carlo based BGCore code system and the results were compared to those obtained with Serpent Monte-Carlo code and deterministic transport code BOXER. One of the major design challenges associated with the SB concept is high power peaking due to the high concentration of fissile material in the seed region. The second objective of this work is to estimate the maximum achievable core power density by evaluation of limiting thermal hydraulic parameters. The analysis showed that both fuel assembly designs have a potential of achieving net breeding. Although hexagonal lattice was found to be somewhat more favorable because it allows achieving higher power density, while having breeding performance comparable to the square lattice case. © Carl Hanser Verlag München.

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.

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.

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.

The potential use of 241Am as proliferation resistant burnable poison in PWRs

This paper discusses the use of 241Am as proliferation resistant burnable poison for light water reactors. Homogeneous addition of small (as little as 0.12%) amounts of 241Am to the conventional light water reactor fuel results in significant increase in 238Pu/Pu ratio in the discharged fuel improving its proliferation resistance. Moreover, 241Am, admixed to the fuel, acts as burnable absorber allowing for substantial reduction in conventional reactivity control means without a notable fuel cycle length penalty. This is possible due to favorable characteristics of 241Am transmutation chain. The fuel cycle length penalty of introducing 241Am into the core is evaluated and discussed, as well as the impact of He production in the fuel pins and degradation of reactivity feedback coefficients. Proliferation resistance and reactivity control features related to the use of 241Am are compared to those of using 237Np, which has also been suggested as an additive to the conventional fuel in order to improve its proliferation resistance. It was found that 241Am admixture is more favorable than 237Np admixture because of the smaller fuel cycle length penalty and higher burnable poison savings. Addition of either 237Np or 241Am would provide substantial but not ultimate protection from misuse of Pu originating in the spent fuel from the commercial power reactors. Therefore, the benefits from application of the concept would have to be carefully evaluated against the additional costs and proliferation risks associated with manufacturing of 237Np or 241Am doped fuel. Although this work concerns specifically with PWRs, the conclusions could also be applied to BWRs and, to some extent, to other thermal spectrum reactor types. © 2009 Elsevier Ltd. All rights reserved.

Decay heat in fast reactors with transuranic fuels

In this work, we performed an evaluation of decay heat power of advanced, fast spectrum, lead and molten salt-cooled reactors, with flexible conversion ratio. The decay heat power was calculated using the BGCore computer code, which explicitly tracks over 1700 isotopes in the fuel throughout its burnup and subsequent decay. In the first stage, the capability of the BGCore code to accurately predict the decay heat power was verified by performing a benchmark calculation for a typical UO2 fuel in a Pressurized Water Reactor environment against the (ANSI/ANS-5.1-2005, "Decay Heat Power in Light Water Reactors," American National Standard) standard. Very good agreement (within 5%) between the two methods was obtained. Once BGCore calculation capabilities were verified, we calculated decay power for fast reactors with different coolants and conversion ratios, for which no standard procedure is currently available. Notable differences were observed for the decay power of the advanced reactor as compared with the conventional UO2 LWR. The importance of the observed differences was demonstrated by performing a simulation of a Station Blackout transient with the RELAP5 computer code for a lead-cooled fast reactor. The simulation was performed twice: using the code-default ANS-79 decay heat curve and using the curve calculated specifically for the studied core by BGCore code. The differences in the decay heat power resulted in failure to meet maximum cladding temperature limit criteria by ∼100 °C in the latter case, while in the transient simulation with the ANS-79 decay heat curve, all safety limits were satisfied. The results of this study show that the design of new reactor safety systems must be based on decay power curves specific to each individual case in order to assure the desired performance of these systems. © 2009 Elsevier B.V. All rights reserved.