Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MULTIPLICATION FACTORS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

Computation of a Canadian SCWR unit cell with deterministic and Monte Carlo codes

The Canadian SCWR has the potential to achieve the goals that the generation IV nuclear reactors must meet. As part of the optimization process for this design concept, lattice cell calculations are routinely performed using deterministic codes. In this study, the first step (self-shielding treatment) of the computation scheme developed with the deterministic code DRAGON for the Canadian SCWR has been validated. Some options available in the module responsible for the resonance self-shielding calculation in DRAGON 3.06 and different microscopic cross section libraries based on the ENDF/B-VII.0 evaluated nuclear data file have been tested and compared to a reference calculation performed with the Monte Carlo code SERPENT under the same conditions. Compared to SERPENT, DRAGON underestimates the infinite multiplication factor in all cases. In general, the original Stammler model with the Livolant-Jeanpierre approximations are the most appropriate self-shielding options to use in this case of study. In addition, the 89 groups WIMS-AECL library for slight enriched uranium and the 172 groups WLUP library for a mixture of plutonium and thorium give the most consistent results with those of SERPENT. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2012

Is Monte Carlo embarrassingly parallel?

Monte Carlo is often stated as being embarrassingly parallel. However, running a Monte Carlo calculation, especially a reactor criticality calculation, in parallel using tens of processors shows a serious limitation in speedup and the execution time may even increase beyond a certain number of processors. In this paper the main causes of the loss of efficiency when using many processors are analyzed using a simple Monte Carlo program for criticality. The basic mechanism for parallel execution is MPI. One of the bottlenecks turn out to be the rendez-vous points in the parallel calculation used for synchronization and exchange of data between processors. This happens at least at the end of each cycle for fission source generation in order to collect the full fission source distribution for the next cycle and to estimate the effective multiplication factor, which is not only part of the requested results, but also input to the next cycle for population control. Basic improvements to overcome this limitation are suggested and tested. Also other time losses in the parallel calculation are identified. Moreover, the threading mechanism, which allows the parallel execution of tasks based on shared memory using OpenMP, is analyzed in detail. Recommendations are given to get the maximum efficiency out of a parallel Monte Carlo calculation. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2013

The impact of fuel particle size distribution on neutron transport in stochastic media

This paper presents a study of the particle size distribution impact on neutron transport in three-dimensional stochastic media. An eigenvalue problem is simulated in a cylindrical container consisting of fissile fuel particles with five different size distributions: constant, uniform, power, exponential and Gaussian. We construct 15 cases by altering the fissile particle volume packing fraction and its optical thickness, but keeping the mean chord length of the spherical fuel particle the same at different size distributions. The tallied effective multiplication factor (k{sub eff}) and flux distribution along axial and radial directions are compared between different size distributions. At low packing fraction and low optical thickness, the size distribution has a significant impact on radiation transport in stochastic media, which can cause as high as {approx}270 pcm difference in k{sub eff} value and {approx}2.6% relative error difference in peak flux. As the packing fraction and optical thickness increase, the impact gradually dissipates. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Criticality calculation of non-ordinary systems

The specific features of calculation of the effective multiplication factor using the Monte Carlo method for weakly coupled and non-asymptotic multiplying systems are discussed. Particular examples are considered and practical recommendations on detection and Monte Carlo calculation of systems typical in numerical substantiation of nuclear safety for VVER fuel management problems are given. In particular, the problems of the choice of parameters for the batch mode and the method for normalization of the neutron batch, as well as finding and interpretation of the eigenvalue spectrum for the integral fission matrix, are discussed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2012

Comparison of XSUSA and "Two-Step" Approaches for Full Core Uncertainty Quantification

While there are multiple sources of error that are introduced into the standard computational regime for simulating reactor cores, rigorous uncertainty analysis methods are available primarily for quantifying the effects of cross section uncertainties. Two methods for propagating cross section uncertainties through core simulators are the XSUSA statistical approach and the 'Two-Step' method. The XSUSA approach, which is based on the SUSA code package, is fundamentally a stochastic sampling method. Alternatively, the Two-Step method utilizes generalized perturbation theory in the first step and stochastic sampling in the second step. The consistency of these two methods in quantifying uncertainties in the multiplication factor and in the core power distribution will be examined in the framework of phase I-3 of the UAM Benchmark. Using the TMI core as a base model for analysis, the XSUSA and Two-Step methods are applied with certain limitations and the results are compared to those produced by other stochastic sampling-based codes. Based on the uncertainty analysis results, conclusions are made for which method is currently a more viable option for computing uncertainties in burnup and transient calculations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2013

Natural equilibria in steady-state neutron diffusion with temperature feedback

The critical diffusion equation with feedback is investigated within the context of steady-state multiphysics. It is proposed that for critical configurations there is no need to include the multiplication factor k in the formulation of the diffusion equation. This is notable because exclusion of k from the coupled system of equations precludes the mathematically tenuous notion of a nonlinear eigenvalue problem. On the other hand, it is shown that if the factor k is retained in the diffusion equation, as is currently common practice, then the resulting problem is equivalent to the constrained minimization of a functional representing the critical equilibrium of neutron and temperature distributions. The unconstrained solution corresponding to k = 1 represents the natural equilibrium of a critical system at steady-state. Computational methods for solving the constrained problem (with k) are briefly reviewed from the literature and a method for the unconstrained problem (without k) is outlined. A numerical example is studied to examine the effects of the constraint in the nonlinear system. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

A demonstration of a whole core neutron transport method in a gas cooled reactor

This paper illustrates a capability of the whole core transport method COMET. Building on previous works which demonstrated the accuracy of the method, this work serves to emphasize the robust capability of the method while also accentuating its efficiency. A set of core configurations is presented based on an operating gas-cooled thermal reactor, Japan's HTTR, and COMET determines the eigenvalue and fission density profile throughout each core configuration. Results for core multiplication factors are compared to MCNP for accuracy and also to compare runtimes. In all cases, the values given by COMET differ by those given by MCNP by less than the uncertainty inherent in the stochastic solution procedure, however, COMET requires runtimes shorter on the order of a few hundred. Figures are provided illustrating the whole core fission density profile, with segments of pins explicitly modeled individually, so that pin-level neutron flux behavior can be seen without any approximation due to simplification strategies such as homogenization. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2012

Safety and power multiplication aspects of mirror fusion-fission hybrids

Recently, in a research project at Uppsala University a simplified neutronic model for a straight field line mirror hybrid has been devised and its most important operation parameters have been calculated under the constraints of a fission power production of 3 GW and that the effective multiplication factor k{sub eff} does not exceed 0.95. The model can be considered as representative for hybrids driven by other types of mirrors too. In order to reduce the demand on the fusion power of the mirror, a modified option of the hybrid has been considered that generates a reduced fission power of 1.5 GW with an increased maximal value k{sub eff}=0.97. The present paper deals with nuclear safety aspects of this type of hybrids. It presents and discusses calculation results of reactivity effects as well as of driver effects.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2012

Comparison of XSUSA and 'two-step' approaches for full-core uncertainty quantification

While there are multiple sources of error that are introduced into the standard computational regime for simulating reactor cores, rigorous uncertainty analysis methods are available primarily for quantifying the effects of cross section uncertainties. Two methods for propagating cross section uncertainties through core simulators are the XSUSA statistical approach and the 'Two-Step' method. The XSUSA approach, which is based on the SUSA code package, is fundamentally a stochastic sampling method. Alternatively, the Two-Step method utilizes generalized perturbation theory in the first step and stochastic sampling in the second step. The consistency of these two methods in quantifying uncertainties in the multiplication factor and in the core power distribution will be examined in the framework of phase 1-3 of the UAM Benchmark. Using the TMI core as a base model for analysis, the XSUSA and Two-Step methods are applied with certain limitations and the results are compared to those produced by other stochastic sampling-based codes. Based on the uncertainty analysis results, conclusions are made for which method is currently a more viable option for computing uncertainties in burnup and transient calculations. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Effect of double heterogeneity treatment on neutronics modeling of HTGR unit cell - 18524

The High Temperature Gas-cooled Reactor (HTGR) utilizes unique coated fuel particle design, which is embedded within a graphite matrix and shaped into spherical or hexagonal fuel element. The feature of the inherent double heterogeneity in the HTGR fuel element therefore raises difficulties in modeling and requires special simulation strategies to address its impact on the self-shielding calculation. In this study, neutronics calculations with various levels of heterogeneity processing were performed for both spherical and hexagonal HTGR fuel elements. The results indicate that ignoring either level of double heterogeneity can lead to significant underestimation of self-shielding effect and consequently yield as much as 8% deviation in the prediction of the multiplication factor. Fission and capture rates were compared for the spherical HTGR fuel model and over 30% deviation was found in the resonance energy range if homogenized model used. Further investigation on the neutron spectrum demonstrated similar tendency for both types of fuel elements when applying the simplified self-shielding models, that is, {approx}20% underestimation of neutron flux in 6.5-6.75 eV energy interval was observed due to the incomplete evaluation of the cross section resonance self-shielding effect. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

Criticality detection method based on FP gamma radiation measurement

A method is proposed to evaluate the extent of subcriticality of an accident-damaged nuclear reactor. In this method the radiation ratio of two FP rare gas nuclides, {sup 88}Kr and {sup 133}Xe, is measured. From the measured value, the value in the core is estimated by correcting for the time required by the gases to diffuse from the core to the measuring point. A simple expression for neutron multiplication factor has been derived, which uses the corrected {sup 88}Kr-to-{sup 133}Xe radiation ratio and the {sup 88}Kr-to-{sup 133}Xe fission yield ratios from {sup 244}Cm and {sup 235}U but no information on the amount or distribution of fissile materials, making the proposed method very simple. The method has the advantage that FP rare gases can easily leak from the reactor core through many openings and gaps, reaching germanium counters without reacting with other materials. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2012

Flux stabilization in neutron problems with fixed sources

Although critical core calculations are the most common in design and safety analysis, fixed source calculations are needed for specific applications, e.g. to compute ex-core detector response functions, to develop new methodologies for dilution and reload error accidents and more in general for all situations involving sub-critical shut-down states. It is well known that the source problem becomes difficult to be solved with core configuration close to criticality, i.e. with the multiplication factor approaching unity, for the occurrence of numerical ill-conditioning and very high number of iterations, possibly leading to failure in the flux convergence. In this work, the Wielandt Eigen-shift technique used in iterative methods of critical problems is developed for source problems too, in order to stabilize the solution. The mathematical basis and the proof of the convergence are discussed. Compared to the existing methods, this technique allows also more control to avoid singular behavior at inner iterations. Numerical tests with a 1D analytical benchmark are reported to prove the robustness of the technique. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2012

Neutron diagnostics for mirror hybrids

Fusion-fission (FuFi) hybrids will need instrumentation to diagnose the deuteriumtritium plasma, whose 14-MeV neutron emission is the driver of the sub-critical fission core. While the fission neutron yield rate (Y{sub fi} and hence power P{sub fi}) can be monitored with standard instrumentation, fusion plasmas in hybrids require special diagnostics where the determination of Y{sub th} ({proportional_to}P{sub fu}) is a challenge. Information on Y{sub fu} is essential for assessing the fusion plasma performance which together with Y{sub fi} allows for the validation of the neutron multiplication factor (k) of the subcritical fission core. Diagnostics for hybrid plasmas are heuristically discussed with special reference to straight field line mirror (SFLM). Relevant DT plasma experience from JET and plans for ITER in the main line of fusion research were used as input. It is shown that essential SFLM plasma information can potentially be obtained with proposed instrumentation, but the state of the hybrid plasma must be predictably robust as derived from fully diagnosed dedicated experiments without interface restrictions of the hybrid application.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2012

Development of HELIOS/CAPP code system for the analysis of block type VHTR cores

In this paper, the HELIOS/CAPP code system developed for the analysis of block type VHTR cores is presented and verified against several VHTR core configurations. Verification results shows that HELIOS code predicts less negative MTC and RTC than McCARD code does and thus HELIOS code overestimates the multiplication factors at the states with high moderator and reflector temperature especially when the B{sub 4}C BP is loaded. In the depletion calculation for the VHTR single cell fuel element, the error of HELIOS code increases as burnup does. It is ascribed to the fact that HELIOS code treats some fission product nuclides with large resonances as non-resonant nuclides. In the 2-D core depletion calculation, a relatively large reactivity error is observed in the case with BP loading while the reactivity error in the case without BP loading is less than 300 pcm. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2013

Benchmark analysis of high temperature engineering test reactor core using McCARD code

A benchmark calculation has been performed for a startup core physics test of Japan's High Temperature Engineering Test Reactor (HTTR). The calculation is carried out by the McCARD code, which adopts the Monte Carlo method. The cross section library is ENDF-B/VII.0. The fuel cell is modeled by the reactivity-equivalent physical transform (RPT) method. Effective multiplication factors with different numbers of fuel columns have been analyzed. The calculation shows that the HTTR becomes critical with 19 fuel columns with an excess reactivity of 0.84% Δk/k. The discrepancies between the measurements and Monte Carlo calculations are 2.2 and 1.4 % Δk/k for 24 and 30 columns, respectively. The reasons for the discrepancy are thought to be the current version of cross section library and the impurity in the graphite which is represented by the boron concentration. In the future, the depletion results will be proposed for further benchmark calculations. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2011

MONTE-CARLO BURNUP CALCULATION UNCERTAINTY QUANTIFICATION AND PROPAGATION DETERMINATION

Reactor burnup or depletion codes are used thoroughly in the fields of nuclear forensics and nuclear safeguards. Two common codes include MONTEBURNS and MCNPX/CINDER. These are Monte-Carlo depletion routines utilizing MCNP for neutron transport calculations and either ORIGEN or CINDER for burnup calculations. Uncertainties exist in the MCNP steps, but this information is not passed to the depletion calculations or saved. To quantify this transport uncertainty and determine how it propagates between burnup steps, a statistical analysis of multiple repeated depletion runs is performed. The reactor model chosen is the Oak Ridge Research Reactor (ORR) in a single assembly, infinite lattice configuration. This model was burned for a 150 day cycle broken down into three steps. The output isotopics as well as effective multiplication factor (k-effective) were tabulated and histograms were created at each burnup step using the Scott Method to determine the bin width. The distributions for each code are a statistical benchmark and comparisons made. It was expected that the gram quantities and k-effective histograms would produce normally distributed results since they were produced from a Monte-Carlo routine, but some of the results appear to not. Statistical analyses are performed using the {chi}{sup 2} test against a normal distribution for the k-effective results and several isotopes including {sup 134}Cs, {sup 137}Cs, {sup 235}U, {sup 238}U, {sup 237}Np, {sup 238}Pu, {sup 239}Pu, and {sup 240}Pu.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2012

Fission Matrix Capability for MCNP Monte Carlo

In a Monte Carlo criticality calculation, before the tallying of quantities can begin, a converged fission source (the fundamental eigenvector of the fission kernel) is required. Tallies of interest may include powers, absorption rates, leakage rates, or the multiplication factor (the fundamental eigenvalue of the fission kernel, k eff ). Just as in the power iteration method of linear algebra, if the dominance ratio (the ratio of the first and zeroth eigenvalues) is high, many iterations of neutron history simulations are required to isolate the fundamental mode of the problem. Optically large systems have large dominance ratios, and systems containing poor neutron communication between regions are also slow to converge. The fission matrix method, implemented into MCNP, addresses these problems.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2015

Joint neutron measurements on the flat-top HEU critical assembly machine - 14388

Flat-Top is a critical assembly with a spherical core (either ∼17 kg of HEU or ∼6 kg of δ-phase Pu) surrounded by a spherical 1000 kg natural uranium reflector. In 2013, joint LANL and CEA measurements were conducted to measure Flat-Top configurations with the HEU core using He-3 detector systems. Similar experiments were previously documented in which the Caliban and Godiva critical assemblies were measured. The core and reflector of Flat-Top have a glory hole in which nuclear material or non-nuclear material may be placed. Five configurations were measured and included subcritical configurations at various reactivity states as well as a configuration at delayed critical. The various reactivity states were achieved by varying the placement of the safety blocks and control rods (composed of natural uranium). A description of the assembly and the specific configurations measured will be presented. The subcritical and delayed critical measurements were analyzed using the Hage-Cifarelli formalism of the Feynman Variance-to-Mean method. This analysis can be used to calculate the multiplication of each configuration, from which the multiplication factor and reactivity are determined. These results are compared to the reactivity states based upon control rod worth curves and MCNP simulations. For this work, a new MCNP model was made to accurately represent the measured configurations. The ICSBEP handbook contains a benchmark of Flat-Top but all models are very simple and do not have the required fidelity to accurately represent the measured configurations. This work also examines the number of neutrons emitted per fission for this assembly (which has a significant contribution of induced fissions in U-238 due to the large natural uranium reflector). (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗