Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “programming”

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,171 records · Page 65

SU-F-J-116: Clinical Experience-Based Verification and Improvement of a 4DCT Program

Purpose: To demonstrate the role of continuous improvement fulfilled by the Medical Physicist in clinical 4DCT and CBCT scanning. Methods: Lung (SABR and Standard) patients’ 4D respiratory motion and image data were reviewed over a 3, 6 and 12 month period following commissioning testing. By identifying trends of clinically relevant parameters and respiratory motions, variables were tested with a programmable motion phantom and assessed. Patient traces were imported to a motion phantom and 4DCT and CBCT imaging were performed. Cos6 surrogate and sup-inf motion was also programmed into the phantom to simulate the long exhale of patients for image contrast tests. Results: Patient surrogate motion amplitudes were 9.9+5.2mm (3–35) at 18+6bpm (6–30). Expiration/Inspiration time ratios of 1.4+0.5second (0.6–2.9) showed image contrast effects evident in the AveCT and 3DCBCT images. Small differences were found for patients with multiple 4DCT data sets. Patient motion assessments were simulated and verified with the phantom within 2mm. Initial image reviews to check for reconstructed artefacts and data loss identified a small number of patients with irregularities in the automatic placement of inspiration and expiration points. Conclusion: The Physicist’s involvement in the continuous improvements of a clinically commissioned technique, processes and workflows continues beyond the commissioning stage of a project. Our experience with our clinical 4DCT program shows that Physics presence is required at the clinical 4DCT scan to assist with technical aspects of the scan and also for clinical image quality assessment prior to voluming. The results of this work enabled the sharing of information from the Medical Physics group with the Radiation Oncologists and Radiation Therapists. This results in an improved awareness of clinical patient respiration variables and how they may affect 4D simulation images and also may also affect the treatment verification images.

60 APPLIED LIFE SCIENCES↗

SU-G-201-01: An Automated Treatment Plan Quality Assurance Program for High-Dose Rate (HDR) Brachytherapy with a VaginalCylinder Applicator

Purpose: Plan specific quality assurance (QA) is an important step in high dose rate (HDR) brachytherapy to ensure the integrity of a treatment plan. The conventional approach is to assemble a set of plan screen-captures in a document and have an independent plan-checker to verify it. Not only is this approach cumbersome and time-consuming, using a document also limits the items that can be verified, hindering plan quality and patient safety. We have initiated efforts to develop a web-based HDR brachytherapy QA system called AutoBrachy QA, for comprehensive and efficient QA. This abstract reports a new plugin in this system for the QA of a cylinder HDR brachytherapy treatment. Methods: A cylinder plan QA module was developed using Python. It was plugged into our AutoBrachy QA system. This module extracted information from CT images and treatment plan. Image processing techniques were employed to obtain geometric parameters, e.g. cylinder diameter. A comprehensive set of eight geometrical and eight dosimetric features of the plan were validated against user specified planning parameter, such as prescription value, treatment depth and length, etc. A PDF document was generated, consisting of a summary QA sheet with all the QA results, as well as images showing plan details. Results: The cylinder QA program has been implemented in our clinic. To date, it has been used in 11 patient cases and was able to successfully perform QA tests in all of them. The QA program reduced the average plan QA time from 7 min using conventional manual approach to 0.5 min. Conclusion: Being a new module in our AutoBrachy QA system, an automated treatment plan QA module for cylinder HDR brachytherapy has been successfully developed and clinically implemented. This module improved clinical workflow and plan integrity compared to the conventional manual approach.

60 APPLIED LIFE SCIENCES↗

MO-DE-BRA-01: Flipped Physics Courses Within a Radiologic Technologist Program: Video Production and Long Term Outcomes

Purpose: To determine if the flipped class format has an effect on examination results for a radiologic technologist (RT) program and discuss benefits from creating video resources. Methods: From 2001–2015, students had taken both a radiological physics and quality control (QC) class as a part of their didactic training. In 2005/2006, the creation of videos of didactic lectures and QC test demonstrations allowed for a flip where content was studied at home while exercises and reviews were done in-class. Final examinations were retrospectively reviewed from this timeframe. 12 multiple choice physics questions (MCP) and 5 short answer QC questions (SAQC) were common to pre and post flip exams. The RT program’s ARRT exam scores were also obtained and compared to national averages. Results: In total, 36 lecture videos and 65 quality control videos were created for the flipped content. Data was ∼2.4GB and distributed to students via USB or CD media. For MCP questions, scores improved by 7.9% with the flipped format and significance (Student’s t-test, p<0.05) was found for 3 of the 12 questions. SAQC questions showed improvement by 14.6% and significance was found for 2 of the 5 questions. Student enrollment increased from ∼14 (2001–2004) to ∼23 students (2005–15). Content was continuously added post-flip due to the efficiency of delivery. The QC class in 2003 covered 45 test setups in-class while 65 were covered with video segments in 2014. Flipped materials are currently being repurposed. In 2015, this video content was restructured into an ARRT exam review guide and in 2016, the content was reorganized for fluoroscopy training for physicians. Conclusion: We believe that flipped classes can improve efficiency of content delivery and improve student performance even with an increase in class size. This format allows for flexibility in learning as well as re-use in multiple applications.

07 ISOTOPE AND RADIATION SOURCES↗

Sci-Fri PM: Radiation Therapy, Planning, Imaging, and Special Techniques - 05: A novel respiratory motion simulation program for VMAT treatment plans: a phantom validation study

Purpose: To develop and validate a computational method to simulate craniocaudal respiratory motion in a VMAT treatment plan. Methods: Three 4DCTs of the QUASAR respiratory motion phantom were acquired with a 2cm water-density spherical tumour embedded in cedar to simulate lung. The phantom was oscillating sinusoidally with an amplitude of 2cm and periods of 3, 4, and 5 seconds. An ITV was contoured and 5mm PTV margin was added. High and a low modulation factor VMAT plans were created for each scan. An in-house program was developed to simulate respiratory motion in the treatment plans by shifting the MLC leaf positions relative to the phantom. Each plan was delivered to the phantom and the dose was measured using Gafchromic film. The measured and calculated plans were compared using an absolute dose gamma analysis (3%/3mm). Results: The average gamma pass rate for the low modulation plan and high modulation plans were 91.1% and 51.4% respectively. The difference between the high and low modulation plans gamma pass rates is likely related to the different sampling frequency of the respiratory curve and the higher MLC leaf speeds in the high modulation plan. A high modulation plan has a slower gantry speed and therefore samples the breathing cycle at a coarser frequency leading to inaccuracies between the measured and planned doses. Conclusion: A simple program, including a novel method for increasing sampling frequency beyond the control point frequency, has been developed to simulate respiratory motion in VMAT plans by shifting the MLC leaf positions.

60 APPLIED LIFE SCIENCES↗

Programming van der Waals interactions with complex symmetries into microparticles using liquid crystallinity

Asymmetric interactions such as entropic (e.g., encoded by nonspherical shapes) or surface forces (e.g., encoded by patterned surface chemistry or DNA hybridization) provide access to functional states of colloidal matter, but versatile approaches for engineering asymmetric van der Waals interactions have the potential to expand further the palette of materials that can be assembled through such bottom-up processes. We show that polymerization of liquid crystal (LC) emulsions leads to compositionally homogeneous and spherical microparticles that encode van der Waals interactions with complex symmetries (e.g., quadrupolar and dipolar) that reflect the internal organization of the LC. Experiments performed using kinetically controlled probe colloid adsorption and complementary calculations support our conclusion that LC ordering can program van der Waals interactions by ~20 k B T across the surfaces of microparticles. Because diverse LC configurations can be engineered by confinement, these results provide fresh ideas for programming van der Waals interactions for assembly of soft matter.

36 MATERIALS SCIENCE↗

ANGEL program: Solution of large systems of linear differential equations describing nonstationary processes using CUDA technology

This paper presents the formulation of the problem and the methodical approach for solving large systems of linear differential equations describing nonstationary processes with the use of CUDA technology; this approach is implemented in the ANGEL program. Results for a test problem on transport of radioactive products over loops of a nuclear power plant are given. The possibilities for the use of the ANGEL program for solving various problems that simulate arbitrary nonstationary processes are discussed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The BURNUP package of applied programs used for computing the isotopic composition of materials of an operating nuclear reactor

This paper described the procedure of implementation and the possibilities of the BURNUP program. The purpose of the program is to predict the change in the nuclear composition of the materials of which a reactor is made in the course of its run and compute the radiation characteristics of the materials after their irradiation.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Program for studying fundamental interactions at the PIK reactor facilities

A research program aimed at studying fundamental interactions by means of ultracold and polarized cold neutrons at the GEK-4-4′ channel of the PIK reactor is presented. The apparatus to be used includes a source of cold neutrons in the heavy-water reflector of the reactor, a source of ultracold neutrons based on superfluid helium and installed in a cold-neutron beam extracted from the GEK-4 channel, and a number of experimental facilities in neutron beams. An experiment devoted to searches for the neutron electric dipole moment and an experiment aimed at a measurement the neutron lifetime with the aid of a large gravitational trap are planned to be performed in a beam of ultracold neutrons. An experiment devoted to measuring neutron-decay asymmetries with the aid of a superconducting solenoid is planned in a beam of cold polarized neutrons from the GEK-4′ channel. The second ultracold-neutron source and an experiment aimed at measuring the neutron lifetime with the aid of a magnetic trap are planned in the neutron-guide system of the GEK-3 channel. In the realms of neutrino physics, an experiment intended for sterile-neutrino searches is designed. The state of affairs around the preparation of the experimental equipment for this program is discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Sequential Quadratic Programming Algorithm for Nonsmooth Problems with Upper- \({\boldsymbol{\mathcal{C}^2}}\) Objective

An optimization algorithm for nonsmooth nonconvex constrained optimization problems with upper- \({\boldsymbol{\mathcal{C}^2}}\) objective functions is proposed and analyzed. Upper- \({\boldsymbol{\mathcal{C}^2}}\) is a weakly concave property that exists in difference of convex (DC) functions and arises naturally in many applications, particularly certain classes of solutions to parametric optimization problems e.g., recourse of stochastic programming and projection onto closed sets. The algorithm can be viewed as an extension of sequential quadratic programming (SQP) to nonsmooth problems with upper- \({\boldsymbol{\mathcal{C}^2}}\) objectives or a simplified bundle method. It is globally convergent with bounded algorithm parameters that are updated with a trust-region criterion. The algorithm handles general smooth constraints through linearization and uses a line search to ensure progress. The potential inconsistencies from the linearization of the constraints are addressed through a penalty method. In conclusion, the capabilities of the algorithm are demonstrated by solving both simple upper- \({\boldsymbol{\mathcal{C}^2}}\) problems and a real-world optimal power flow problem used in current power grid industry practices.

97 MATHEMATICS AND COMPUTING↗

Nucleon form factors program with SBS at JLAB

The physics of the nucleon form factors is the basic part of the Jefferson Laboratory program. We review the achievements of the 6-GeV era and the program with the 12- GeV beam with the SBS spectrometer in Hall A, with a focus on the nucleon ground state properties.

Wojtsekhowski, Bogdan B.↗

Fault Tolerance for OpenSHMEM. In: PGAS '14 Proceedings of the 8th International Conference on Partitioned Global Address Space Programming Models. Article No. 23

On today's supercomputing systems, faults are becoming a norm rather than an exception. Given the complexity required for achieving expected scalability and performance on future systems, this situation is expected to become worse. The systems are expected to function in a nearly constant presence of faults. To be productive on these systems, programming models will require both hardware and software to be resilient to faults. With the growing importance of PGAS programming model and OpenSHMEM, as a part of HPC software stack, a lack of a fault tolerance model may become a liability for its users. Towards this end, in this paper, we discuss the viability of using checkpoint/restart as a fault-tolerance method for OpenSHMEM, propose a selective checkpoint/restart fault-tolerance model, and discuss challenges associated with implementing the proposed model.

Hao, Pengfei↗

Efficient Execution of Recursive Programs on Commodity Vector Hardware

The pursuit of computational efficiency has led to the proliferation of throughput-oriented hardware, from GPUs to increasingly-wide vector units on commodity processors and accelerators. This hardware is designed to efficiently execute data-parallel computations in a vectorized manner. However, many algorithms are more naturally expressed as divide-and-conquer, recursive, task-parallel computations; in the absence of data parallelism, it seems that such algorithms are not well-suited to throughput-oriented architectures. This paper presents a set of novel code transformations that expose the data-parallelism latent in recursive, task-parallel programs. These transformations facilitate straightforward vectorization of task-parallel programs on commodity hardware. We also present scheduling policies that maintain high utilization of vector resources while limiting space usage. Across several task-parallel benchmarks, we demonstrate both efficient vector resource utilization and substantial speedup on chips using Intel's SSE4.2 vector units as well as accelerators using Intel's AVX512 units.

resursion↗

Exploiting Vector and Multicore Parallelsim for Recursive, Data- and Task-Parallel Programs

Modern hardware contains parallel execution resources that are well-suited for data-parallelism-vector units-and task parallelism-multicores. However, most work on parallel scheduling focuses on one type of hardware or the other. In this work, we present a scheduling framework that allows for a unified treatment of task- and data-parallelism. Our key insight is an abstraction, task blocks, that uniformly handles data-parallel iterations and task-parallel tasks, allowing them to be scheduled on vector units or executed independently as multicores. Our framework allows us to define schedulers that can dynamically select between executing task- blocks on vector units or multicores. We show that these schedulers are asymptotically optimal, and deliver the maximum amount of parallelism available in computation trees. To evaluate our schedulers, we develop program transformations that can convert mixed data- and task-parallel pro- grams into task block-based programs. Using a prototype instantiation of our scheduling framework, we show that, on an 8-core system, we can simultaneously exploit vector and multicore parallelism to achieve 14×-108× speedup over sequential baselines.

Ren, Bin↗

Python in the NERSC Exascale Science Applications Program for Data

We describe a new effort at the National Energy Re- search Scientific Computing Center (NERSC) in performance analysis and optimization of scientific Python applications targeting the Intel Xeon Phi (Knights Landing, KNL) many- core architecture. The Python-centered work outlined here is part of a larger effort called the NERSC Exascale Science Applications Program (NESAP) for Data. NESAP for Data focuses on applications that process and analyze high-volume, high-velocity data sets from experimental/observational science (EOS) facilities supported by the US Department of Energy Office of Science. We present three case study applications from NESAP for Data that use Python. These codes vary in terms of “Python purity” from applications developed in pure Python to ones that use Python mainly as a convenience layer for scientists without expertise in lower level programming lan- guages like C, C++ or Fortran. The science case, requirements, constraints, algorithms, and initial performance optimizations for each code are discussed. Our goal with this paper is to contribute to the larger conversation around the role of Python in high-performance computing today and tomorrow, highlighting areas for future work and emerging best practices

Ronaghi, Zahra↗

Design Considerations for GPU-based Mixed Integer Programming on Parallel Computing Platforms

Mixed Integer Programming (MIP) is a powerful abstraction in combinatorial optimization that finds real-life application across many significant sectors. The recent proliferation of graphical processing unit (GPU)-based accelerated computing architectures in large-scale parallel computing or supercomputing presents new opportunities as well as challenges in the advancement of MIP solver technology to effectively use the new accelerated computing platforms and scale to large parallel systems. Here, we recount the conventional processor-based strategies and focus on configurations where the most promising intersection lies between parallel MIP solver approaches and the specific strengths of accelerated parallel platforms. We note that the best potential lies in solving problems whose individual matrix sizes (of the linear program relaxation) fit entirely within one accelerator's memory and whose branch-and-bound (or branch-and-cut) trees cannot be fully contained within a small number of computational nodes. Additionally, we identify ideal features of computational linear algebra support on GPU accelerators that would help advance this direction of scalable parallel solution of MIP problems on GPU-based accelerated computing architectures.

Perumalla, Kalyan↗

Performance Evaluation of Heterogeneous GPU Programming Frameworks for Hemodynamic Simulations

Preparing for the deployment of large scientific and engineering codes on upcoming exascale systems with GPU-dense nodes is made challenging by the unprecedented diversity of device architectures and heterogeneous programming models. In this work, we evaluate the process of porting a massively parallel, fluid dynamics code written in CUDA to SYCL, HIP, and Kokkos with a range of backends, using a combination of automated tools and manual tuning. We use a proxy application along with a custom performance model to inform the results and identify additional optimization strategies. At scale performance of the programming model implementations are evaluated on pre-production GPU node architectures for Frontier and Aurora, as well as on current NVIDIA device-based systems Summit and Polaris. Real-world workloads representing 3D blood flow calculations in complex vasculature are assessed. Our analysis highlights critical trade-offs between code performance, portability, and development time.

Martin, Aristotle↗

Jefferson Lab at 12 GeV: The Science Program

Jefferson Lab’s upgrade of its Continuous Electron Beam Accelerator Facility (CEBAF) has recently been completed. The project involved an upgrade of the accelerator to achieve a maximum beam energy of 12 GeV and the construction of a fourth end station, Experimental Hall D, as well as new detector equipment for two of the three existing halls (A, B, and C). A broad experimental program has been developed to map the nucleon’s intrinsic quark distributions in transverse space and in longitudinal momentum through measurements of deeply exclusive and semi-inclusive processes, and to probe color confinement by studying the spectrum of hadrons with active gluon degrees of freedom in the wave function. Other programs include the forward parton distribution function at large quark momentum fraction χ, the quark and gluon polarized distribution functions, measurements of electromagnetic form factors of the nucleon ground state and of nucleon resonance transitions at short distances, and the exploration of physics beyond the Standard Model in high-precision parity-violating processes and in the search for signals of dark matter. The higher beam energy is also suitable for exploration of quark hadronization properties using the nucleus as a laboratory. This review highlights major areas of hadron and nuclear science that will be the focus of the first 5 years of operation.

3D imaging↗

Characterizing and Mitigating Work Time Inflation in Task Parallel Programs

Task parallelism raises the level of abstraction in shared memory parallel programming to simplify the development of complex applications. However, task parallel applications can exhibit poor performance due to thread idleness, scheduling overheads, and work time inflation – additional time spent by threads in a multithreaded computation beyond the time required to perform the same work in a sequential computation. We identify the contributions of each factor to lost efficiency in various task parallel OpenMP applications and diagnose the causes of work time inflation in those applications. Increased data access latency can cause significant work time inflation in NUMA systems. Our locality framework for task parallel OpenMP programs mitigates this cause of work time inflation. Our extensions to the Qthreads library demonstrate that locality-aware scheduling can improve performance up to 3X compared to the Intel OpenMP task scheduler.

Olivier, Stephen L.↗