Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distributed computation”

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,279 records · Page 71

Methods and apparatus for capture and storage of semantic information with sub-files in a parallel computing system

Techniques are provided for storing files in a parallel computing system using sub-files with semantically meaningful boundaries. A method is provided for storing at least one file generated by a distributed application in a parallel computing system. The file comprises one or more of a complete file and a plurality of sub-files. The method comprises the steps of obtaining a user specification of semantic information related to the file; providing the semantic information as a data structure description to a data formatting library write function; and storing the semantic information related to the file with one or more of the sub-files in one or more storage nodes of the parallel computing system. The semantic information provides a description of data in the file. The sub-files can be replicated based on semantically meaningful boundaries.

97 MATHEMATICS AND COMPUTING↗

Power system effects and mitigation recommendations for DER cyberattacks

Extensive deployment of interoperable distributed energy resources (DER) is increasing the power system cyber security attack surface. National and jurisdictional interconnection standards require DER to include a range of autonomous and commanded grid‐support functions, which can drastically influence power quality, voltage, and bulk system frequency. Here, the authors investigate the impact to the cyber‐physical power system in scenarios where communications and operations of DER are controlled by an adversary. The findings show that each grid‐support function exposes the power system to distinct types and magnitudes of risk. The physical impact from cyber actions was analysed in cases of DER providing distribution system voltage regulation and transmission system support. Finally, recommendations are presented for minimising the risk using engineered parameter limits and segmenting the control network to minimise common‐mode vulnerabilities.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A feasibility study of using couch-based real time dosimetric device in external beam radiotherapy

Purpose: Measurement of actual dose delivered during radiotherapy treatment aids in checking the accuracy of dose delivered to the patient. In this study, a couch-based real time dosimetric device has been proposed to measure the exit or entrance dose to a patient during external beam radiotherapy. The utility and feasibility of such a device using a 2D array of diodes has been demonstrated. Methods: Two MAPCHECK devices: MAPCHECK (1175) and MAPCHECK 2 (both SunNuclear) were embedded in a foam block in the treatment couch of a Varian 21iX linear accelerator. The angular dependence of the detector response for both devices was studied before implementing the MAPCHECKs for experimental purposes. An Alderson Rando head phantom was scanned with the MAPCHECK and MAPCHECK 2 devices separately and four different treatment plans were generated with target volumes at three different positions simulating typical clinical situations. The analytical anisotropic algorithm (AAA) was used to compute the doses in an Eclipse treatment planning system (Varian Medical Systems). The Rando phantom with the MAPCHECK device was exposed in Clinac 21iX linear accelerator. The measured dose distribution was compared with the calculated dose distribution to check for the accuracy in dose delivery. Results: Measured and computed dose distribution were found to agree with more than 93% of pixels passing at 3% and 3 mm gamma criteria for all the treatment plans. The couch-based real time dosimetry system may also be applied for noncoplanar beams where electronic portal imaging device (EPID) is not practical to measure the dose. Other advantages include checking the beam stability during the patient treatment, performing routine morning quality assurance (QA) tests in the linear accelerator, and to perform pretreatment verification of intensity modulated radiation therapy (IMRT). One of the drawbacks of this system is that it cannot be used for measuring the dose at 90 deg. or 270 deg. gantry angles. Conclusions: This preliminary study shows that a 2D array of detectors may be used as part of the treatment couch for real time patient dosimetry in studying the dose delivered to the patient in real time and also for performing routine quality assurance.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Practical Method of Adaptive Radiotherapy for Prostate Cancer Using Real-Time Electromagnetic Tracking

Purpose: We have created an automated process using real-time tracking data to evaluate the adequacy of planning target volume (PTV) margins in prostate cancer, allowing a process of adaptive radiotherapy with minimal physician workload. We present an analysis of PTV adequacy and a proposed adaptive process. Methods and Materials: Tracking data were analyzed for 15 patients who underwent step-and-shoot multi-leaf collimation (SMLC) intensity-modulated radiation therapy (IMRT) with uniform 5-mm PTV margins for prostate cancer using the Calypso Registered-Sign Localization System. Additional plans were generated with 0- and 3-mm margins. A custom software application using the planned dose distribution and structure location from computed tomography (CT) simulation was developed to evaluate the dosimetric impact to the target due to motion. The dose delivered to the prostate was calculated for the initial three, five, and 10 fractions, and for the entire treatment. Treatment was accepted as adequate if the minimum delivered prostate dose (D{sub min}) was at least 98% of the planned D{sub min}. Results: For 0-, 3-, and 5-mm PTV margins, adequate treatment was obtained in 3 of 15, 12 of 15, and 15 of 15 patients, and the delivered D{sub min} ranged from 78% to 99%, 96% to 100%, and 99% to 100% of the planned D{sub min}. Changes in D{sub min} did not correlate with magnitude of prostate motion. Treatment adequacy during the first 10 fractions predicted sufficient dose delivery for the entire treatment for all patients and margins. Conclusions: Our adaptive process successfully used real-time tracking data to predict the need for PTV modifications, without the added burden of physician contouring and image analysis. Our methods are applicable to other uses of real-time tracking, including hypofractionated treatment.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Constructing probabilistic scenarios for wide-area solar power generation

Optimizing thermal generation commitments and dispatch in the presence of high penetrations of renewable resources such as solar energy requires a characterization of their stochastic properties. In this study, we describe novel methods designed to create day-ahead, wide-area probabilistic solar power scenarios based only on historical forecasts and associated observations of solar power production. Each scenario represents a possible trajectory for solar power in next-day operations with an associated probability computed by algorithms that use historical forecast errors. Scenarios are created by segmentation of historic data, fitting non-parametric error distributions using epi-splines, and then computing specific quantiles from these distributions. Additionally, we address the challenge of establishing an upper bound on solar power output. Our specific application driver is for use in stochastic variants of core power systems operations optimization problems, e.g., unit commitment and economic dispatch. These problems require as input a range of possible future realizations of renewables production. However, the utility of such probabilistic scenarios extends to other contexts, e.g., operator and trader situational awareness. Finally, we compare the performance of our approach to a recently proposed method based on quantile regression, and demonstrate that our method performs comparably to this approach in terms of two widely used methods for assessing the quality of probabilistic scenarios: the Energy score and the Variogram score.

14 SOLAR ENERGY↗

Diaspora: Resilience-enabling services for science from HPC to edge

Scientific applications of interest to DOE must increasingly engage distributed resources (e.g., instruments, remote computers, data stores, edge devices) and deliver more stringent levels of service (e.g., uninterrupted processing of experiment data streams). In such systems, state is distributed and components can fail in many ways, often silently, making application resilience a major concern. Addressing the resilience needs of such applications requires methods for gaining knowledge of resources and applications and for translating that knowledge into action. We are working on addressing these needs in the context of multi-messenger astronomy, where detecting and responding to unusual transient events in multiple cosmic messengers (gravitational wave, electromagnetic, high- energy particles) from different instruments leads to a federated learning problem.

47 OTHER INSTRUMENTATION↗

Probing water distribution in compressed fuel-cell gas-diffusion layers using X-ray computed tomography

X-ray computed tomography was used to investigate geometrical land and channel effects on spatial liquid-water distribution in gas-diffusion layers (GDLs) of polymer-electrolyte fuel cells under different levels of compression. At low compression, a uniform liquid-water front was observed due to water redistribution and uniform porosity; however, at high compression, the water predominantly advanced at locations under the channel for higher liquid pressures. Lastly, at low compression, no apparent correlation between the spatial liquid water and porosity distributions was observed, whereas at high compression, a strong correlation was shown, indicating a potential for smart GDL architecture design with modulated porosity.

25 ENERGY STORAGE↗

Quantum Key Distribution for Critical Infrastructures: Towards Cyber-Physical Security for Hydropower and Dams

Hydropower facilities are often remotely monitored or controlled from a centralized remote control room. Additionally, major component manufacturers monitor the performance of installed components, increasingly via public communication infrastructures. While these communications enable efficiencies and increased reliability, they also expand the cyber-attack surface. Communications may use the internet to remote control a facility’s control systems, or it may involve sending control commands over a network from a control room to a machine. The content could be encrypted and decrypted using a public key to protect the communicated information. These cryptographic encoding and decoding schemes become vulnerable as more advances are made in computer technologies, such as quantum computing. In contrast, quantum key distribution (QKD) and other quantum cryptographic protocols are not based upon a computational problem, and offer an alternative to symmetric cryptography in some scenarios. Although the underlying mechanism of quantum cryptogrpahic protocols such as QKD ensure that any attempt by an adversary to observe the quantum part of the protocol will result in a detectable signature as an increased error rate, potentially even preventing key generation, it serves as a warning for further investigation. In QKD, when the error rate is low enough and enough photons have been detected, a shared private key can be generated known only to the sender and receiver. We describe how this novel technology and its several modalities could benefit the critical infrastructures of dams or hydropower facilities. The presented discussions may be viewed as a precursor to a quantum cybersecurity roadmap for the identification of relevant threats and mitigation.

97 MATHEMATICS AND COMPUTING↗

Using Parameter Sweep in WaterTAP to Analyze New Water Treatment Technologies

We describe a powerful and generalized parameter sweep tool in this report that was originally developed to analyze the performance of existing and novel water treatment models being developed in WaterTAP. Since WaterTAP is built upon IDAES and Pyomo, the parameter sweep tool can be used to systematically explore and debug the behavior of most Pyomo and IDAES numerical models. In order to enable meaningful analyses, the parameter sweep tool has been designed with the following features: 1) Model flexibility: The parameter sweep tool does not enforce any restrictions on the types of models that can be used with it. As long as a Pyomo model can be solved and the parameter is active and mutable, the tool only needs functions that describe how to run the model, the sweep parameters, and the output quantities of interest. 2) Flexible sampling: The parameter sweep tool has inbuilt functions to generate samples from a random distribution or a multidimensional Euclidean space. Furthermore, the users have to ability to supply samples generated from a tool of their choice. 3) Multiple sweep types: A user can choose from one of 3 types of parameter sweeps depending on their needs. 4) Detailed outputs: Outputs generated by the parameter sweep tool can be stored in detailed H5 file or user-friendly CSV files for post processing. 5) Parallel computing: The parameter sweep supports shared and distributed memory parallel computing to enable the use of high performance computers (HPC) for large-scale analyses. 6) Modular: The parameter sweep tool is self-contained and can easily be integrated within an outer-loop analysis or as desired by the user. 7) Ease of use: The tool is well documented and a simple sweep can be easily executed by following the online documentation in a few lines of code. We demonstrate the use of the parameter sweep tool on a simple water treatment system from the WaterTAP repository and show its parallel scaling performance on an Apple laptop and NREL's Eagle HPC. The parameter sweep tool is actively being used with models currently being developed within WaterTAP and we expect its use to grow beyond it to other IDAES and Pyomo models.

97 MATHEMATICS AND COMPUTING↗

Accelerated Distributed Energy Management for Microgrids

This paper considers the optimal energy management for a microgrid. The goal of optimal energy management is to minimize the total generation cost while meeting total demand and satisfying individual generator output limits. To achieve this goal, we develop a distributed algorithm based on the consensus theory and gradient estimation technique. Compared to the existing algorithms with diminishing step-sizes, the proposed algorithm enjoys a faster convergence speed due to the fixed step-size, and thus increases the computation speed, and reduces the computation and communication burden. The proposed distributed algorithm is applicable to both fixed and time-varying directed communication networks. We show that the convergence of the proposed distributed algorithm is achieved if the fixed communication network is strongly connected, or if the time-varying communication network is uniformly jointly strongly connected. Numerical simulation results are used to illustrate and demonstrate the effectiveness of the proposed algorithm.

Du, Wen↗

Variable Time-Step Implementation for Rapid Quasi-Static Time-Series (QSTS) Simulations of Distributed PV

Distribution system analysis with high penetrations of distributed PV require quasi-static time-series (QSTS) analysis to model the variability introduced on the distribution system, but current QSTS algorithms are prohibitively burdensome and computationally intensive. This paper proposes a variable timestep algorithm to calculate the critical time periods when QSTS simulations should be solved at higher or lower time-resolution and to backtrack for any critical periods that were missed. This variable time-step solver is a new method of performing timeseries simulations with high accuracy while performing the simulation more than 50 times faster. The scalability of the algorithm is demonstrated using a real utility distribution system model with thousands of buses.

14 SOLAR ENERGY↗

Targeting accuracy of single-isocenter intensity-modulated radiosurgery for multiple lesions

To investigate the targeting accuracy of intensity-modulated SRS (IMRS) plans designed to simultaneously treat multiple brain metastases with a single isocenter. A home-made acrylic phantom able to support a film (EBT3) in its coronal plane was used. The phantom was CT scanned and three coplanar small targets (a central and two peripheral) were outlined in the Eclipse system. Peripheral targets were 6 cm apart from the central one. A reference IMRS plan was designed to simultaneously treat the three targets, but only a single isocenter located at the center of the central target was used. After positioning the phantom on the linac using the room lasers, a CBCT scan was acquired and the reference plan were mapped on it, by placing the planned isocenter at the intersection of the landmarks used in the film showing the linac isocenter. The mapped plan was then recalculated and delivered. The film dose distribution was derived using a cloud computing application ( (www.radiochromic.com)) that uses a triple-channel dosimetry algorithm. Comparison of dose distributions using the gamma index (5%/1 mm) were performed over a 5 × 5 cm{sup 2} region centered over each target. 2D shifts required to get the best gamma passing rates on the peripheral target regions were compared with the reported ones for the central target. The experiment was repeated ten times in different sessions. Average 2D shifts required to achieve optimal gamma passing rates (99%, 97%, 99%) were 0.7 mm (SD: 0.3 mm), 0.8 mm (SD: 0.4 mm) and 0.8 mm (SD: 0.3 mm), for the central and the two peripheral targets, respectively. No statistical differences (p > 0.05) were found for targeting accuracy between the central and the two peripheral targets. The study revealed a targeting accuracy within 1 mm for off-isocenter targets within 6 cm of the linac isocenter, when a single-isocenter IMRS plan is designed.

61 RADIATION PROTECTION AND DOSIMETRY↗

Lattice-QCD Computable Quark Correlation Functions at Three-Loop Order and Extraction of Splitting Functions

We present the first complete next-to-next-to-next-to-leading-order calculation of the matching coefficients that link unpolarized flavor nonsinglet parton distribution functions with lattice QCD computable correlation functions. By using this high-order result, we notice a reduction in theoretical uncertainties compared to relying solely on previously known lower-order matching coefficients. Furthermore, based on this result we have extracted the three-loop unpolarized flavor nonsinglet splitting function, which is in agreement with the state-of-the-art result. Because of the simplicity of our method, it has the potential to advance the calculation of splitting functions to the desired four-loop order.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

SU-F-T-122: 4Dand 5D Proton Dose Evaluation with Monte Carlo

Purpose: We evaluated uncertainties in therapeutic proton doses of a lung treatment, taking into account intra-fractional geometry changes, such as breathing, and inter-fractional changes, such as tumor shrinkage and weight loss. Methods: A Monte Carlo study was performed using four dimensional CT image sets (4DCTs) and weekly repeat imaging (5DCTs) to compute fixed RBE (1.1) and variable RBE weighted dose in an actual lung treatment geometry. The MC2 Monte Carlo system was employed to simulate proton energy deposition and LET distributions according to a thoracic cancer treatment plan developed with a 3D-CT in a commercial treatment planning system, as well as in each of the phases of 4DCT sets which were recorded weekly throughout the course of the treatment. A cumulative dose distribution in relevant structures was computed and compared to the predictions of the treatment planning system. Results: Using the Monte Carlo method, dose deposition estimates with the lowest possible uncertainties were produced. Comparison with treatment planning predictions indicates that significant uncertainties may be associated with therapeutic lung dose prediction from treatment planning systems, depending on the magnitude of inter- and intra-fractional geometry changes. Conclusion: As this is just a case study, a more systematic investigation accounting for a cohort of patients is warranted; however, this is less practical because Monte Carlo simulations of such cases require enormous computational resources. Hence our study and any future case studies may serve as validation/benchmarking data for faster dose prediction engines, such as the track repeating algorithm, FDC.

60 APPLIED LIFE SCIENCES↗

Elucidating the geometric and electronic structure of a fully sulfided analog of an Anderson polyoxomolybdate cluster

The catalytic activity of transition metal sulfide (TMS) clusters in small molecule activation, redox transformations, and charge transfer has inspired the design of novel TMS-based materials for energy-related catalysis and chemical applications. Polyoxometalates (POMs), known for their structural diversity, can in principle be transformed into TMS clusters; however, fully sulfided analogs are rarely isolated, likely due to the strong tendency of uncapped TMS clusters to agglomerate. Here, we report the geometric and electronic structure of a capping ligand-free fully sulfided analog of heptamolybdate Anderson POM [Mo VI 7 O 24 ] 6− , synthesized through the sulfidation of a nanoconfined POM secured within a porous Zr-metal organic framework (NU-1000). A combined computational and experimental analysis indicates that the sulfided counterpart of the Anderson POM is geometrically and electronically more sophisticated than the parent POM. Comparison of experimental pair distribution function (PDF) data with computational simulations confirms that, unlike the oxygen-only [Mo VI 7 O 24 ] 6− cluster, the [Mo IV 7 (μ 3 -S) 6 (μ 2 -SH) 6 (S 2 ) 6 ] 2− polythiometalate (PTM) exhibits diverse sulfur anions (S 2− , HS − , S 2 2− ). DFT calculations indicate that H 2 S acts as a reducing agent, and together with terminal disulfide (S 2 2− ) ligands in the PTM structure, facilitates the complete reduction of all seven Mo VI centers in the parent POM to Mo IV . These findings are supported by X-ray photoelectron spectroscopy (XPS), which confirms exclusive Mo IV , and elemental analysis, which shows quantitative sulfur incorporation. Difference envelope density (DED) mapping further reveals that the PTM clusters are spatially confined within the MOF pores, preventing agglomeration and preserving molecular integrity.

Rabbani, S. M. Gulam [The Ohio State University, C↗

A high-performance parallel algorithm for nonnegative matrix factorization

Non-negative matrix factorization (NMF) is the problem of determining two non-negative low rank factors W and H, for the given input matrix A, such that A ≈ WH. NMF is a useful tool for many applications in different domains such as topic modeling in text mining, background separation in video analysis, and community detection in social networks. Despite its popularity in the data mining community, there is a lack of efficient distributed algorithms to solve the problem for big data sets. We propose a high-performance distributed-memory parallel algorithm that computes the factorization by iteratively solving alternating non-negative least squares (NLS) subproblems for W and H. It maintains the data and factor matrices in memory (distributed across processors), uses MPI for interprocessor communication, and, in the dense case, provably minimizes communication costs (under mild assumptions). As opposed to previous implementations, our algorithm is also flexible: (1) it performs well for both dense and sparse matrices, and (2) it allows the user to choose any one of the multiple algorithms for solving the updates to low rank factors W and H within the alternating iterations. We demonstrate the scalability of our algorithm and compare it with baseline implementations, showing significant performance improvements.

Kannan, Ramakrishnan↗

HPC-NMF: A High-Performance Parallel Algorithm for Nonnegative Matrix Factorization

NMF is a useful tool for many applications in different domains such as topic modeling in text mining, background separation in video analysis, and community detection in social networks. Despite its popularity in the data mining community, there is a lack of efficient distributed algorithms to solve the problem for big data sets. We propose a high-performance distributed-memory parallel algorithm that computes the factorization by iteratively solving alternating non-negative least squares (NLS) subproblems for $\WW$ and $\HH$. It maintains the data and factor matrices in memory (distributed across processors), uses MPI for interprocessor communication, and, in the dense case, provably minimizes communication costs (under mild assumptions). As opposed to previous implementation, our algorithm is also flexible: It performs well for both dense and sparse matrices, and allows the user to choose any one of the multiple algorithms for solving the updates to low rank factors $\WW$ and $\HH$ within the alternating iterations.

Kannan, Ramakrishnan↗