Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer Based Procedures”

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,135 records · Page 63

METHODOLOGY AND APPLICATION OF PHYSICAL SECURITY EFFECTIVENESS BASED ON DYNAMIC FORCE-ON-FORCE MODELING

This paper describes ongoing work within the Light Water Reactor Sustainability (LWRS) Program at Idaho National Laboratory (INL) to optimize security and cost of nuclear power plants (NPPs). It reviews the conservatisms in conventional physical security posture and regulations. It introduces the dynamic risk assessment tool developed at INL, Event Modeling Risk Assessment using Linked Diagrams (EMRALD). The dynamic assessment methodology leverages EMRALD to process results of force-on-force (FOF) simulations and crediting safety mitigation actions from probabilistic risk assessment (PRA) models as well as diverse and flexible coping strategies (FLEX) mitigation strategies. Timing information from these simulations are compared against the available time to perform mitigations obtained from Reactor Excursion and Leak Analysis Program (RELAP5) simulations. To illustrate the methodology, a station blackout (SBO) attack scenario was modeled in commercially available FOF simulation tools. The simulation results provide valuable insights into possible attack outcomes and as the probabilistic risk of a core damage event given these outcomes. Safety mitigation procedures were modeled in EMRALD, and were dependent on the attack outcomes by considering human operator uncertainties. RELAP5 simulations incorporating human and hardware uncertainties were performed to estimate the distribution of time-to-core damage. The results demonstrate that, even in the extreme case of a successful adversarial attack, plant mitigation strategies provide significantly high-likelihood of preventing radiological release. The proposed modeling and simulation framework of integrating FLEX equipment with FOF models enables the NPPs to credit FLEX portable equipment in the plant security posture, resulting in an efficient and optimized physical security.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Para-aortic Lymphadenectomy for Gynecologic Cancers: Introducing the “Trans-Retro-Peritoneal (TRP) Single-Port Access”

Background: In gynecological surgery to date, two distinct types of endoscopic accesses have been used to perform para-aortic lymphadenectomies: transperitoneal and extraperitoneal, each with advantages and disadvantages.{sup 1} We propose to develop a new mini-invasive access to perform an endoscopic extraperitoneal para-aortic lymphadenectomy via a single-port umbilical device that, to our knowledge, has never been described. Methods: This innovative approach combines both an extraperitoneal and intraperitoneal procedure via the same umbilical incision using one single trocar. A 3–4 cm trans-umbilical incision is performed and a three-channel single-port device (Gelpoint Mini°; Applied Medical, Rancho Santa Margarita, CA, USA) is introduced in the peritoneal cavity. After peritoneal exploration, the peritoneum overlying the aorta bifurcation is grabbed, raised to the umbilicus, opened, and the single-port device is then re-introduced into the retroperitoneal space. Results: The intervention has been successful in three patients with locally advanced cervical cancer (two International Federation of Gynecology and Obstetrics [FIGO] stage IB2, and one FIGO stage IVA) scheduled for concomitant radiochemotherapy after exclusion of any suspicious lymph nodes by 18F-fluorodeoxyglucose-positron emission tomography (FDG-PET), according to our standards of practice based on the European Society of Gynaecological Oncology (ESGO) and National Comprehensive Cancer Network (NCCN) guidelines.{sup 2,3} Due to the type of cancer, lymphadenectomy was limited to the infra-mesenteric nodes,{sup 2,4} although the dissection went up to the left renal vein in each case. We retrieved 13, 20, and 25 lymph nodes. No complications occurred. Conclusions: We describe a promising technique that combines all the advantages of the two previously described accesses without their disadvantages, and with the cosmetic benefit of one almost invisible single trans-umbilical scar.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Qualification of the First ICS-3000 ION Chromatograph for use at the Defense Waste Processing Facility

The ICS-3000 Ion Chromatography (IC) system installed in 221-S M-13 has been qualified for use. The qualification was a head to head comparison of the ICS-3000 with the currently used DX-500 IC system. The crosscheck work included standards for instrument calibration and calibration verifications and standards for individual anion analysis, where the standards were traceable back to the National Institute of Standards and Technology (NIST). In addition the crosscheck work included the analysis of simulated Sludge Receipt Adjustment Tank (SRAT) Receipt, SRAT Product, and Slurry Mix Evaporator (SME) samples, along with radioactive Sludge Batch 5 material from the SRAT and SME tanks. Based upon the successful qualification of the ICS-3000 in M-13, it is recommended that this task proceed in developing the data to qualify, by a head to head comparison of the two ICS-3000 instruments, a second ICS-3000 to be installed in M-14. The Defense Waste Processing Facility (DWPF) requires the analysis of specific anions at various stages of its processing of high level waste (HLW). The anions of interest to the DWPF are fluoride, formate, chloride, nitrite, nitrate, sulfate, oxalate, and phosphate. The anion analysis is used to evaluate process chemistry including formic acid/nitric acid additions to establish optimum conditions for mercury stripping, reduction-oxidation (REDOX) chemistry for the melter, nitrite destruction, organic acid constituents, etc. The DWPF Laboratory (Lab) has been using Dionex DX-500 ion chromatography (IC) systems since 1998. The vendor informed DWPF in 2006 that the instruments would no longer be supported by service contracts after 2008. DWPF purchased three new ICS-3000 systems in September of 2006. The ICS-3000 instruments are (a) designed to be more stable using an eluent generator to make eluent, (b) require virtually no daily chemical handling by the analysts, (c) require less line breaks in the hood, and (d) generally require less maintenance due to the pump configuration only using water versus the current system where the pump uses various hydroxide concentrations. The ICS-3000 instruments also allow the DWPF to maintain current service contracts, which support routine preventive maintenance and emergency support for larger problems such as component failure. One of the three new systems was set up in the DWPF Lab trailers in January of 2007 to be used for the development of methods and procedures. This system will continue to be used for training, new method development and potential improvements to current methods. The qualification of the other two ICS-3000 instruments is to be a phased effort. This effort is to be supported by the Applied Computational Engineering and Statistical (ACES) group of the Savannah River National Laboratory (SRNL) as authorized by the Technical Task Request (TTR) and as directed by the corresponding Task Technical and Quality Assurance (TT&QA) plan. The installation of the first 'rad' system into the M-13 Lab module required modifications to both the Lab module and to the radiohood. The installation was completed in July 2008. The testing of this system was conducted as directed by the TT&QA plan. The purpose of this technical report is to provide a review of the data generated by these tests that will lead to the recommendation for the qualification of the M-13 ICS-3000 instrument. With the successful qualification of this first ICS-3000, plans will be developed for the installation of the second 'rad' system in the M-14 Lab module later in fiscal year 2009. When the second 'rad' ICS-3000 system is installed, the DX-500 systems will be removed and retired from service.

47 OTHER INSTRUMENTATION↗

TH-AB-207A-04: Assessment of Patients’ Cumulative Effective Dose From CT Examinations

Purpose: The Joint Commission requires institutions to consider patient’s age and recent imaging exams when deciding on the most appropriate type of imaging exam. Additionally, knowing patient’s imaging history can help prevent duplicate scans. Radiation dose management software affords new opportunities to identify and utilize patients with high cumulative doses as one proxy for subsequent review of imaging history and opportunities in avoiding redundant exams. Methods: Using dose management software (Radimetrics, Bayer Healthcare) a total of 72073 CT examinations performed from Jan 2015 to Jan 2016 were examined to categorize patients with a cumulative effective dose of 100 mSv and above. This threshold was selected based on epidemiological studies on populations exposed to radiation, which demonstrate a statistical increase of cancer risk at doses above 100 mSv. Histories of patients with highest cumulative dose and highest number of exams were further investigated by a Radiologist for appropriateness of recurrent studies and potential opportunities for reduction. Results: Out of 34762 patients, 927 (2.7%) were identified with a cumulative dose of 100 mSv and above. The highest cumulative dose (842 mSv) belonged to an oncology patient who underwent 2 diagnostic exams and 9 interventional ablative CT guided procedures. The patient with highest number of exams (56 counts) and cumulative dose of 170 mSv was a 17 year old trauma patient. An imaging history review of these two patients did not suggest any superfluous scans. Conclusion: Our limited pilot study suggests that recurrent CT exams for patients with oncologic or severe trauma history may be warranted and appropriate. As a result, for future studies we will be focusing on high dose patient cohorts not associated with oncology or severe trauma. Additionally, the review process itself has suggested areas for potential improvement in patient care, including improved documentation and Radiologist involvement in patient management. Dr. McNitt-Gray’s disclosures: Institutional research agreement, Siemens Healthcare; Past recipient, research grant support, Siemens Healthcare; Consultant, Toshiba America Medical Systems; Consultant, Samsung Electronics.

60 APPLIED LIFE SCIENCES↗

First Human Experience with Directly Image-able Iodinated Embolization Microbeads

PurposeTo describe first clinical experience with a directly image-able, inherently radio-opaque microspherical embolic agent for transarterial embolization of liver tumors.MethodologyLC Bead LUMI™ is a new product based upon sulfonate-modified polyvinyl alcohol hydrogel microbeads with covalently bound iodine (~260 mg I/ml). 70–150 μ LC Bead LUMI™ iodinated microbeads were injected selectively via a 2.8 Fr microcatheter to near complete flow stasis into hepatic arteries in three patients with hepatocellular carcinoma, carcinoid, or neuroendocrine tumor. A custom imaging platform tuned for LC LUMI™ microbead conspicuity using a cone beam CT (CBCT)/angiographic C-arm system (Allura Clarity FD20, Philips) was used along with CBCT embolization treatment planning software (EmboGuide, Philips).ResultsLC Bead LUMI™ image-able microbeads were easily delivered and monitored during the procedure using fluoroscopy, single-shot radiography (SSD), digital subtraction angiography (DSA), dual-phase enhanced and unenhanced CBCT, and unenhanced conventional CT obtained 48 h after the procedure. Intra-procedural imaging demonstrated tumor at risk for potential under-treatment, defined as paucity of image-able microbeads within a portion of the tumor which was confirmed at 48 h CT imaging. Fusion of pre- and post-embolization CBCT identified vessels without beads that corresponded to enhancing tumor tissue in the same location on follow-up imaging (48 h post).ConclusionLC Bead LUMI™ image-able microbeads provide real-time feedback and geographic localization of treatment in real time during treatment. The distribution and density of image-able beads within a tumor need further evaluation as an additional endpoint for embolization.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

iCycle: Integrated, multicriterial beam angle, and profile optimization for generation of coplanar and noncoplanar IMRT plans

Purpose: To introduce iCycle, a novel algorithm for integrated, multicriterial optimization of beam angles, and intensity modulated radiotherapy (IMRT) profiles. Methods: A multicriterial plan optimization with iCycle is based on a prescription called wish-list, containing hard constraints and objectives with ascribed priorities. Priorities are ordinal parameters used for relative importance ranking of the objectives. The higher an objective priority is, the higher the probability that the corresponding objective will be met. Beam directions are selected from an input set of candidate directions. Input sets can be restricted, e.g., to allow only generation of coplanar plans, or to avoid collisions between patient/couch and the gantry in a noncoplanar setup. Obtaining clinically feasible calculation times was an important design criterium for development of iCycle. This could be realized by sequentially adding beams to the treatment plan in an iterative procedure. Each iteration loop starts with selection of the optimal direction to be added. Then, a Pareto-optimal IMRT plan is generated for the (fixed) beam setup that includes all so far selected directions, using a previously published algorithm for multicriterial optimization of fluence profiles for a fixed beam arrangement Breedveld et al.[Phys. Med. Biol. 54, 7199-7209 (2009)]. To select the next direction, each not yet selected candidate direction is temporarily added to the plan and an optimization problem, derived from the Lagrangian obtained from the just performed optimization for establishing the Pareto-optimal plan, is solved. For each patient, a single one-beam, two-beam, three-beam, etc. Pareto-optimal plan is generated until addition of beams does no longer result in significant plan quality improvement. Plan generation with iCycle is fully automated. Results: Performance and characteristics of iCycle are demonstrated by generating plans for a maxillary sinus case, a cervical cancer patient, and a liver patient treated with SBRT. Plans generated with beam angle optimization did better meet the clinical goals than equiangular or manually selected configurations. For the maxillary sinus and liver cases, significant improvements for noncoplanar setups were seen. The cervix case showed that also in IMRT with coplanar setups, beam angle optimization with iCycle may improve plan quality. Computation times for coplanar plans were around 1-2 h and for noncoplanar plans 4-7 h, depending on the number of beams and the complexity of the site. Conclusions: Integrated beam angle and profile optimization with iCycle may result in significant improvements in treatment plan quality. Due to automation, the plan generation workload is minimal. Clinical application has started.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Unsupervised machine learning based on non-negative tensor factorization for analyzing reactive-mixing

Analysis of reactive-diffusion simulations representing complex mixing processes requires a large number of independent model runs. For each high-fidelity model simulation, the model inputs are varied and the predicted mixing behavior is represented by temporal and spatial changes in species concentration. It is then required to discern how the model inputs (such as diffusivity, dispersion, anisotropy, and velocity field properties) impact the mixing process. This task is challenging and typically involves interpretation of large model outputs representing temporal and spatial changes of species concentration within the model domain. Yet, the task can be automated and substantially simplified by applying Machine Learning (ML) methods. Here, we present an application of an unsupervised ML method (called NTF k ) using Non-negative Tensor Factorization (NTF) coupled with a custom clustering procedure based on k-means to reveal the temporal and spatial features in product concentrations. An attractive and unique aspect of the proposed ML method is that it ensures the extracted features are non-negative, which is important to obtain a meaningful deconstruction of the mixing processes. The ML methodology is applied to a large set of high-resolution finite-element model simulations representing anisotropic reaction-diffusion processes in perturbed vortex-based velocity fields. The applied finite-element method ensures that spatial and temporal species concentration are always non-negative, even in the case of high anisotropic contrasts. The simulated reaction is a fast irreversible bimolecular reaction A+B→C , where species A and B react to form species C. The reactive-diffusion model input parameters that control mixing include properties of the velocity field (such as vortex structures), anisotropic dispersion, and molecular diffusion. We demonstrate the applicability of the ML feature extraction method to produce a meaningful deconstruction of model outputs to discriminate between different physical processes impacting the reactants, their mixing, and the spatial distribution of the product C. The introduced ML analysis allowed us to identify additive temporal and spatial features that characterize mixing behavior. The application of the proposed NTF k approach is not limited to reactive-mixing. NTFk can be readily applied to any observed or simulated datasets that can be represented as tensors (multi-dimensional arrays) and have separable latent signatures or features.

97 MATHEMATICS AND COMPUTING↗

Population of anatomically variable 4D XCAT adult phantoms for imaging research and optimization

Purpose: The authors previously developed the 4D extended cardiac-torso (XCAT) phantom for multimodality imaging research. The XCAT consisted of highly detailed whole-body models for the standard male and female adult, including the cardiac and respiratory motions. In this work, the authors extend the XCAT beyond these reference anatomies by developing a series of anatomically variable 4D XCAT adult phantoms for imaging research, the first library of 4D computational phantoms. Methods: The initial anatomy of each phantom was based on chest-abdomen-pelvis computed tomography data from normal patients obtained from the Duke University database. The major organs and structures for each phantom were segmented from the corresponding data and defined using nonuniform rational B-spline surfaces. To complete the body, the authors manually added on the head, arms, and legs using the original XCAT adult male and female anatomies. The structures were scaled to best match the age and anatomy of the patient. A multichannel large deformation diffeomorphic metric mapping algorithm was then used to calculate the transform from the template XCAT phantom (male or female) to the target patient model. The transform was applied to the template XCAT to fill in any unsegmented structures within the target phantom and to implement the 4D cardiac and respiratory models in the new anatomy. Each new phantom was refined by checking for anatomical accuracy via inspection of the models. Results: Using these methods, the authors created a series of computerized phantoms with thousands of anatomical structures and modeling cardiac and respiratory motions. The database consists of 58 (35 male and 23 female) anatomically variable phantoms in total. Like the original XCAT, these phantoms can be combined with existing simulation packages to simulate realistic imaging data. Each new phantom contains parameterized models for the anatomy and the cardiac and respiratory motions and can, therefore, serve as a jumping point from which to create an unlimited number of 3D and 4D variations for imaging research. Conclusions: A population of phantoms that includes a range of anatomical variations representative of the public at large is needed to more closely mimic a clinical study or trial. The series of anatomically variable phantoms developed in this work provide a valuable resource for investigating 3D and 4D imaging devices and the effects of anatomy and motion in imaging. Combined with Monte Carlo simulation programs, the phantoms also provide a valuable tool to investigate patient-specific dose and image quality, and optimization for adults undergoing imaging procedures.

60 APPLIED LIFE SCIENCES↗

SU-F-T-231: Improving the Efficiency of a Radiotherapy Peer-Review System for Quality Assurance

Purpose: To improve the efficiency of a radiotherapy peer-review system using a commercially available software application for plan quality evaluation and documentation. Methods: A commercial application, FullAccess (Radialogica LLC, Version 1.4.4), was implemented in a Citrix platform for peer-review process and patient documentation. This application can display images, isodose lines, and dose-volume histograms and create plan reports for peer-review process. Dose metrics in the report can also be benchmarked for plan quality evaluation. Site-specific templates were generated based on departmental treatment planning policies and procedures for each disease site, which generally follow RTOG protocols as well as published prospective clinical trial data, including both conventional fractionation and hypo-fractionation schema. Once a plan is ready for review, the planner exports the plan to FullAccess, applies the site-specific template, and presents the report for plan review. The plan is still reviewed in the treatment planning system, as that is the legal record. Upon physician’s approval of a plan, the plan is packaged for peer review with the plan report and dose metrics are saved to the database. Results: The reports show dose metrics of PTVs and critical organs for the plans and also indicate whether or not the metrics are within tolerance. Graphical results with green, yellow, and red lights are displayed of whether planning objectives have been met. In addition, benchmarking statistics are collected to see where the current plan falls compared to all historical plans on each metric. All physicians in peer review can easily verify constraints by these reports. Conclusion: We have demonstrated the improvement in a radiotherapy peer-review system, which allows physicians to easily verify planning constraints for different disease sites and fractionation schema, allows for standardization in the clinic to ensure that departmental policies are maintained, and builds a comprehensive database for potential clinical outcome evaluation.

60 APPLIED LIFE SCIENCES↗

SU-E-T-469: Implementation of VAs Web-Based Radiotherapy Incident Reporting and Analysis System (RIRAS)

Purpose: This Web-based Radiotherapy Incident Reporting and Analysis System (RIRAS) is a tool to improve quality of care for radiation therapy patients. This system is an important facet of continuing effort by our community to maintain and improve safety of radiotherapy.Material and Methods: VA’s National Radiation Oncology Program office has embarked on a program to electronically collect adverse events and good-catch data of radiation treatment of over 25,000 veterans treated with radiotherapy annually. This VA-Intranet based software design has made use of dataset taxonomies and data dictionaries defined in AAPM/ASTRO reports on error reporting. We used proven industrial and medical event reporting techniques to avoid several common problems faced in effective data collection such as incomplete data due to data entry fatigue by the reporters, missing data due to data difficult to obtain or not familiar to most reporters, missing reports due to fear of reprisal etc. This system encompasses the entire feedback loop of reporting an incident, analyzing it for salient details, and developing interventions to prevent it from happening again. The analysis reports with corrective, learning actions are shared with the reporter/facility and made public to the community (after deidentification) as part of the learning process. Results: Till date 50 incident/good catches have been reported in RIRAS and we have completed analysis on 100% of these reports. This is done due to the fact that each reported incidents is investigated and a complete analysis/patient-safety-work-product report is generated by radiation oncology domain-experts. Conclusions Because of the completeness of the data, the system has enabled us to analyze process steps and track trends of major errors which in the future will lead to implementing system wide process improvement steps and safe standard operating procedures for each radiotherapy treatment modality/technique and fulfills our goal of “Effecting Quality While Treating Safely”. RIRAS developed and copyrighted by TSG Innovations Inc.

60 APPLIED LIFE SCIENCES↗

SU-E-T-195: Gantry Angle Dependency of MLC Leaf Position Error

Purpose: The aim of this study was to investigate the gantry angle dependency of the multileaf collimator (MLC) leaf position error. Methods: An automatic MLC quality assurance system (AutoMLCQA) was developed to evaluate the gantry angle dependency of the MLC leaf position error using an electronic portal imaging device (EPID). To eliminate the EPID position error due to gantry rotation, we designed a reference maker (RM) that could be inserted into the wedge mount. After setting up the EPID, a reference image was taken of the RM using an open field. Next, an EPID-based picket-fence test (PFT) was performed without the RM. These procedures were repeated at every 45° intervals of the gantry angle. A total of eight reference images and PFT image sets were analyzed using in-house software. The average MLC leaf position error was calculated at five pickets (-10, -5, 0, 5, and 10 cm) in accordance with general PFT guidelines using in-house software. This test was carried out for four linear accelerators. Results: The average MLC leaf position errors were within the set criterion of <1 mm (actual errors ranged from -0.7 to 0.8 mm) for all gantry angles, but significant gantry angle dependency was observed in all machines. The error was smaller at a gantry angle of 0° but increased toward the positive direction with gantry angle increments in the clockwise direction. The error reached a maximum value at a gantry angle of 90° and then gradually decreased until 180°. In the counter-clockwise rotation of the gantry, the same pattern of error was observed but the error increased in the negative direction. Conclusion: The AutoMLCQA system was useful to evaluate the MLC leaf position error for various gantry angles without the EPID position error. The Gantry angle dependency should be considered during MLC leaf position error analysis.

07 ISOTOPE AND RADIATION SOURCES↗

SU-F-T-252: An Investigation of Gamma Knife Frame Definition Error When Using a Pre-Planning Workflow

Purpose: To determine causal factors related to high frame definition error when treating GK patients using a pre-planning workflow. Methods: 160 cases were retrospectively reviewed. All patients received treatment using a pre-planning workflow whereby stereotactic coordinates are determined from a CT scan acquired after framing using a fiducial box. The planning software automatically detects the fiducials and compares their location to expected values based on the rigid design of the fiducial system. Any difference is reported as mean and maximum frame definition error. The manufacturer recommends these values be less than 1.0 mm and 1.5 mm. In this study, frame definition error was analyzed in comparison with a variety of factors including which neurosurgeon/oncologist/physicist was involved with the procedure, number of post used during framing (3 or 4), type of lesion, and which CT scanner was utilized for acquisition. An analysis of variance (ANOVA) approach was used to statistically evaluate the data and determine causal factors related to instances of high frame definition error. Results: Two factors were identified as significant: number of post (p=0.0003) and CT scanner (p=0.0001). Further analysis showed that one of the four scanners was significantly different than the others. This diagnostic scanner was identified as an older model with localization lasers not tightly calibrated. The average value for maximum frame definition error using this scanner was 1.48 mm (4 posts) and 1.75 mm (3 posts). For the other scanners this value was 1.13 mm (4 posts) and 1.40 mm (3 posts). Conclusion: In utilizing a pre-planning workflow the choice of CT scanner matters. Any scanner utilized for GK should undergo routine QA at a level appropriate for radiation oncology. In terms of 3 vs 4 post, it is hypothesized that three posts provide less stability during CT acquisition. This will be tested in future work.

60 APPLIED LIFE SCIENCES↗

Development of a Reference Governor-based Control Scheme for Integrated Energy Systems

The concept of an Integrated Energy System (IES) is meant to combine different energy technologies in synergistic ways to achieve a more secure and economical energy supply. The RAVEN-based HYBRID framework is used to find the optimal installed capacity and the optimal economical dispatch of each component of the IES. In practice what is achievable is limited by not only the inherent physical response characteristics of individual plant units but also by equipment operating limits. The inherent response characteristics are presently accounted for in the RAVEN-based HYBRID framework by dynamic models for plant units. The equipment operating limits, which are not yet addressed in HYBRID, are the subject of this report. These limits are set to achieve a design service life for a plant unit and are based on material degradation in the presence of mechanical fatigue induced by thermal cycling and by structure creep resulting from elevated temperature conditions. In general, the plant unit designer determines and specifies the allowable operating range of key process variables. It is then up to the integrated system design engineer to ensure those conditions are met. This can be accomplished through a combination of limiting the demand rates made on the plant unit and designing a control system that limits variations of key process variables. The typical control system design procedure, however, is not well suited for achieving this for complex dynamic systems such as encountered in integrated energy systems where multiple plant units must be simultaneously coordinated. This report describes how equipment operating limitations in an IES can be formally met through a control system design procedure know as the Reference Governor. The approach augments the existing controllers in a plant unit with an advance controller that modifies the issued setpoints so that selected process variables remain within limits. The design procedure is applied to the reference multi-unit system that has been the focus of Modeling and Simulation task in the Integrated Energy System program. Results describing application to the problem of limiting key process variables during operational transients for two of the plant units, the balance of plant and the gas turbine, are then presented. It is demonstrated that the Reference Governor method is an effective means for ensuring equipment operating constraints can be met in an IES. The software that implements the method has been made accessible to others on the Modeling and Simulation upload site.

42 ENGINEERING↗

Metrics for Assessment of Smart Grid Data Integrity Attacks

There is an emerging consensus that the nation’s electricity grid is vulnerable to cyber attacks. This vulnerability arises from the increasing reliance on using remote measurements, transmitting them over legacy data networks to system operators who make critical decisions based on available data. Data integrity attacks are a class of cyber attacks that involve a compromise of information that is processed by the grid operator. This information can include meter readings of injected power at remote generators, power flows on transmission lines, and relay states. These data integrity attacks have consequences only when the system operator responds to compromised data by redispatching generation under normal or contingency protocols. These consequences include (a) financial losses from sub-optimal economic dispatch to service loads, (b) robustness/resiliency losses from placing the grid at operating points that are at greater risk from contingencies, and (c) systemic losses resulting from cascading failures induced by poor operational choices. This paper is focused on understanding the connections between grid operational procedures and cyber attacks. We first offer two examples to illustrate how data integrity attacks can cause economic and physical damage by misleading operators into taking inappropriate decisions. We then focus on unobservable data integrity attacks involving power meter data. These are coordinated attacks where the compromised data are consistent with the physics of power flow, and are therefore passed by any bad data detection algorithm. We develop metrics to assess the economic impact of these attacks under re-dispatch decisions using optimal power flow methods. These metrics can be use to prioritize the adoption of appropriate countermeasures including PMU placement, encryption, hardware upgrades, and advance attack detection algorithms.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Correction for FDG PET dose extravasations: Monte Carlo validation and quantitative evaluation of patient studies

Purpose: Current procedure guidelines for whole body [18F]fluoro-2-deoxy-D-glucose (FDG)-positron emission tomography (PET) state that studies with visible dose extravasations should be rejected for quantification protocols. Our work is focused on the development and validation of methods for estimating extravasated doses in order to correct standard uptake value (SUV) values for this effect in clinical routine. Methods: One thousand three hundred sixty-seven consecutive whole body FDG-PET studies were visually inspected looking for extravasation cases. Two methods for estimating the extravasated dose were proposed and validated in different scenarios using Monte Carlo simulations. All visible extravasations were retrospectively evaluated using a manual ROI based method. In addition, the 50 patients with higher extravasated doses were also evaluated using a threshold-based method. Results: Simulation studies showed that the proposed methods for estimating extravasated doses allow us to compensate the impact of extravasations on SUV values with an error below 5%. The quantitative evaluation of patient studies revealed that paravenous injection is a relatively frequent effect (18%) with a small fraction of patients presenting considerable extravasations ranging from 1% to a maximum of 22% of the injected dose. A criterion based on the extravasated volume and maximum concentration was established in order to identify this fraction of patients that might be corrected for paravenous injection effect. Conclusions: The authors propose the use of a manual ROI based method for estimating the effectively administered FDG dose and then correct SUV quantification in those patients fulfilling the proposed criterion.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

A Stochastic Maximum Principle for a Stochastic Differential Game of a Mean-Field Type

We construct a stochastic maximum principle (SMP) which provides necessary conditions for the existence of Nash equilibria in a certain form of N-agent stochastic differential game (SDG) of a mean-field type. The information structure considered for the SDG is of a possible asymmetric and partial type. To prove our SMP we take an approach based on spike-variations and adjoint representation techniques, analogous to that of S. Peng (SIAM J. Control Optim. 28(4):966-979, 1990) in the optimal stochastic control context. In our proof we apply adjoint representation procedures at three points. The first-order adjoint processes are defined as solutions to certain mean-field backward stochastic differential equations, and second-order adjoint processes of a first type are defined as solutions to certain backward stochastic differential equations. Second-order adjoint processes of a second type are defined as solutions of certain backward stochastic equations of a type that we introduce in this paper, and which we term conditional mean-field backward stochastic differential equations. From the resulting representations, we show that the terms relating to these second-order adjoint processes of the second type are of an order such that they do not appear in our final SMP equations. A comparable situation exists in an article by R. Buckdahn, B. Djehiche, and J. Li (Appl. Math. Optim. 64(2):197-216, 2011) that constructs a SMP for a mean-field type optimal stochastic control problem; however, the approach we take of using these second-order adjoint processes of a second type to deal with the type of terms that we refer to as the second form of quadratic-type terms represents an alternative to a development, to our setting, of the approach used in their article for their analogous type of term.

97 MATHEMATICS AND COMPUTING↗

Petascale solvers for anisotropic PDEs in atmospheric modelling on GPU clusters

Memory bound applications such as solvers for large sparse systems of equations remain a challenge for GPUs. Fast solvers should be based on numerically efficient algorithms and implemented such that global memory access is minimised. To solve systems with trillions ( O ( 10 12 ) ) unknowns the code has to make efficient use of several million individual processor cores on large GPU clusters. We describe the multi-GPU implementation of two algorithmically optimal iterative solvers for anisotropic PDEs which are encountered in (semi-) implicit time stepping procedures in atmospheric modelling. In this application the condition number is large but independent of the grid resolution and both methods are asymptotically optimal, albeit with different absolute performance. In particular, an important constant in the discretisation is the CFL number; only the multigrid solver is robust to changes in this constant. We parallelise the solvers and adapt them to the specific features of GPU architectures, paying particular attention to efficient global memory access. We achieve a performance of up to 0.78 PFLOPs when solving an equation with 0.55 · 10 12 unknowns on 16384 GPUs; this corresponds to about 3% of the theoretical peak performance of the machine and we use more than 40% of the peak memory bandwidth with a Conjugate Gradient (CG) solver. Although the other solver, a geometric multigrid algorithm, has a slightly worse performance in terms of FLOPs per second, overall it is faster as it needs less iterations to converge; the multigrid algorithm can solve a linear PDE with half a trillion unknowns in about one second.

Computer Science↗

A hardware investigation of robotic SPECT for functional and molecular imaging onboard radiation therapy systems

Purpose: To construct a robotic SPECT system and to demonstrate its capability to image a thorax phantom on a radiation therapy flat-top couch, as a step toward onboard functional and molecular imaging in radiation therapy. Methods: A robotic SPECT imaging system was constructed utilizing a gamma camera detector (Digirad 2020tc) and a robot (KUKA KR150 L110 robot). An imaging study was performed with a phantom (PET CT Phantom{sup TM}), which includes five spheres of 10, 13, 17, 22, and 28 mm diameters. The phantom was placed on a flat-top couch. SPECT projections were acquired either with a parallel-hole collimator or a single-pinhole collimator, both without background in the phantom and with background at 1/10th the sphere activity concentration. The imaging trajectories of parallel-hole and pinhole collimated detectors spanned 180° and 228°, respectively. The pinhole detector viewed an off-centered spherical common volume which encompassed the 28 and 22 mm spheres. The common volume for parallel-hole system was centered at the phantom which encompassed all five spheres in the phantom. The maneuverability of the robotic system was tested by navigating the detector to trace the phantom and flat-top table while avoiding collision and maintaining the closest possible proximity to the common volume. The robot base and tool coordinates were used for image reconstruction. Results: The robotic SPECT system was able to maneuver parallel-hole and pinhole collimated SPECT detectors in close proximity to the phantom, minimizing impact of the flat-top couch on detector radius of rotation. Without background, all five spheres were visible in the reconstructed parallel-hole image, while four spheres, all except the smallest one, were visible in the reconstructed pinhole image. With background, three spheres of 17, 22, and 28 mm diameters were readily observed with the parallel-hole imaging, and the targeted spheres (22 and 28 mm diameters) were readily observed in the pinhole region-of-interest imaging. Conclusions: Onboard SPECT could be achieved by a robot maneuvering a SPECT detector about patients in position for radiation therapy on a flat-top couch. The robot inherent coordinate frames could be an effective means to estimate detector pose for use in SPECT image reconstruction.

62 RADIOLOGY AND NUCLEAR MEDICINE↗