Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computer Graphics”

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

The state of Monte Carlo neutron transport on GPU accelerators - A review - 25609

A review of the current state of Monte Carlo neutron transport of Graphical Processing Units (GPUs) has been performed. Up-to-date data has been solicited from researchers in the field and combined with data available in literature to provide the current state of Monte Carlo neutron transport on GPU hardware. Nine implementations of neutron transport on GPUs have been summarized and the performance results have been presented to allow the reader to compare performance. A comparison of the results has shown that researchers obtained similar performance when scaled to the performance of NVIDIA Tesla K40 GPUs. A strong scaling equation will be presented that can be used to describe the GPU performance against different numbers of CPU cores. One implementation of a production Monte Carlo code on GPU hardware, the code Shift, showed that newer NVIDIA Pascal and Volta based GPUs provide a significant performance improvement. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2023

MATBOX, an Open-Source Microstructure Analysis Toolbox for Meshing, Generation, Segmentation, and Characterization of 3D Heterogenous Volumes

Battery performance is strongly correlated with electrode microstructural properties. To account for its impact, lithium-ion battery (LIB) models either abstract the microstructural heterogeneity of composite electrodes using effective macroscopic properties (macro- or meso- scale models) or directly solve the system of equations on the microstructure geometry or mesh (microstructure-scale models). Therefore, to be adequate, both families of models require information from the microstructure geometry, which can be provided by the numerical tool presented in this work. MATBOX is a MATLAB open-source application [1] developed by NREL for performing various microstructure-related tasks including microstructure numerical generation, image filtering and microstructure segmentation, microstructure characterization and correlation, visualization, and microstructure meshing. MATBOX was originally developed for the analysis of LIB electrode microstructures; however, the algorithms provided by the toolbox are widely applicable to other heterogeneous materials. The toolbox provides a user-friendly experience thanks to a Graphical-User Interface, requires no coding by the user, and is well documented. This presentation will illustrate various MATBOX features for the characterization of a LIB electrode, including a fully automated Representative Volume Element (RVE) analysis, the numerical generation of complex 'virtual' microstructure, including dual-layer electrodes and carbon-binder additive phase, and the meshing of a complex NMC/graphite full cell microstructure suitable for 3D finite-element modeling. Other modules (segmentation, visualization, and correlation) will be briefly presented. Thanks to its modular, open-source approach, MATBOX can easily incorporate third-party algorithms to eventually build a standard in the field that will benefit the whole scientific community. Effective diffusion coefficient [2], additive phase numerical generation [3], and meshing [4] third-party algorithms have been already integrated in the toolbox with more to come.

DIRECT ENERGY CONVERSION,MATHEMATICS AND COMPUTING↗
DOE OSTI2015

A case study of CUDA FORTRAN and OpenACC for an atmospheric climate kernel

The porting of a key kernel in the tracer advection routines of the Community Atmosphere Model – Spectral Element (CAM-SE) to use Graphics Processing Units (GPUs) using OpenACC is considered in comparison to an existing CUDA FORTRAN port. The development of the OpenACC kernel for GPUs was substantially simpler than that of the CUDA port. Also, OpenACC performance was about 1.5× slower than the optimized CUDA version. Particular focus is given to compiler maturity regarding OpenACC implementation for modern FORTRAN, and it is found that the Cray implementation is currently more mature than the PGI implementation. Still, for the case that ran successfully on PGI, the PGI OpenACC runtime was slightly faster than Cray. The results show encouraging performance for OpenACC implementation compared to CUDA while also exposing some issues that may be necessary before the implementations are suitable for porting all of CAM-SE. Furthermore, most notable are that GPU shared memory should be used by future OpenACC implementations and that derived type support should be expanded.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

SU-E-T-36: A GPU-Accelerated Monte-Carlo Dose Calculation Platform and Its Application Toward Validating a ViewRay Beam Model

Purpose: To build a fast, accurate and easily-deployable research platform for Monte-Carlo dose calculations. We port the dose calculation engine PENELOPE to C++, and accelerate calculations using GPU acceleration. Simulations of a Co-60 beam model provided by ViewRay demonstrate the capabilities of the platform. Methods: We built software that incorporates a beam model interface, CT-phantom model, GPU-accelerated PENELOPE engine, and GUI front-end. We rewrote the PENELOPE kernel in C++ (from Fortran) and accelerated the code on a GPU. We seamlessly integrated a Co-60 beam model (obtained from ViewRay) into our platform. Simulations of various field sizes and SSDs using a homogeneous water phantom generated PDDs, dose profiles, and output factors that were compared to experiment data. Results: With GPU acceleration using a dated graphics card (Nvidia Tesla C2050), a highly accurate simulation – including 100*100*100 grid, 3×3×3 mm3 voxels, <1% uncertainty, and 4.2×4.2 cm2 field size – runs 24 times faster (20 minutes versus 8 hours) than when parallelizing on 8 threads across a new CPU (Intel i7-4770). Simulated PDDs, profiles and output ratios for the commercial system agree well with experiment data measured using radiographic film or ionization chamber. Based on our analysis, this beam model is precise enough for general applications. Conclusions: Using a beam model for a Co-60 system provided by ViewRay, we evaluate a dose calculation platform that we developed. Comparison to measurements demonstrates the promise of our software for use as a research platform for dose calculations, with applications including quality assurance and treatment plan verification.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Bi{sub 1−x}Nb{sub x}O{sub 1.5+x} (x=0.0625, 0.12) fast ion conductors: Structures, stability and oxide ion migration pathways

A combined experimental and computational study of Bi{sub 1−x}Nb{sub x}O{sub 1.5+x} (x=0.0625 and 0.12) has been carried out using laboratory X-ray, neutron and electron diffraction, impedance measurements and ab-initio molecular dynamics. We demonstrate that Bi{sub 0.9375}Nb{sub 0.0625}O{sub 1.5625}, previously reported to adopt a cubic fluorite-type superstructure, can form two different polymorphs depending on the synthetic method: a metastable cubic phase is produced by quenching; while slower cooling yields a stable material with a tetragonal √2×√2×1 superstructure, which undergoes a reversible phase transition into the cubic form at ~680 °C on subsequent reheating. Neutron diffraction reveals that the tetragonal superstructure arises mainly from ordering in the oxygen sublattice, with Bi and Nb remaining disordered, although structured diffuse scattering observed in the electron diffraction patterns suggests a degree of short-range ordering. Both materials are oxide ion conductors. On thermal cycling, Bi{sub 0.88}Nb{sub 0.12}O{sub 1.62} exhibits a decrease in conductivity of approximately an order of magnitude due to partial transformation into the tetragonal phase, but still exhibits conductivity comparable to yttria-stabilised zirconia (YSZ). Ab-initio molecular dynamics simulations performed on Bi{sub 0.9375}Nb{sub 0.0625}O{sub 1.5625} show that oxide ion diffusion occurs by O{sup 2−} jumps between edge- and corner-sharing OM{sub 4} groups (M=Bi, Nb) via tetrahedral □M{sub 4} and octahedral □M{sub 6} vacancies. - Graphical abstract: Oxide ion migration in tetragonal Bi{sub 0.9375}Nb{sub 0.0625}O{sub 1.5625} occurs by O{sup 2−} jumps between edge- and corner-sharing OM{sub 4} groups (M=Bi, Nb) via tetrahedral M{sub 4} and octahedral M{sub 6} vacancies. - Highlights: • Bi{sub 0.9375}Nb{sub 0.0625}O{sub 1.5625} adopts a tetragonal √2×√2×1 fluorite superstructure. • Superstructure is due to ordering in the O-sublattice, with Bi/Nb disordered. • Bi{sub 0.9375}Nb{sub 0.0625}O{sub 1.5625} is a good oxide ion conductor. • O{sup 2−} jump between OM{sub 4} groups (M=Bi, Nb) via tetrahedral and octahedral vacancies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Li{sub 2}Ca{sub 2}Si{sub 2}O{sub 7}: Structural, spectroscopic and computational studies on a sorosilicate

Synthesis experiments in the system Li{sub 2}O–CaO–SiO{sub 2} resulted in the formation of single-crystals of Li{sub 2}Ca{sub 2}Si{sub 2}O{sub 7}. Structural investigations were based on single-crystal diffraction. At ambient conditions the compound has the following basic crystallographic data: hexagonal symmetry, space group P6{sub 1}22, a=5.0961(2) Å, c=41.264(2) Å, V=928.07(6) Å{sup 3}, Z=6. Structure solution was performed using direct methods. The final least-squares refinement calculations converged at a residual of R(|F|)=0.0260. From a structural point the lithium calcium silicate belongs to the group of pyrosilicates containing [Si{sub 2}O{sub 7}]-groups. Additional lithium and calcium cations are incorporated between the silicate dimers and are coordinated by four and six nearest oxygen neighbours, respectively. Each [LiO{sub 4}]-tetrahedron shares two common corners with directly neighboring tetrahedra forming zweier single-chains which are running parallel to 〈1 0 0〉 in z-levels defined by the presence of the 6{sub 1}{sup [0} {sup 0} {sup 1]}-screw axes. From the corner-sharing [LiO{sub 4}]- and [SiO{sub 4}]-moieties a three dimensional framework can be constructed. An interesting feature of this framework is the presence of an O{sup [3]}-type bridging oxygen linking three tetrahedra (one [LiO{sub 4}]- and two [SiO{sub 4}]-units). Structural similarities with other silicates are discussed in detail. The high-temperature behavior of the Si–O, Ca–O and Li–O bond distances in Li{sub 2}Ca{sub 2}Si{sub 2}O{sub 7} was investigated by in{sub -}situ single-crystal X-ray diffraction in the range between 65 and 700 °C. From the evolution of the lattice parameters, the thermal expansion tensor α{sub ij} has been determined. The structural characterization has been supplemented by micro-Raman spectroscopy. Interpretation of the spectroscopic data including the allocation of the bands to certain vibrational species has been aided by DFT-calculations. - Graphical abstract: Framework of [SiO{sub 4}]- and [LiO{sub 4}]-tetrahedra in the crystal structure of Li{sub 2}Ca{sub 2}Si{sub 2}O{sub 7}. - Highlights: • We report the crystal structure of the sorosilicate Li{sub 2}Ca{sub 2}Si{sub 2}O{sub 7} at ambient conditions. • The thermal expansion tensor as well as the evolution of the structure between 25 and 700 °C was studied. • A topological analysis of the tetrahedral framework based on [SiO{sub 4}]- and [LiO{sub 4}]-units is presented. • The crystal structure of Li{sub 2}Ca{sub 2}Si{sub 2}O{sub 7} is discussed with respect to related compounds. • Interpretation of the Raman spectra of the crystals has been aided by DFT-calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

CPMC-Lab: A Matlab package for Constrained Path Monte Carlo calculations

We describe CPMC-Lab, a Matlab program for the constrained-path and phaseless auxiliary-field Monte Carlo methods. These methods have allowed applications ranging from the study of strongly correlated models, such as the Hubbard model, to ab initio calculations in molecules and solids. The present package implements the full ground-state constrained-path Monte Carlo (CPMC) method in Matlab with a graphical interface, using the Hubbard model as an example. The package can perform calculations in finite supercells in any dimensions, under periodic or twist boundary conditions. Importance sampling and all other algorithmic details of a total energy calculation are included and illustrated. This open-source tool allows users to experiment with various model and run parameters and visualize the results. It provides a direct and interactive environment to learn the method and study the code with minimal overhead for setup. Furthermore, the package can be easily generalized for auxiliary-field quantum Monte Carlo (AFQMC) calculations in many other models for correlated electron systems, and can serve as a template for developing a production code for AFQMC total energy calculations in real materials. Several illustrative studies are carried out in one- and two-dimensional lattices on total energy, kinetic energy, potential energy, and charge- and spin-gaps.

Computer Science↗
DOE OSTI2014

Achieving Robustness to Uncertainty for Financial Decision-making

This report investigates the concept of robustness analysis to support financial decision-making. Financial models, that forecast future stock returns or market conditions, depend on assumptions that might be unwarranted and variables that might exhibit large fluctuations from their last-known values. The analysis of robustness explores these sources of uncertainty, and recommends model settings such that the forecasts used for decision-making are as insensitive as possible to the uncertainty. A proof-of-concept is presented with the Capital Asset Pricing Model. The robustness of model predictions is assessed using info-gap decision theory. Info-gaps are models of uncertainty that express the “distance,” or gap of information, between what is known and what needs to be known in order to support the decision. The analysis yields a description of worst-case stock returns as a function of increasing gaps in our knowledge. The analyst can then decide on the best course of action by trading-off worst-case performance with “risk”, which is how much uncertainty they think needs to be accommodated in the future. The report also discusses the Graphical User Interface, developed using the MATLAB® programming environment, such that the user can control the analysis through an easy-to-navigate interface. Three directions of future work are identified to enhance the present software. First, the code should be re-written using the Python scientific programming software. This change will achieve greater cross-platform compatibility, better portability, allow for a more professional appearance, and render it independent from a commercial license, which MATLAB® requires. Second, a capability should be developed to allow users to quickly implement and analyze their own models. This will facilitate application of the software to the evaluation of proprietary financial models. The third enhancement proposed is to add the ability to evaluate multiple models simultaneously. When two models reflect past data with similar accuracy, the more robust of the two is preferable for decision-making because its predictions are, by definition, less sensitive to the uncertainty.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Radiological Performance Modeling: Platform Comparison

Zapadnoe Uranium Mill Tailings Site, Pridneprovsky Chemical Plant Ukraine: Two exposure scenarios were utilized to develop performance assessment models in GoldSim and NORMALYSA. The first scenario is an onsite worker exposed to radioactivity through external ground irradiation as well as radon and radioactive dust inhalation. The second being an offsite Resident 300 years after institutional control exposed through the ingestion of garden vegetables irrigated by contaminated ground water, ingestion of fish from Konoplyanka River, and agricultural activities in contaminated garden. NORMALYSA utilizes a modular interface to connect extensive pre-created modules to develop contaminant transport and risk assessment models for safety assessments and performance assessments at NORM and legacy sites. Developed by Facilia AB, the program was tested with the creation of the Zapadnoe model under the International Atomic Energy Agency MODARIA Il Working Group 1 (WG1). Default values are provided for input parameters where site specific data and distributions lack. Probabilistic simulations and inputs are limited to parameters with assigned Gaussian probability functions. Similar results for transient worker inhalation dose were found in comparison with the Zapadnoe GoldSim model. However, the NORMALYSA model has results larger by an order of magnitude for external radiation in the transient worker and for aquatic food ingestion, inhalation, and garden food ingestion to the resident. Performance Assessment Modeling: Radiological performance models promote good management and support in remediation and environmental decision making by optimizing radioactive waste disposal site performance. a performance assessment model of the Zapadnoe uranium mill tailings site and a model for calculating authorized limits at Los Alamos National Laboratory are created to demonstrate and highlight the respective capabilities of three modeling platforms: GoldSim, NORMALYSA, and RESRAD-Onsite. GoldSim is a general-purpose modeling platform with the ability to simulate site-specific FT radionuclide transport and risk assessment using a graphical interface to create and link elements that act as indivisible building blocks of the model. The Zapadnoe performance assessment model created in GoldSim is composed of input parameters collected from site specific parameters from WG1 as well as standard assumptions seen in the Generic PA Model from Neptune and Company. The multivariate Zapadnoe simulation contains both deterministic and probabilistic simulation capabilities with the ability to quantitatively represent the uncertainty and associated risks of the site via the Monte Carlo simulation method. Model parameters, documentation, and assumptions are built directly into the model creating a transparent, reproducible, and more purposeful site specific model. To demonstrate this, an ingestion of water pathway was included in addition to the exposure pathways included in the NORMALYSA model. GoldSim offers a wide variety of probability distributions to choose from as well as many tools to visualize the uncertainty of the models and explore which pathways and radionuclides dominate dose. In contrast to the NORMALYSA model, the GoldSim model uses dynamic simulation to calculate tritium concentrations. Authorized limits support the evaluation of radiological health effects from radionuclides present in bulk contaminated soil in order to assist in the process of clearing land for release under different uses. Calculations of authorized limits in RESRAD-Onsite and GoldSim models support the evaluation for 15 radionuclides using a residential exposure scenario, indoor and outdoor inhalation of dust and vapor phase radionuclides (tritium), incidental soil ingestion, and ingestion grown fruits and Vegetables. RESRAD-Onsite, developed by Argonne National Laboratory is a dedicated computer code for calculating radiation doses and cancer risks to an individual located on top of radioactively contaminated soils. It can also be used to provide guidelines for deriving radionuclide soil levels corresponding to a specific dose criterion. The menu driven user interface facilitates selection of radionuclides and soil concentrations, as well as defines parameter values related to materials, environmental factors, and exposure pathways. The program contains radionuclide-specific data such as dose conversion factors and risk coefficients, as well as soil-water partition coefficients and transfer factors for biotic media. The GoldSim model successfully emulates the RESRAD-Onsite code using identical radiological decay and ingrowth factors, radiation dose conversion factors, and other input parameters. DE EEDE ETS Each model supports dynamic modeling of radionuclide ingrowth and decay over time, and integrate this capability with the transport of radionuclides in the environment via processes such as advection and/or diffusion in fluids, or within porous solid media, The applicability of each program is heavily dependent on project context. In the context of developing Authorized Limits, RESRAD-Onsite may be a more desirable option as it was developed for this particular application. For simple disposal or legacy sites or sites that do not need site-specific model structure, NORMALYSA may be a desirable modeling platform, GoldSim can be used to model any problem for which RESRAD or NORMALYSA is applicable, but its strength lies in its flexibility, analytical, and probabilistic capabilities.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

An open source fast fluid dynamics model for data center thermal management

Although computational fluid dynamics (CFD) has been widely adopted to improve data center thermal management, the high computational demand limits its applications, such as multivariate optimal design and operation. Fast fluid dynamics (FFD), which has been applied for fast airflow simulation, shows great potential. However, few research applied FFD for optimal design and operation of data center thermal management. This research improves the FFD model for data centers and conducts a comprehensive evaluation and demonstration. First, the FFD model is improved by solving the advection and diffusion equations together using an upwind scheme instead of a semi-Lagrangian advection solver in the conventional FFD model. Second, new features for data centers are added, such as a pressure correction method to simulate plenum airflow and dynamic boundary conditions for IT racks. The new FFD model is first validated with two indoor environment cases and the results show that the new FFD model has slightly better overall prediction accuracy and faster speed compared to the conventional FFD model. It is also observed that both FFD models achieve acceptable accuracy, except for a few localized disparities with experimental data, which might be due to simplified handling of turbulence viscosity near the boundaries. Furthermore, validation with a real data center shows that the FFD model achieves a similar level of accuracy as CFD when compared to the experimental measurements with some level of uncertainties. It is then demonstrated for data center optimal design and operation, which saves 53.4–58.8% of annual energy while still meeting the thermal requirements. In conclusion, with a much faster speed and comparable accuracy compared to CFD, the FFD model parallelized on a graphics processing unit is promising for practical model-based data center early design and operation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2014

Analysis of the Space Propulsion System Problem Using RAVEN

This paper presents the solution of the space propulsion problem using a PRA code currently under development at Idaho National Laboratory (INL). RAVEN (Reactor Analysis and Virtual control ENviroment) is a multi-purpose Probabilistic Risk Assessment (PRA) software framework that allows dispatching different functionalities. It is designed to derive and actuate the control logic required to simulate the plant control system and operator actions (guided procedures) and to perform both Monte- Carlo sampling of random distributed events and Event Tree based analysis. In order to facilitate the input/output handling, a Graphical User Interface (GUI) and a post-processing data-mining module are available. RAVEN allows also to interface with several numerical codes such as RELAP5 and RELAP-7 and ad-hoc system simulators. For the space propulsion system problem, an ad-hoc simulator has been developed and written in python language and then interfaced to RAVEN. Such simulator fully models both deterministic (e.g., system dynamics and interactions between system components) and stochastic behaviors (i.e., failures of components/systems such as distribution lines and thrusters). Stochastic analysis is performed using random sampling based methodologies (i.e., Monte-Carlo). Such analysis is accomplished to determine both the reliability of the space propulsion system and to propagate the uncertainties associated to a specific set of parameters. As also indicated in the scope of the benchmark problem, the results generated by the stochastic analysis are used to generate risk-informed insights such as conditions under witch different strategy can be followed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Inducing multiple functionalities in ZnS nanoparticles by doping Ni{sup +2} ions

Graphical abstract: Figure shows spherical ZnS nanoparticle (light red) containing Ni as dopant ions (light green inside) exhibiting both optical (straight lines) as well as magnetic (elliptical field lines) activities simultaneously. Highlights: ► A simple method has been adopted to synthesize Ni{sup +2} doped ZnS nanoparticles. ► The enhancing and quenching of the PL intensity depended upon ‘purity of the phase’. ► The change in packing of the molecules in ZnS NPs directly depended on concentration of dopant. ► Magnetic and optical properties were found to be dopant dependent (e.g. Ni{sup +2} ions). ► Proper choice of solvent for synthesis shows a marked effect on magnetism. - Abstract: In this study, we propose that it is possible to induce multiple functionalities such as optical activity, electrical conductance and magnetism in single ZnS/Ni nanoparticles and exploit the same by only changing the external stimuli such as magnetic field, wavelength of light, electric field etc. Such type of material finds great significance in the field of electronics as well as in bioimaging. For the purpose, we have synthesized cubic ZnS:Ni{sup 2+} nanoparticles (NPs) using a simple wet-chemical method. Synthesized ZnS:Ni{sup 2+} NPs had been characterized by X-ray diffraction pattern (XRD) and energy dispersive X-ray (EDX) analysis. The amounts of Ni and Zn in the material were determined by inductively coupled plasma-atomic emission spectroscopy (ICP-AES). Morphology of the NPs is studied by transmission electron microscopy (TEM). Optical studies are carried out using UV–visible (UV–Vis) and Photoluminescence (PL) spectroscopies. Quantum efficiency of the material was also computed. The enhancement and quenching of the PL intensity are correlated with the ‘purity of the phase’. It is observed that change in packing of the molecules of the host material directly depended on concentration of dopant ion. The thermal study of the material was carried out using differential scanning calorimetry (DSC) to support the earlier proposition. Magnetic hysteresis behavior was studied using vibrating sample magnetometer (VSM). An interesting phenomenon of defect-induced magnetism was observed, which was enhanced on doping with Ni{sup +2}.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Distributed memory, GPU accelerated Fock construction for hybrid, Gaussian basis density functional theory

With the growing reliance of modern supercomputers on accelerator-based architecture such a graphics processing units (GPUs), the development and optimization of electronic structure methods to exploit these massively parallel resources has become a recent priority. While significant strides have been made in the development GPU accelerated, distributed memory algorithms for many modern electronic structure methods, the primary focus of GPU development for Gaussian basis atomic orbital methods has been for shared memory systems with only a handful of examples pursing massive parallelism. Here in this work, we present a set of distributed memory algorithms for the evaluation of the Coulomb and exact exchange matrices for hybrid Kohn–Sham DFT with Gaussian basis sets via direct density-fitted (DF-J-Engine) and seminumerical (sn-K) methods, respectively. The absolute performance and strong scalability of the developed methods are demonstrated on systems ranging from a few hundred to over one thousand atoms using up to 128 NVIDIA A100 GPUs on the Perlmutter supercomputer.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

A Decision-Based Iterative Modeling Framework for Probabilistic Performance Assessments - 18441

Radiological probabilistic performance assessments (PPAs) employ highly complex computer models to simulate the release of radionuclides from wastes in an engineered disposal facility, migration of these radionuclides over time, and resulting radiological dose to potentially exposed individuals. An iterative approach to PPA modeling is described to help define the appropriate complexity of the computer model. An initial model is developed that addresses all relevant processes and media, but that employs preliminary statistical parameter distributions to minimize the initial effort expended. Where appropriate, transport and exposure submodels with simplifying assumptions are also used in the initial model. Different assumptions and boundaries may be explored with the initial model to assess their importance in constraining model results. Exploration of the results of the initial model are used to determine the benefits of model refinements, and to guide such refinements to focus on specific areas of the model that most influence the model outputs of interest. This approach promotes efficiency and model parsimony, but does not compromise defensibility when performed in a transparent manner. Key model boundaries for a PPA include the time frame of the modeling (temporal boundaries), where future human receptors may be located (spatial boundaries), and what types of land uses and activities may exist in the future (exposure scenarios). Simplified transport and exposure models can be useful in helping decision makers and stakeholders establish such boundaries by analyzing the consequences of selecting different bounds and understanding which assumptions and boundaries have the greatest influence on the magnitude and distribution of model results. Pathway analysis and parameter sensitivity analysis performed on PPA model results are key components of the iterative modeling approach. The objective of pathway analysis is to identify the critical radionuclides and pathways contributing to radiation dose during the modeling performance period. An environmental 'pathway' is described by the connections from a waste package or cell in the disposal facility to a future human exposure event. These connections are commonly visualized using graphical means. A simple example of a potentially complete pathway for radon-222 is: radium-containing waste → radon emanation → upwards diffusion in cap material → diffusion/advection into a future building → inhalation exposure to a building occupant Parameter sensitivity analysis is conducted on the results of a PPA to identify those specific model input parameters whose distribution of values contributes the most to the distribution of target model output, such as annual radiation dose. Pathway and sensitivity analyses on initial modeling results provide information for guiding subsequent model development, and possibly interim decisions, by informing questions tied to the modeling objectives and performance criteria. For example, questions related to radiological doses from existing disposals may include: 1. Is radiation dose from radiological exposures either obviously negligible or obviously unacceptable based on initial modeling results? 2. If dose is not obviously negligible or unacceptable relative to applicable performance criteria, what model refinements would reduce uncertainty in the results to better inform decisions? (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

TU-FG-201-07: Development of SRS Conical Collimator Collision Prediction Software for Radiation Treatment Safety

Purpose: Stereotactic radiosurgery (SRS) treatments using conical collimators can potentially result in gantry collision with treatment table due to limited collision-clear spaces. An in-house software was developed to help the SRS treatment planner mitigate potential SRS conical collimator (Varian Medical System, Palo Alto, CA) collisions with the treatment table. This software was designed to remove treatment re-planning secondary to unexpected collisions. Methods: A BrainLAB SRS ICT Frameless Extension used for SRS treatments in our clinic was mathematically modelled using surface points registered to the 3D co-ordinate space of the couch extension. The surface points are transformed based on the treatment isocenter point and potential collisions are determined in 3D space for couch and gantry angle combinations. The distance between the SRS conical collimators and LINAC isocenter is known. The collision detection model was programmed in MATLAB (Mathwork, Natick, MA) to display graphical plots of the calculations, and the plotted data is used to avoid the gantry and couch angle combinations that would likely result in a collision. We have utilized the cone collision tool for 23 SRS cone treatment plans (8 retrospective and 15 prospective for 10 patients). Results: Twenty one plans strongly agreed with the software tool prediction for collision. However, in two plans, a collision was observed with a 0.5 cm margin when the software predicted no collision. Therefore, additional margins were added to the clearance criteria in the program to achieve a lower risk of actual collisions. Conclusion: Our in-house developed collision check software successfully avoided SRS cone re-planning by 91.3% due to a reduction in cone collisions with the treatment table. Future developments to our software will include a CT image data set based collision prediction model as well as a beam angle optimization tool to avoid normal critical tissues as well as previously treated lesions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

A GPU-accelerated continuous and discontinuous Galerkin non-hydrostatic atmospheric model

In the article, we present a Graphics Processing Unit (GPU)-accelerated nodal discontinuous Galerkin method for the solution of the three-dimensional Euler equations that govern the motion and thermodynamic state of the atmosphere. Acceleration of the dynamical core of atmospheric models plays an important practical role in not only getting daily forecasts faster, but also in obtaining more accurate (high resolution) results within a given simulation time limit. We use algorithms suitable for the single instruction multiple thread architecture of GPUs to accelerate our model by two orders of magnitude relative to one core of a CPU. Tests on one node of the Titan supercomputer show a speedup of up to 15 times using the K20X GPU as compared to that on the 16-core AMD Opteron CPU. The scalability of the multi-GPU implementation is tested using 16,384 GPUs, which resulted in a weak scaling efficiency of about 90%. Finally, the accuracy and performance of our GPU implementation is verified using several benchmark problems representative of different scales of atmospheric dynamics.

Computer Science↗
DOE OSTI2013

Implementation of Information Management System for Radiation Safety of Personnel at the Russian Northwest Center for Radioactive Waste Management 'SevRAO' - 13131

The report is an overview of the information-analytical system designed to assure radiation safety of workers. The system was implemented in the Northwest Radioactive Waste Management Center 'SevRAO' (which is a branch of the Federal State Unitary Enterprise 'Radioactive Waste Management Enterprise RosRAO'). The center is located in the Northwest Russia. In respect to 'SevRAO', the Federal Medical-Biological Agency is the regulatory body, which deals with issues of radiation control. The main document to regulate radiation control is 'Reference levels of radiation factors in radioactive wastes management center'. This document contains about 250 parameters. We have developed a software tool to simplify control of these parameters. The software includes: input interface, the database, dose calculating module and analytical block. Input interface is used to enter radiation environment data. Dose calculating module calculates the dose on the route. Analytical block optimizes and analyzes radiation situation maps. Much attention is paid to the GUI and graphical representation of results. The operator can enter the route at the industrial site or watch the fluctuations of the dose rate field on the map. Most of the results are presented in a visual form. Here we present some analytical tasks, such as comparison of the dose rate in some point with control levels at this point, to be solved for the purpose of radiation safety control. The program helps to identify points making the largest contribution to the collective dose of the personnel. The tool can automatically calculate the route with the lowest dose, compare and choose the best route. The program uses several options to visualize the radiation environment at the industrial site. This system will be useful for radiation monitoring services during the operation, planning of works and development of scenarios. The paper presents some applications of this system on real data over three years - from March 2009 to February 2012. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2015

Scaling soft matter physics to thousands of graphics processing units in parallel

We describe a multi-graphics processing unit (GPU) implementation of the Ludwig application, which specialises in simulating a variety of complex fluids via lattice Boltzmann fluid dynamics coupled to additional physics describing complex fluid constituents. We describe our methodology in augmenting the original central processing unit (CPU) version with GPU functionality in a maintainable fashion. We present several optimisations that maximise performance on the GPU architecture through tuning for the GPU memory hierarchy. We describe how we implement particles within the fluid in such a way to avoid a major diversion of the CPU and GPU codebases, whilst minimising data transfer at each time step. We detail our halo-exchange communication phase for the code, which exploits overlapping to allow efficient parallel scaling to many GPUs. We present results showing that the application demonstrates excellent scaling to at least 8192 GPUs in parallel, the largest system tested at the time of writing. The GPU version (on NVIDIA K20X GPUs) is around 3.5–5 times faster that the CPU version (on fully utilised AMD Opteron 6274 16-core CPUs), comparing equal numbers of CPUs and GPUs.

Computer Science↗