Improving resilience of scientific software through a domain-specific approach
Not Available
Reguly, I. Z. (ORCID:0000000243854204)↗
SEARCH · Search NASA
Results for “Parallel applications”
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.
Accelerating science: The usage of commercial clouds in ATLAS Distributed Computing
The ATLAS experiment at CERN is one of the largest scientific machines built to date and will have ever growing computing needs as the Large Hadron Collider collects an increasingly larger volume of data over the next 20 years. ATLAS is conducting R&D projects on Amazon Web Services and Google Cloud as complementary resources for distributed computing, focusing on some of the key features of commercial clouds: lightweight operation, elasticity and availability of multiple chip architectures. The proof of concept phases have concluded with the cloud-native, vendoragnostic integration with the experiment’s data and workload management frameworks. Google Cloud has been used to evaluate elastic batch computing, ramping up ephemeral clusters of up to O(100k) cores to process tasks requiring quick turnaround. Amazon Web Services has been exploited for the successful physics validation of the Athena simulation software on ARM processors. We have also set up an interactive facility for physics analysis allowing endusers to spin up private, on-demand clusters for parallel computing with up to 4 000 cores, or run GPU enabled notebooks and jobs for machine learning applications. The success of the proof of concept phases has led to the extension of the Google Cloud project, where ATLAS will study the total cost of ownership of a production cloud site during 15 months with 10k cores on average, fully integrated with distributed grid computing resources and continue the R&D projects.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Rosette-scan video-rate atomic force microscopy: Trajectory patterning and control design
We introduce an analysis and a systematic design methodology for a novel nonraster scan method based on a rosette pattern and demonstrate its application in video-rate atomic force microscopy. This pattern is traced when the lateral axes of a parallel kinematic scanner are commanded to follow a combination of two sinusoids with identical amplitudes and different frequencies. We create an internal-model-based controller to enhance the tracking performance of this pattern and implement the scheme on a microelectromechanical system scanner. The conclusions reveal high-precision tracking of the rosette pattern in order to acquire time-lapsed atomic force microscope images at the rate of 10 frames/s.
42 ENGINEERING↗
A Unit Testing Framework for Scientific Legacy Code
Large-scale scientific applications play important roles in supporting research. However, it is often very expensive and time-consuming to make changes to, maintain and evolve the scientific code due to its complexity and poor programming skills of researchers. Therefore, in order to visualize scientific code architecture to optimize software design, understand undocumented source code, and analyze software flow and functionality, we first introduce a unit testing framework (UTF). Then, because such infrastructure's performance is very crucial in practical use since the scientific legacy applications simulate instances in a long period of time, we improve the UTF by applying Message Passing based Parallelization and parallel I/O operations. Furthermore, due to the scientific code has enormous state data and the I/O capacity on the server is limited, we apply in situ data analysis method to encounter fewer resource limitations, and adopt signal processing to greatly reduce data transfer. Last, we demonstrated the correctness and high-efficiency of our framework for legacy Earth model on Titan supercomputer.
97 MATHEMATICS AND COMPUTING↗
Inverse Normalization Method for Cross-Sectional Imaging and Velocimetry of Two-Phase Flows Based on Electrical Capacitance Tomography
Proper normalization of the measured capacitance data is a prerequisite in electrical capacitance tomography (ECT). Conventional methods, such as parallel and series normalization, assume a higher permittivity perturbation in a homogeneous (lower) permittivity background. Although this is applicable to a number of cases, some two-phase flows are better modeled instead as a lower permittivity perturbation in a homogeneous (higher) permittivity background. A different normalization method is proposed for such cases, which provides similar image quality as conventional, but is particularly beneficial for velocimetry applications based on sensitivity gradient of ECT sensor. Furthermore, simulation and experimental results are provided to illustrate the advantages for the proposed normalization technique.
47 OTHER INSTRUMENTATION↗
New Horizons for High-Performance Computing
Here we provide an overview of the past, present, and a diverse collection of future computer architecture alternatives for HPC. The end of Moore’s Law influenced the current HPC architecture focus on accelerated compute nodes composed of CPU and GPU computing components integrated into massively parallel processor architecture systems. There are many alternatives for future HPC directions, with different technologies, computing ecosystems, opportunities for lead user application-driven customization, and the role of open innovation business models. This paper provides an overview of these different new horizons for HPC, an organizing principle to focus future computing research, different public-private partnership models, and the critical role of workforce development.
97 MATHEMATICS AND COMPUTING↗
srMO-BO-3GP: A sequential regularized multi-objective Bayesian optimization for constrained design applications using an uncertain Pareto classifier
Bayesian optimization (BO) is an efficient and flexible global optimization framework that is applicable to a very wide range of engineering applications. To leverage the capability of the classical BO, many extensions, including multi-objective, multi-fidelity, parallelization, and latent-variable modeling, have been proposed to address the limitations of the classical BO framework. In this work, we propose a novel multi-objective BO formalism, called srMO-BO-3GP, to solve multi-objective optimization problems in a sequential setting. Three different Gaussian processes (GPs) are stacked together, where each of the GPs is assigned with a different task. The first GP is used to approximate a single-objective computed from the multi-objective definition, the second GP is used to learn the unknown constraints, and the third one is used to learn the uncertain Pareto frontier. At each iteration, a multi-objective augmented Tchebycheff function is adopted to convert multi-objective to single-objective, where the regularization with a regularized ridge term is also introduced to smooth the single-objective function. Finally, we couple the third GP along with the classical BO framework to explore the convergence and diversity of the Pareto frontier by the acquisition function for exploitation and exploration. The proposed framework is demonstrated using several numerical benchmark functions, as well as a thermomechanical finite element model for flip-chip package design optimization.
42 ENGINEERING↗
Communication Lower Bounds and Optimal Algorithms for Symmetric Matrix Computations
In this article, we focus on the communication costs of three symmetric matrix computations: (i) multiplying a matrix with its transpose, known as a symmetric rank-k update (SYRK) (ii) adding the result of the multiplication of a matrix with the transpose of another matrix and the transpose of that result, known as a symmetric rank-2k update (SYR2K) (iii) performing matrix multiplication with a symmetric input matrix (SYMM). All three computations appear in the Level 3 Basic Linear Algebra Subroutines (BLAS) and have wide use in applications involving symmetric matrices. We establish communication lower bounds for these kernels using sequential and distributed-memory parallel computational models, and we show that our bounds are tight by presenting communication-optimal algorithms for each setting. Our lower bound proofs rely on applying a geometric inequality for symmetric computations and analytically solving constrained nonlinear optimization problems. As a result, the symmetric matrix and its corresponding computations are accessed and performed according to a triangular block partitioning scheme in the optimal algorithms.
Al Daas, Hussam [Rutherford Appleton Laboratory, D↗
Parameter Optimization Toolbox for NS-3 network optimization, NS-3 Parameter Optimization Framework [SWR-18-60]
This simulation-based parameter optimization framework is proposed to tune parameters of different types of communication networks using ns-3 to achieve the optimal network performance. It consists of three main components: an ns-3 packet reporting module; a sampler running simulations with all possible parameter sets for the input parameter variables by using a parallel executor at each generation; and a hybrid optimization algorithm for tuning configurable parameters of hybrid designs and application parameter variables. The proposed hybrid metaheuristic optimization algorithm combines an evolutionary algorithm with a gradient descent function to quickly achieve an approximate globally optimum solution. This software is designed to be used in a multi-core processing Linux environment and run over a long duration of time. The execution time varies depending mainly upon the nature of the ns-3 configuration being simulated. This software includes a custom ns-3 QoS measurement application which must be included with the ns-3 source code during installation of the software.
Hasandka, Adarsh↗
The Need for Nb3Sn Coated Cu Cavities for Future Accelerators
Based on current efforts in the U.S. on the novel concept of parallel-feed RF accelerator structures, and in the U.S. and abroad in producing Nb3Sn films on either Cu or bronze, we rec-ommend that the Particle Physics community foster R&D in Superconducting Nb3Sn coated Cu RF Cavities instead of costly bulk Nb. The paper includes methods to process the coated cavi-ties at temperatures consistent with Cu retaining its shape. A devoted global effort in develop-ing Cu cavity structures coated with Nb3Sn would make the ILC or Higgs factories more afforda-ble and more likely to be built. Not only do parallel-feed RF structures enable both higher ac-celerating gradients and higher efficiencies, but they would be applicable to both Cu and Nb3Sn coated Cu cells. Increased effort on these two techniques would synergize expenditures to-wards progress, which will converge on the choice of technology for the RF of an ILC or any fu-ture accelerator. The current methods of Nb3Sn coatings on Cu or bronze can be geared also towards standard cavity cells. In conclusion, the use of distributed coupling structure topology within improved performance parameters together with Nb3Sn coating technology can lead to a paradigm shift for superconducting linacs, with higher gradient, higher temperature of opera-tion, and reduced overall costs for any future collider.
43 PARTICLE ACCELERATORS↗
Computational and experimental study of laminar flames
During the past three years, our research has centered on an investigation of the effects of complex chemistry and detailed transport on the structure and extinction of hydrocarbon flames in coflowing axisymmetric configurations. We have pursued both computational and experimental aspects of the research in parallel on both steady-state and time-dependent systems. The computational work has focused on the application of accurate and efficient numerical methods for the solution of the steady-state and time-dependent boundary value problems describing the various reacting systems. Detailed experimental measurements were performed on axisymmetric coflow flames using two-dimensional imaging techniques. Previously, spontaneous Raman scattering, chemiluminescence, and laser-induced fluorescence were used to measure the temperature, major and minor species profiles. Particle image velocimetry (PIV) has been used to investigate velocity distributions and for calibration of time-varying flames. Laser-induced incandescence (LII) with an extinction calibration was used to determine soot volume fractions, while soot surface temperatures were measured with three-color optical pyrometry using a color digital camera. A blackbody calibration of the camera allows for determination of soot volume fraction as well, which can be compared with the LII measurements. More recently, we have concentrated on a detailed characterization of soot using a variety of techniques including time-resolved LII (TiRe-LII) for soot primary particles sizes, multi-angle light scattering (MALS) for soot radius of gyration, and spectrally-resolved line of sight attenuation (spec-LOSA). Combining the information from all of these soot measurements can be used to determine the soot optical properties, which are observed to vary significantly depending on spatial location and fuel dilution. Our goal has been to obtain a more fundamental understanding of the important fluid dynamic and chemical interactions in these flames so that this information can be used effectively in combustion modeling.
02 PETROLEUM↗
Methods and Tools to allow molecular flow simulations to be coupled to higher level continuum descriptions of flows in porous/fractured media and aerosol/dust dynamics
The purpose of this project was to develop methods and tools that will aid in safety evaluation of nuclear fuels and licensing of nuclear reactors relating to accidents.The objectives were to develop more detailed and faster computations of fission product transport and aerosol evolution as they generally relate to nuclear fuel and/or nuclear reactor accidents. The two tasks in the project related to molecular transport in nuclear fuel and aerosol transport in reactor vessel and containment. For both the tasks, explorations of coupling of Direct Simulation Monte Carlo with Navier-Stokes solvers or the Sectional method were not successful. However, Mesh free methods for the Direct Simulation Monte Carlo method were successfully explored.These explorations permit applications to porous and fractured media, and arbitrary geometries.The computations were carried out in Mathematica and are fully parallelized. The project has resulted in new computational tools (algorithms and programs) that will improve the fidelity of computations to actual physics, chemistry and transport of fission products in the nuclear fuel and aerosol in reactor primary and secondary containments.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
FleCSPH notes
FleCSI is a compile-time con gurable framework designed to support multi-physics application development. As such, FleCSI provides a very general set of infrastructure design patterns that can be specialized and extended to suit the needs of a broad variety of solver and data requirements. FleCSI currently supports multi-dimensional mesh topology, geometry, and adjacency information, as well as n-dimensional hashed-tree data structures, graph partitioning interfaces, and dependency closures. FleCSI introduces a functional programming model with control, execution, and data abstractions that are consistent both with MPI and with state-of-the-art, task-based runtimes such as Legion and Charm++. The abstraction layer insulates developers from the underlying runtime, while allowing support for multiple runtime systems including conventional models like asynchronous MPI. The intent is to provide developers with a concrete set of user-friendly programming tools that can be used now, while allowing exibility in choosing runtime implementations and optimization that can be applied to future architectures and runtimes. FleCSI's control and execution models provide formal nomenclature for describing poorly understood concepts such as kernels and tasks. FleCSI's data model provides a low-buy-in approach that makes it an attractive option for many application projects, as developers are not locked into particular layouts or data structure representations. FleCSI currently provides a parallel but not distributed implementation of Binary, Quad and Oct-tree topology. This implementation is base on space lling curves domain decomposition, the Morton order. The current FleCSI version requires the implementation of a driver and a specialization driver. The role of the specialization driver is to provide the data distribution. This feature is not complete in FleCSI code and we provide it. The next step will be to incorporate it directly from FleCSPH to FleCSI as we reach a good level of performance. Then the driver represent the general execution of the resolution without worrying of the data locality and communications. As FleCSI is an On-Development code the structure may change in the future and we keep track of these changes in FleCSPH.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
SPARC: Demonstrate burst-buffer-based checkpoint/restart on ATS-1.
Recent high-performance computing (HPC) platforms such as the Trinity Advanced Technology System (ATS-1) feature burst buffer resources that can have a dramatic impact on an application’s I/O performance. While these non-volatile memory (NVM) resources provide a new tier in the storage hierarchy, developers must find the right way to incorporate the technology into their applications in order to reap the benefits. Similar to other laboratories, Sandia is actively investigating ways in which these resources can be incorporated into our existing libraries and workflows without burdening our application developers with excessive, platform-specific details. This FY18Q1 milestone summaries our progress in adapting the Sandia Parallel Aerodynamics and Reentry Code (SPARC) in Sandia’s ATDM program to leverage Trinity’s burst buffers for checkpoint/restart operations. We investigated four different approaches with varying tradeoffs in this work: (1) simply updating job script to use stage-in/stage out burst buffer directives, (2) modifying SPARC to use LANL’s hierarchical I/O (HIO) library to store/retrieve checkpoints, (3) updating Sandia’s IOSS library to incorporate the burst buffer in all meshing I/O operations, and (4) modifying SPARC to use our Kelpie distributed memory library to store/retrieve checkpoints. Team members were successful in generating initial implementation for all four approaches, but were unable to obtain performance numbers in time for this report (reasons: initial problem sizes were not large enough to stress I/O, and SPARC refactor will require changes to our code). When we presented our work to the SPARC team, they expressed the most interest in the second and third approaches. The HIO work was favored because it is lightweight, unobtrusive, and should be portable to ATS-2. The IOSS work is seen as a long-term solution, and is favored because all I/O work (including checkpoints) can be deferred to a single library.
97 MATHEMATICS AND COMPUTING↗
A Call to Arms Control: Synergies between Nonproliferation Applications of Neutrino Detectors and Large-Scale Fundamental Neutrino Physics Experiments (A Snowmass White Paper)
The High Energy Physics community can benefit from a natural synergy in research activities into next-generation large-scale water and scintillator neutrino detectors, now being studied for remote reactor monitoring, discovery and exclusion applications in cooperative nonproliferation contexts. Since approximately 2010, US nonproliferation researchers, supported by the National Nuclear Security Administration (NNSA), have been studying a range of possible applications of relatively large (100 ton) to very large (hundreds of kiloton) water and scintillator neutrino detectors. In parallel, the fundamental physics community has been developing detectors at similar scales and with similar design features for a range of high-priority physics topics, primarily in fundamental neutrino physics. These topics include neutrino oscillation studies at beams and reactors, solar, and geological neutrino measurements, supernova studies, and others. Examples of ongoing synergistic work at U.S. national laboratories and universities include prototype gadolinium-doped water and water-based and opaque scintillator test-beds and demonstrators, extensive testing and industry partnerships related to large area fast position-sensitive photomultiplier tubes, and the development of concepts for a possible underground kiloton-scale water-based detector for reactor monitoring and technology demonstrations. Some opportunities for engagement between the two communities include bi-annual Applied Antineutrino Physics conferences, collaboration with U.S. National Laboratories engaging in this research, and occasional NNSA funding opportunities supporting a blend of nonproliferation and basic science R&D, directed at the U.S. academic community.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Effects of Hybridization on Selective Catalytic Reduction (SCR) Thermal Management of a Medium Heavy-Duty Hybrid Work Truck
The increased market penetration of hybrid electric powertrains in medium heavy-duty (MHD) applications has provided a novel platform for vehicle research. One example of such a platform is the MHD parallel hybrid truck developed by Odyne Systems, LLC. In collaboration with Odyne Systems, LLC and the Department of Energy (DOE), Oak Ridge National Laboratory (ORNL) developed a validated vehicle plant model for this truck and tested the Odyne powertrain in a hardware-in-the-loop (HIL) environment. While testing in the HIL environment, the effects of reduced engine load, and thus catalyst heating, on the selective catalytic reduction (SCR) catalyst produced diminished hybrid improvement as the level of energy storage usage increased. This article will discuss these results and the potentially unforeseen interactions with modern aftertreatment systems when hybridizing conventional powertrains.
33 ADVANCED PROPULSION SYSTEMS↗
Study of solid rocket motors for a space shuttle booster. Volume 4: Mass properties report
Mass properties data for the 156 inch diameter, parallel burn, solid propellant rocket engine for the space shuttle booster are presented. Design ground rules and assumptions applicable to generation of the mass properties data are described, together with pertinent data sources.
Vonderesch, A. H.↗