Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed algorithms”

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,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

A simplified approach to characterizing a kilovoltage source spectrum for accurate dose computation

Purpose: To investigate and validate the clinical feasibility of using half-value layer (HVL) and peak tube potential (kVp) for characterizing a kilovoltage (kV) source spectrum for the purpose of computing kV x-ray dose accrued from imaging procedures. To use this approach to characterize a Varian Registered-Sign On-Board Imager Registered-Sign (OBI) source and perform experimental validation of a novel in-house hybrid dose computation algorithm for kV x-rays. Methods: We characterized the spectrum of an imaging kV x-ray source using the HVL and the kVp as the sole beam quality identifiers using third-party freeware Spektr to generate the spectra. We studied the sensitivity of our dose computation algorithm to uncertainties in the beam's HVL and kVp by systematically varying these spectral parameters. To validate our approach experimentally, we characterized the spectrum of a Varian Registered-Sign OBI system by measuring the HVL using a Farmer-type Capintec ion chamber (0.06 cc) in air and compared dose calculations using our computationally validated in-house kV dose calculation code to measured percent depth-dose and transverse dose profiles for 80, 100, and 125 kVp open beams in a homogeneous phantom and a heterogeneous phantom comprising tissue, lung, and bone equivalent materials. Results: The sensitivity analysis of the beam quality parameters (i.e., HVL, kVp, and field size) on dose computation accuracy shows that typical measurement uncertainties in the HVL and kVp ({+-}0.2 mm Al and {+-}2 kVp, respectively) source characterization parameters lead to dose computation errors of less than 2%. Furthermore, for an open beam with no added filtration, HVL variations affect dose computation accuracy by less than 1% for a 125 kVp beam when field size is varied from 5 Multiplication-Sign 5 cm{sup 2} to 40 Multiplication-Sign 40 cm{sup 2}. The central axis depth dose calculations and experimental measurements for the 80, 100, and 125 kVp energies agreed within 2% for the homogeneous and heterogeneous block phantoms, and agreement for the transverse dose profiles was within 6%. Conclusions: The HVL and kVp are sufficient for characterizing a kV x-ray source spectrum for accurate dose computation. As these parameters can be easily and accurately measured, they provide for a clinically feasible approach to characterizing a kV energy spectrum to be used for patient specific x-ray dose computations. Furthermore, these results provide experimental validation of our novel hybrid dose computation algorithm.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

SU-E-T-27: A Tool for Routine Quality Assurance of Radiotherapy Dose Calculation Software

Purpose: Dose calculation software is thoroughly evaluated when it is commissioned; however, evaluation of periodic software updates is typically limited in scope due to staffing constraints and the need to quickly return the treatment planning system to clinical service. We developed a tool for quickly and comprehensively testing and documenting dose calculation software against measured data. Methods: A tool was developed using MatLab (The MathWorks, Natick, MA) for evaluation of dose calculation algorithms against measured data. Inputs to the tool are measured data, reference DICOM RT PLAN files describing the measurements, and dose calculations in DICOM format. The tool consists of a collection of extensible modules that can perform analysis of point dose, depth dose curves, and profiles using dose difference, distance-to-agreement, and the gamma-index. Each module generates a report subsection that is incorporated into a master template, which is converted to final form in portable document format (PDF). Results: After each change to the treatment planning system, a report can be generated in approximately 90 minutes. The tool has been in use for more than 5 years, spanning 5 versions of the eMC and 4 versions of the AAA. We have detected changes to the algorithms that affected clinical practice once during this period. Conclusion: Our tool provides an efficient method for quality assurance of dose calculation software, providing a complete set of tests for an update. Future work includes the addition of plan level tests, allowing incorporation of, for example, the TG-119 test suite for IMRT, and integration with the treatment planning system via an application programming interface. Integration with the planning system will permit fully-automated testing and reporting at scheduled intervals.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2014

SU-E-T-577: Commissioning of a Deterministic Algorithm for External Photon Beams

Purpose: We report commissioning results for a deterministic algorithm for external photon beam treatment planning. A deterministic algorithm solves the radiation transport equations directly using a finite difference method, thus improve the accuracy of dose calculation, particularly under heterogeneous conditions with results similar to that of Monte Carlo (MC) simulation. Methods: Commissioning data for photon energies 6 – 15 MV includes the percentage depth dose (PDD) measured at SSD = 90 cm and output ratio in water (Spc), both normalized to 10 cm depth, for field sizes between 2 and 40 cm and depths between 0 and 40 cm. Off-axis ratio (OAR) for the same set of field sizes was used at 5 depths (dmax, 5, 10, 20, 30 cm). The final model was compared with the commissioning data as well as additional benchmark data. The benchmark data includes dose per MU determined for 17 points for SSD between 80 and 110 cm, depth between 5 and 20 cm, and lateral offset of up to 16.5 cm. Relative comparisons were made in a heterogeneous phantom made of cork and solid water. Results: Compared to the commissioning beam data, the agreement are generally better than 2% with large errors (up to 13%) observed in the buildup regions of the FDD and penumbra regions of the OAR profiles. The overall mean standard deviation is 0.04% when all data are taken into account. Compared to the benchmark data, the agreements are generally better than 2%. Relative comparison in heterogeneous phantom is in general better than 4%. Conclusion: A commercial deterministic algorithm was commissioned for megavoltage photon beams. In a homogeneous medium, the agreement between the algorithm and measurement at the benchmark points is generally better than 2%. The dose accuracy for a deterministic algorithm is better than a convolution algorithm in heterogeneous medium.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2014

SU-E-CAMPUS-T-06: Radiochromic Film Analysis Based On Principal Components

Purpose: An algorithm to convert the color image of scanned EBT2 radiochromic film [Ashland, Covington KY] into a dose map was developed based upon a principal component analysis. The sensitive layer of the EBT2 film is colored so that the background streaks arising from variations in thickness and scanning imperfections may be distinguished by color from the dose in the exposed film. Methods: Doses of 0, 0.94, 1.9, 3.8, 7.8, 16, 32 and 64 Gy were delivered to radiochromic films by contact with a calibrated Sr-90/Y-90 source. They were digitized by a transparency scanner. Optical density images were calculated and analyzed by the method of principal components. The eigenimages of the 0.94 Gy film contained predominantly noise, predominantly background streaking, and background streaking plus the source, respectively, in order from the smallest to the largest eigenvalue. Weighting the second and third eigenimages by −0.574 and 0.819 respectively and summing them plus the constant 0.012 yielded a processed optical density image with negligible background streaking. This same weighted sum was transformed to the red, green and blue space of the scanned images and applied to all of the doses. The curve of processed density in the middle of the source versus applied dose was fit by a twophase association curve. A film was sandwiched between two polystyrene blocks and exposed edge-on to a different Y-90 source. This measurement was modeled with the GATE simulation toolkit [Version 6.2, OpenGATE Collaboration], and the on-axis depth-dose curves were compared. Results: The transformation defined using the principal component analysis of the 0.94 Gy film minimized streaking in the backgrounds of all of the films. The depth-dose curves from the film measurement and simulation are indistinguishable. Conclusion: This algorithm accurately converts EBT2 film images to dose images while reducing noise and minimizing background streaking. Supported by a sponsored research agreement from IsoTherapeutics Group, LLC.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2015

Monte Carlo validation of the TrueBeam 10XFFF phase–space files for applications in lung SABR

Purpose: To establish the clinical acceptability of universal Monte Carlo phase–space data for the 10XFFF (flattening filter free) photon beam on the Varian TrueBeam Linac, including previously unreported data for small fields, output factors, and inhomogeneous media. The study was particularly aimed at confirming the suitability for use in simulations of lung stereotactic ablative radiotherapy treatment plans. Methods: Monte Carlo calculated percent depth doses (PDDs), transverse profiles, and output factors for the TrueBeam 10 MV FFF beam using generic phase–space data that have been released by the Varian MC research team were compared with in-house measurements and published data from multiple institutions (ten Linacs from eight different institutions). BEAMnrc was used to create field size specific phase–spaces located underneath the jaws. Doses were calculated with DOSXYZnrc in a water phantom for fields ranging from 1 × 1 to 40 × 40 cm{sup 2}. Particular attention was paid to small fields (down to 1 × 1 cm{sup 2}) and dose per pulse effects on dosimeter response for high dose rate 10XFFF beams. Ion chamber measurements were corrected for changes in ion collection efficiency (P{sub ion}) with increasing dose per pulse. MC and ECLIPSE ANISOTROPIC ANALYTICAL ALGORITHM (AAA) calculated PDDs were compared to Gafchromic film measurement in inhomogeneous media (water, bone, lung). Results: Measured data from all machines agreed with Monte Carlo simulations within 1.0% and 1.5% for PDDs and in-field transverse profiles, respectively, for field sizes >1 × 1 cm{sup 2} in a homogeneous water phantom. Agreements in the 80%–20% penumbra widths were better than 2 mm for all the fields that were compared. For all the field sizes considered, the agreement between their measured and calculated output factors was within 1.1%. Monte Carlo results for dose to water at water/bone, bone/lung, and lung/water interfaces as well as within lung agree with film measurements to within 2.8% for 10 × 10 and 3 × 3 cm{sup 2} field sizes. This represents a significant improvement over the performance of the ECLIPSE AAA. Conclusions: The 10XFFF phase–space data offered by the Varian Monte Carlo research team have been validated for clinical use using measured, interinstitutional beam data in water and with film dosimetry in inhomogeneous media.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-155: Validation of a Commercial Monte Carlo Dose Calculation Algorithm for Proton Therapy

Purpose: Compare proton pencil beam scanning dose measurements to GATE/GEANT4 (GMC) and RayStation™ Monte Carlo (RMC) simulations. Methods: Proton pencil beam models of the IBA gantry at the Seattle Proton Therapy Center were developed in the GMC code system and a research build of the RMC. For RMC, a preliminary beam model that does not account for upstream halo was used. Depth dose and lateral profiles are compared for the RMC, GMC and a RayStation™ pencil beam dose (RPB) model for three spread out Bragg peaks (SOBPs) in homogenous water phantom. SOBP comparisons were also made among the three models for a phantom with a (i) 2 cm bone and a (ii) 0.5 cm titanium insert. Results: Measurements and GMC estimates of R80 range agree to within 1 mm, and the mean point-to-point dose difference is within 1.2% for all integrated depth dose (IDD) profiles. The dose differences at the peak are 1 to 2%. All of the simulated spot sigmas are within 0.15 mm of the measured values. For the three SOBPs considered, the maximum R80 deviation from measurement for GMC was −0.35 mm, RMC 0.5 mm, and RPB −0.1 mm. The minimum gamma pass using the 3%/3mm criterion for all the profiles was 94%. The dose comparison for heterogeneous inserts in low dose gradient regions showed dose differences greater than 10% at the distal edge of interface between RPB and GMC. The RMC showed improvement and agreed with GMC to within 7%. Conclusion: The RPB dosimetry show clinically significant differences (> 10%) from GMC and RMC estimates. The RMC algorithm is superior to the RPB dosimetry in heterogeneous media. We suspect modelling of the beam’s halo may be responsible for a portion of the remaining discrepancy and that RayStation will reduce this discrepancy as they finalize the release. Erik Traneus is employed as a Research Scientist at RaySearch Laboratories. The research build of the RayStation TPS used in the study was made available to the SCCA free of charge. RaySearch did not provide any monetary support other than a license to use the research build of the TPS.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Online distributed price-based control of DR resources with competitive guarantees

Demand response (DR) of building HVAC load can provide crucial demand-side flexibility for the future smart grid. Compared to direct load control, price-based control can respect the customers’ autonomy and privacy. However, it is challenging for price-based control to attain provable performance guarantees under future uncertainty. In this paper, we propose a framework for a utility to perform price-based control of flexible building load within the utility’s service area, in order to attain competitive performance guarantees in terms of controlling the system peak demand under future uncertainty. By adopting a two-step approach, our online price-based control solution can attain a provable competitive ratio for all possible realizations within a given uncertainty set. Simulation experiments demonstrate that, with a robustification procedure, our solution can perform well not only for worst-case inputs, but also for average-case inputs.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2013

Polytope

Polytope is a library that wraps several algorithms for generating the Voronoi diagram in both serial and distributed parallel.

Source record↗
DOE OSTI2020

WindWatts Computational Framework and Web UI [SWR-20-100]

This software provides a collection of algorithms, an API, and a functional Web UI for the Distributed Wind’s WindWatts project. DW WindWatts is a DOE WETO-funded project aimed at the development of tools to supporting the distributed wind industry, particularly with respect to siting and resource assessment. The computational framework and the back-end API are powering the easy to use front-end services available at: https://windwatts.nrel.gov. https://github.com/NREL/dw-tap-api; https://github.com/NREL/dw-tap; https://github.com/NREL/windwatts-data For reference, this software was previously known as DW TAP Computational Framework.

Phillips, Caleb↗
DOE OSTI2016

GridPACK™ : A Framework for Developing Power Grid Simulations on High-Performance Computing Platforms

This paper describes the GridPACK™ framework, which is designed to help power grid engineers develop modeling software capable of running on high performance computers. The framework makes extensive use of software templates to provide high level functionality while at the same time allowing developers the freedom to express whatever models and algorithms they are using. GridPACK™ contains modules for setting up distributed power grid networks, assigning buses and branches with arbitrary behaviors to the network, creating distributed matrices and vectors and using parallel linear and non-linear solvers to solve algebraic equations. It also provides mappers to create matrices and vectors based on properties of the network and functionality to support IO and to mana

Palmer, Bruce J.↗
DOE OSTI2014

Anomaly Detection in Dynamic Networks

Anomaly detection in dynamic communication networks has many important security applications. These networks can be extremely large and so detecting any changes in their structure can be computationally challenging; hence, computationally fast, parallelisable methods for monitoring the network are paramount. For this reason the methods presented here use independent node and edge based models to detect locally anomalous substructures within communication networks. As a first stage, the aim is to detect changes in the data streams arising from node or edge communications. Throughout the thesis simple, conjugate Bayesian models for counting processes are used to model these data streams. A second stage of analysis can then be performed on a much reduced subset of the network comprising nodes and edges which have been identified as potentially anomalous in the first stage. The first method assumes communications in a network arise from an inhomogeneous Poisson process with piecewise constant intensity. Anomaly detection is then treated as a changepoint problem on the intensities. The changepoint model is extended to incorporate seasonal behavior inherent in communication networks. This seasonal behavior is also viewed as a changepoint problem acting on a piecewise constant Poisson process. In a static time frame, inference is made on this extended model via a Gibbs sampling strategy. In a sequential time frame, where the data arrive as a stream, a novel, fast Sequential Monte Carlo (SMC) algorithm is introduced to sample from the sequence of posterior distributions of the change points over time. A second method is considered for monitoring communications in a large scale computer network. The usage patterns in these types of networks are very bursty in nature and don’t fit a Poisson process model. For tractable inference, discrete time models are considered, where the data are aggregated into discrete time periods and probability models are fitted to the communication counts. In a sequential analysis, anomalous behavior is then identified from outlying behavior with respect to the fitted predictive probability models. Seasonality is again incorporated into the model and is treated as a changepoint model on the transition probabilities of a discrete time Markov process. Second stage analytics are then developed which combine anomalous edges to identify anomalous substructures in the network.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

IR-thermography-based investigation of critical heat flux in subcooled flow boiling of water at atmospheric and high pressure conditions

Here we report on MIT’s THM work in Q4 2016 and Q1 2017. The goal of this project is to design, construct and execute tests of flow boiling critical heat flux (CHF) at high-pressure using high-resolution and high-speed video and infrared (IR) thermometry, to generate unique data to inform the development of and validate mechanistic boiling heat transfer and CHF models. In FY2016, a new test section was designed and fabricated. Data was collected at atmospheric conditions at 10, 25 and 50 K subcoolings, and three mass fluxes, i.e. 500, 750 and 1000 kg/m2/s. Starting in Q4 2016 and continuing forward, new post-processing techniques have been developed to analyze the data collected. These new algorithms analyze the time-dependent temperature and heat flux distributions to calculate nucleation site density, nucleation frequency, growth and wait time, dry area fraction, and the complete heat flux partitioning. In Q1 2017 a new flow boiling loop was designed and constructed to support flow boiling tests up 10 bar pressure and 180 °C. Initial shakedown and testing has been completed. The flow loop and test section are now ready to begin high-pressure flow boiling testing.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2018

Interatomic Potentials Models for Cu-Ni and Cu-Zr Alloys

This study explores a Bayesian calibration framework for the RAMPAGE alloy potential model for Cu-Ni and Cu-Zr systems, respectively. In RAMPAGE potentials, it is proposed that once calibrated potentials for individual elements are available, the inter-species interactions can be described by fitting a Morse potential for pair interactions with three parameters, while densities for the embedding function can be scaled by two parameters from the elemental densities. Global sensitivity analysis tools were employed to understand the impact each parameter has on the MD simulation results. A transitional Markov Chain Monte Carlo algorithm was used to generate samples from the multimodal posterior distribution consistent with the discrepancy between MD simulation results and DFT data. For the Cu-Ni system the posterior predictive tests indicate that the fitted interatomic potential model agrees well with the DFT data, justifying the basic RAMPAGE assumptions. For the Cu-Zr system, where the phase diagram suggests more complicated atomic interactions than in the case of Cu-Ni, the RAMPAGE potential captured only a subset of the DFT data. The resulting posterior distribution for the 5 model parameters exhibited several modes, with each mode corresponding to specific simulation data and a suboptimal agreement with the DFT results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Gamma-ray raster imaging with robotic data collection. FY18 report

This project has developed an alternative imaging capability for measuring radioactive materials in situ. This approach uses a robotic-mounted gamma-ray detector which can move around an area of interest, sampling the space at a high frequency. By rastering across the area, an image can be created with no collimation and a high efficiency. A regression model was developed to map out material distribution as the robot moves through the environment. An informative path planning algorithm was also developed to maximize the information collected throughout the movement path. The use of a robotic mount allows data collection for long periods of time unattended, and it will also eliminate uncertainties in positioning typically introduced by personnel. This approach will be particularly relevant for gloveboxes, shielded cells, or process piping which may have complex, non-uniform distributions of material.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2017

The Bird project: Using Big Data tools to support Search Analytics

The Bird project explored the use of big data analytics tool to improve the findability of information within the Sandia internal network. We were able to perform query classification utilizing the supervised learning algorithms in the Apache Spark library. By relying on the distributed processing capabilities provided by the Apache Hadoop framework, we successfully processed the large query log files needed to train the models in this effort. The capabilities developed in this project are being used to enhance the effectiveness of the enterprise search engine.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2019

SGEMP Algorithm Development and Demonstration for ICEPIC

Algorithms designed to facilitate calculation of gas-modified System-Generated Electromagnetic Pulse (SGEMP) phenomena have been developed for use in ICEPIC. Improvements and extensions of these and related models are considered for helping them be more widely applicable for SGEMP modeling. Confluent Sciences has developed a Buneman-Villasenor unit test tool for assessing particle merge algorithms and has applied that to several candidate algorithms. Statistical analysis of Maxwellian and near-Maxwellian velocity distribution functions has been used to understand limits and tradeoffs of particle merging. Collisional phenomena and plasma models that relieve the burden of a fully kinetic approach have also been developed and assessed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

NRC Multiphysics Analysis Capability Deployment (FY2021--Part 1)

This report details progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Support for Advanced Non-Light Water Reactors.” The tasks completed for this report are as follows: First, Task 1d: The net radiation transfer method was implemented into MOOSE for modeling reactor cavity cooling system geometries. RCCS models for two experiments were created: (1) Natural Convection Shutdown Heat Removal Test Facility (NSTF) experiment R022, and (2) HTTR VCS mockup. For validation, computed temperature distributions were compared to measured temperatures. Next, Task 4c: An algorithm for computing the pebble bed reactor equilibrium core isotopic com-position was developed and an initial version is implemented into the reactor multi-physics code Griffin. Initial results for a simplified axisymmetric pebble bed reactor are presented. Finally, Task 7: generation of a reference plant model for molten salt cooled pebble bed reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Conceptual Machine Learning Strategy for Maintenance of Molten Salt Heat Exchanger

This report addresses the requirements of Milestone M3.2 Define maintenance scenarios and benchmarks. Development of digital twins (DT) for molten salt reactor (MSR) components is crucial for reducing operating and maintenance costs (O&M) and ensuring commercial viability of these reactors. Our focus is on development of DT for MSR primary system heat exchanger (HX), a critical component, the fault in which can reduce operating efficiency and force reactor shutdown. We are investigating feasibility of a conceptual DT of HX consisting of internal distributed temperature sensing with fiber optics and machine learning (ML) algorithms to detect and localize faults.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗