Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “source 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 1,567 records · Page 87

Branson: A Mini-App for Studying Parallel IMC, Version 1.0

This code solves the gray thermal radiative transfer (TRT) equations in parallel using simple opacities and Cartesian meshes. Although Branson solves the TRT equations it is not designed to model radiation transport: Branson contains simple physics and does not have a multigroup treatment, nor can it use physical material data. The opacities have are simple polynomials in temperature there is a limited ability to specify complex geometries and sources. Branson was designed only to capture the computational demands of production IMC codes, especially in large parallel runs. It was also intended to foster collaboration with vendors, universities and other DOE partners. Branson is similar in character to the neutron transport proxy-app Quicksilver from LLNL, which was recently open-sourced.

Long, Alex↗

Quantification of LEU Holdup using gamma ray imaging and inverse transport solver

Holdup is the residual amount of special nuclear material (SNM) remaining in a processing facility after the bulk materials have been cleaned out. In commercial uranium processing facilities, quantification of holdup is a major challenge because of the highly variable shapes and sizes of the deposits. Any method that attempts to generalize and calibrate deposit shapes in order to quantify holdup will be prone to high uncertainties. Uncertainties on the order of ±50% are typical in holdup results. In international safeguards applications, a ±50% uncertainty can result in a large amount of material unaccounted for (MUF) thereby increasing the difficulty of detecting material diversion and facility misuse. An imaging-based methodology has been developed with the objective of significantly reducing this uncertainty by using the true deposit shape, instead of relying on oversimplified geometric assumptions. The project is a collaboration between ORNL, Y-12, and the University of Tennessee, Knoxville, TN. Uranium sources of known masses were measured using the Germanium Gamma-ray Imager (GeGI), a high-resolution imaging spectrometer, creating a pixelated map for each spectral bin. Two different gamma imaging methods are employed in this work: coded aperture imaging and Compton imaging. A validated MonteCarlo model of the detector has been developed using the GEANT4 code for determining the intrinsic response of the detector, its enclosure, and the coded aperture mask. An inverse transport solver based on the Markov Chain Monte-Carlo approach known as Differential Evolution Adaptive Metropolis (DREAM) is employed to use the measurement data from the image pixels (coded aperture or Compton) to solve for the mass of 235 U in the deposit. A reliable method based on the DREAM solver has been developed to flag the infinite thickness condition of a uranium deposit. The project team is working towards improving the image reconstruction for Compton imaging so that a better localization of the source can be achieved. Besides treating the coded aperture and Compton imaging methods independently, the project is also evaluating a combined method that uses the Compton scatter data from a coded aperture measurement. GEANT4 simulations are being performed to evaluate the combined approach. The impact on the DREAM optimization as the source thickness progressively approaches infinite thickness is being evaluated. A number of uranium sources available at ORNL have been measured, and the DREAM results have been tested and validated for the coded aperture imaging. A similar effort will be carried out to validate the Compton based method once the development of algorithms for better localization are complete. The imaging based quantification is very amenable to unattended monitoring of holdup accumulation at key measurement points. A proof of concept measurement has been completed to demonstrate this capability The current work used the high energy resolution imager GeGI. However, the approach and methodologies are applicable to other imagers such as the cadmium zin telluride (CZT) based imager manufactured by H3D, Inc.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Design and performance testing of a novel far field Gamma system to assay radioactive waste in 400 and 200 litre drums - 15339

A Far Field Gamma ray Measurement system has been developed for measuring and assaying both 400 and 200 litre waste drums with a wide range of activity at nuclear power stations. These include radioactive waste drums with surface dose-rates of less than 2 mSv/h (200 mrem/h) and, at the higher range of activity, 400 litre drums with surface dose-rates of up to 50 mSv/h (5 rem/h). The assay system incorporates an electro-mechanically cooled high purity Germanium profile detector and advanced digital spectroscopic analysis electronics. A novel optimised conical-trapezoidal collimator has been incorporated to reduce the effect of the radioactive background in addition to the conventional detector lead shield with tin and copper graded lining to reduce the effect of lead X-rays. Two features have been incorporated to reduce detector dead time and maintain detector resolution when measuring high activity drums. The first is the inclusion of two automated tungsten filters of different thickness as an integral part of the detector collimator. The filters are used in conjunction with an automated rail system on which the detector platform is mounted. The detector can be positioned at different distances from the drum surface in order to reduce detector dead time. The deployment of the filters and the adjustment of the drum - detector position are automatically controlled based on a user defined dose-rate limit table. The data is provided by two Geiger-Muller dose-rate detectors, which are used to measure the drum surface dose-rate. This data is subsequently used to adjust the detector position and in cases of higher dose-rate deploy one or other of the tungsten filters. The drum rotation platform incorporates a load cell to determine drum weight and hence the drum density. This information is used to calculate an attenuation correction based on the assumption of uniform drum density and uniform distribution of activity within the drum. Where the waste drums have a more complicated and known regular internal structure, such as a small radioactive region surrounded by an annulus of shielding material, the analysis algorithms incorporated into the far field geometry and attenuation correction spectroscopic data analysis code are able to make a more accurate determination of drum activity than would be possible assuming a completely uniform drum. The detector energy and efficiency calibration is achieved using a point source with multiple gamma ray energy peaks. The response of the high purity Germanium detector is also modelled using the Monte Carlo Neutron Photon code and the model is benchmarked using the measured point source calibration data. The performance of the instrument has been determined for a range of uniform density matrix drums containing a set of Eu-152 line or rod sources located in re-entrant tubes positioned on an equal volume basis. When rotated, the test drums with volume distributed line sources simulate waste drums of different density each with a uniformly distributed radioactive source. The measurements have been validated using Monte Carlo simulations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

BISON Simulated and Experimental Fission Product Release Comparisons from Reradiated AGR-3/4 Compacts During High Temperature Heating Tests

The fuel performance modeling code BISON was used to predict the release of fission products iodine-131 (131I), xenon-133 (133Xe), and krypton-85 (85Kr) from four re-irradiated AGR-3/4 fuel compacts containing tristructural isotropic (TRISO) coated particles during high-temperature isothermal heating tests. The AGR-3/4 fuel compacts were irradiated in the Advanced Test Reactor (ATR) as part of the third and fourth series of planned experiments to support the Advanced Gas Reactor (AGR) Program. They were subsequently stored and re-irradiated in the Neutron Radiography (NRAD) reactor for approximately five days and then stored for another five to eight days before being subjected to isothermal heating tests in the Fuel Accident Condition Simulation (FACS furnace) for 200 to 300 hours at temperatures between 1000°C and 1600°C to evaluate fission product release at elevated temperatures. New nuclide-specific fission product source term models for the three nuclides of interest were developed using the reactor multiphysics code Griffin and implemented into BISON to support this work. The new source term models were incorporated into coupled compact- and particle-scale BISON simulations, which predict spatially- and temporally-resolved radionuclide generation, radioactive decay, transport, and release throughout the entire irradiation history, including the initial ATR irradiation, NRAD re-irradiations, FACS heating tests, and intermediate periods spent in storage. The experimentally measured fission product release from the heating tests were compared to modeling release predictions calculated by BISON to evaluate how well the code compares to experimental results. Overall, the experimental measured and BISON predicted comparative release results varied but generally agreed to within 5 particle equivalents. Comparative release results identified general observations to take into consideration to help refine future models and reduce uncertainties associated with both the measurement results and predictive results. This includes developing new uranium oxycarbide (UCO) specific kernel diffusivities for the three isotopes examined to more accurately reflect the material properties of the fuel form. Deriving new diffusivities will aid in producing a more informed BISON model

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Beam extraction and high stability operation of high current electron cyclotron resonance proton ion source

A high current electron cyclotron resonance proton ion source is designed and developed for the low energy high intensity proton accelerator at Bhabha Atomic Research Centre. The plasma discharge in the ion source is stabilized by minimizing the reflected microwave power using four stub auto tuner and magnetic field. The optimization of extraction geometry is performed using PBGUNS code by varying the aperture, shape, accelerating gap, and the potential on the electrodes. While operating the source, it was found that the two layered microwave window (6 mm quartz plate and 2 mm boron nitride plate) was damaged (a fine hole was drilled) by the back-streaming electrons after continuous operation of the source for 3 h at beam current of 20–40 mA. The microwave window was then shifted from the line of sight of the back-streaming electrons and located after the water-cooled H-plane bend. In this configuration the stable operation of the high current ion source for several hours is achieved. The ion beam is extracted from the source by biasing plasma electrode, puller electrode, and ground electrode to +10 to +50 kV, −2 to −4 kV, and 0 kV, respectively. The total ion beam current of 30–40 mA is recorded on Faraday cup at 40 keV of beam energy at 600–1000 W of microwave power, 800–1000 G axial magnetic field and (1.2–3.9) × 10{sup −3} mbar of neutral hydrogen gas pressure in the plasma chamber. The dependence of beam current on extraction voltage, microwave power, and gas pressure is investigated in the range of operation of the ion source.

43 PARTICLE ACCELERATORS↗

The Latest Code Development Progress of JSPEC

The JLab Simulation Package on Electron Cooling (JSPEC) is an open source software developed at Jefferson Lab for electron cooling and intrabeam scattering (IBS) simulations. IBS is an important factor that leads to the growth of the beam emittance and hence the reduction of the luminosity in a high density ion collider ring. Electron cooling is an effective measure to overcome the IBS effect. Although JSPEC is initiated to fulfill the simulation needs in JLab Electron Ion Collider project, it can be used as a general design tool for other accelerators. JSPEC provides various models of the ion beam and the electron beam and it calculates the expansion rate and simulates the evolution of the ion beam under the IBS and/or electron cooling effect. In this report, we will give a brief introduction of JSPEC and then present the latest code development progress of JSPEC, including new models, algorithms, and the user interface.

Zhang, Y.↗

Conceptual moderator studies for the Spallation Neutron Source short-pulse second target station

Candidate moderator configurations for a short-pulse second target station (STS) at the Oak Ridge National Laboratory Spallation Neutron Source (SNS) have been identified using a global optimizer framework built around the MCNPX particle transport code. Neutron brightness metrics were selected as the figure-of-merit. We assumed that STS would use one out of six proton pulses produced by an SNS accelerator upgraded to operate at 1.3 GeV proton energy, 2.8 MW power and 60 Hz repetition rate. The simulations indicate that the peak brightness can be increased by a factor of 5 and 2.5 on a per proton pulse basis compared to the SNS first target station for both coupled and decoupled para-hydrogen moderators, respectively. Additional increases by factors of 3 and 2 were demonstrated for coupled and decoupled moderators, respectively, by reducing the area of neutron emission from 100 × 100 mm{sup 2} to 20 × 20 mm{sup 2}. This increase in brightness has the potential to translate to an increase of beam intensity at the instruments’ sample positions even though the total neutron emission of the smaller moderator is less than that of the larger. This is especially true for instruments with small samples (beam dimensions). The increased fluxes in the STS moderators come at accelerated poison and de-coupler burnout and higher radiation-induced material damage rates per unit power, which overall translate into lower moderator lifetimes. A first effort was undertaken to group decoupled moderators into a cluster collectively positioning them at the peak neutron production zone in the target and having a three-port neutron emission scheme that complements that of a cylindrical coupled moderator.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

RectifHydPlus Data Pipeline

The RectifHydPlus Data Pipeline is an open source and fully reproducible data processing pipeline for creating RectifHydPlus—a dataset of historical monthly net electricity generation for all US hydropower plants (>10MW). The pipeline is coded in R, applying tidyverse libraries and code principles, and using the targets data pipeline framework. All data inputs to the RectifHydPlus Data Pipeline are available from public sources. References to all data inputs, as well as instructions for running the RectifHydPlus Data Pipeline, are available on the GitLab code repository: https://code.ornl.gov/turnersw/rectifhydplus

Turner, SeanWilliam Donald [Oak Ridge National Lab↗

RectifHydPlus Data Pipeline v1.1.0

The RectifHydPlus Data Pipeline is an open source and fully reproducible data processing pipeline for creating RectifHydPlus—a dataset of historical monthly net electricity generation for all US hydropower plants (>10MW). The pipeline is coded in R, applying tidyverse libraries and code principles, and using the targets data pipeline framework. All data inputs to the RectifHydPlus Data Pipeline are available from public sources. References to all data inputs, as well as instructions for running the RectifHydPlus Data Pipeline, are available on the GitLab code repository: https://code.ornl.gov/turnersw/rectifhydplus

Turner, SeanWilliam Donald [Oak Ridge National Lab↗

TU-AB-BRC-08: Egs-brachy, a Fast and Versatile Monte Carlo Code for Brachytherapy Applications

Purpose: To introduce egs-brachy, a new, fast, and versatile Monte Carlo code for brachytherapy applications. Methods: egs-brachy is an EGSnrc user-code based on the EGSnrc C++ class library (egs++). Complex phantom, applicator, and source model geometries are built using the egs++ geometry module. egs-brachy uses a tracklength estimator to score collision kerma in voxels. Interaction, spectrum, energy fluence, and phase space scoring are also implemented. Phase space sources and particle recycling may be used to improve simulation efficiency. HDR treatments (e.g. stepping source through dwell positions) can be simulated. Standard brachytherapy seeds, as well as electron and miniature x-ray tube sources are fully modelled. Variance reduction techniques for electron source simulations are implemented (Bremsstrahlung cross section enhancement, uniform Bremsstrahlung splitting, and Russian Roulette). TG-43 parameters of seeds are computed and compared to published values. Example simulations of various treatments are carried out on a single 2.5 GHz Intel Xeon E5-2680 v3 processor core. Results: TG-43 parameters calculated with egs-brachy show excellent agreement with published values. Using a phase space source, 2% average statistical uncertainty in the PTV ((2mm){sup 3} voxels) can be achieved in 10 s for 100 {sup 125}I or {sup 103}Pd seeds in a 36.2 cm{sup 3} prostate PTV, 31 s for 64 {sup 103}Pd seeds in a 64 cm{sup 3} breast PTV, and 56 s for a miniature x-ray tube in a 27 cm{sup 3} breast PTV. Comparable uncertainty is reached in 12 s in a (1 mm){sup 3} water voxel 5 mm away from a COMS 16mm eye plaque with 13 {sup 103}Pd seeds. Conclusion: The accuracy of egs-brachy has been demonstrated through benchmarking calculations. Calculation times are sufficiently fast to allow full MC simulations for routine treatment planning for diverse brachytherapy treatments (LDR, HDR, miniature x-ray tube). egs-brachy will be available as free and open-source software to the medical physics research community. This work is partially funded by the Canada Research Chairs program, the Natural Sciences and Engineering Research Council of Canada, and the Ontario Ministry of Research and Innovation (Ontario Early Researcher Award).

60 APPLIED LIFE SCIENCES↗

Superdense coding for quantum networking environments

Quantum networks provide a versatile infrastructure for communication, computing, and sensing with quantum information. Novel sources and detectors for transmitting and receiving quantum states are critical elements in the development and eventual deployment of robust quantum networks. Alongside performance, the compatibility of quantum network devices with modern networking infrastructure is an important requirement for deployment. We present results on the integration of quantum communication using superdense coding transmitted over optical fiber links into network environments. Our approach takes advantage of a novel complete Bell-state measurement setup that relies on hyper-entanglement in the temporal and polarization degrees of freedom for a two-photon state emitted from a quantum light source. Using linear optics and common single-photon detectors, we record a single-qubit channel capacity of 1.665±0.018. We then demonstrate a full experimental implementation of hybrid, quantum-classical communication protocol for image transfer applications. Our devices integrate with existing fiber optical network and software-defined transmitters and receivers as part of a modular design to provide an extensible quantum communication system that can adapt to future quantum technology goals.

Humble, Travis↗

Initial results of a full kinetic simulation of RF H{sup −} source including Coulomb collision process

In order to evaluate Electron Energy Distribution Function (EEDF) more correctly for radio frequency inductively coupled plasma (RF-ICP) in hydrogen negative ion sources, the Electromagnetic Particle-In-Cell (EM-PIC) simulation code has been improved by taking into account electron-electron Coulomb collision. Binary collision model is employed to model Coulomb collision process and we have successfully modeled it. The preliminary calculation including Coulomb collision has been done and it is shown that Coulomb collision doesn’t have significant effects under the condition: electron density n{sub e} ∼ 10{sup 18} m{sup −3} and high gas pressure p{sub H{sub 2}} = 3 Pa, while it is necessary to include Coulomb collision under high electron density and low gas pressure conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

sourcePy

Pollutant source identification techniques (of which there are many variations) are either locked behind researchers writing their own code for each use case or GUI platforms that are easy to use but inflexible and opaque. The Python package sourcePy brings together many of the pollutant source identification algorithms, giving the user full control out of the box. It aims to create a platform for source identification experiments where the full analysis from beginning to end can be done in Python, with a level of specificity in design that isn't available in the GUI options. sourcePy provides users with a few key features: -A Python interface with HYSPLIT, which can be used to generate trajectories and concentration plumes -Several Python classes which standardize the preparation and processing of data related to source identification experiments -Example scripts and notebooks that allow even new python users to get started with their own experiments quickly -Visualization methods

Arseneau, Isaac [Oak Ridge National Laboratory (OR↗

ADETS User Manual

ADETS (Automated DEpletion Transport System) calculates coupled neutronic/isotopic results for nuclear systems and produces a large number of criticality and burnup results based on various material feed/removal specifications, power(s), and time intervals. ADETS is a fully automated tool that links the LANL MCNP Monte Carlo transport code with the SCALE system, specifically ORIGEN-ARP (radioactive decay and burnup code). Additionally, ADETS can compute the dose at various distances of the reactor using a point gamma source, whose characteristics are computed by ORIGEN-ARP. In addition, ADETS is fully coupled with the Uncertanty Quantification, Probabilistic Risk Assessment, Parameter Optimization and Data Analysis code RAVEN (developed at the Idaho National Laboratory as well). RAVEN is released with ADETS but MCNP and SCALE (ORIGEN-ARP) are not included.

42 ENGINEERING↗

Optical properties of Peaked Spectrum radio sources

ABSTRACT In this work, we study the optical properties of compact radio sources selected from the literature in order to determine the impact of the radio-jet in their circumnuclear environment. Our sample includes 58 Compact Steep Spectrum (CSS) and GigaHertz Peaked Spectrum (GPS) and 14 Megahertz-Peaked spectrum (MPS) radio sources located at z ≤ 1. The radio luminosity (LR) of the sample varies between Log LR ∼ 23.2 and 27.7 W Hz−1. We obtained optical spectra for all sources from SDSS-DR12 and performed a stellar population synthesis using the starlight code. We derived stellar masses (M⋆), ages 〈t⋆〉, star formation rates (SFR), metallicities 〈Z⋆〉 and internal reddening AV for all young AGNs of our sample. A visual inspection of the SDSS images was made to assign a morphological class for each source. Our results indicate that the sample is dominated by intermediate to old stellar populations and there is no strong correlation between optical and radio properties of these sources. Also, we found that young AGNs can be hosted by elliptical, spiral and interacting galaxies, confirming recent findings. When comparing the optical properties of CSS/GPS and MPS sources, we do not find any significant difference. Finally, the Mid-Infrared WISE colours analysis suggests that the compact radio sources defined as powerful AGNs are, in general, gas-rich systems.

79 ASTRONOMY AND ASTROPHYSICS↗

Open Source High Fidelity Modeling Of A Type-5 Wind Turbine Drivetrain

The purpose of this code is to provide a high-fidelity model of Type-5 wind turbine drivetrain for grid integration and transient stability studies. The increasing integration of renewable energy resources in evolving bulk power system (BPS) is impacting the system’s inertia. Type-5 wind turbine generation has the potential to behave like a traditional synchronous generator and can help mitigate the aforementioned impact. A hydraulic torque converter (TC) and gearbox with torque limiting feature are integral parts of a Type-5 wind turbine unit. A generic and high-fidelity model of a Type-5 wind turbine drivetrain is not openly and widely available for grid integration and transient stability studies. This hinders appropriate assessment of Type-5 wind power plant's contribution to the BPS resilience. This work develops and validates a TC model based on those generally used in automobile's transmission system.

Hussain, Tanveer↗

Modelling a Heaving Point-Absorber with a Closed-Loop Control System Using the DualSPHysics Code

The present work addresses the need for an efficient, versatile, accurate and open-source numerical tool to be used during the design stage of wave energy converters (WECs). The device considered here is the heaving point-absorber developed and tested by Sandia National Laboratories. The smoothed particle hydrodynamics (SPH) method, as implemented in DualSPHysics, is proposed since its meshless approach presents some important advantages when simulating floating devices. The dynamics of the power take-off system are also modelled by coupling DualSPHysics with the multi-physics library Project Chrono. A satisfactory matching between experimental and numerical results is obtained for: (i) the heave response of the device when forced via its actuator; (ii) the vertical forces acting on the fixed device under regular waves and; (iii) the heave response of the WEC under the action of both regular waves and the actuator force. This proves the ability of the numerical approach proposed to simulate accurately the fluid–structure interaction along with the WEC’s closed-loop control system. In addition, radiation models built from the experimental and WAMIT results are compared with DualSPHysics by plotting the intrinsic impedance in the frequency domain, showing that the SPH method can be also employed for system identification.

42 ENGINEERING↗