Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “high performance visualization”

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,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

ORNL Campus Sustainability and Decarbonization using Waste Heat Recovery from the Oak Ridge Leadership Computing Facility’s High-Performance Computing Data Center

Heat pumps are a clean and efficient technology that can be powered by renewable electricity to transfer heat using a refrigerant from one place to another by different heat sources, making buildings clean and environmentally friendly. With the support of the ORNL Laboratory Modernization Division, this project explored and evaluated an innovative solution that uses water-water cost-effective midtemperature heat pump (MTHP) technology to leverage the low-grade waste heat from ORNL Frontier and the data center to deliver 85°C hot water, which replaces hot steam generated using natural gas combustion boilers for water heating or space heating in the buildings of ORNL campus. Two scenarios were studied. In the first scenario, which considered the 5600-5700-5800 complex only, Carrier’s commercial 1,000 kW MTHP technology achieves more than 6,640 MWh/year energy savings, an emission reduction of 858 TCO2e/year CO2, and a payback time of 4.85 years. In the second scenario, which considered the 5600-5700-5800 complex and Buildings 5100, 5200, and 5300, the CO2 emission reduction is 1,483 TCO2e/year, the operating cost savings are $0.21 million annually, and the payback time is 3.74 years. Additionally, a comprehensive HP ShowCase Tool was developed for evaluating the optimal solution to improve sustainability and decarbonization of the buildings on the ORNL campus. The tool is an Excel-based tool integrated with VBA (Visual Basic for Applications) coding. The tool includes collected ORNL campus building information and an MTHP library, which comprises collected commercial and ORNL-defined MTHPs. The tool was used to evaluate the sustainability and decarbonization of the ORNL campus. The tool can be widely used or referenced for heat pump solutions and building decarbonization renovation strategies to modernize ORNL facilities and energy use–intensive equipment to enable efficient, sustainable, and resilient operations in the future.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1992

Mixing in straight shear layers

Planar laser-induced fluorescence measurements were performed in a liquid plane mixing layer to extract the probability density function (pdf) of the mixture fraction of a passive scalar across the layer. Three Reynolds number (Re) cases were studied, 10,000, 33,000 and 90,000, with Re based on velocity difference and visual thickness. The results show that a non-marching pdf (central hump invariant from edge to edge of the layer) exists for Re = 10,000 but that a marching type pdf characterizes the Re = 33,000 and Re = 90,000 cases. For all cases, a broad range of mixture fraction values is found at each location across the layer. Streamwise and spanwise ramps across the layer, and structure-to-structure variation were observed and are believed to be responsible for the above behavior of the composition field. Tripping the boundary layer on the high-speed side of the splitter plate for each of the above three cases resulted in increased three-dimensionality and a change in the composition field. Average and average mixed fluid compositions are reported for all cases.

Karasso, P. S.↗
NASA NTRS2014

NASA Tech Briefs, June 2014

Topics include: Real-Time Minimization of Tracking Error for Aircraft Systems; Detecting an Extreme Minority Class in Hyperspectral Data Using Machine Learning; KSC Spaceport Weather Data Archive; Visualizing Acquisition, Processing, and Network Statistics Through Database Queries; Simulating Data Flow via Multiple Secure Connections; Systems and Services for Near-Real-Time Web Access to NPP Data; CCSDS Telemetry Decoder VHDL Core; Thermal Response of a High-Power Switch to Short Pulses; Solar Panel and System Design to Reduce Heating and Optimize Corridors for Lower-Risk Planetary Aerobraking; Low-Cost, Very Large Diamond-Turned Metal Mirror; Very-High-Load-Capacity Air Bearing Spindle for Large Diamond Turning Machines; Elevated-Temperature, Highly Emissive Coating for Energy Dissipation of Large Surfaces; Catalyst for Treatment and Control of Post-Combustion Emissions; Thermally Activated Crack Healing Mechanism for Metallic Materials; Subsurface Imaging of Nanocomposites; Self-Healing Glass Sealants for Solid Oxide Fuel Cells and Electrolyzer Cells; Micromachined Thermopile Arrays with Novel Thermo - electric Materials; Low-Cost, High-Performance MMOD Shielding; Head-Mounted Display Latency Measurement Rig; Workspace-Safe Operation of a Force- or Impedance-Controlled Robot; Cryogenic Mixing Pump with No Moving Parts; Seal Design Feature for Redundancy Verification; Dexterous Humanoid Robot; Tethered Vehicle Control and Tracking System; Lunar Organic Waste Reformer; Digital Laser Frequency Stabilization via Cavity Locking Employing Low-Frequency Direct Modulation; Deep UV Discharge Lamps in Capillary Quartz Tubes with Light Output Coupled to an Optical Fiber; Speech Acquisition and Automatic Speech Recognition for Integrated Spacesuit Audio Systems, Version II; Advanced Sensor Technology for Algal Biotechnology; High-Speed Spectral Mapper; "Ascent - Commemorating Shuttle" - A NASA Film and Multimedia Project DVD; High-Pressure, Reduced-Kinetics Mechanism for N-Hexadecane Oxidation; Method of Error Floor Mitigation in Low-Density Parity-Check Codes; X-Ray Flaw Size Parameter for POD Studies; Large Eddy Simulation Composition Equations for Two-Phase Fully Multicomponent Turbulent Flows; Scheduling Targeted and Mapping Observations with State, Resource, and Timing Constraints;

Source record↗
DOE OSTI2022

Optimizing Magnetic Force Microscopy Resolution and Sensitivity to Visualize Nanoscale Magnetic Domains

Magnetic force microscopy (MFM) enables mapping local magnetic fields across a sample surface with nanoscale resolution. To perform MFM, an atomic force microscopy (AFM) probe whose tip has been magnetized vertically (i.e., perpendicular to the probe cantilever) is oscillated at a fixed height above the sample surface. The resultant shifts in the oscillation phase or frequency, which are proportional to the magnitude and sign of the vertical magnetic force gradient at each pixel location, are then tracked and mapped. Although the spatial resolution and sensitivity of the technique increases with decreasing lift height above the surface, this seemingly straightforward path to improved MFM images is complicated by considerations such as minimizing topographical artifacts due to shorter range van der Waals forces, increasing the oscillation amplitude to further improve sensitivity, and the presence of surface contaminants (in particular water due to humidity under ambient conditions). In addition, due to the orientation of the probe's magnetic dipole moment, MFM is intrinsically more sensitive to samples with an out-of-plane magnetization vector. Here, high-resolution topographical and magnetic phase images of single and bicomponent nanomagnet artificial spin-ice (ASI) arrays obtained in an inert (argon) atmosphere glovebox with <0.1 ppm O 2 and H 2 O are reported. Further, optimization of lift height and drive amplitude for high resolution and sensitivity while simultaneously avoiding the introduction of topographical artifacts is discussed, and detection of the stray magnetic fields emanating from either end of the nanoscale bar magnets (~250 nm long and <100 nm wide) aligned in the plane of the ASI sample surface is shown. Likewise, using the example of a Ni-Mn-Ga magnetic shape memory alloy (MSMA), MFM is demonstrated in an inert atmosphere with magnetic phase sensitivity capable of resolving a series of adjacent magnetic domains each ~200 nm wide.

47 OTHER INSTRUMENTATION↗
DOE OSTI2018

Dynamic Data Repartitioning for Load-Balanced Parallel Particle Tracing

We present a novel dynamic load-balancing algorithm based on data repartitioning for parallel particle tracing in flow visualization. Instead of static data assignment, we dynamically repartition the data into blocks and reassign the blocks to processes to balance the workload distribution among the processes. Block repartitioning is performed based on a dynamic workload estimation method that predicts the workload in the flow field on the fly as the input. In our approach, we allow data duplication in the repartitioning, enabling the same data blocks to be assigned to multiple processes. Load balance is achieved by regularly exchanging the blocks (together with the particles in the blocks) among processes according to the output of the data repartitioning. Compared with other load-balancing algorithms, our approach does not need any preprocessing on the raw data and does not require any dedicated process for work scheduling, while it has the capability to balance uneven workload efficiently. Results show improved load balance and high efficiency of our method on tracing particles in both steady and unsteady flow.

data repartitioning↗
NASA NTRS2011

Space Radiation Effects on Human Cells: Modeling DNA Breakage, DNA Damage Foci Distribution, Chromosomal Aberrations and Tissue Effects

Future long-tem space travel will face challenges from radiation concerns as the space environment poses health risk to humans in space from radiations with high biological efficiency and adverse post-flight long-term effects. Solar particles events may dramatically affect the crew performance, while Galactic Cosmic Rays will induce a chronic exposure to high-linear-energy-transfer (LET) particles. These types of radiation, not present on the ground level, can increase the probability of a fatal cancer later in astronaut life. No feasible shielding is possible from radiation in space, especially for the heavy ion component, as suggested solutions will require a dramatic increase in the mass of the mission. Our research group focuses on fundamental research and strategic analysis leading to better shielding design and to better understanding of the biological mechanisms of radiation damage. We present our recent effort to model DNA damage and tissue damage using computational models based on the physics of heavy ion radiation, DNA structure and DNA damage and repair in human cells. Our particular area of expertise include the clustered DNA damage from high-LET radiation, the visualization of DSBs (DNA double strand breaks) via DNA damage foci, image analysis and the statistics of the foci for different experimental situations, chromosomal aberration formation through DSB misrepair, the kinetics of DSB repair leading to a model-derived spectrum of chromosomal aberrations, and, finally, the simulation of human tissue and the pattern of apoptotic cell damage. This compendium of theoretical and experimental data sheds light on the complex nature of radiation interacting with human DNA, cells and tissues, which can lead to mutagenesis and carcinogenesis later in human life after the space mission.

Ponomarev, A. L.↗
DOE OSTI2019

Imaging doxorubicin and polymer-drug conjugates of doxorubicin-induced cardiotoxicity with bispecific anti-myosin-anti-DTPA antibody and Tc-99m-labeled polymers

Background: Radiolabeled anti-myosin imaging is well-established for imaging doxorubicin-induced cardiotoxicity. However, to enable imaging of drug-induced cardiotoxicity in small experimental animals, pretargeting with bispecific anti-myosin-anti-DTPA-Fab-Fab’ and targeting with high-specific radioactivity Tc-99m-DTPA-succinylated-polylysine (DSPL) was developed. Methods: Mice were injected biweekly with 10 mg/kg Dox or its equivalent as D-Dox-PGA. Tc-99m-DSPL myocardial activity after pretargeting with bsAb-Fab-Fab’ was determined after gamma imaging performed at day 7 for Dox-treated mice and day 39 for all others. Results: Mice treated with 10 mg/kg Dox lost 10% total body weight in 1 week and 20% after a second dose. Pretargeted mice treated with 30 mg/kg cumulative D-Dox-PGA dose showed no loss of body weight for the duration of the study. Cardiotoxicity was confirmed by gamma imaging and scintillation counting (1.9 ± 0.25 [mean% ID/g ± SD]) after 1 dose of Dox. Mice injected with 3 × 10 mg/kg Dox equivalent as D-Dox-PGA (0.4 ± 0.04, P < .01) and untreated 2 control groups (0.20 ± 0.05 and 0.19 ± 0.04, P < .01) showed significantly lower myocardial anti-myosin radioactivity relative to the 10 mg/kg Dox group. Conclusion: Pretargeting with bsAb-Fab-Fab’ and targeting with Tc-99m labeled high-specific activity polymers enabled early visualization of doxorubicin induce cardiotoxicity in mice. Tolerated dose of D-Dox-PGA was greater than to 30 mg/kg Dox-equivalent dose with minimal cardiotoxicity.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

In-vitro and in-vivo imaging of MMP activity in cartilage and joint injury

Non-destructive detection of cartilage-degrading activities represents an advance in osteoarthritis (OA) research, with implications in studies of OA pathogenesis, progression, and intervention strategies. Matrix metalloproteinases (MMPs) are principal cartilage degrading enzymes that contribute to OA pathogenesis. MMPSense750 is an in-vivo fluorimetric imaging probe with the potential to continuously and non-invasively trace real-time MMP activities, but its use in OA-related research has not been reported. Our objective is to detect and characterize the early degradation activities shortly after cartilage or joint injury with MMPSense750. We determined the appropriate concentration, assay time, and linear range using various concentrations of recombinant MMPs as standards. We then quantified MMP activity from cartilage explants subjected to either mechanical injury or inflammatory cytokine treatment in-vitro. Finally, we performed in-vivo MMP imaging of a mouse model of post-traumatic OA. Our in-vitro results showed that the optimal assay time was highly dependent on the MMP enzyme. In cartilage explant culture media, mechanical impact or cytokine treatment increased MMP activity. Injured knees of mice showed significantly higher fluorescent signal than uninjured knees. We conclude that MMPSense750 detects human MMP activities and can be used for in-vitro study with cartilage, as well as in-vivo studies of knee injury, and can offering real-time insight into the degradative processes that occurring within the joint before structural changes become evident radiographically. - Highlights: • MMPSense750 is near-infrared fluorescent probe which can detect MMP activity. • MMPSense750 can detect human MMP-3, -9, and -13. • The reaction kinetics with MMPSense750 were different for the three MMPs. • MMPSense750 can visualized real time MMP activity in mouse injured knees. • MMPSense750 is convenient tool to evaluate real-time MMP activity non-invasively.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1963

THIRTY DAY ENDURANCE TEST OF A 30KW ARC JET ENGINE

Details are presented on a 30-day run of a radiation cooled thermal arc jet engine. The engine performed satisfactorily throughout the test period, maintaining a thrust of 250 grams and a propellant specific impulse greater than 1000 seconds a t an input power level of 30 kilowatts. The resulting propulsive efficiency was greater than 40 percent. Visual examination of the arc discharge through a quartz constrictor indicates that there is a smooth, constant diameter, bright core extending from the cathode completely through the constrictor to the cathode. Photographic observations on the arc discharge, using a red filter, indicate that the center of the column is darker than the outside: and this, in turn, suggests that the maximum core temperature is greater than 20,000°K. Calculations on the energy flow rate through a fully developed laminar column, combined with the experimental observation that the longitudinal energy flux gradient along the column is small, suggest that there may be a region of high volumetric energy release in the vicinity of the cathode tip.

Arc jet↗
NASA NTRSDate not supplied

SIRIUS: Simulation Infrastructure for Research on Interoperating Unmanned Systems

This paper presents the Simulation Infrastructure for Research on Interoperating Unmanned Systems (SIRIUS), a research framework for simulation and analysis of future conceptual Urban Air Mobility (UAM) operations. SIRIUS is being developed under the auspices of the NASA Air Traffic Management eXploration project, UAM subproject (ATM-X UAM). SIRIUS provides an intuitive, highly configurable graphical user interface to design complex traffic scenarios and airspace configurations representative of conceptual UAM operations. Aircraft simulated with SIRIUS can be equipped with flight-tested capabilities for detect and avoid (DAA), geofencing, distributed merging and spacing, path conformance, and path planning while executing time-constrained, 4D trajectories generated by a UAM ground operations system. Central to the design of the SIRIUS simulation framework is the capability to evaluate the integration and interoperability of ground-based separation services (e.g., strategic separation) with extended DAA functionality (e.g., path monitoring, separation provision, merging and spacing, etc.) The simulation environment also supports modelling of wind, navigation, and sensor uncertainties, as well as communication delays. SIRIUS enables distributed simulation of large-scale scenarios. An interactive graphical analysis capability helps isolate, visualize, and compare relevant vehicle state data and widely used measures of performance metrics across multiple scenarios.

Andrew Peters↗
DOE OSTI2025

Collaborative: in situ visual analytics technologies for extreme scale combustion simulations

This project aims to drastically enhance the usability of in situ analysis and visualization for extreme-scale scientific simulations. Current exascale computing capabilities promise to offer greater predictive ability of simulations and to further push the frontiers of science and technology. However, to validate the simulation output at extreme scale, examine the modeled phenomena, and discover previously unknowns from the output data, the output must be reduced or transformed in situ as it is being generated during the simulation such that the amount of data to examine and store is kept to a minimum. Such in situ approaches allow us to process and analyze the data and any embedded geometry to an extent that would be prohibitively expensive, if not impossible, to perform as a post hoc task. While in situ processing has been demonstrated to be a feasible and promising approach, its full potential has not yet been leveraged. In this project, we have developed comprehensive enhancements to in situ technology based on probability distributions in data. Our research focuses on jointly developing new ways of interacting with massive statistical samples while creatively utilizing new state-of-the-art computational resources to push the boundaries of in situ exploration. Moreover, we have developed new time-dependent techniques to enable previously unattainable capabilities in areas such as intelligent simulation steering and precise feature identification. We have experimentally studied our design and implementation at NERSC and OLCF, and are able to leverage existing in situ infrastructures whenever possible. While the exemplar in this project is combustion, many other fields for which turbulent transport is important, e.g., fusion, climate, astrophysics among others, encounter similar issues as simulations scale up to the exascale. This project shows its potential to generate high impact on DOE missions since the resulting technology promises to improve scientists’ ability to rapidly and correctly interpret and tune extreme-scale simulations, leading to new scientific understanding and advancements.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1998

A Geometry Based Infra-Structure for Computational Analysis and Design

The computational steps traditionally taken for most engineering analysis suites (computational fluid dynamics (CFD), structural analysis, heat transfer and etc.) are: (1) Surface Generation -- usually by employing a Computer Assisted Design (CAD) system; (2) Grid Generation -- preparing the volume for the simulation; (3) Flow Solver -- producing the results at the specified operational point; (4) Post-processing Visualization -- interactively attempting to understand the results. For structural analysis, integrated systems can be obtained from a number of commercial vendors. These vendors couple directly to a number of CAD systems and are executed from within the CAD Graphical User Interface (GUI). It should be noted that the structural analysis problem is more tractable than CFD; there are fewer mesh topologies used and the grids are not as fine (this problem space does not have the length scaling issues of fluids). For CFD, these steps have worked well in the past for simple steady-state simulations at the expense of much user interaction. The data was transmitted between phases via files. In most cases, the output from a CAD system could go to Initial Graphics Exchange Specification (IGES) or Standard Exchange Program (STEP) files. The output from Grid Generators and Solvers do not really have standards though there are a couple of file formats that can be used for a subset of the gridding (i.e. PLOT3D data formats). The user would have to patch up the data or translate from one format to another to move to the next step. Sometimes this could take days. Specifically the problems with this procedure are:(1) File based -- Information flows from one step to the next via data files with formats specified for that procedure. File standards, when they exist, are wholly inadequate. For example, geometry from CAD systems (transmitted via IGES files) is defined as disjoint surfaces and curves (as well as masses of other information of no interest for the Grid Generator). This is particularly onerous for modern CAD systems based on solid modeling. The part was a proper solid and in the translation to IGES has lost this important characteristic. STEP is another standard for CAD data that exists and supports the concept of a solid. The problem with STEP is that a solid modeling geometry kernel is required to query and manipulate the data within this type of file. (2) 'Good' Geometry. A bottleneck in getting results from a solver is the construction of proper geometry to be fed to the grid generator. With 'good' geometry a grid can be constructed in tens of minutes (even with a complex configuration) using unstructured techniques. Adroit multi-block methods are not far behind. This means that a million node steady-state solution can be computed on the order of hours (using current high performance computers) starting from this 'good' geometry. Unfortunately, the geometry usually transmitted from the CAD system is not 'good' in the grid generator sense. The grid generator needs smooth closed solid geometry. It can take a week (or more) of interaction with the CAD output (sometimes by hand) before the process can begin. One way Communication. (3) One-way Communication -- All information travels on from one phase to the next. This makes procedures like node adaptation difficult when attempting to add or move nodes that sit on bounding surfaces (when the actual surface data has been lost after the grid generation phase). Until this process can be automated, more complex problems such as multi-disciplinary analysis or using the above procedure for design becomes prohibitive. There is also no way to easily deal with this system in a modular manner. One can only replace the grid generator, for example, if the software reads and writes the same files. Instead of the serial approach to analysis as described above, CAPRI takes a geometry centric approach. This makes the actual geometry (not a discretized version) accessible to all phases of the analysis. The connection to the geometry is made through an Application Programming Interface (API) and NOT a file system. This API isolates the top-level applications (grid generators, solvers and visualization components) from the geometry engine. Also this allows the replacement of one geometry kernel with another, without effecting these top-level applications. For example, if UniGraphics is used as the CAD package then Parasolid (UG's own geometry engine) can be used for all geometric queries so that no solid geometry information is lost in a translation. This is much better than STEP because when the data is queried, the same software is executed as used in the CAD system. Therefore, one analyzes the exact part that is in the CAD system. CAPRI uses the same idea as the commercial structural analysis codes but does not specify control. Software components of the CAD system are used, but the analysis suite, not the CAD operator, specifies the control of the software session. This also means that the license issues (may be) minimized and individuals need not have to know how to operate a CAD system in order to run the suite.

Haimes, Robert↗
DOE OSTI2026

Holistic energy analysis method for thermal management architectures of data centers

Modern high-performance computing (HPC) data centers (DCs), particularly those supporting energy-intensive artificial intelligence (AI) workloads, face escalating thermal management challenges that degrade performance through thermal throttling and drive up cooling power consumption and operational costs. To address this challenge, many have developed a wide variety of thermal management solutions (single-phase, two-phase, direct, indirect, hybrid, and more) which attempt to cool HPC DCs effectively while attempting to minimize overall system power consumption. However, the analysis of these solutions and methods to effectively compare one with another is lacking. Overall power usage effectiveness (PUE) and total-power usage effectiveness (TUE) provide a metric to quantify power consumption but fail to identify components in the system which require further optimization. To address this, we propose a holistic analytical framework – the waterfall diagram (WFD) – which leverages a waterfall chart methodology, offering a comprehensive visualization of both the thermal management system loop and heat flow pathways from individual server components to the outdoor ambient. Use of the WFD enables graphical estimations of power efficiency and cooling performance across each component of a DC cooling system and complements Sankey-style energy flow visualizations by additionally resolving stage-wise temperature changes and incremental TUE contributions. The framework is used in conjunction with simulation-based approaches, to conduct a detailed pressure drop and flow distribution analysis aimed at identifying the optimal coolant distribution architecture for a single-phase direct-to-chip water-cooled DC, which serves as the baseline for subsequent WFD analysis. Among the evaluated architectures, the 3 U modular coolant distribution architecture is found to demonstrate the best performance, considering minimal pressure drop and uniform flow distribution. In addition, TUE is calculated for each cooling loop component based on its associated pressure drop and corresponding pumping power, which are integrated into the WFD. This correlation between TUE and local temperature offers immediate insight into the power efficiency and thermal performance contributions of individual components, facilitating further development and optimization. Examples of WFD applications are presented under varying thermal loads and ambient conditions, demonstrating reasonable cooling strategies. Notably, the 3 U modular architecture maintains a consistent chip case temperature of 85°C, achieving a TUE of 1.016 at ambient temperature of 47°C, and a TUE of 1.026 at ambient temperature of 52°C. The WFD methodology provides an efficient, holistic, and streamlined framework for DC thermal management architecture assessment and enables design optimization which is important for addressing the thermal-fluidic energy challenges of current and next-generation DCs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2018

Structures of the intermediates of Kok’s photosynthetic water oxidation clock

Inspired by the period-four oscillation in flash-induced oxygen evolution of photosystem II discovered by Joliot in 1969, Kok performed additional experiments and proposed a five-state kinetic model for photosynthetic oxygen evolution, known as Kok’s S-state clock or cycle. The model comprises four (meta)stable intermediates (S 0 , S 1 , S 2 and S 3 ) and one transient S 4 state, which precedes dioxygen formation occurring in a concerted reaction from two water-derived oxygens bound at an oxo-bridged tetra manganese calcium (Mn 4 CaO 5 ) cluster in the oxygen-evolving complex. This reaction is coupled to the two-step reduction and protonation of the mobile plastoquinone QB at the acceptor side of PSII. Here in this paper, using serial femtosecond X-ray crystallography and simultaneous X-ray emission spectroscopy with multi-flash visible laser excitation at room temperature, we visualize all (meta)stable states of Kok’s cycle as high-resolution structures (2.04–2.08 Å). In addition, we report structures of two transient states at 150 and 400 µs, revealing notable structural changes including the binding of one additional ‘water’, Ox, during the S 2 →S 3 state transition. Our results suggest that one water ligand to calcium (W3) is directly involved in substrate delivery. The binding of the additional oxygen Ox in the S 3 state between Ca and Mn1 supports O–O bond formation mechanisms involving O5 as one substrate, where Ox is either the other substrate oxygen or is perfectly positioned to refill the O5 position during O 2 release. Thus, our results exclude peroxo-bond formation in the S 3 state, and the nucleophilic attack of W3 onto W2 is unlikely.

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

Accelerating Multigrid-based Hierarchical Scientific Data Refactoring on GPUs

Rapid growth in scientific data and a widening gap between computational speed and I/O bandwidth make it increasingly infeasible to store and share all data produced by scientific simulations. Instead, we need methods for reducing data volumes: ideally, methods that can scale data volumes adaptively so as to enable negotiation of performance and fidelity tradeoffs in different situations. Multigrid-based hierarchical data representations hold promise as a solution to this problem, allowing for flexible conversion between different fidelities so that, for example, data can be created at high fidelity and then transferred or stored at lower fidelity via logically simple and mathematically sound operations. However, the effective use of such representations has been hindered until now by the relatively high costs of creating, accessing, reducing, and otherwise operating on such representations. We describe here highly optimized data refactoring kernels for GPU accelerators that enable efficient creation and manipulation of data in multigrid-based hierarchical forms. We demonstrate that our optimized design can achieve up to 250 TB/s aggregated data refactoring throughput—83% of theoretical peak—on 1024 nodes of the Summit supercomputer. We showcase our optimized design by applying it to a large-scale scientific visualization workflow and the MGARD lossy compression software.

Data refactoring↗
DOE OSTI2013

Quantitative Monte Carlo-based holmium-166 SPECT reconstruction

Purpose: Quantitative imaging of the radionuclide distribution is of increasing interest for microsphere radioembolization (RE) of liver malignancies, to aid treatment planning and dosimetry. For this purpose, holmium-166 ({sup 166}Ho) microspheres have been developed, which can be visualized with a gamma camera. The objective of this work is to develop and evaluate a new reconstruction method for quantitative {sup 166}Ho SPECT, including Monte Carlo-based modeling of photon contributions from the full energy spectrum.Methods: A fast Monte Carlo (MC) simulator was developed for simulation of {sup 166}Ho projection images and incorporated in a statistical reconstruction algorithm (SPECT-fMC). Photon scatter and attenuation for all photons sampled from the full {sup 166}Ho energy spectrum were modeled during reconstruction by Monte Carlo simulations. The energy- and distance-dependent collimator-detector response was modeled using precalculated convolution kernels. Phantom experiments were performed to quantitatively evaluate image contrast, image noise, count errors, and activity recovery coefficients (ARCs) of SPECT-fMC in comparison with those of an energy window-based method for correction of down-scattered high-energy photons (SPECT-DSW) and a previously presented hybrid method that combines MC simulation of photopeak scatter with energy window-based estimation of down-scattered high-energy contributions (SPECT-ppMC+DSW). Additionally, the impact of SPECT-fMC on whole-body recovered activities (A{sup est}) and estimated radiation absorbed doses was evaluated using clinical SPECT data of six {sup 166}Ho RE patients.Results: At the same noise level, SPECT-fMC images showed substantially higher contrast than SPECT-DSW and SPECT-ppMC+DSW in spheres ≥17 mm in diameter. The count error was reduced from 29% (SPECT-DSW) and 25% (SPECT-ppMC+DSW) to 12% (SPECT-fMC). ARCs in five spherical volumes of 1.96–106.21 ml were improved from 32%–63% (SPECT-DSW) and 50%–80% (SPECT-ppMC+DSW) to 76%–103% (SPECT-fMC). Furthermore, SPECT-fMC recovered whole-body activities were most accurate (A{sup est}= 1.06 × A − 5.90 MBq, R{sup 2}= 0.97) and SPECT-fMC tumor absorbed doses were significantly higher than with SPECT-DSW (p = 0.031) and SPECT-ppMC+DSW (p = 0.031).Conclusions: The quantitative accuracy of {sup 166}Ho SPECT is improved by Monte Carlo-based modeling of the image degrading factors. Consequently, the proposed reconstruction method enables accurate estimation of the radiation absorbed dose in clinical practice.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2015

SU-E-J-275: Review - Computerized PET/CT Image Analysis in the Evaluation of Tumor Response to Therapy

Purpose: To review the literature in using computerized PET/CT image analysis for the evaluation of tumor response to therapy. Methods: We reviewed and summarized more than 100 papers that used computerized image analysis techniques for the evaluation of tumor response with PET/CT. This review mainly covered four aspects: image registration, tumor segmentation, image feature extraction, and response evaluation. Results: Although rigid image registration is straightforward, it has been shown to achieve good alignment between baseline and evaluation scans. Deformable image registration has been shown to improve the alignment when complex deformable distortions occur due to tumor shrinkage, weight loss or gain, and motion. Many semi-automatic tumor segmentation methods have been developed on PET. A comparative study revealed benefits of high levels of user interaction with simultaneous visualization of CT images and PET gradients. On CT, semi-automatic methods have been developed for only tumors that show marked difference in CT attenuation between the tumor and the surrounding normal tissues. Quite a few multi-modality segmentation methods have been shown to improve accuracy compared to single-modality algorithms. Advanced PET image features considering spatial information, such as tumor volume, tumor shape, total glycolytic volume, histogram distance, and texture features have been found more informative than the traditional SUVmax for the prediction of tumor response. Advanced CT features, including volumetric, attenuation, morphologic, structure, and texture descriptors, have also been found advantage over the traditional RECIST and WHO criteria in certain tumor types. Predictive models based on machine learning technique have been constructed for correlating selected image features to response. These models showed improved performance compared to current methods using cutoff value of a single measurement for tumor response. Conclusion: This review showed that computerized PET/CT image analysis holds great potential to improve the accuracy in evaluation of tumor response. This work was supported in part by the National Cancer Institute Grant R01CA172638.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TU-H-206-01: An Automated Approach for Identifying Geometric Distortions in Gamma Cameras

Purpose: To develop a clinically-deployable, automated process for detecting artifacts in routine nuclear medicine (NM) quality assurance (QA) bar phantom images. Methods: An artifact detection algorithm was created to analyze bar phantom images as part of an ongoing QA program. A low noise, high resolution reference image was acquired from an x-ray of the bar phantom with a Philips Digital Diagnost system utilizing image stitching. NM bar images, acquired for 5 million counts over a 512×512 matrix, were registered to the template image by maximizing mutual information (MI). The MI index was used as an initial test for artifacts; low values indicate an overall presence of distortions regardless of their spatial location. Images with low MI scores were further analyzed for bar linearity, periodicity, alignment, and compression to locate differences with respect to the template. Findings from each test were spatially correlated and locations failing multiple tests were flagged as potential artifacts requiring additional visual analysis. The algorithm was initially deployed for GE Discovery 670 and Infinia Hawkeye gamma cameras. Results: The algorithm successfully identified clinically relevant artifacts from both systems previously unnoticed by technologists performing the QA. Average MI indices for artifact-free images are 0.55. Images with MI indices < 0.50 have shown 100% sensitivity and specificity for artifact detection when compared with a thorough visual analysis. Correlation of geometric tests confirms the ability to spatially locate the most likely image regions containing an artifact regardless of initial phantom orientation. Conclusion: The algorithm shows the potential to detect gamma camera artifacts that may be missed by routine technologist inspections. Detection and subsequent correction of artifacts ensures maximum image quality and may help to identify failing hardware before it impacts clinical workflow. Going forward, the algorithm is being deployed to monitor data from all gamma cameras within our health system.

60 APPLIED LIFE SCIENCES↗