Monte Carlo methods for radiation transport analysis on vector computers

The development of advanced computers with special capabilities for vectorized or parallel calculations demands the development of new calculational methods. The very nature of the Monte Carlo process precludes direct conversion of old (scalar) codes to the new machines. Instead, major changes in global algorithms and careful selection of compatible physics treatments are required. Recent results for Monte Carlo in multigroup shielding applications and in continuous-energy reactor lattice analysis have demonstrated that Monte Carlo methods can be successfully vectorized. The significant effort required for stylized coding and major algorithmic changes is worthwhile, and significant gains in computational efficiency are realized. Speedups of at least twenty to forty times faster than CDC-7600 scalar calculations have been achieved on the CYBER-205 without sacrificing the accuracy of standard Monte Carlo methods. Speedups of this magnitude provide reductions in statistical uncertainties for a given amount of computing time, permit more detailed and realistic problems to be analyzed, and make the Monte Carlo method more accessible to nuclear analysts. Following overviews of the Monte Carlo method for particle transport analysis and of vector computer hardware and software characteristics, both general and specific aspects of the vectorization of Monte Carlo are discussed. Finally, numerical results obtained from vectorized Monte Carlo codes run on the CYBER-205 are presented.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/24996/1/0000423.pd

On-the-fly doppler broadening of unresolved resonance region cross sections via probability band interpolation

In this work we present a scheme for computing temperature-dependent unresolved resonance region cross sections in Monte Carlo neutron transport simulations. This approach relies on the generation of equiprobable cross section magnitude bands on an energy-temperature mesh. The bands are then interpolated in energy and temperature to obtain a cross section value. This is in contrast to the typical procedure of pre-generating probability tables at all temperatures present in the simulation. As part of this work, a flexible probability table generation capability is integrated into the continuous-energy neutron transport code OpenMC [1]. Both single-level and multi-level Breit-Wigner formalisms are supported, as is modeling of the resonance structure of competitive reactions. A user-specified cross section band tolerance is enabled with batch statistics. Probability tables are generated for all 268 ENDF/B-VII.1 [2] isotopes that have an unresolved resonance region evaluation. Integral benchmark simulations of the Big Ten critical assembly show that, for a system that is sensitive to the unresolved resonance region, a temperature interval of ∼200 K around 293.6 K is sufficient to reproduce the keff value that is obtained with probability tables generated exactly at room temperature. A finer mesh of < 50 K is required to reproduce some cross section values at the common target relative difference of 0.1

Direct, on-the-fly calculation of unresolved resonance region cross sections in Monte Carlo simulations

The theory, implementation, and testing of a method for on-the-fly unresolved resonance region cross section calculations in continuous-energy Monte Carlo neutron transport codes are presented. With this method, each time that a cross section value is needed within the simulation, a realization of unresolved resonance parameters is generated about the desired energy and temperature-dependent single-level Breit-Wigner resonance cross sections are computed directly via use of the analytical ψ − χ Doppler integrals. Results indicate that, in room-temperature simulations of a system that is known to be highly sensitive to the effects of resonance structure in unresolved region cross sections, the on-the-fly treatment produces results that are in excellent agreement with those produced with the well-established probability table method. Additionally, similar agreement is observed between results obtained from the on-the-fly and probability table methods for another intermediate spectrum system at temperatures of 293.6 K and 2500 K. With relatively tight statistical uncertainties at the ∼ 10 pcm level, all on-the-fly and probability table keff eigenvalues agree to within 2σ. Also, we use the on-the-fly approach to show that accounting for the resonance structure of competitive reaction cross sections can have non-negligible effects for intermediate/fast spectrum systems. Biases of up to 90 pcm are observed. Finally, the consequences of the on-the-fly method with respect to simulation runtime and memory requirements are briefly discussed.United States. Department of Energy (Consortium for Advanced Simulation of Light Water Reactors. Contract DE-AC05-00OR22725

Improved method for implicit Monte Carlo

The Implicit Monte Carlo (IMC) method has been used for over 30 years to analyze radiative transfer problems, such as those encountered in stellar atmospheres or inertial confinement fusion. Reference [2] provided an exact error analysis of IMC for 0-D problems and demonstrated that IMC can exhibit substantial errors when timesteps are large. These temporal errors are inherent in the method and are in addition to spatial discretization errors and approximations that address nonlinearities (due to variation of physical constants). In Reference [3], IMC and four other methods were analyzed in detail and compared on both theoretical grounds and the accuracy of numerical tests. As discussed in, two alternative schemes for solving the radiative transfer equations, the Carter-Forest (C-F) method and the Ahrens-Larsen (A-L) method, do not exhibit the errors found in IMC; for 0-D, both of these methods are exact for all time, while for 3-D, A-L is exact for all time and C-F is exact within a timestep. These methods can yield substantially superior results to IMC

Transit times and mean ages for nonautonomous and autonomous compartmental systems

We develop a theory for transit times and mean ages for nonautonomous compartmental systems. Using the McKendrick-von F\"orster equation, we show that the mean ages of mass in a compartmental system satisfy a linear nonautonomous ordinary differential equation that is exponentially stable. We then define a nonautonomous version of transit time as the mean age of mass leaving the compartmental system at a particular time and show that our nonautonomous theory generalises the autonomous case. We apply these results to study a nine-dimensional nonautonomous compartmental system modeling the terrestrial carbon cycle, which is a modification of the Carnegie-Ames-Stanford approach (CASA) model, and we demonstrate that the nonautonomous versions of transit time and mean age differ significantly from the autonomous quantities when calculated for that model

A Mixed-Method Multiple Case Study of Three Business Models for Local Healthy Food Delivery Systems in Underprivileged Urban Areas

Define a USDA food desert Analyze the three business models for benefits and deficiencies Identify methods for starting a sustainable mobile farmers market program. Over 40% of Berrien County Michigan land use is agricultural. Many products are fruits and vegetables. Yet, the county has six identified USDA defined Food Deserts. Past research, based on a mobile farmers market, confirmed local trends and deficits. The purpose for this research is to define a sustainable business model that delivers healthy local food options to USDA Food Deserts combating food inequity. A mixed method multiple case study was created to test three business models in one of the local Food Deserts. Results identify a sustainable model that provides local healthy food options. Model One developed lunchtime stops serving robust local businesses and community epicenters, with a goal of creating lunchtime profit, focusing evenings on service stops at a potential loss. Model Two required local farmers to provide produce at no cost to the market. Market locations were service stops embedded inside the residential community and stops increased to 15 per week. Model Three purchased local produce and focused stops on work and shopping locals, targeting the ALICE poverty segment. Models were tested for two weeks each with data collected for correlational comparison. Results show a high demand for lunch hour food options. Labor at such stops was increased over service stops. Model Two proved willingness from local farmers to support neighboring Food Deserts, yet sales didn’t equal market costs. Model Three demonstrated increased demand from ALICE segments. Conclusions indicate that sustainability could be reached with a hybrid model with limited volunteer intern positions. A three to five-year plan should be built with funding support ebbing with market growth

11.1: Invited Paper : Advances in Blue Phosphorescent Organic Light‐Emitting Devices

This paper discusses the latest developments towards a commercial blue phosphorescent organic light emitting device (PHOLED™) technology. Progress towards achieving a high efficiency, long‐lived saturated blue PHOLED is discussed. First, a high efficiency (20% EQE, 45 cd/A), light blue (0.17, 0.39) PHOLED is presented. Next, long‐lived blue PHOLEDs having chromaticity co‐ordinates (0.17, 0.38) and (0.16, 0.29) are estimated to degrade to half their initial luminance of 200cd/m 2 after >100,000 hrs and 17,500 hrs, respectively. Finally, results from PHOLEDs designed to increase blue color saturation and lifetime are presented.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/92011/1/1.2433213.pd

THE EDINBURGH RANDOMIZED TRIAL OF BREAST-CANCER SCREENING - RESULTS AFTER 10 YEARS OF FOLLOW-UP

The Edinburgh Randomised Trial of Breast Cancer Screening recruited 44,288 women aged 45-64 years into the initial cohort of the trial during 1978-81, and 10 years of follow-up is now complete. A total of 22,944 women were randomised into the study group and were offered screening for 7 years; the remaining women formed the control group. After 10 years, breast cancer mortality is 14-21% lower in the study group than in the controls depending on the precise definition of the end point. These differences are not statistically significant; for breast cancer as the underlying cause of death the relative risk is 0.82 (95% confidence interval 0.61-1.11). Rates of locally advanced and metastatic cancer were substantially lower in the study group, but screening has failed to achieve marked reductions in rates of small node-positive cancers. Those women who accepted the final invitation to screening have been monitored over the 3 year period prior to their first screen under the UK service screening programme. Interval cases, expressed as a proportion of the control incidence, increased from 12% in the first year to 67% in the third year. The reduction in breast cancer mortality for older women (aged at least 50 years) is the same as that for the total study group for this duration of follow-up. For analyses of breast cancer mortality in younger women updates recruited to the trial from 1982 to 1985 (10,383 women with 6-8 years' follow-up) have been included. The reduction in breast cancer mortality for women aged 45-49 years at entry was 22% (relative risk = 0.78, 95% confidence interval = 0.46-1.31)

47.4: Blue Phosphorescent Organic Light Emitting Device Stability Analysis

A model based on defect generation by exciton‐polaron annihilation interactions between the emitter and host molecules, in a blue phosphorescent OLED, is shown to fit well with experimental data. A blue PHOLED with (0.15, 0.25) chromaticity is shown to have a half‐life, from 1,000 nits, of 690 hrs.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/92134/1/1.3069766.pd