Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “production simulation codes”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Reactive flow modeling of initial density effect on divergence JB-9014 detonation driving

A serious of experiments were designed and the results were represented in this paper, in which 2 mm thickness cooper shells were impacted by explosives named JB-9014 with different densities, and the surface velocities of the OFHC shells were measured. The comparison of experimental data shows the free surface velocity of the OFHC shell increase with the IHE density. Numerical modeling, which occupied phenomenological reactive flow rate model using the two-dimensional Lagrange hydrodynamic code, were carried out to simulate the above experiments, and empirical adjustments on detonation velocity and pressure and Pier Tang’s adjustments on EOS of detonation products were both introduced in our numerical simulation work. The computational results agree well with that of experiments, and the numerical results with original parameters of products and the adjusted ones of JB-9014 could describe the density effect distinctly.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Validation of Actuator Line and Actuator Disk Models with Filtered Lifting Line Corrections Implemented in Nalu-Wind Large Eddy Simulations of the Atmospheric Boundary Layer

Turbine generator power from simulations using Actuator Line Models and Actuator Disk Models with a Filtered Lifting Line Correction are compared to field data of a V27 turbine. Preliminary results of the wake characteristics are also presented. Turbine quantities of interest from traditional ALM and ADM with the Gaussian kernel (..epsilon..) set at the optimum value for matching power production and that resolve the kernel at all mesh sizes are also presented. The atmospheric boundary layer is simulated using Nalu-Wind, a Large Eddy Simulation code which is part of the ExaWind code suite. The effect of mesh resolution on quantities of interest is also examined.

Actuator Disk Models↗
NASA NTRS1994

Nonlinear evolution of tidally distorted accretion disks: Two-dimensional simulations

According to a previously published linear analysis, the tidal distortion of accretion disks in binary star systems produces a local hydrodynamic instability to m = 1 internal waves, which may have arbitrarily small wavelengths in the absence of viscosity. The instability is three-dimensional and approximately incrompressible. To explore the nonlinear outcome of this instability, we develop a shearing-sheet approximation on scales comparable to the disk thickness. The large-scale azimuthal variation of the disk is represented by varying the local metric with time (local orbital phase). The hydrodynamic equations can then be posed two-dimensionally on local meridional planes. We solve these equations with a second-order gasdynamical code based on the Total-Variation-Diminishing scheme. Our simulations confirm the predicted linear growth rate. The modes saturate chaotically at velocities scaling as the product of the linear growth rate and the wavelength. If the wavelength is small compared with the disk thickness, the modes remain nearly incompressible even when nonlinear. The two-dimensional power spectrum of velocities after saturation is roughly isotropic and extends over a broad range of scales in an approximately power-law fashion. We measure the heating rate associated with the nonlinear dissipation of the modes. The dissipation implies a secular torque on the disk and a return of angular momentum to the secondary star via the tidal potential. The estimated torque is somewhat larger than the tidal torque produced by maximal disk viscosity (alpha approximately 1). At least in these two-dimensional simulations, however, there is no significant angular momentum flux within the disk.

Ryu, Dongsu↗
DOE OSTI2020

Temperature response and brine availability to heated boreholes in bedded salt

Abstract There is a growing need for disposal of high‐level nuclear waste. To reduce uncertainty associated with brine availability to repository excavations in salt formations, a collaboration between Sandia, Los Alamos, and Lawrence Berkeley National Laboratories is performing a series of borehole‐scale coupled process tests. Here, we report on the first round of the Brine Availability Test in Salt (BATS) project, a “shakedown” experiment called Phase 1s. Experimental testing included placing a resistive heater, a 260‐W radiative heater, and a 750‐W radiative heater within previously drilled horizontal boreholes at the Waste Isolation Pilot Plant (WIPP) while monitoring temperature and water inflow. The experiments successfully achieved the targeted temperature of 120 °C when using the 750‐W radiative heater. Simulations using FEHM (Finite Element Heat and Mass transfer code) and TOUGH‐FLAC (Transport Of Unsaturated Groundwater and Heat–Fast Lagrangian Analysis of Continua) were able to accurately predict the coupled thermo–hydro–mechanical–chemical response of salt, matching the observed temperature and brine production. Due to the extremely low permeability of salt, these systems take many years to reach steady state when perturbed by mining activities. Long‐term numerical simulations are used to develop the initial pressure and saturation conditions. The inclusion of a damaged rock zone with higher permeability around the borehole also affects the saturation and pressure distributions and plays an important role in dissipating the potential for thermal pressurization. Knowledge gained from this round of experimentation and simulation will be used to conduct the next BATS project experiment in new boreholes at WIPP.

58 GEOSCIENCES↗
DOE OSTI2021

Design Study for Bio-Oil Production from Biomass Using a Dual Fluidized-Bed Reactor

To evaluate the re-design and reconfiguration of a dual-fluidized bed (DFB) gasification system into a recirculating pyrolysis reactor, Computation Fluid Dynamic (CFD) simulations of the system were conducted. The Barracuda Virtual Reactor® computational particle fluid dynamic code was used to perform simulations of the pyrolysis process. Modeling of the chemical reaction kinetics for both gas phase and solid particle phase were included. The recirculating pyrolysis reactor shown in Fig. 1a is based on a bubbling-bed biomass pyrolyzer and a riser combustor to convert the remaining char. The operational differences in the re-configuration of the DFB gasification system into a recirculating pyrolysis system for the production of bio-oil are (1) replacement of a low-surface area inert bed material with a high-surface-area bed material that has acidic properties to provide catalytic activity for the production of bio-oil with reduced oxygen content, (2) lower temperature and residence time for bio-oil production from pyrolysis, (3) replacement of the fluidization gas in the bubbling bed pyrolyzer from steam to nitrogen, and (4) the reduction of pyrolyzer freeboard volume. The bed material used for catalytic pyrolysis is Sasol 300 (300-micron dia., bulk density 0.94 kg/l, and surface area 130 m2/g) and is a theta-alumina with mild acidity. This is in comparison with previous standard bed material Carbo HSP (430-micron dia., bulk density 2.01 kg/l, and surface area 0.03 m2/g) used for gasification. For bio-oil production, pyrolysis in the bubbling bed requires temperatures in the range of 550 C in comparison with gasification temperatures near 850 C. To attain this lower temperature requires management of the energy mass balances, with control of the bed material recirculation rate between bubbling bed pyrolyzer and riser combustor, the introduction of a nitrogen purge in the pyrolyzer, and adjusting the pressure balance between the two vessels. To extract bio-oil from the pyrolysis reactor with a snorkel, two different freeboard configurations were evaluated. In Fig. 1 b the existing high freeboard configuration is shown and in Fig. 1 c the reduced freeboard design is presented. The introduction of a nitrogen purge for the high free board configuration provided the highest bio-oil production from the CFD simulations. To decrease the bed-material circulation rate, primary and secondary air on the combustor side were reduced, and the pressure on the combustor was slightly increased. A portion of the biomass and bio-oil was observed to be transported to the combustor, leading to a smaller pyrolysis yield. The control of the pyrolyzer temperature is performed by controlling the circulation rate. Good fluidization of the bubbling bed and cascade PID control are required to keep the temperature from oscillating due to large time delay experienced when changing primary and secondary air. (a) (b) (c) Figure 1. (a) Dual Fluidized Bed pyrolysis system configuration, (b) mole fraction of nonpolar biooil in high freeboard configuration, and (c) mole fraction of nonpolar bio-oil in low freeboard configuration with nitrogen purge introduced in both configurations.

BASIC BIOLOGICAL SCIENCES,BIOMASS FUELS↗
DOE OSTI2016

CASL Verification and Validation Plan

This report documents the Consortium for Advanced Simulation of LWRs (CASL) verification and validation plan. The document builds upon input from CASL subject matter experts, most notably the CASL Challenge Problem Product Integrators, CASL Focus Area leaders, and CASL code development and assessment teams. This document will be a living document that will track progress on CASL to do verification and validation for both the CASL codes (including MPACT, CTF, BISON, MAMBA) and for the CASL challenge problems (CIPS, PCI, DNB). The CASL codes and the CASL challenge problems are at differing levels of maturity with respect to validation and verification. The gap analysis will summarize additional work that needs to be done. Additional VVUQ work will be done as resources permit. This report is prepared for the Department of Energy’s (DOE’s) CASL program in support of milestone CASL.P13.02.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2016

Multi-Sensor Cloud and Aerosol Retrieval Simulator and Remote Sensing from Model Parameters : Aerosols - Part 2

The Multi-sensor Cloud Retrieval Simulator (MCRS) produces a simulated radiance product from any high-resolution general circulation model with interactive aerosol as if a specific sensor such as the Moderate Resolution Imaging Spectroradiometer (MODIS) were viewing a combination of the atmospheric column and land ocean surface at a specific location. Previously the MCRS code only included contributions from atmosphere and clouds in its radiance calculations and did not incorporate properties of aerosols. In this paper we added a new aerosol properties module to the MCRS code that allows users to insert a mixture of up to 15 different aerosol species in any of 36 vertical layers. This new MCRS code is now known as MCARS (Multi-sensor Cloud and Aerosol Retrieval Simulator). Inclusion of an aerosol module into MCARS not only allows for extensive, tightly controlled testing of various aspects of satellite operational cloud and aerosol properties retrieval algorithms, but also provides a platform for comparing cloud and aerosol models against satellite measurements. This kind of two-way platform can improve the efficacy of model parameterizations of measured satellite radiances, allowing the assessment of model skill consistently with the retrieval algorithm. The MCARS code provides dynamic controls for appearance of cloud and aerosol layers. Thereby detailed quantitative studies of the impacts of various atmospheric components can be controlled. In this paper we illustrate the operation of MCARS by deriving simulated radiances from various data field output by the Goddard Earth Observing System version 5 (GEOS-5) model. The model aerosol fields are prepared for translation to simulated radiance using the same model sub grid variability parameterizations as are used for cloud and atmospheric properties profiles, namely the ICA technique. After MCARS computes modeled sensor radiances equivalent to their observed counterparts, these radiances are presented as input to operational remote-sensing algorithms. Specifically, the MCARS-computed radiances are input into the processing chain used to produce the MODIS Data Collection 6 aerosol product (MOYD04). TheMOYD04 product is of course normally produced from MOYD021KM MODIS Level-1B radiance product directly acquired by the MODIS instrument. MCARS matches the format and metadata of a MOYD021KM product. The resulting MCARS output can be directly provided to MODAPS (MODIS Adaptive Processing System) as input to various operational atmospheric retrieval algorithms. Thus the operational algorithms can be tested directly without needing to make any software changes to accommodate an alternative input source. We show direct application of this synthetic product in analysis of the performance of the MOD04 operational algorithm. We use biomass-burning case studies over Amazonia employed in a recent Working Group on Numerical Experimentation (WGNE)-sponsored study of aerosol impacts on numerical weather prediction (Freitas et al., 2015). We demonstrate that a known low bias in retrieved MODIS aerosol optical depth appears to be due to a disconnect between actual column relative humidity and the value assumed by the MODIS aerosol product.

aerosol retrieval↗
DOE OSTI2020

Wind Farm Layout Optimization with Loads Considerations

The objective of this paper is to improve the annual energy production of a wind farm by optimizing the layout of a wind farm, while considering fatigue loads on turbines. In this paper, the loads are estimated using the edgewise bending moment computed using CCBlade, a steady-state blade element momentum code. The edgewise bending moment is then used to calculate fatigue damage using Miner's rule. The fatigue damage is used to constrain the layout optimization problem. We show that our method can predict blade root damage with similar trends to damage calculated with other methods, such as a complex, computationally expensive large-eddy simulation and unsteady aeroelastic code. We also optimize wind farm layouts and show that, for a simple problem with two wind directions and ten turbines, the fatigue damage can be constrained without sacrificing wind farm power production.

17 WIND ENERGY↗
DOE OSTI2014

Beam dynamics study of a 30 MeV electron linear accelerator to drive a neutron source

An experimental neutron facility based on 32 MeV/18.47 kW electron linac has been studied by means of PARMELA simulation code. Beam dynamics study for a traveling wave constant gradient electron accelerator is carried out to reach the preferential operation parameters (E = 30 MeV, P = 18 kW, dE/E < 12.47% for 99% particles). The whole linac comprises mainly E-gun, pre-buncher, buncher, and 2 accelerating columns. A disk-loaded, on-axis-coupled, 2π/3-mode type accelerating rf cavity is considered for this linac. After numerous optimizations of linac parameters, 32 MeV beam energy is obtained at the end of the linac. As high electron energy is required to produce acceptable neutron flux. The final neutron flux is estimated to be 5 × 10{sup 11} n/cm{sup 2}/s/mA. Future development will be the real design of a 30 MeV electron linac based on S band traveling wave.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

Multi-Physics Coupling to Model a RbCl-RbCl-Ga Target Stack

A computational model has been developed to simulate the performance of the Los Alamos Isotope Production Facility's (IPF) routine production target configuration, which consists of two Inconel-encapsulated RbCl salt targets and one niobium-encapsulated gallium target separated by water cooling channels. The methodology involves iteratively coupling particle transport and thermal hydraulic codes to determine the time dependent behavior of the target. Phase and density changes seen in targets are known to impact beam penetration and distribution. In addition, evidence of energy distributions seen in experimental radionuclide yield measurements at IPF do not match those seen in previous MCNP predictions. Density reduction as a result of heat deposition affects proton energy loss along the particle's track, impacting both radionuclide yields and target design. Accurate modeling of the thermal behavior of these multi-phase targets and the corresponding impact on beam behavior will significantly improve our understanding of incident and exit proton energy distributions and resulting radionuclide yields in downstream targets. Thermal hydraulic modeling, including computational fluid dynamics (CFD) and heat transfer, was performed with ANSYS CFX while particle transport was performed with the Monte Carlo N-Particle code (MCNP). (authors)

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2019

Criticality Accident Alarm System Shielding Benchmark: Integral Experiment Request 498, Critical Engineering Decision 1 Report

This report supports Critical/Subcritical Experiment Decision 1 (CED-1) for Integral Experiment Request (IER) 498, a shielding benchmark produced using the Godiva IV critical assembly at the National Critical Experiments Research Center (NCERC). Shielding benchmarks are datasets collected to provide a comparison for computer models that simulate the effect of shielding on radiation transmission from a source to a detector. The purpose of a CED-1 report is to produce a conceptual design for an integral experiment and to identify the long lead-time items. The design presented in this report favors certain detailed design options such as material selections. However, alternative options can be selected as the project evolves. Dimensions are specified only as needed to evaluate feasibility and performance. Development of Criticality Accident Alarm System (CAAS) shielding benchmarks supports two main programmatic needs. The first need is to improve confidence in models of facilities that house nuclear material handling operations. Better models can support effective CAAS detector placement when criticality may result in dose, and it enables more accurate analysis of criticality accident impacts at various locations in a facility. A key quantity that can be verified is the secondary gamma production caused by neutron interactions with various materials. The second programmatic need supported by shielding benchmarks is the ongoing development of methods being integrated into nuclear simulation codes such as the Monte Carlo N-Particle (MCNP) transport code and the SCALE code system used in efforts funded by the Nuclear Criticality Safety Program (NCSP). Previous work performed under the International Criticality Safety Benchmark Experimental Program (ICSBEP) to produce shielding benchmark data is primarily limited to a set of measurements at the SILENE critical assembly at the Valduc facility in France. The experiment used an open configuration and a burst mode operation, which limited detector options to thermoluminescent dosimeters (TLDs) and activation foils. Because TLDs do not record gamma spectra, it was difficult to verify the gamma energy spectrum from the experiment. In this work, the Godiva IV assembly will be used because it is a fast system capable of producing a wide range of source strengths, including super-prompt critical bursts. The room return for the facility is expected to be substantial and challenging to characterize, so a more closed, shielded experimental configuration is considered. Furthermore, because the physics of radiation transport outside the assembly does not change between the burst mode and the delayed critical mode, it is anticipated that most measurements will be made with the Godiva IV assembly in delayed critical mode, operating within a range that will allow electronic counters to collect accurate data. Scoping calculations concluded that a workable room return shield is comprised of 40 cm of 5% borated polyethylene backed by 7 cm of lead. This configuration was shown to produce a ratio of signal–to–room return of over 100 for up to 2 mean free paths of shield sample material, and over 10 for up to 4 mean free paths of shield sample material, ensuring that useful results can be obtained. With no shielding sample in the room return shield, it was shown that the neutron and gamma spectra do not appreciably change with a change in room size. This scoping calculation does not preclude the use of tungsten powder in a binder, some other effective gamma shielding material, or a hybrid design. Dose rates for various source intensities at an unshielded detector are listed herein to inform detector selection during CED-2, and calibration chains are identified. Preliminary estimates of costs and schedules are outlined but remain to be verified. Values could double or triple as estimates become more accurate and are principally included here to demonstrate feasibility. Preliminary estimates indicate that the apparatus can be installed within a feasible schedule and budget. Two measurements per day are assumed. Producing a shielding benchmark at NCERC using the Godiva IV assembly has been shown to be feasible, with manageable uncertainties and risks. A clear, continuing need for shielding benchmark data remains present, and this report identifies and evaluates a method for creating integrated shielding benchmark datasets.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2003

Team Software Development for Aerothermodynamic and Aerodynamic Analysis and Design

A collaborative approach to software development is described. The approach employs the agile development techniques: project retrospectives, Scrum status meetings, and elements of Extreme Programming to efficiently develop a cohesive and extensible software suite. The software product under development is a fluid dynamics simulator for performing aerodynamic and aerothermodynamic analysis and design. The functionality of the software product is achieved both through the merging, with substantial rewrite, of separate legacy codes and the authorship of new routines. Examples of rapid implementation of new functionality demonstrate the benefits obtained with this agile software development process. The appendix contains a discussion of coding issues encountered while porting legacy Fortran 77 code to Fortran 95, software design principles, and a Fortran 95 coding standard.

Alexandrov, N.↗
NASA NTRS2012

New Modeling Approaches to Study DNA Damage by the Direct and Indirect Effects of Ionizing Radiation

DNA is damaged both by the direct and indirect effects of radiation. In the direct effect, the DNA itself is ionized, whereas the indirect effect involves the radiolysis of the water molecules surrounding the DNA and the subsequent reaction of the DNA with radical products. While this problem has been studied for many years, many unknowns still exist. To study this problem, we have developed the computer code RITRACKS [1], which simulates the radiation track structure for heavy ions and electrons, calculating all energy deposition events and the coordinates of all species produced by the water radiolysis. In this work, we plan to simulate DNA damage by using the crystal structure of a nucleosome and calculations performed by RITRACKS. The energy deposition events are used to calculate the dose deposited in nanovolumes [2] and therefore can be used to simulate the direct effect of the radiation. Using the positions of the radiolytic species with a radiation chemistry code [3] it will be possible to simulate DNA damage by indirect effect. The simulation results can be compared with results from previous calculations such as the frequencies of simple and complex strand breaks [4] and with newer experimental data using surrogate markers of DNA double ]strand breaks such as . ]H2AX foci [5].

Plante, Ianik↗
NASA NTRS2019

Operational Stability Limits in Rotating Detonation Engine Numerical Simulations

An instability is described which arises in computational fluid dynamic (CFD) simulations of semi-idealized rotating detonation engines (RDE) configured with a throat at the exit. Its existence is verified by examining output from two independently developed CFD codes simulating the same configuration and producing solutions that agree well. The instability is shown to be thermo-acoustic in that a spatial integral of the product of pressure and heat release fluctuations develops a regular oscillation which grows in time. The instability can become severe enough to cause detonation failure. Its onset is shown to be closely linked to the size of the exit throat and the size of the inlet restriction; both parameters that strongly influence RDE performance. It is shown that the instability places a cap on ideal RDEperformance, but that an optimized exhaust throat and inlet restriction combination still yields substantial pressure gain. Other parametric sensitivities are also examined in terms ofinstability growth. These include axial length, inlet manifold pressure, and air-fuel ratio.

Paxson, Daniel E.↗
NASA NTRS2016

Validation of High-Fidelity CFD/CAA Framework for Launch Vehicle Acoustic Environment Simulation against Scale Model Test Data

A hybrid Computational Fluid Dynamics and Computational Aero-Acoustics (CFD/CAA) modeling framework has been developed for launch vehicle liftoff acoustic environment predictions. The framework couples the existing highly-scalable NASA production CFD code, Loci/CHEM, with a high-order accurate discontinuous Galerkin solver developed in the same production framework, Loci/THRUST, to accurately resolve and propagate acoustic physics across the entire launch environment. Time-accurate, Hybrid RANS/LES CFD modeling is applied for predicting the acoustic generation physics at the plume source, and a high-order accurate unstructured discontinuous Galerkin (DG) method is employed to propagate acoustic waves away from the source across large distances using high-order accurate schemes. The DG solver is capable of solving 2nd, 3rd, and 4th order Euler solutions for non-linear, conservative acoustic field propagation. Initial application testing and validation has been carried out against high resolution acoustic data from the Ares Scale Model Acoustic Test (ASMAT) series to evaluate the capabilities and production readiness of the CFD/CAA system to resolve the observed spectrum of acoustic frequency content. This paper presents results from this validation and outlines efforts to mature and improve the computational simulation framework.

Liever, Peter A.↗
NASA NTRS2016

Validation of High-Fidelity CFD/CAA Framework for Launch Vehicle Acoustic Environment Simulation against Scale Model Test Data

A hybrid Computational Fluid Dynamics and Computational Aero-Acoustics (CFD/CAA) modeling framework has been developed for launch vehicle liftoff acoustic environment predictions. The framework couples the existing highly-scalable NASA production CFD code, Loci/CHEM, with a high-order accurate Discontinuous Galerkin solver developed in the same production framework, Loci/THRUST, to accurately resolve and propagate acoustic physics across the entire launch environment. Time-accurate, Hybrid RANS/LES CFD modeling is applied for predicting the acoustic generation physics at the plume source, and a high-order accurate unstructured mesh Discontinuous Galerkin (DG) method is employed to propagate acoustic waves away from the source across large distances using high-order accurate schemes. The DG solver is capable of solving 2nd, 3rd, and 4th order Euler solutions for non-linear, conservative acoustic field propagation. Initial application testing and validation has been carried out against high resolution acoustic data from the Ares Scale Model Acoustic Test (ASMAT) series to evaluate the capabilities and production readiness of the CFD/CAA system to resolve the observed spectrum of acoustic frequency content. This paper presents results from this validation and outlines efforts to mature and improve the computational simulation framework.

Liever, Peter A.↗
DOE OSTI2019

NOBLE METAL MASS TRANSPORT MODEL FOR MOLTEN SALT REACTOR ANALYSIS IN VERA-CS

Molten Salt Reactor designs require robust multiphysics modelling capabilities to further the development of this next generation technology. Due to the unique liquid fuel design, traditional reactor physics analysis codes do not capture all of the important physics of these reactors. This work builds upon a general species transport model implemented into the multiphysics core simulator suite VERA-CS, by adding source and sink models for non-soluble fission products, specifically the noble metals. Noble metal fission products are born in the fuel-salt, but do not form stable fluorides, and instead plate out on various surfaces in the reactor fuel loop. A concentration gradient diffusion-convection based mass transfer model is implemented into the general species transport model within CTF-the subchannel thermal hydraulics code in VERA-CS-and verified. Two noble metal decay chains are analyzed, and their steady state bulk and surface concentrations are calculated for a simple flow loop in CTF that is representative of the Oak Ridge National Laboratory Molten Salt Reactor Experiment (MSRE). Non-soluble species 99Mo to 99Tc captures wall plate out of decaying species, and non-soluble 135Te to soluble 135I and gaseous 135Xe captures wall loss of decaying species. Results are compared with MSRE data and the possible multiphysics effects of this phenomena are identified. Future work to enhance the robustness of the physical model and couple it with neutronics to evaluate these effects of Noble Metal mass transport is discussed.

Walker, Samuel↗
DOE OSTI2018

An improved target velocity sampling algorithm for free gas elastic scattering

We present an improved algorithm for sampling the target velocity when simulating elastic scattering in a Monte Carlo neutron transport code that correctly accounts for the energy dependence of the scattering cross section. The algorithm samples the relative velocity directly, thereby avoiding a potentially inefficient rejection step based on the ratio of cross sections. Here, we have shown that this algorithm requires only one rejection step, whereas other methods of similar accuracy require two rejection steps. The method was verified against stochastic and deterministic reference results for upscattering percentages in 238 U. Simulations of a light water reactor pin cell problem demonstrate that using this algorithm results in a 3% or less penalty in performance when compared with an approximate method that is used in most production Monte Carlo codes

22 GENERAL STUDIES OF NUCLEAR REACTORS↗