Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Code of Accounts”

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 1,603 records · Page 89

A Modular Approach to Model Oscillating Control Surfaces Using Navier Stokes Equations

The use of active controls for rotorcraft is becoming more important for modern aerospace configurations. Efforts to reduce the vibrations of helicopter blades with use of active-controls are in progress. Modeling oscillating control surfaces using the linear aerodynamics theory is well established. However, higher-fidelity methods are needed to account for nonlinear effects, such as those that occur in transonic flow. The aeroelastic responses of a wing with an oscillating control surface, computed using the transonic small perturbation (TSP) theory, have been shown to cause important transonic flow effects such as a reversal of control surface effectiveness that occurs as the shock wave crosses the hinge line. In order to account for flow complexities such as blade-vortex interactions of rotor blades higher-fidelity methods based on the Navier-Stokes equations are used. Reference 6 presents a procedure that uses the Navier-Stokes equations with moving-sheared grids and demonstrates up to 8 degrees of control-surface amplitude, using a single grid. Later, this procedure was extended to accommodate larger amplitudes, based on sliding grid zones. The sheared grid method implemented in EulerlNavier-Stokes-based aeroelastic code ENS AERO was successfully applied to active control design by industry. Recently there are several papers that present results for oscillating control surface using Reynolds Averaged Navier-Stokes (RANS) equations. References 9 and 10 report 2-D cases by filling gaps with overset grids. Reference 9 compares integrated forces with the experiment at low oscillating frequencies whereas Ref. 10 reports parametric studies but with no validation. Reference II reports results for a 3D case by modeling the gap region with a deformed grid and compares force results with the experiment only at the mid-span of flap. In Ref. II grid is deformed to match the control surface deflections at the section where the measurements are made. However, there is no indication in Ref. II that the gaps are explicitly modeled as in Ref. 6. Computations using overset grids are reported in Ref. 12 for a case by adding moving control surface to an existing blade but with no validation either with an experiment or another computation.

Navier↗

In service inspection and repair of sodium cooled ASTRID prototype

In the frame of the large R and D work which is performed for the future ASTRID sodium cooled prototype, In Service Inspection and Repair (ISI and R) has been identified as a major issue to be taken into account in order to enlarge the plant safety, to consolidate its availability and to protect the associated investment. After the first part of pre-conceptual design phase (2008-2012), the running second part of pre-conceptual phase (2013-2015) allows to increase the ISI and R tool ability for immersed sodium structures of ASTRID, at about 200 deg. C, on the basis of consolidated specifications and thanks to their qualification through more and more realistic laboratory tests and simulation with CIVA code. ISI and R items are being developed and qualified during a pluri-annual program which mainly deals with the reactor block structures, the primary components and circuit, and the Power Conversion System. It ensures a strong connection between the reactor designers and inspection specialists, as the optimization of inspectability and repairability is looked at: this already induced specific rules for design, in order to shorten and ease the ISI and R operations, which have been merged into RCC-MRx rules. In the frame of increasing technology readiness level with corresponding performance demonstration, this paper presents R and D dealing with the ISI and R items: it highlights the sensor development (both ultrasonic and electromagnetic concepts, compatible with sodium at 200 deg. C), then their applications for ASTRID structure control (under sodium telemetry, imaging and NDE). Activity for repair is also presented (a single laser tool for sodium sweeping, machining and welding), and finally the effort for associated robotic (generic program for ASTRID applications, specific technological tools for sodium medium, tight immersed bell). The main results of testing and simulation are given for telemetry, vision, NDE applications, laser process repair and under sodium sealing. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Disclination-dislocation based model for grain boundary stress field evolution due to slip transmission history and influence on subsequent dislocation transmission

This work demonstrates how the structure of a grain boundary (GB) and its evolution due to slip transmission history influences subsequent dislocation transmission. First, a model for the evolution of stress fields within grain boundaries that accounts for the effects of coherent dislocation transmission is introduced. Starting with a disclination-based construct of GBs at minimum energy (equilibrium), the model describes the evolution of the GB stress field to a state characteristic of excess energy (non-equilibrium) due to the incorporation of residual Burgers vector content following sequential slip transmission events. Several essential features of this model are verified via molecular dynamics simulations of lattice dislocation absorption. Second, this model is implemented into a discrete dislocation dynamics (DDD) code and simulations are performed to understand the influence of Non-equilibrium GB stress fields, conditioned by the slip transmission history, on subsequent dislocation transmission. DDD simulations reveal that the critical resolved stress necessary for slip propagation can be reduced by with continued absorption of residual dislocation content. Moreover, DDD simulations prove that a comprehensive consideration of both the binding and driving stresses, and the evolution of the transmission configuration is necessary to quantify the influence of the mechanical state of the GB on slip propagation. In general, this work provides important insights into the role of GB structure evolution, conditioned by prior deformation history, on intergranular plasticity.

42 ENGINEERING↗

Feedback control design for non-inductively sustained scenarios in NSTX-U using TRANSP

This paper examines a method for real-time control of non-inductively sustained scenarios in NSTX-U by using TRANSP, a time-dependent integrated modeling code for prediction and interpretive analysis of tokamak experimental data, as a simulator. The actuators considered for control in this work are the six neutral beam sources and the plasma boundary shape. To understand the response of the plasma current, stored energy, and central safety factor to these actuators and to enable systematic design of control algorithms, simulations were run in which the actuators were modulated and a linearized dynamic response model was generated. A multi-variable model-based control scheme that accounts for the coupling and slow dynamics of the system while mitigating the effect of actuator limitations was designed and simulated. Simulations show that modest changes in the outer gap and heating power can improve the response time of the system, reject perturbations, and track target values of the controlled values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SU-E-T-438: Commissioning of An In-Vivo Quality Assurance Method Using the Electronic Portal Imaging Device

Purpose: Patient specific pre-treatment quality assurance (QA) using arrays of detectors or film have been the standard approach to assure the correct treatment is delivered to the patient. This QA approach is expensive, labor intensive and does not guarantee or document that all remaining fractions were treated properly. The purpose of this abstract is to commission and evaluate the performance of a commercially available in-vivo QA software using the electronic portal imaging device (EPID) to record the daily treatments. Methods: The platform EPIgray V2.0.2 (Dosisoft), which machine model compares ratios of TMR with EPID signal to predict dose was commissioned for an Artiste (Siemens Oncology Care Systems) and a Truebeam (Varian medical systems) linear accelerator following the given instructions. The systems were then tested on three different phantoms (homogeneous stack of solid water, anthropomorphic head and pelvis) and on a library of patient cases. Simple and complex fields were delivered at different exposures and for different gantry angles. The effects of the table attenuation and the EPID sagging were evaluated. Gamma analysis of the measured dose was compared to the predicted dose for complex clinical IMRT cases. Results: Commissioning of the EPIgray system for two photon energies took 8 hours. The difference between the dose planned and the dose measured with EPIgray was better than 3% for all phantom scenarios tested. Preliminary results on patients demonstrate an accuracy of 5% is achievable in high dose regions for both 3DCRT and IMRT. Large discrepancies (>5%) were observed due to metallic structures or air cavities and in low dose areas. Flat panel sagging was visible and accounted for in the EPIgray model. Conclusion: The accuracy achieved by EPIgray is sufficient to document the safe delivery of complex IMRT treatments. Future work will evaluate EPIgray for VMAT and high dose rate deliveries. This work is supported by Dosisoft, Cachan, France.

60 APPLIED LIFE SCIENCES↗

Confidence intervals with a priori parameter bounds

We review the methods of constructing confidence intervals that account for a priori information about one-sided constraints on the parameter being estimated. We show that the so-called method of sensitivity limit yields a correct solution of the problem. Derived are the solutions for the cases of a continuous distribution with non-negative estimated parameter and a discrete distribution, specifically a Poisson process with background. For both cases, the best upper limit is constructed that accounts for the a priori information. A table is provided with the confidence intervals for the parameter of Poisson distribution that correctly accounts for the information on the known value of the background along with the software for calculating the confidence intervals for any confidence levels and magnitudes of the background (the software is freely available for download via Internet).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Fire Protection Engineering Design Brief Template: Hydrogen Refueling Station

Building a hydrogen infrastructure system is critical to supporting the development of alternate- fuel vehicles. This report provides a methodology for implementing a performance-based design of an outdoor hydrogen refueling station that does not meet specific prescriptive requirements in NFPA 2, The Hydrogen Technologies Code . Performance-based designs are a code-compliant alternative to meeting prescriptive requirements. Compliance is demonstrated by comparing a prescriptive-based fueling station design with a performance-based design approach using Quantitative Risk Assessment (QRA) methods and hydrogen risk assessment tools. This template utilizes the Sandia-developed QRA tool, Hydrogen Risk Analysis Models (HyRAM), which combines reduced-order deterministic models that characterize hydrogen release and flame behavior with probabilistic risk models to quantify risk values. Each project is unique and this template is not intended to account for site-specific characteristics. Instead, example content and a methodology are provided for a representative hydrogen refueling site which can be built upon for new hydrogen applications.

08 HYDROGEN↗

Crystal Plasticity Model of Reactor Pressure Vessel Embrittlement in GRIZZLY

The integrity of reactor pressure vessels (RPVs) is of utmost importance to ensure safe operation of nuclear reactors under extended lifetime. Microstructure-scale models at various length and time scales, coupled concurrently or through homogenization methods, can play a crucial role in understanding and quantifying irradiation-induced defect production, growth and their influence on mechanical behavior of RPV steels. A multi-scale approach, involving atomistic, meso- and engineering-scale models, is currently being pursued within the GRIZZLY project to understand and quantify irradiation-induced embrittlement of RPV steels. Within this framework, a dislocation-density based crystal plasticity model has been developed in GRIZZLY that captures the effect of irradiation-induced defects on the flow stress behavior and is presented in this report. The present formulation accounts for the interaction between self-interstitial loops and matrix dislocations. The model predictions have been validated with experiments and dislocation dynamics simulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

MCRadar: A Monte Carlo Solver for Cloud and Precipitation Radar

Multiple scattering produces anomalous echoes in observed radar profiles that cannot be explained by other phenomena. These effects are most obvious for spaceborne platforms, at shorter wavelengths, and in convection, as multiple scattering is governed by antenna beam width, optical depth, and albe do; however, multiple scattering has been observed in a range of precipitating conditions. To account for the effects of multiple scattering, Monte Carlo integration is employed in a flexible framework to enable arbitrary radar configurations with finite Gaussian beams for three-dimensional atmospheric scenarios. The three-dimensional nature of the model, coupled with the finite antenna response, also allows for consideration of nonuniform beam filling effects that often coincide with multiple scattering. Examples of airborne and spaceborne radars at various wavelengths are used to illustrate the effects of multiple scattering and nonuniform beam filling. The MCRadar code is currently available in the development version of ARTS.

Adams, Ian Stuart↗

Study on Modeling of Nuclear Fuel End Drop during a Hypothetical Accident

It is well known that shell elements have been used to model the cladding of single nuclear fuel pins. These have been used, for example, to analyze the response of a fuel assembly by assuming that all fuel pins in a fuel assembly deform in the same pattern. It is increasingly significant to model the entire fuel assembly in order to capture more realistically the behavior of each individual nuclear fuel during a hypothetical accident in a manner that can account for interactions among single pins. However, modeling each individual fuel pin in an entire fuel assembly using shell elements will result in models with large amount of elements and degrees of freedom that take correspondingly larger amount of computational times to execute. Therefore, modeling the nuclear fuel assembly with beam elements may become an option. This paper presents the results of a set of analysis with the LS-DYNA finite element software that investigated the differences and similarities of modeling nuclear fuel assemblies using either shell elements or beam elements. This study provides additional insights on the appropriateness of using beam element models of the fuel pins in the simulation of a hypothetical accident. In this study, two single pin models, a single pin model using shell elements and a single beam model using beam elements, are constructed to analyze end drop scenarios with the LS-DYNA code. Each model is then analyzed for a 19 ft end drop scenario (linear behavior) and 30 ft end drop scenario (non-linear behavior). Comparisons are made between the results obtained with the single pin shell model and with the single pin beam model to study the differences and similarities resulting from the two modeling approaches. In this study, a B and W 15x15 PWR fuel is analyzed for both 30 ft and 19 ft end hypothetical accident drops using both single pin finite element model with shell elements and beam elements for the fuel cladding. It is found that for the 19 ft end drop, the difference in the Max/Min total axial strains is about 12% between the two models. For the 30 ft end drop, the maximum total axial strains calculated with both models is similar, and a difference of about 12% is observed for the minimum total axial strains calculated by the two models. Advantages of the single pin shell model are that it permits consideration of internal pressures and may capture local buckling. Therefore, use of beam elements to model the fuel cladding needs to take into account the possible significance of local buckling and internal pressurization for the problem being analyzed. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Study of electron transport in a Hall thruster by axial–radial fully kinetic particle simulation

Electron transport across a magnetic field in a magnetic-layer-type Hall thruster was numerically investigated for the future predictive modeling of Hall thrusters. The discharge of a 1-kW-class magnetic-layer-type Hall thruster designed for high-specific-impulse operation was modeled using an r-z two-dimensional fully kinetic particle code with and without artificial electron-diffusion models. The thruster performance results showed that both electron transport models captured the experimental result within discrepancies less than 20% in thrust and discharge current for all the simulated operation conditions. The electron cross-field transport mechanism of the so-called anomalous diffusion was self-consistently observed in the simulation without artificial diffusion models; the effective electron mobility was two orders of magnitude higher than the value obtained using the classical diffusion theory. To account for the self-consistently observed anomalous transport, the oscillation of plasma properties was speculated. It was suggested that the enhanced random-walk diffusion due to the velocity oscillation of low-frequency electron flow could explain the observed anomalous diffusion within an order of magnitude. The dominant oscillation mode of the electron flow velocity was found to be 20 kHz, which was coupled to electrostatic oscillation excited by global ionization instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Soft X-ray emission from laser-irradiated gold foils

Here, this paper reports measurements of soft-x-ray emission from gold foils irradiated by 6 ns laser pulses, and analysis and simulations of the observations. These foils can be used as x-ray sources to drive a wide range of experiments. A multichannel, photodiode array measured the time-resolved, soft-x-ray emission. A soft-x-ray framing camera imaged the emission in selected energy bands. Foil thicknesses were from 0.5 to 1.5 μm. The imaging data show that the region emitting soft x-rays grows throughout the laser drive, on both the front and rear surfaces. Analysis of the emitted radiation flux from the rear surface, taking the time-dependent spot size into account, showed that the peak effective temperature of 0.5-μm-thick foils is near 88 eV, while that of 0.75-μm-thick foils is near 78 eV. A Monte Carlo method was used to evaluate the component of the uncertainty in the effective temperature introduced by variations in signal voltages and by uncertainty in the size of the emitting spot. This was found to be near ±2 eV in most cases. Simple theoretical considerations explain the main features of the observations. The Florence code, working with atomic physics from sophisticated models, proved able to reproduce the main features of the observed spectra with 1D simulations in which the laser energy flux was adjusted on the basis of the observed lateral spreading of energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Flight experiment of thermal energy storage

Thermal energy storage (TES) enables a solar dynamic system to deliver constant electric power through periods of sun and shade. Brayton and Stirling power systems under current considerations for missions in the near future require working fluid temperatures in the 1100 to 1300+ K range. TES materials that meet these requirements fall into the fluoride family of salts. These salts store energy as a heat of fusion, thereby transferring heat to the fluid at constant temperature during shade. The principal feature of fluorides that must be taken into account is the change in volume that occurs with melting and freezing. Salts shrink as they solidify, a change reaching 30 percent for some salts. The location of voids that form as result of the shrinkage is critical when the solar dynamic system reemerges into the sun. Hot spots can develop in the TES container or the container can become distorted if the melting salt cannot expand elsewhere. Analysis of the transient, two-phase phenomenon is being incorporated into a three-dimensional computer code. The code is capable of analysis under microgravity as well as 1 g. The objective of the flight program is to verify the predictions of the code, particularly of the void location and its effect on containment temperature. The four experimental packages comprising the program will be the first tests of melting and freezing conducted under microgravity. Each test package will be installed in a Getaway Special container to be carried by the shuttle. The package will be self-contained and independent of shuttle operations other than the initial opening of the container lid and the final closing of the lid. Upon the return of the test package from flight, the TES container will be radiographed and finally partitioned to examine the exact location and shape of the void. Visual inspection of the void and the temperature data during flight will constitute the bases for code verification.

Namkoong, David↗

Comparison of FDMA and CDMA for second generation land-mobile satellite communications

Code Division Multiple Access (CDMA) and Frequency Division Multiple Access (FDMA) (both analog and digital) systems capacities are compared on the basis of identical link availabilities and physical propagation models. Parameters are optimized for a bandwidth limited, multibeam environment. For CDMA, the benefits of voice activated carriers, antenna discrimination, polarization reuse, return link power control and multipath suppression are included in the analysis. For FDMA, the advantages of bandwidth efficient modulation/coding combinations, voice activated carriers, polarization reuse, beam placement, and frequency staggering were taken into account.

Yongacoglu, A.↗

Numerical Modeling of HgCdTe Solidification: Effects of Phase Diagram, Double-Diffusion Convection and Microgravity Level

Melt convection, along with species diffusion and segregation on the solidification interface are the primary factors responsible for species redistribution during HgCdTe crystal growth from the melt. As no direct information about convection velocity is available, numerical modeling is a logical approach to estimate convection. Furthermore influence of microgravity level, double-diffusion and material properties should be taken into account. In the present study, HgCdTe is considered as a binary alloy with melting temperature available from a phase diagram. The numerical model of convection and solidification of binary alloy is based on the general equations of heat and mass transfer in two-dimensional region. Mathematical modeling of binary alloy solidification is still a challenging numericial problem. A Rigorous mathematical approach to this problem is available only when convection is not considered at all. The proposed numerical model was developed using the finite element code FIDAP. In the present study, the numerical model is used to consider thermal, solutal convection and a double diffusion source of mass transport.

Bune, Andris V.↗

CMC Property Variability and Life Prediction Methods for Turbine Engine Component Application

The ever increasing need for lower density and higher temperature-capable materials for aircraft engines has led to the development of Ceramic Matrix Composites (CMCs). Today's aircraft engines operate with >3000"F gas temperatures at the entrance to the turbine section, but unless heavily cooled, metallic components cannot operate above approx.2000 F. CMCs attempt to push component capability to nearly 2700 F with much less cooling, which can help improve engine efficiency and performance in terms of better fuel efficiency, higher thrust, and reduced emissions. The NASA Glenn Research Center has been researching the benefits of the SiC/SiC CMC for engine applications. A CMC is made up of a matrix material, fibers, and an interphase, which is a protective coating over the fibers. There are several methods or architectures in which the orientation of the fibers can be manipulated to achieve a particular material property objective as well as a particular component geometric shape and size. The required shape manipulation can be a limiting factor in the design and performance of the component if there is a lack of bending capability of the fiber as making the fiber more flexible typically sacrifices strength and other fiber properties. Various analysis codes are available (pcGINA, CEMCAN) that can predict the effective Young's Moduli, thermal conductivities, coefficients of thermal expansion (CTE), and various other properties of a CMC. There are also various analysis codes (NASAlife) that can be used to predict the life of CMCs under expected engine service conditions. The objective of this summer study is to utilize and optimize these codes for examining the tradeoffs between CMC properties and the complex fiber architectures that will be needed for several different component designs. For example, for the pcGINA code, there are six variations of architecture available. Depending on which architecture is analyzed, the user is able to specify the fiber tow size, tow spacing, weave parameter, and angle of orientation of fibers. By holding the volume fraction of the fibers constant, variations in tow spacing can be explored for different architectures. The CMC material properties are usually calculated assuming the component is manufactured perfectly. However, this is typically not the case so that a quantification of the material property variability is needed to account for processing and/or manufacturing imperfections. The overall inputs and outputs are presented using a regression software to rapidly investigate the tradeoffs associated with fiber architecture, material properties, and ultimately cost. This information is then propagated through lifing models and Larson-Miller data to assess timehemperature-dependent CMC strength. In addition, a first order cost estimation will be quantified from a current qualitative perspective. This cost estimation includes the manufacturing challenges, such as tooling, as well as the component cost for a particular application. Ultimately, a cost to performance ratio should be established that compares the effectiveness of CMCs to their current rival, nickel superalloys.

Cheplak, Matthew L.↗

Toward the development of intrafraction tumor deformation tracking using a dynamic multi-leaf collimator

Purpose: Intrafraction deformation limits targeting accuracy in radiotherapy. Studies show tumor deformation of over 10 mm for both single tumor deformation and system deformation (due to differential motion between primary tumors and involved lymph nodes). Such deformation cannot be adapted to with current radiotherapy methods. The objective of this study was to develop and experimentally investigate the ability of a dynamic multi-leaf collimator (DMLC) tracking system to account for tumor deformation. Methods: To compensate for tumor deformation, the DMLC tracking strategy is to warp the planned beam aperture directly to conform to the new tumor shape based on real time tumor deformation input. Two deformable phantoms that correspond to a single tumor and a tumor system were developed. The planar deformations derived from the phantom images in beam's eye view were used to guide the aperture warping. An in-house deformable image registration software was developed to automatically trigger the registration once new target image was acquired and send the computed deformation to the DMLC tracking software. Because the registration speed is not fast enough to implement the experiment in real-time manner, the phantom deformation only proceeded to the next position until registration of the current deformation position was completed. The deformation tracking accuracy was evaluated by a geometric target coverage metric defined as the sum of the area incorrectly outside and inside the ideal aperture. The individual contributions from the deformable registration algorithm and the finite leaf width to the tracking uncertainty were analyzed. Clinical proof-of-principle experiment of deformation tracking using previously acquired MR images of a lung cancer patient was implemented to represent the MRI-Linac environment. Intensity-modulated radiation therapy (IMRT) treatment delivered with enabled deformation tracking was simulated and demonstrated. Results: The first experimental investigation of adapting to tumor deformation has been performed using simple deformable phantoms. For the single tumor deformation, the A{sub u}+A{sub o} was reduced over 56% when deformation was larger than 2 mm. Overall, the total improvement was 82%. For the tumor system deformation, the A{sub u}+A{sub o} reductions were all above 75% and the total A{sub u}+A{sub o} improvement was 86%. Similar coverage improvement was also found in simulating deformation tracking during IMRT delivery. The deformable image registration algorithm was identified as the dominant contributor to the tracking error rather than the finite leaf width. The discrepancy between the warped beam shape and the ideal beam shape due to the deformable registration was observed to be partially compensated during leaf fitting due to the finite leaf width. The clinical proof-of-principle experiment demonstrated the feasibility of intrafraction deformable tracking for clinical scenarios. Conclusions: For the first time, we developed and demonstrated an experimental system that is capable of adapting the MLC aperture to account for tumor deformation. This work provides a potentially widely available management method to effectively account for intrafractional tumor deformation. This proof-of-principle study is the first experimental step toward the development of an image-guided radiotherapy system to treat deforming tumors in real-time.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

On-the-fly Doppler Broadening with Probability Table Interpolation for Unresolved Resonance Region in RMC

With the rapid development of computational power and parallel algorithms, Monte Carlo method has been widely investigated and used in neutron transport simulations due to its advantages in geometry modelling and usage of continuous energy point-wise cross sections. Nowadays the road has been paved to couple Monte Carlo codes with thermal-hydraulics codes to make multi-physics calculations. One of the most important aspects of thermal-hydraulics feedback is the detailed temperature distribution, which demands temperature dependent cross sections in neutron transport simulations. The conventional pre-generated cross section libraries at each temperature point are no longer feasible for detailed temperature distributions because of tremendous memory footprint. Hence, different on-the-fly Doppler broadening methods have been proposed for resolved resonance energy region and thermal energy region, such as the interpolation, regress model fitting, multipole representation (MPR), target motion sampling (TMS) and different integration algorithms based on Doppler broadening equation. Compared to resolved resonance energy region and thermal energy region, the representation of cross sections in unresolved resonance region is slightly different. As a result of the limitation of experiment measurements, the overlapped resonances cannot be distinguished from each other and the statistic average resonance parameters are given in the evaluated nuclear data file. In early applications, infinite-dilute pointwise cross sections were obtained from average parameters and were directly used in neutron simulations. However, it was proven that these infinite-dilute pointwise cross sections cannot capture the energy self-shielding effect, thus leading to an overestimate of reaction rates and usually an underestimate of k{sub eff} owing to the increased absorption rates and increased scattering probability to energy below unresolved resonance region where the capture-to-fission ratio is typically higher. This phenomenon is of great importance especially for problems where there is a considerable flux in unresolved resonance region. Then the probability table method was put forward by Levitt to take energy self-shielding effect into account, which provides probability distribution function of cross sections in unresolved resonance region instead of specific cross section values. The probability distribution function is presented in the form of tables and the mean values are equal to the infinite dilution cross sections. To obtain the temperature dependent cross sections in unresolved energy region, there are two different methods. One is to calculate cross sections directly from average resonance parameters and their statistical distributions , which is effective and precise but time-consuming. The other is to interpolate between pre-generated probability table files, which is much more efficient. In this paper, on-the-fly Doppler broadening of cross sections in unresolved resonance regions based on probability tables interpolation was implemented in RMC codes as a supplement to on-the-fly Doppler broadening for resolved resonance energy region and thermal region, which is based on the integration of Doppler broadening equation. The results of criticality calculation and depletion calculation show that on-the-fly Doppler broadening for the whole energy region is accurate and efficient. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗