Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo Code”

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

(U) LNK3DNT Capability in SENSMG Extended

The SENSMG multigroup neutron sensitivity code can read materials and geometry from redoin and lnk3dnt files written by the PARTISN multigroup discrete ordinates code or the MCNP6 Monte Carlo code. This report describes the user interface for this capability and the code verification that was done. Reference 5 described an initial version of this capability. Two important extensions have been made: 1) Reaction-rate ratios are now allowed as responses with redoin/lnk3dnt inputs; and 2) Sensitivities of reaction-rate ratios to mass densities and interface locations are now calculated (for all problems, not just redoin/lnk3dnt problems).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

(U) LNK3DNT Capability in SENSMG Tested with Transport Correction

The SENSMG multigroup neutron sensitivity code can read materials and geometry from redoin and lnk3dnt files written by the PARTISN multigroup discrete ordinates code or the MCNP6 Monte Carlo code. Reference 5 described an initial version of this capability. Reference 6 introduced two important extensions: 1) Reaction-rate ratios are now allowed as responses with redoin/lnk3dnt inputs; and 2) Sensitivities of reaction-rate ratios to mass densities and interface locations are now calculated (for all problems, not just redoin/lnk3dnt problems). In this report, we document the testing of the redoin/lnk3dnt capability when used with transport correction. PARTISN’s transport corrections are described in Ref. 3, and Refs. 7 and 8 discuss the calculation of cross-section sensitivities when transport corrections are applied. References 7 and 8 verify SENSMG for problems with transport correction using central differences. This report verifies the capability with redoin/lnk3dnt inputs by comparison with regular SENSMG inputs. Unlike Refs. 5 and 6, in this report we use a spherical test problem because the sensitivity of the response with respect to the highest-order scattering is zero in cylinders due to the use of flux moments in the inner products.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Modeling the Ex-core Detector Response into a Production Option in MPACT

In the past few years, the Consortium for Advanced Simulation of Light Water Reactors (CASL) has overseen a successful activity to extend the radiation transport capability from the reactor in-core geometry to the ex-core geometry by coupling the deterministic neutronics code MPACT, with the Monte Carlo code Shift. However, transient problems such as RIAs require an efficient prediction of the ex-core detector response to determine the timeframe for a reactor trip. Therefore, in addition to the high-fidelity MPACT-Shift coupled calculation, separate efforts have been made to explore simplified methods for ex-core detector response with MPACT standalone calculation. In a previous milestone L3_RTM.MCH.P17.06, several simplified methods for ex-core detector response were developed in MPACT. In particular, the double kernel method is based on the default MPACT core geometry, so extensive modeling or transport calculation with full reflector and ex-core details are not needed in this method. Since the 1-D diffusion and neutron streaming kernels are used for extended ex-core calculations, the extra computational efforts are very small. The work implemented in the current milestone provides a fully integrated option for MPACT to approximately calculate the ex-core detector response using the double kernel method. To overcome the difficulty of obtaining the consistent diffusion coefficients (these were obtained through separate CMFD calculation and edits), the 1-D diffusion solver is replaced by a 1-D Sn transport solver to compute the neutron flux at the vessel outer surface. The new VERA inputs for ex-core detector are processed in MPACT. The detector signals are linked to the SCRAM logic as an option to trip the core during transient calculation. Two sets of 3-D quarter core problems are run with MPACT standalone and VeraShift to verify the radial and axial effects on detector response by varying the moderator density and control rod position. For the cases with moderator density changes, the detector responses from MPACT are within 2 sigma error of VeraShift. For the cases at different control rod positions, the axial offsets computed from MPACT are within 2% error of VeraShift results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A highly optimized vectorized code for Monte Carlo simulations of SU(3) lattice gauge theories

New methods are introduced for improving the performance of the vectorized Monte Carlo SU(3) lattice gauge theory algorithm using the CDC CYBER 205. Structure, algorithm and programming considerations are discussed. The performance achieved for a 16(4) lattice on a 2-pipe system may be phrased in terms of the link update time or overall MFLOPS rates. For 32-bit arithmetic, it is 36.3 microsecond/link for 8 hits per iteration (40.9 microsecond for 10 hits) or 101.5 MFLOPS.

Barkai, D.↗

Analysis of MERCI decay heat measurement for PWR UO{sub 2} fuel rod

Decay heat measurements, called the MERCI experiment, were conducted at Commissariat a l'Energie Atomique (CEA)/Saclay to characterize accurately residual power at short cooling time and verify its prediction by decay code and nuclear data. The MOSAIC calorimeter, developed and patented by CEA/Grenoble (DTN/SE2T), enables measurement of the decay heat released by a pressurized water reactor (PWR) fuel rod sample between 200 and 4 W within a precision of 1%. The MERCI experiment included three phases. At first, a UO{sub 2} fuel rod sample was irradiated in the CEA/Saclay experimental reactor OSIRIS. The burnup achieved at the end of irradiation was similar to 3.5 GWd/tonne. The second phase was the transfer of the fuel rod sample from its irradiation location to a hot cell, to be inserted inside the MOSAIC calorimeter. It took 26 min to carry out the transfer. Finally, decay heat released by the PWR sample was measured from 27 min to 42 days after shutdown. Post irradiation examinations were performed to measure concentrations of some heavy nuclei (U, Pu) and fission products (Cs, Nd). The decay heat was predicted using a calculation scheme based on the PEPIN2 depletion code, the TRIPOLI-4 Monte Carlo code, and the JEFF3.1.1 nuclear data file. The MERCI experiment analysis shows that the discrepancy between the calculated and the experimental decay heat values is included between -10% at 27 min and +6% at 12 h, 30 min otter shutdown. From 4 up to 42 days of cooling time, the difference between calculation and measurement is about ± 1%, i.e., experimental uncertainty. The MERCI experiment represents a significant contribution for code validation; the time range above 10{sup 5} s has not been validated previously. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Geant4 and beyond: recent progress in precision physics modeling - Paper 100

This extended abstract briefly summarizes ongoing research activity on the evaluation and experimental validation of physics methods for photon and electron transport. The analysis includes physics models currently implemented in Geant4 as well as modeling methods used in other Monte Carlo codes, or not yet considered in general purpose Monte Carlo simulation systems. The validation of simulation models is performed with the support of rigorous statistical methods, which involve goodness-of-fit tests followed by categorical analysis. All results are quantitative, and are fully documented. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

Commissioning of a Varian Clinac iX 6 MV photon beam using Monte Carlo simulation

Monte Carlo modelling of a linear accelerator is the first and most important step in Monte Carlo dose calculations in radiotherapy. Monte Carlo is considered today to be the most accurate and detailed calculation method in different fields of medical physics. In this research, we developed a photon beam model for Varian Clinac iX 6 MV equipped with MilleniumMLC120 for dose calculation purposes using BEAMnrc/DOSXYZnrc Monte Carlo system based on the underlying EGSnrc particle transport code. Monte Carlo simulation for this commissioning head LINAC divided in two stages are design head Linac model using BEAMnrc, characterize this model using BEAMDP and analyze the difference between simulation and measurement data using DOSXYZnrc. In the first step, to reduce simulation time, a virtual treatment head LINAC was built in two parts (patient-dependent component and patient-independent component). The incident electron energy varied 6.1 MeV, 6.2 MeV and 6.3 MeV, 6.4 MeV, and 6.6 MeV and the FWHM (full width at half maximum) of source is 1 mm. Phase-space file from the virtual model characterized using BEAMDP. The results of MC calculations using DOSXYZnrc in water phantom are percent depth doses (PDDs) and beam profiles at depths 10 cm were compared with measurements. This process has been completed if the dose difference of measured and calculated relative depth-dose data along the central-axis and dose profile at depths 10 cm is ≤ 5%. The effect of beam width on percentage depth doses and beam profiles was studied. Results of the virtual model were in close agreement with measurements in incident energy electron 6.4 MeV. Our results showed that photon beam width could be tuned using large field beam profile at the depth of maximum dose. The Monte Carlo model developed in this study accurately represents the Varian Clinac iX with millennium MLC 120 leaf and can be used for reliable patient dose calculations. In this commissioning process, the good criteria of dose difference in PDD and dose profiles were achieve using incident electron energy 6.4 MeV.

61 RADIATION PROTECTION AND DOSIMETRY↗

TNSL Overview

Thermal neutron scattering law (TNSL) data describe low-energy neutrons scattering off of bound materials, and can have a significant impact on modeling any system with slow neutrons, including nuclear reactors. Previous work to introduce TNSL data to neutron transport codes at LLNL focused on COG and TART [1], with the limitation that these codes require highly specialized data processing and formatting. We have recently increased efforts to process TNSL data with the central LLNL nuclear data processing code FUDGE, to be stored in the generalized nuclear database structure (GNDS) for use in any general transport code with the ability to read GNDS data. The first step in this effort is to verify that the TNSL processing with FUDGE yields results comparable to results obtained using the LANL nuclear data processing code NJOY. The next step is to verify the transport of thermal neutrons in Mercury (a Monte Carlo code) and Ardra (a deterministic code) against one another, as well as against the LANL Monte Carlo neutron transport code MCNP. This verification step has not been completed, due to a number of discrepancies between results obtained using differently processed data. There is ongoing effort to understand differences between FUDGE and NJOY. Finally, we map out our current capability to validate TNSL data against benchmark systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Application of the FW-CADIS variance reduction method to calculate a precise N-flux distribution for the FRJ-2 research reactor - 15603

In order to avoid unnecessary radiation exposure and to minimize the radioactive waste and respectively cost by the decommissioning of nuclear facilities, the radiological characteristics and the activity atlases are required. Therefore, a significant step is the accurate determination of the neutron fluence distribution resulting from the operation of a reactor. For this purpose, computer codes based on the Monte-Carlo method are employed because of the flexibility in the modelling of the complex structures and physical processes. In the scope of the research project founded by the German Federal Ministry of Education and Research (BMBF), and implemented by the Institute of Nuclear Engineering and Technology Transfer (NET) -RWTH Aachen, sophisticated simulation tools are developed, which have the capability to cover the range of demands for activity and dose rate distribution using the example of the German research reactor FRJ-2. Due to the complex geometry of FRJ-2, which consists of different structures and configuration, the MCNP5 Monte Carlo Code is the optimal method for a precise neutron transport simulation. This work presents a method to achieve a precise N-Flux distribution for the FRJ-2 research reactor using hybrid deterministic Monte Carlo method. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The History of Monte Carlo and MCNP at Los Alamos [Slides]

The Monte Carlo method for radiation particle transport has its origins at LANL dating back to the 1940’s. The creators of these methods were Drs. Stanislaw Ulam, John von Neumann, Robert Richtmyer, and Nicholas Metropolis. Monte Carlo methods for particle transport have been driving computational developments since the beginning of modern computers; this continues today. In the 1950’s and 1960’s, these new methods were organized into a series of special-purpose Monte Carlo codes, including MCS, MCN, MCP, and MCG. These codes were able to transport neutrons and photons for specialized LANL applications. In 1977, these separate codes were combined to create the LANL Monte Carlo N-Particle (MCNP) radiation particle transport code. In 1983, MCNP3 was released for public distribution to the Radiation Safety Information Computational Center (RSICC). MCNP6.3 will be released late in 2021. Each year, LANL has ~100 MCNP new users and RSICC ~1250 new licenses distributed. This talk will review the Los Alamos history of the development of the modern Monte Carlo method and MCNP.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Study on neutronics calculation method of SCWR with annular fuel assembly - 26804

Many Supercritical Water Reactor designs have been proposed all over the world. Annular fuel was employed in SCWR by Xi'an JiaoTong University to achieve lower fuel center temperature and lower cladding surface temperature. Different analysis codes have been developed for these designs. In this article, a neutronics calculation code for SCWR with annular fuel is developed. This code is compared with Monte-Carlo code doing both assembly and core calculation. The results showed that the spectra difference of PWRs and SCWRs is the main reason for assembly calculation errors and conventional assembly homogenization method is not suitable for SCWRs. Thus the group-collapsing and spatial-averaging method were studied to improve the accuracy. The few-group structure which contains 14 energy groups is used. The 'super-assembly' model is employed to consider effect of reflector. The results showed that the accuracy of effective multiplication factor and radial power was improved. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗