Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Program Optimization”

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,675 records · Page 93

Branch Flow Model: Relaxations and Convexification-Part I

We propose a branch flow model for the analysis and optimization of mesh as well as radial networks. The model leads to a new approach to solving optimal power flow (OPF) that consists of two relaxation steps. The first step eliminates the voltage and current angles and the second step approximates the resulting problem by a conic program that can be solved efficiently. For radial networks, we prove that both relaxation steps are always exact, provided there are no upper bounds on loads. For mesh networks, the conic relaxation is always exact but the angle relaxation may not be exact, and we provide a simple way to determine if a relaxed solution is globally optimal. We propose convexification of mesh networks using phase shifters so that OPF for the convexified network can always be solved efficiently for an optimal solution. We prove that convexification requires phase shifters only outside a spanning tree of the network and their placement depends only on network topology, not on power flows, generation, loads, or operating constraints. Part I introduces our branch flow model, explains the two relaxation steps, and proves the conditions for exact relaxation. Part II describes convexification of mesh networks, and presents simulation results.

Farivar, M↗

Branch Flow Model: Relaxations and Convexification-Part II

We propose a branch flow model for the analysis and optimization of mesh as well as radial networks. The model leads to a new approach to solving optimal power flow (OPF) that consists of two relaxation steps. The first step eliminates the voltage and current angles and the second step approximates the resulting problem by a conic program that can be solved efficiently. For radial networks, we prove that both relaxation steps are always exact, provided there are no upper bounds on loads. For mesh networks, the conic relaxation is always exact but the angle relaxation may not be exact, and we provide a simple way to determine if a relaxed solution is globally optimal. We propose convexification of mesh networks using phase shifters so that OPF for the convexified network can always be solved efficiently for an optimal solution. We prove that convexification requires phase shifters only outside a spanning tree of the network and their placement depends only on network topology, not on power flows, generation, loads, or operating constraints. Part I introduces our branch flow model, explains the two relaxation steps, and proves the conditions for exact relaxation. Part II describes convexification of mesh networks, and presents simulation results.

Farivar, M↗

Optimizing Power Line Undergrounding Decisions under Varying Wildfire Risk and Weather Scenarios

Abstract—The threat of wildfire ignitions from electric power equipment has led utilities to increasingly turn to preemptive power shutoffs, which, while effective in reducing grid-induced wildfire risk, can cause significant load loss. Undergrounding power lines is an alternative strategy for preventing grid-induced wildfires. However, undergrounding lines is costly, so an efficient undergrounding plan must balance reductions in wildfire risk and load loss with the cost of undergrounding lines. We propose a robust optimization model to identify which power lines to underground to maximize load served while limiting wildfire risk across a range of wildfire risk and weather scenarios. Since solving this problem may be computationally heavy for large power grids and many operating scenarios, we present a delayed constraint generation algorithm to iteratively add scenarios until an optimal solution is found. We evaluate the performance of this framework on the RTS-GMLC with scenarios representing a year of operating conditions and compare it with a stochastic programming formulation. Our results indicate that our undergrounding model is successful in reducing load shed and risk compared to baseline cases in which no mitigation action is taken and only power shutoffs are implemented (no undergrounding). The robust formulation also reduces more load shed than the stochastic formulation in the most extreme scenarios. Index Terms—grid resilience, optimization, transmission systems, underground power lines, wildfire risk.

Taylor, S. [Department of Electrical and Computer ↗

Fully Distributed Acoustic and Magnetic Field Monitoring Via a Single Fiber Line for Optimized Production of Unconventional Resource Plays

This is the final technical report for the Virginia Tech Center for Photonics Technology (VT-CPT) research project entitled “Fully Distributed Acoustic and Magnetic Field Monitoring Via a Single Fiber Line for Optimized Production of Unconventional Resource Plays”. In the the four-year effort, VT-CPT collaborated with the Departments of Mathmatics at VT and Sentek Instrument, LLC (Sentek) to develop a fiber-optic sensing system capable of real-time simultaneous distributed measurement of multiple subsurface, drilling, and production parameters. The ultra-sensitive fiber optic distributed acoustic sensing technology, picoDAS, developed by Sentek was integrated with a novel multi-material optical sensing fiber fabricated by CPT-VT to obtain distributed acoustic and magnetic field measurements with ultrahigh sensitivity and high spatial resolution. The sensing technology successfully developed and demonstrated under this research program is truly unique has application beyond subsurface imaging to include carbon storage and electrical grid monitoring, as well as healthcare and nuclear fusion. The project, sponsored by the Advanced Technology Program, began in 2019 and had the original goal of developing the next generation of harsh environment sensing systems. Exhaustive theoretical modeling was performed to optimize the multi-material optical fiber design. VT developed the processing techniques employed to successfully fabricate single mode fibers with metal cladding wires. Sensing fibers with nickel and Metglas cladding wires were fully characterized, inscribed with sensors, and fully integrated with the picoDAS system. Laboratory scale testing was performed to demonstrate performance upon exposure to lateral and transverse magnetic fields. Preliminary field trials were performed for prototype magnetic and acoustic sensing systems upon near surface deployment at VT. The distributed acoustic and magnetic field sensing system was advanced from a (Technology Readiness Level) TRL=2 to a TRL=5 via laboratory scale system validation in relevant environments. The truly one-of-its kind distributed magnetic and acoustic sensing system is expected to find immediate applications that require subsurface monitoring, to include carbon storage site monitoring, and has the potential to be a distributive technology in healthcare and electric grid markets. The technologies developed by Virginia Tech’s Center for Photonics Technology will support the mission of the National Energy Technology Laboratory (NETL) to drive innovation and deliver solutions for an environmentally sustainable and prosperous energy future by ensuring affordable, abundant, and reliable energy that drives a robust economy and national security. Technical accomplishments during the program are summarized briefly here and described in detail in the remainder of the report.

02 PETROLEUM↗

Progress towards developing neutron tolerant magnetostrictive and piezoelectric transducers

Current generation light water reactors (LWRs), sodium cooled fast reactors (SFRs), small modular reactors (SMRs), and next generation nuclear plants (NGNPs) produce harsh environments in and near the reactor core that can severely tax material performance and limit component operational life. To address this issue, several Department of Energy Office of Nuclear Energy (DOE-NE) research programs are evaluating the long duration irradiation performance of fuel and structural materials used in existing and new reactors. In order to maximize the amount of information obtained from Material Testing Reactor (MTR) irradiations, DOE is also funding development of enhanced instrumentation that will be able to obtain in-situ, real-time data on key material characteristics and properties, with unprecedented accuracy and resolution. Such data are required to validate new multi-scale, multi-physics modeling tools under development as part of a science-based, engineering driven approach to reactor development. It is not feasible to obtain high resolution/microscale data with the current state of instrumentation technology. However, ultrasound-based sensors offer the ability to obtain such data if it is demonstrated that these sensors and their associated transducers are resistant to high neutron flux, high gamma radiation, and high temperature. To address this need, the Advanced Test Reactor National Scientific User Facility (ATR-NSUF) is funding an irradiation, led by PSU, at the Massachusetts Institute of Technology Research Reactor to test the survivability of ultrasound transducers. As part of this effort, PSU and collaborators have designed, fabricated, and provided piezoelectric and magnetostrictive transducers that are optimized to perform in harsh, high flux, environments. Four piezoelectric transducers were fabricated with either aluminum nitride, zinc oxide, or bismuth titanate as the active element that were coupled to either Kovar or aluminum waveguides and two magnetostrictive transducers were fabricated with Remendur or Galfenol as the active elements. Pulse-echo ultrasonic measurements of these transducers are made in-situ. This paper will present an overview of the test design including selection criteria for candidate materials and optimization of test assembly parameters, data obtained from both out-of-pile and in-pile testing at elevated temperatures, and an assessment based on initial data of the expected performance of ultrasonic devices in irradiation conditions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

TU-CD-304-01: FEATURED PRESENTATION and BEST IN PHYSICS (THERAPY): Trajectory Modulated Arc Therapy: Development of Novel Arc Delivery Techniques Integrating Dynamic Table Motion for Extended Volume Treatments

Purpose: Integration of coordinated robotic table motion with inversely-planned arc delivery has the potential to resolve table-top delivery limitations of large-field treatments such as Total Body Irradiation (TBI), Total Lymphoid Irradiation (TLI), and Cranial-Spinal Irradiation (CSI). We formulate the foundation for Trajectory Modulated Arc Therapy (TMAT), and using Varian Developer Mode capabilities, experimentally investigate its practical implementation for such techniques. Methods: A MATLAB algorithm was developed for inverse planning optimization of the table motion, MLC positions, and gantry motion under extended-SSD geometry. To maximize the effective field size, delivery trajectories for TMAT TBI were formed with the table rotated at 270° IEC and dropped vertically to 152.5cm SSD. Preliminary testing of algorithm parameters was done through retrospective planning analysis. Robotic delivery was programmed using custom XML scripting on the TrueBeam Developer Mode platform. Final dose was calculated using the Eclipse AAA algorithm. Initial verification of delivery accuracy was measured using OSLDs on a solid water phantom of varying thickness. Results: A comparison of DVH curves demonstrated that dynamic couch motion irradiation was sufficiently approximated by static control points spaced in intervals of less than 2cm. Optimized MLC motion decreased the average lung dose to 68.5% of the prescription dose. The programmed irradiation integrating coordinated table motion was deliverable on a TrueBeam STx linac in 6.7 min. With the couch translating under an open 10cmx20cm field angled at 10°, OSLD measurements along the midline of a solid water phantom at depths of 3, 5, and 9cm were within 3% of the TPS AAA algorithm with an average deviation of 1.2%. Conclusion: A treatment planning and delivery system for Trajectory Modulated Arc Therapy of extended volumes has been established and experimentally demonstrated for TBI. Extension to other treatment techniques such as TLI and CSI is readily achievable through the developed platform. Grant Funding by Varian Medical Systems.

60 APPLIED LIFE SCIENCES↗

Modeling and optimization of steady flow of natural gas and hydrogen mixtures in pipeline networks

Here, we extend the canonical problems of simulation and optimization of steady-state gas flows in pipeline networks with compressors to the transport of mixtures of highly heterogeneous gases injected throughout a network. Our study is motivated by proposed projects to blend hydrogen generated using clean energy into existing natural gas pipeline systems as part of efforts to reduce the reliance of energy systems on fossil fuels. Flow in a pipe is related to endpoint pressures by a basic Weymouth equation model, with an ideal gas equation of state, where the wave speed depends on the hydrogen concentration. At vertices, in addition to mass balance, we also consider mixing of incoming flows of varying hydrogen concentrations. The problems of interest are the heterogeneous gas flow simulation (HGFS), which determines system pressures and flows given fixed boundary conditions and compressor settings, as well as the heterogeneous gas flow optimization (HGFO), which extremizes an objective by determining optimal boundary conditions and compressor settings. We examine conditions for uniqueness of solutions to the HGFS, as well as compare and contrast mixed-integer and continuous nonlinear programming formulations for the HGFO. We develop computational methods to solve both problems, and examine their performance using four test networks of increasing complexity.

08 HYDROGEN↗

Smart Composite Pressure Vessels (SCPV) with Integrated Health Monitoring

The U.S. Department of Energy (DOE) is promoting and developing more energy efficient and environmentally friendly technologies that will enable America to use less petroleum. Hydrogen fuel cells – which directly convert the chemical energy in hydrogen to electricity with only water and heat as byproducts – are a very attractive solution that can enable this to happen. To date, DOE's efforts have culminated in commercial demonstration of on-board vehicular hydrogen storage systems that can allow for a driving range of greater than 300 miles. This requires storing 5 Kgs of hydrogen onboard a light vehicle. In order to store this quantity of hydrogen, the hydrogen gas needs to be stored in a composite overwrapped pressure vessel (COPV) at a very high internal pressure. The objective of the program was to demonstrate continuous and predictable health-monitoring of composite pressure vessels. A higher confidence in the operational safety of the vessels will lead to reduced burst factor of safety imposed by regulatory standards and hence will help reduce the structural wall thickness that drives the cost of the vessel. Highly optimized Type III (metal lined) and Type IV (polymer lined) COPVs that are manufactured using filament winding process have been designed and qualified for the above mentioned application. However, the main structural component of the COPVs is carbon fiber and the high cost of carbon fiber composite in a pressure vessel is a primary challenge in reducing the cost of gaseous hydrogen storage. Continuous and remote monitoring of structural health of the COPVs as well as optimization of the strength translation of carbon fibers have the potential to allow for reduced factor of safety, thereby reduced amount of carbon and associated cost of the vessel. In the current program, the project team including Steelhead Composites (SHC), University of Tennessee, Knoxville (UTK), Teijin Carbon America, Oak Ridge National Laboratory (ORNL) and LUNA designed, fabricated, and tested smart composite pressure vessels with integrated sensors. Iterative loops of design and testing using coupons, subscale STEB vessels and full size vessels proved that a new generation of high performance carbon fiber can be used for efficient design of Type 3 vessels for H2 storage. Novel analysis techniques were developed to predict the initiation and propagation of interlaminar damage inside the composite shell due to impact damage, a real threat in practical operation. Such an analysis scheme is typically not used in tank design but is an essential tool for health monitoring of composite vessels. Remote sensing of the key signatures of the tank’s operating parameters such as pressure, acceleration and humidity was demonstrated using a unique device that can wirelessly transmit and stream the data to a remote server. Fiber optic sensors were successfully integrated during fabrication of the tanks. These sensors provide a wealth of information regarding the structural health of the vessel when it is subsequently pressurized or subjected to impact damage. Excellent correlation was demonstrated between the fiber optic sensor and mechanical strain gage data, and between the measurements and analytical predictions.

08 HYDROGEN↗

MO-F-CAMPUS-I-03: Tissue Equivalent Material Phantom to Test and Optimize Coherent Scatter Imaging for Tumor Classification

Purpose: To accurately model different breast geometries using a tissue equivalent phantom, and to classify these tissues in a coherent x-ray scatter imaging system. Methods: A breast phantom has been designed to assess the capability of coded aperture coherent x-ray scatter imaging system to classify different types of breast tissue (adipose, fibroglandular, tumor). The tissue-equivalent phantom was modeled as a hollow plastic cylinder containing multiple cylindrical and spherical inserts that can be positioned, rearranged, or removed to model different breast geometries. Each enclosure can be filled with a tissue-equivalent material and excised human tumors. In this study, beef and lard, placed inside 2-mm diameter plastic Nalgene containers, were used as surrogates for fibroglandular and adipose tissue, respectively. The phantom was imaged at 125 kVp, 40 mA for 10 seconds each with a 1-mm pencil beam. The raw data were reconstructed using a model-based reconstruction algorithm and yielded the location and form factor, or momentum transfer (q) spectrum of the materials that were imaged. The measured material form factors were then compared to the ground truth measurements acquired by x-ray diffraction (XRD) imaging. Results: The tissue equivalent phantom was found to accurately model different types of breast tissue by qualitatively comparing our measured form factors to those of adipose and fibroglandular tissue from literature. Our imaging system has been able to define the location and composition of the various materials in the phantom. Conclusion: This work introduces a new tissue equivalent phantom for testing and optimization of our coherent scatter imaging system for material classification. In future studies, the phantom will enable the use of a variety of materials including excised human tissue specimens in evaluating and optimizing our imaging system using pencil- and fan-beam geometries. United States Department of Homeland Security Duke University Medical Center - Department of Radiology Carl E Ravin Advanced Imaging Laboratories Duke University Medical Physics Graduate Program.

60 APPLIED LIFE SCIENCES↗

Multi‐dimensional lock‐free arrays for multithreaded mode‐directed tabling in Prolog

Summary This work proposes a new design for the supporting data structures used to implement multithreaded tabling in Prolog systems. Tabling is an implementation technique that improves the expressiveness of traditional Prolog systems in dealing with recursion and redundant computations. Mode‐directed tabling is an extension to the tabling technique that supports the definition of alternative criteria for specifying how answers are aggregated, thus being very suitable for problems where the goal is to dynamically calculate optimal or selective answers. In this work, we leverage the intrinsic potential that mode‐directed tabling has to express dynamic programming problems by creating a new design that improves the representation of multi‐dimensional arrays in the context of multithreaded tabling. To do so, we introduce a new mode for indexing arguments in mode‐directed tabled evaluations, named dim , where each dim argument features a uni‐dimensional lock‐free array. Experimental results using well‐known dynamic programming problems on a 32‐core machine show that the new design introduces less overheads and clearly improves the execution time for sequential and multithreaded tabled evaluations.

Areias, Miguel↗

Constrained Bayesian optimization of criticality experiments

The design of criticality experiments is typically an iterative process that employs a Monte Carlo transport code. The goal is to find a design that optimizes some variable, like the sensitivity of a response to a cross section, while simultaneously ensuring criticality. The high fidelity of the Monte Carlo code is a great asset, but it makes exploring the design space computationally expensive. Herein, we present how a constrained Bayesian optimization algorithm can be used to efficiently design a criticality experiment. It uses Gaussian processes as a surrogate model to probe the design space and to reduce the number of code executions that are needed to find the optimum. Furthermore, we demonstrate constrained Bayesian optimization with a Pu-239/polyethylene solution system and a TEX experiment that is designed for criticality safety validation of a nuclear waste model at the Hanford Site. For both systems, a global optimum was found within 75 Monte Carlo simulations.

42 ENGINEERING↗

Ten questions concerning occupant-centric control and operations

Occupant-Centric Control and Operation (OCC) represents a transformative approach to building management, integrating sensing of indoor environmental quality, occupant presence, and occupant-building interactions. These data are then utilized to optimize both operational efficiency and occupant comfort. This paper summarizes the findings from the IEA-EBC Annex 79 research program's subtask on real world implementations of OCC during the past 5 years. First, in Q1 and Q2, we provide a definition and categorization of OCC. Q3 addresses the role of building operators for OCC, while Q4 describes the implications for designers. Then, Q5 and Q6 discuss the role and possibilities of OCC for load flexibility, and for pandemic induced paradigm shifts in the built environment, respectively. In Q7, we provide a taxonomy and selection process of OCC, while Q8 details real world implementation case studies. Finally, Q9 explains the limits of OCC, and Q10 provides a vision for future research opportunities. Furthermore, our findings offer valuable insights for researchers, practitioners, and policy makers, contributing to the ongoing discourse on the future of building operations management.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A multi-criteria decision analysis-based approach for dispatch of electric microgrids

This paper presents a decision support system (DSS) for multi-objective dispatch of an electric microgrid con-sidering cost of operation, peak load reduction, and emissions. Discrete compromise programming (DCP) is used as the multi-criteria decision analysis (MCDA) technique for providing the decision support to the distribution system operator (DSO) by ranking various alternatives based on preference of the objec-tives. Simulations are performed on a notional electric microgrid based on the IEEE 13-node test feeder and the results are pre-sented with appropriate visualization.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Are We Appropriately Selecting Therapy For Patients With Cervical Cancer? Longitudinal Patterns-of-Care Analysis for Stage IB-IIB Cervical Cancer

Purpose: We performed a patterns-of-care analysis evaluating the effects of newer technology and recent research findings on treatment decisions over 26 years to determine whether patients with cervical cancer are being appropriately selected for treatment to optimize the therapeutic ratio. Methods and Materials: A retrospective analysis was conducted using the Surveillance, Epidemiology and End Results (SEER) program from 1983 to 2009. We identified 10,933 women with stage IB-IIB cervical carcinoma. Results: Of the 10,933 subjects identified, 40.1% received surgery, 26.8% received radiation (RT), and 33.1% received surgery plus RT. RT use increased after 2000 compared to prior to 2000, with a corresponding decrease in surgery and surgery plus RT. Among patients with risk factors including tumor size >4 cm, positive parametria, and positive lymph nodes, declining use of surgery plus RT was observed. However, 23% of patients with tumors >4 cm, 20% of patients with positive parametria, and 55% of node-positive patients continued to receive surgery plus RT as of 2009. Factors associated with increased use of surgery plus RT included patient age <50 and node-positive status. Conclusions: In this largest patterns-of-care analysis to date for patients with locally advanced cervical cancer, we found a substantial proportion of patients continue to undergo surgery followed by radiation, despite randomized data supporting the use of definitive radiation therapy, with lower morbidity than surgery and radiation.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Memory-optimized shift operator alternating direction implicit finite difference time domain method for plasma

Through introducing the alternating direction implicit (ADI) technique and the memory-optimized algorithm to the shift operator (SO) finite difference time domain (FDTD) method, the memory-optimized SO-ADI FDTD for nonmagnetized collisional plasma is proposed and the corresponding formulae of the proposed method for programming are deduced. In order to further the computational efficiency, the iteration method rather than Gauss elimination method is employed to solve the equation set in the derivation of the formulae. Complicated transformations and convolutions are avoided in the proposed method compared with the Z transforms (ZT) ADI FDTD method and the piecewise linear JE recursive convolution (PLJERC) ADI FDTD method. The numerical dispersion of the SO-ADI FDTD method with different plasma frequencies and electron collision frequencies is analyzed and the appropriate ratio of grid size to the minimum wavelength is given. The accuracy of the proposed method is validated by the reflection coefficient test on a nonmagnetized collisional plasma sheet. The testing results show that the proposed method is advantageous for improving computational efficiency and saving computer memory. The reflection coefficient of a perfect electric conductor (PEC) sheet covered by multilayer plasma and the RCS of the objects coated by plasma are calculated by the proposed method and the simulation results are analyzed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Toward performance-portable PETSc for GPU-based exascale systems

The Portable Extensible Toolkit for Scientific computation (PETSc) library delivers scalable solvers for nonlinear time-dependent differential and algebraic equations and for numerical optimization. The PETSc design for performance portability addresses fundamental GPU accelerator challenges and stresses flexibility and extensibility by separating the programming model used by the application from that used by the library, and it enables application developers to use their preferred programming model, such as Kokkos, RAJA, SYCL, HIP, CUDA, or OpenCL, on upcoming exascale systems. Furthermore, a blueprint for using GPUs from PETSc-based codes is provided, and case studies emphasize the flexibility and high performance achieved on current GPU-based systems.

97 MATHEMATICS AND COMPUTING↗

The role of rotation-vibration coupling in symmetric and asymmetric isotopomers of ozone

A theoretical framework and a computer code (SpectrumSDT) are developed for accurate calculations of coupled rotational–vibrational states in triatomic molecules using hyper-spherical coordinates and taking into account the Coriolis coupling effect. Concise final formulas are derived for the construction of the Hamiltonian matrix using an efficient combination of the variational basis representation and discrete variable representation methods with locally optimized basis sets and grids. First, the new code is tested by comparing its results with those of the APH3D program of Kendrick et al. [Kendrick, Pack, Walker, and Hayes, J. Chem. Phys. 110, 6673 (1999)]. Then, accurate calculations of the rovibrational spectra are carried out for doubly substituted symmetric ( 18 O 16 O 18 O) and asymmetric ( 18 O 18 O 16 O) ozone isotopomers for the total angular momentum up to J = 5. Together with similar data recently reported for the singly substituted symmetric ( 16 O 18 O 16 O) and asymmetric ( 16 O 16 O 18 O) ozone isotopomers, these calculations quantify the role of the Coriolis coupling effect in the large mass-independent isotopic enrichment of ozone, observed in both laboratory experiments and the atmosphere of the Earth. It is found that the Coriolis effect in ozone is relatively small, as evidenced by deviations of its rotational constants from the symmetric-top-rotor behavior, magnitudes of parity splittings (Λ-doubling), and ratios of rovibrational partition functions for asymmetric vs symmetric ozone molecules. It is concluded that all of these characteristics are influenced by the isotopic masses as much as they are influenced by the overall symmetry of the molecule. It is therefore unlikely that the Coriolis coupling effect could be responsible for symmetry-driven mass-independent fractionation of oxygen isotopes in ozone.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

GPU-Based Static Data-Flow Analysis for Fast and Scalable Android App Vetting

Many popular vetting tools for Android applications use static code analysis techniques. In particular, Interprocedural Data-Flow Graph (IDFG) construction is the computation at the core of Android static data-flow analysis and consumes most of the analysis time. Many analysis tools use a worklist algorithm, an iterative fixed-point approach, to construct the IDFG. In this paper, we observe that a straightforward GPU parallelization of the worklist algorithm leads to significant underutilization of the GPU resources. We identify four performance bottlenecks, namely, frequent dynamic memory allocations, high branch divergence, workload imbalance, and irregular memory access patterns. Accordingly, we propose GDroid, a GPU-based worklist algorithm implementation with multiple fine-grained optimizations tailored to common characteristics of Android applications. The optimizations considered are: matrix-based data structure, memory access-based node grouping, and worklist merging. Our experimental evaluation, performed on 1000 Android applications, shows that the proposed optimizations are beneficial to performance, and GDroid can achieve up to 128X speedups against a plain GPU implementation.

Android security↗