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,225 records · Page 68

Cybersecurity Value-at-Risk Framework

As more variable renewable energy sources are added to the grid, the role of hydropower as a reliable baseline and firming resource is growing more critical. However, the U.S hydropower fleet is not fully prepared to face modern issues such as cybersecurity threats. Hydropower accounts for 37% of U.S. utility-scale renewable electricity but is challenged by diverse infrastructure and legacy devices that predate modern security practices. While new cybersecurity solutions cannot simply be added to current hydropower generation and operation technologies, custom cybersecurity assessments can reveal system-specific threats and risk probabilities and identify mitigating enhancements.

cybersecurity valuation methodology↗

Diagnosis of Performance Faults in LargeScale MPI Applications via Probabilistic Progress-Dependence Inference

Debugging large-scale parallel applications is challenging. Most existing techniques provide little information about failure root causes. Further, most debuggers significantly slow down program execution, and run sluggishly with massively parallel applications. This paper presents a novel technique that scalably infers the tasks in a parallel program on which a failure occurred, as well as the code in which it originated. Our technique combines scalable runtime analysis with static analysis to determine the least-progressed task(s) and to identify the code lines at which the failure arose. We present a novel algorithm that infers probabilistically progress dependence among MPI tasks using a globally constructed Markov model that represents tasks' control-flow behavior. In comparison to previous work, our algorithm infers more precisely the least-progressed task. Further, we combine this technique with static backward slicing analysis, further isolating the code responsible for the current state. A blind study demonstrates that our technique isolates the root cause of a concurrency bug in a molecular dynamics simulation, which only manifests itself at 7,996 tasks or more. We extensively evaluate fault coverage of our technique via fault injections in 10 HPC benchmarks and show that our analysis takes less than a few seconds on thousands of parallel tasks.

97 MATHEMATICS AND COMPUTING↗

Reduced Order Model of Transactive Bidding Loads

Transactive energy (TE) has been identified to provide better grid efficiency and reliability by market-based transactive exchanges between energy producers and energy consumers. Simulations of TE systems are crucial to evaluate the benefits and impacts of different transactive mechanisms. However, such simulations can be time consuming due to the information exchange between various participants and complex co-simulation environments. In this paper, we develop a reduced order model to speed up the simulation of transactive systems in TE simulation platform (TESP) while achieving very low error between the reduced order and full model results. Specifically, the developed reduced order model consists of an aggregate responsive load agent which utilizes two Recurrent Neural Networks (RNNs) with Long Short-Term Memory units (LSTMs) to enable transactive elements to collectively participate in the TE system. The proposed aggregate responsive load (ARL) agent is able to produce similar transactive behaviors to the full simulation model while achieving significant simulation time reduction. Finally, we also show that the developed model enables generalization of simulation results across different dates and across different number of loads included in the simulations.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Multifractal Characterization of Distribution Synchrophasors for Cybersecurity Defense of Smart Grids

“Source ID Mix” spoofing emerged as a new type of cyber-attack on Distribution Synchrophasors (DS) where adversaries have the capability to swap the source information of DS without changing the measurement values. Accurate detection of such a highly-deceptive attack is a challenging task especially when the spoofing attack happens on short fragments of DS recorded within a relatively small geographical scale. Herein this letter proposes an effective approach to detect this cyber-attack by realizing the multifractal characteristics of DS measurements. First, the multifractal cross-correlation of DS measured at multiple intra-state locations is revealed. Then the derived correlation is integrated with weighted two-dimensional multifractal surface interpolation to reconstruct quasi high-resolution signals. Finally, informative location-specific signatures are extracted from the high-resolution DS and they are integrated with advanced machine learning techniques for source authentication. Experiments using the real-life DS are performed to verify the proposed method.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Integrated Framework of Multisource Data Fusion for Outage Location in Looped Distribution Systems

Accurate outage location is essential for expediting post-outage power restoration, minimizing outage duration, and enhancing the resilience of distribution networks. With the advent of advanced metering infrastructure, data-driven outage location methods have significantly advanced beyond traditional approaches that rely on manual inspections. However, existing methods still face critical challenges, like reliance on single-source data, limited ability to handle partially observable systems or difficulties with loop networks. To the best of our knowledge, no single approach has comprehensively addressed all of these challenges at once. To this end, this paper proposes a comprehensive multisource data fusion framework for outage locations via probabilistic graph networks. The framework consists of three key phases. First, a novel method for reconstituting distribution networks with loops is developed, transforming looped networks into multiple radial subnetworks that retain all outage causalities of the original network. Second, Bayesian network (BN) models are established for each subnetwork, integrating multiple data sources and network structures. Finally, a joint Gibbs sampling mechanism, featuring forward and backward information flow, is designed to merge data from separate BN models and maximize the utilization of limited evidence, ensuring accurate outage location identification. In conclusion, the framework was validated on two modified public test systems, and comparative studies confirmed its effectiveness.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Fault-Tolerant Decentralized Control for Large-Scale Inverter-Based Resources for Active Power Tracking

Integration of inverter-based resources (IBRs) which lack the intrinsic characteristics such as the inertial response of the traditional synchronous-generator (SG)-based sources presents a new challenge in the form of analyzing the grid stability under their presence. While the dynamic composition of IBRs differs from that of the SGs, the control objective remains similar in terms of tracking the desired active power. This letter presents a decentralized primal-dual-based fault-tolerant control framework for the power allocation in IBRs. Overall, a hierarchical control algorithm is developed with a lower level addressing the current control and the parameter estimation for the IBRs and the higher level acting as the reference power generator to the low level based on the desired active power profile. The decentralized network-based algorithm adaptively splits the desired power between the IBRs taking into consideration the health of the IBRs transmission lines. The proposed framework is tested through a simulation on the network of IBRs and the high-level controller performance is compared against the existing framework in the literature. The proposed algorithm shows significant performance improvement in the magnitude of power deviation and settling time to the nominal value under faulty conditions as compared to the algorithm in the literature.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Comparison of Gafchromic EBT2 and EBT3 for patient-specific quality assurance: Cranial stereotactic radiosurgery using volumetric modulated arc therapy with multiple noncoplanar arcs

Purpose: Patient-specific quality assurance in volumetric modulated arc therapy (VMAT) brain stereotactic radiosurgery raises specific issues on dosimetric procedures, mainly represented by the small radiation fields associated with the lack of lateral electronic equilibrium, the need of small detectors and the high dose delivered (up to 30 Gy). Gafchromic{sup TM} EBT2 and EBT3 films may be considered the dosimeter of choice, and the authors here provide some additional data about uniformity correction for this new generation of radiochromic films.Methods: A new analysis method using blue channel for marker dye correction was proposed for uniformity correction both for EBT2 and EBT3 films. Symmetry, flatness, and field-width of a reference field were analyzed to provide an evaluation in a high-spatial resolution of the film uniformity for EBT3. Absolute doses were compared with thermoluminescent dosimeters (TLD) as baseline. VMAT plans with multiple noncoplanar arcs were generated with a treatment planning system on a selected pool of eleven patients with cranial lesions and then recalculated on a water-equivalent plastic phantom by Monte Carlo algorithm for patient-specific QA. 2D quantitative dose comparison parameters were calculated, for the computed and measured dose distributions, and tested for statistically significant differences.Results: Sensitometric curves showed a different behavior above dose of 5 Gy for EBT2 and EBT3 films; with the use of inhouse marker-dye correction method, the authors obtained values of 2.5% for flatness, 1.5% of symmetry, and a field width of 4.8 cm for a 5 × 5 cm{sup 2} reference field. Compared with TLD and selecting a 5% dose tolerance, the percentage of points with ICRU index below 1 was 100% for EBT2 and 83% for EBT3. Patients analysis revealed statistically significant differences (p < 0.05) between EBT2 and EBT3 in the percentage of points with gamma values <1 (p= 0.009 and p= 0.016); the percent difference as well as the mean difference between calculated and measured isodoses (20% and 80%) were found not to be significant (p= 0.074, p= 0.185, and p= 0.57).Conclusions: Excellent performances in terms of dose homogeneity were obtained using a new blue channel method for marker-dye correction on both EBT2 and EBT3 Gafchromic{sup TM} films. In comparison with TLD, the passing rates for the EBT2 film were higher than for EBT3; a good agreement with estimated data by Monte Carlo algorithm was found for both films, with some statistically significant differences again in favor of EBT2. These results suggest that the use of Gafchromic{sup TM} EBT2 and EBT3 films is appropriate for dose verification measurements in VMAT stereotactic radiosurgery; taking into account the uncertainty associated with Gafchromic film dosimetry, the use of adequate action levels is strongly advised, in particular, for EBT3.

61 RADIATION PROTECTION AND DOSIMETRY↗

TU-F-18C-05: Evaluation of a Method to Calculate Patient-Oriented MGD Coefficients Using Estimates of Glandular Tissue Distribution

Purpose: Evaluate mammographic mean glandular dose (MGD) coefficients for particular known tissue distributions using a novel formalism that incorporates the effect of the heterogeneous glandular tissue distribution, by comparing them with MGD coefficients derived from the corresponding anthropomorphic computer breast phantom. Methods: MGD coefficients were obtained using MCNP5 simulations with the currently used homogeneous assumption and the heterogeneously-layered breast (HLB) geometry and compared against those from the computer phantom (ground truth). The tissue distribution for the HLB geometry was estimated using glandularity map image pairs corrected for the presence of non-glandular fibrous tissue. Heterogeneity of tissue distribution was quantified using the glandular tissue distribution index, Idist. The phantom had 5 cm compressed breast thickness (MLO and CC views) and 29% whole breast glandular percentage. Results: Differences as high as 116% were found between the MGD coefficients with the homogeneous breast core assumption and those from the corresponding ground truth. Higher differences were found for cases with more heterogeneous distribution of glandular tissue. The Idist for all cases was in the [−0.8{sup −}+0.3] range. The use of the methods presented in this work results in better agreement with ground truth with an improvement as high as 105 pp. The decrease in difference across all phantom cases was in the [9{sup −}105] pp range, dependent on the distribution of glandular tissue and was larger for the cases with the highest Idist values. Conclusion: Our results suggest that the use of corrected glandularity image pairs, as well as the HLB geometry, improves the estimates of MGD conversion coefficients by accounting for the distribution of glandular tissue within the breast. The accuracy of this approach with respect to ground truth is highly dependent on the particular glandular tissue distribution studied. Predrag Bakic discloses current funding from NIH, NSF, and DoD, former funding from Real Time Tomography, LLC and a current research collaboration with Barco and Hologic.

60 APPLIED LIFE SCIENCES↗

Radiation dose calculations for CT scans with tube current modulation using the approach to equilibrium function

Purpose: The approach to equilibrium function has been used previously to calculate the radiation dose to a shift-invariant medium undergoing CT scans with constant tube current [Li, Zhang, and Liu, Med. Phys. 39, 5347–5352 (2012)]. The authors have adapted this method to CT scans with tube current modulation (TCM). Methods: For a scan with variable tube current, the scan range was divided into multiple subscan ranges, each with a nearly constant tube current. Then the dose calculation algorithm presented previously was applied. For a clinical CT scan series that presented tube current per slice, the authors adopted an efficient approach that computed the longitudinal dose distribution for one scan length equal to the slice thickness, which center was at z = 0. The cumulative dose at a specific point was a summation of the contributions from all slices and the overscan. Results: The dose calculations performed for a total of four constant and variable tube current distributions agreed with the published results of Dixon and Boone [Med. Phys. 40, 111920 (14pp.) (2013)]. For an abdomen/pelvis scan of an anthropomorphic phantom (model ATOM 701-B, CIRS, Inc., VA) on a GE Lightspeed Pro 16 scanner with 120 kV, N × T = 20 mm, pitch = 1.375, z axis current modulation (auto mA), and angular current modulation (smart mA), dose measurements were performed using two lines of optically stimulated luminescence dosimeters, one of which was placed near the phantom center and the other on the surface. Dose calculations were performed on the central and peripheral axes of a cylinder containing water, whose cross-sectional mass was about equal to that of the ATOM phantom in its abdominal region, and the results agreed with the measurements within 28.4%. Conclusions: The described method provides an effective approach that takes into account subject size, scan length, and constant or variable tube current to evaluate CT dose to a shift-invariant medium. For a clinical CT scan, dose calculations may be performed with a water-containing cylinder whose cross-sectional mass is equal to that of the subject. This method has the potential to substantially improve evaluations of patient dose from clinical CT scans, compared to CTDI{sub vol}, size-specific dose estimate (SSDE), or the dose evaluated for a TCM scan with a constant tube current equal to the average tube current of the TCM scan.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Impact of SPECT corrections on 3D-dosimetry for liver transarterial radioembolization using the patient relative calibration methodology

Purpose: Many centers aim to plan liver transarterial radioembolization (TARE) with dosimetry, even without CT-based attenuation correction (AC), or with unoptimized scatter correction (SC) methods. This work investigates the impact of presence vs absence of such corrections, and limited spatial resolution, on 3D dosimetry for TARE. Methods: Three voxelized phantoms were derived from CT images of real patients with different body sizes. Simulations of {sup 99m}Tc-SPECT projections were performed with the SIMIND code, assuming three activity distributions in the liver: uniform, inside a “liver’s segment,” or distributing multiple uptaking nodules (“nonuniform liver”), with a tumoral liver/healthy parenchyma ratio of 5:1. Projection data were reconstructed by a commercial workstation, with OSEM protocol not specifically optimized for dosimetry (spatial resolution of 12.6 mm), with/without SC (optimized, or with parameters predefined by the manufacturer; dual energy window), and with/without AC. Activity in voxels was calculated by a relative calibration, assuming identical microspheres and {sup 99m}Tc-SPECT counts spatial distribution. 3D dose distributions were calculated by convolution with {sup 90}Y voxel S-values, assuming permanent trapping of microspheres. Cumulative dose-volume histograms in lesions and healthy parenchyma from different reconstructions were compared with those obtained from the reference biodistribution (the “gold standard,” GS), assessing differences for D95%, D70%, and D50% (i.e., minimum value of the absorbed dose to a percentage of the irradiated volume). γ tool analysis with tolerance of 3%/13 mm was used to evaluate the agreement between GS and simulated cases. The influence of deep-breathing was studied, blurring the reference biodistributions with a 3D anisotropic gaussian kernel, and performing the simulations once again. Results: Differences of the dosimetric indicators were noticeable in some cases, always negative for lesions and distributed around zero for parenchyma. Application of AC and SC reduced systematically the differences for lesions by 5%–14% for a liver segment, and by 7%–12% for a nonuniform liver. For parenchyma, the data trend was less clear, but the overall range of variability passed from −10%/40% for a liver segment, and −10%/20% for a nonuniform liver, to −13%/6% in both cases. Applying AC, SC with preset parameters gave similar results to optimized SC, as confirmed by γ tool analysis. Moreover, γ analysis confirmed that solely AC and SC are not sufficient to obtain accurate 3D dose distribution. With breathing, the accuracy worsened severely for all dosimetric indicators, above all for lesions: with AC and optimized SC, −38%/−13% in liver’s segment, −61%/−40% in the nonuniform liver. For parenchyma, D50% resulted always less sensitive to breathing and sub-optimal correction methods (difference overall range: −7%/13%). Conclusions: Reconstruction protocol optimization, AC, SC, PVE and respiratory motion corrections should be implemented to obtain the best possible dosimetric accuracy. On the other side, thanks to the relative calibration, D50% inaccuracy for the healthy parenchyma from absence of AC was less than expected, while the optimization of SC was scarcely influent. The relative calibration therefore allows to perform TARE planning, basing on D50% for the healthy parenchyma, even without AC or with suboptimal corrections, rather than rely on nondosimetric methods.

60 APPLIED LIFE SCIENCES↗

Streaming Data Reorganization at Scale with DeltaFS Indexed Massive Directories

We report complex storage stacks providing data compression, indexing, and analytics help leverage the massive amounts of data generated today to derive insights. It is challenging to perform this computation, however, while fully utilizing the underlying storage media. This is because, while storage servers with large core counts are widely available, single-core performance and memory bandwidth per core grow slower than the core count per die. Computational storage offers a promising solution to this problem by utilizing dedicated compute resources along the storage processing path. We present DeltaFS Indexed Massive Directories (IMDs), a new approach to computational storage. DeltaFS IMDs harvest available (i.e., not dedicated) compute, memory, and network resources on the compute nodes of an application to perform computation on data. We demonstrate the efficiency of DeltaFS IMDs by using them to dynamically reorganize the output of a real-world simulation application across 131,072 CPU cores. DeltaFS IMDs speed up reads by 1,740x while only slightly slowing down the writing of data during simulation I/O for in situ data processing.

97 MATHEMATICS AND COMPUTING↗

E4D_RT

Time-lapse ERT imaging for monitoring both natural and engineered subsurface processes has advanced rapidly over the past 15-20 years. However, imaging results generally required a significant amount of manual and computational effort, and therefore were not available while the process was occurring. Although the value of real-time imaging was recognized, several obstacles prevented it's implementation. E4D_RT addresses these obstacles by 1) providing specialized algorithms that negate the need for user intervention, thereby automating the time-lapse data processing steps, 2) linking field data collection systems with parallel supercomputing systems via wireless data transfer link, and 3) addressing the computational burdens by utilizing distributed memory supercomputing resources, thereby enabling rapid data processing and imaging results.

Source record↗

NOODLES [SWR-22-78]

NOODLES is a protocol specification for collaborative visualization. It allows software tools of any type to participate in an analysis or visualization session. Use cases include, but are not limited to, distributed analysis, workflow interoperability, computational steering, etc.

Brunhart-Lupo, Nicholas↗

pnnl/frequency_sensitivity

This research code base includes functions to compute statistics of image distributions in (spatial DFT) frequency space, train deep learning image classifiers on these distributions with variable depth and weight decay, and finally measure the sensitivity of the trained models to perturbations along (spatial DFT) frequency components.

Central, PNNL Developer↗

fp-tools

The output of a molecular dynamics simulation is a data file containing properties of the motion of the simulated particles as a function of time (e.g., positions, velocities). From this output, scientists can extract various properties which give insight into the physics of the system; for example, the radial distribution function can be computed from the particle positions, which gives insight into how close the system is to melting or freezing. Many publicly available codes to perform molecular dynamics simulations exist and have been well adopted by the scientific community (e.g., LAMMPS, VASP). However, it remains standard practice for scientists to write their own post-processing scripts to extract various properties of interest from the molecular dynamics simulation data. Since in many cases scientists are interested in the same set of “textbook” properties, there is a large duplication of effort in writing these post-processing scripts. Our proposed code, FP-Tools, is a C++ toolkit designed to help scientists extract and analyze quantities of interest from the output of a molecular dynamics simulation, eliminating the need for them to write these codes themselves and thus reducing this duplication of effort. The quantities which are computed by our code are well documented in the literature (either in publications or textbooks), and the algorithms we implement to compute these properties are also well known in the field. We are not introducing new science or methods here; rather, our primary goal is to provide a useful tool to the community in the form of a well-documented, easy to use package.

Hartman, Leah↗

Simulation as a sustainable trading zone: Aiming at intergenerational justice

The paper, drawing on the example of simulation codes used in nuclear physics and high-energy physics, seeks to call attention to the ethical implications of discontinuing support for simulation codes and the loss of knowledge embodied in them. Predicated on the concept of trading zones and actor network models, the paper addresses the problem of extinction of simulation codes and attempts to understand their evolution and development within those frameworks. We show that simulation codes of closed type develop to the level of creoles, becoming local languages and standards of scientific centers and disappearing as their few main developers leave, whereas codes of open types become universal languages, imposing problemsolving patterns on the entire community and crowding out other codes. The paper suggests that because of simulations’ reliance on tacit knowledge, practices entrenched in codes cannot be exhaustively explicated or transmitted through writing alone; on the contrary, the life cycle of a simulation code is determined by the life cycle of its trading zone. We examine the extent to which both of these phenomena pose a risk to the preservation of knowledge. Bearing upon intergenerational ethics, we draw analogies between the pure intergenerational problem (PIP) and the problem of preserving the knowledge implemented in simulation codes and transmitting it to future generations. We argue that for the complete transfer of knowledge, it is necessary to develop and maintain inhabitability and sustainability of simulation trading zones in a controllable way, at least until the demand for these codes is warranted to cease in the future.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Evaluation of camera-based freehand SPECT in preoperative sentinel lymph node mapping for melanoma patients

Assessment of lymphatic status via sentinel lymph node (SLN) biopsy is an integral and crucial part of melanoma surgical oncology. The most common technique for sentinel node mapping is preoperative planar scintigraphy of an injected gamma-emitting lymphatic tracer followed by intraoperative node localization using a non-imaging gamma probe with auditory feedback. In recent years, intraoperative visualization of SLNs in 3D has become possible by coupling the probe to an external system capable of tracking its location and orientation as it is read out, thereby enabling computation of the 3D distribution of the tracer (freehand SPECT). In this project, the non-imaging probe of the fhSPECT system was replaced by a unique handheld gamma camera containing an array of sodium iodide crystals optically coupled to an array of silicon photomultipliers (SiPMs). A feasibility study was performed in which preoperative SLN mapping was performed using camera fhSPECT and the number of detected nodes was compared to that visualized by lymphoscintigraphy, probe fhSPECT, and to the number ultimately excised under non-imaging probe guidance. Among five subjects, SLNs were detected in nine lymphatic basins, with one to five SLNs detected per basin. A basin-by-basin comparison showed that the number of SLNs detected using camera fhSPECT exceeded that using lymphoscintigraphy and probe fhSPECT in seven of nine basins and five of five basins, respectively. (Probe fhSPECT scans were not performed for four basins.) It exceeded the number excised under non-imaging probe guidance for seven of nine basins and equaled the number excised for the other two basins. Freehand SPECT using a prototype SiPM-based gamma camera demonstrates high sensitivity for detection of SLNs in a preoperative setting. Camera fhSPECT is a potential means for efficiently obtaining real-time 3D activity distribution maps in applications such as image-guided percutaneous biopsy, and surgical SLN biopsy or radioguided tumor excision.

60 APPLIED LIFE SCIENCES↗

Overcoming the fundamental limit of quantum transduction via intraband entanglement

A quantum transducer converts an input signal to an output probe at a distant frequency band while maintaining the quantum information with high fidelity, which is crucial for quantum networking and distributed quantum sensing and computing. In terms of microwave–optical quantum transduction, the state-of-the-art quantum transducers suffer low transduction efficiency from weak nonlinear coupling, wherein increasing pump power to enhance efficiency inevitably leads to thermal noise from heating. Moreover, we reveal that the efficiency-bandwidth product of a cavity electro-optical or electro-optomechanical transducer is fundamentally limited by pump power and nonlinear coupling coefficient, irrespective of cavity engineering efforts. To overcome this fundamental limit, we propose to noiselessly boost the transduction efficiency by consuming intraband entanglement (e.g., microwave–microwave or optical–optical entanglement in the case of microwave–optical transduction). Via a squeezer–coupler–antisqueezer sandwich structure, the protocol enhances the transduction efficiency to unity in the ideal lossless case, given an arbitrarily weak pump and nonlinear coupling. In practical cavity systems, our entanglement-assisted protocol surpasses the non-assisted fundamental limit of the efficiency-bandwidth product and reduces the threshold cooperativity for positive quantum capacity by a factor proportional to two-mode squeezing gain. Given a fixed cooperativity, our approach increases the broadband quantum capacity by orders of magnitude. The entanglement-assisted advantage is robust to ancilla loss and cavity detuning.

Shi, Haowei [Southern California U.] (ORCID:000000↗