Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Software Library”

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 937 records · Page 52

Rapid Optimization Library

ROL provides interfaces to and implementations of algorithms for gradient-based unconstrained and constrained optimization. ROL can be used to optimize the response of any client simulation code that evaluates scalar-valued response functions. If the client code can provide gradient information for the response function, ROL will take advantage of it, resulting in faster runtimes. ROL's interfaces are matrix-free, in other words ROL only uses evaluations of scalar-valued and vector-valued functions. ROL can be used to solve optimal design problems and inverse problems based on a variety of simulation software.

Kouri, DrewKouri↗

Rapid Optimization Library

ROL provides interfaces to and implementations of algorithms for gradient-based unconstrained and constrained optimization. ROL can be used to optimize the response of any client simulation code that evaluates scalar-valued response functions. If the client code can provide gradient information for the response function, ROL will take advantage of it, resulting in faster runtimes. ROL's interfaces are matrix-free, in other words ROL only uses evaluations of scalar-valued and vector-valued functions. ROL can be used to solve optimal design problems and inverse problems based on a variety of simulation software.

Denis Rldzal, Drew Kouri↗

The deal.II library, Version 9.3

This paper provides an overview of the new features of the finite element library deal.II, version 9.3.

97 MATHEMATICS AND COMPUTING↗

Memory forensic analysis of a programmable logic controller in industrial control systems

In industrial control systems (ICS), programmable logic controllers (PLCs) are used to automate physical processes such as nuclear plants and power grid stations, and are often subject to cyber attacks. As in conventional IT domain, the memory analysis of the PLCs can help answer important forensic questions about the attack, such as the presence of malicious firmware, injection of modified control logic (the program running on the PLC), and manipulation of I/O devices (e.g., sensors and actuators). Unlike conventional IT domain, PLCs have heterogeneous hardware architecture, proprietary firmware and control software, making it challenging to employ a unified framework for their memory forensics. For merely extracting artifacts of forensic importance, reverse-engineering the firmware is a tedious task, and the effort needs to be repeated for every PLC model. As a community, a step-wise approach to tackle this challenge is to analyze the memory of specific PLCs, and subsequently find a generic framework applicable to all PLCs. Our work is a step forward in this direction. By following a methodology that focuses on the functional layer of PLCs instead of reverse engineering the firmware, we analyze the digital forensic artifacts available in a common PLC, Allen-Bradley ControlLogix 1756-L61. Before diving into the memory dump, we analyze the PLC control software to create a list of important artifacts that are sure to exist in the PLC memory dump. The approach employs a setup where PLC control software RSLogix-5000 is connected to the PLC, and the memory dump can be obtained as and when needed. We create test cases that sequentially highlight each category of artifacts, followed by an examination of the resultant impact on memory. After attaining the listed artifacts, we employ conventional string and known data searches to extract interesting information present in this PLC's memory. The memory analysis profile, presented as a Python library and shared with the community, can help a forensic investigator to readily extract forensic artifacts from the same model's controller. The adopted approach may help researchers in creating memory profile of other PLCs, and ultimately formulating a generic PLC memory analysis framework.

Rais, Muhammad Haris↗

Computational Modeling of a 3D Printed Recuperator and Subsequent Experimental Loop for Supercritical Carbon Dioxide Cycles

Oak Ridge National Laboratory (ORNL), in collaboration with mechanical-thermal energy storage (mTES) provider EarthEn, a US Department of Energy (DOE) Lab-Embedded Entrepreneurship Program (LEEP) recipient at ORNL’s Innovation Crossroads 2023, is utilizing a state-of-the-art patented 3D printing technique to design an additively manufactured (AM) supercritical CO2 (sCO2) recuperator (REC) for EarthEn’s charge/discharge cycle. The AM REC will be printed at ORNL’s Manufacturing Demonstration Facility using Inconel Alloy 718 and tested on a closed-loop, ∼100 kW scale experimental facility that is under construction. The testing will compare the printed design against a commercial-off-the-shelf Printed Circuit Heat Exchanger (PCHE) REC. The design of the sCO2 facility is guided by a Modelica-based system model which is primarily dependent on the open-source TRANSFORM library developed at ORNL and uses the open-source CoolProp library for thermophysical properties of sCO2 via the External Media library. It is envisioned that an iterative process will be followed between the physical loop and the system model wherein the initial experimental data will be used to tune the model, which in turn will be used to guide future loop operation. Simultaneously, the AM REC is being designed using computer-aided design models, and it is also being analyzed for hydraulic and thermomechanical response using commercial computational fluid dynamics software, Simcenter STAR-CCM+, on highperformance computing resources.1

See, Nate [ORNL] (ORCID:0000000178581202)↗

Electronic journals: Their use by teachers/researchers of engineering and social sciences

Libraries must attend the needs of their different users. Academics are usually a particular kind of users with specific needs. Universities are environments where scientific communication is essential and where electronic format of journals is becoming more and more frequently used. This way it becomes increasingly important to understand how academics from different scientific areas use the available electronic resources. The aim of this study is to better understand the existing differences among the users of electronic journals in Engineering and Social Sciences. The research undertaken was mainly focused on the study of the use of electronic journals by teachers/researchers from the Faculties of Engineering and of Arts from the University of Porto, Portugal. In this study an international survey was used in order to characterize the levels of use and access of electronic journals by these communities. The ways of seeking and using scientific information, namely in terms frequency of access, the number of articles consulted, the use of databases and the preference of publishing in electronic journals were analyzed. A set of comparisons were established and results indicate an extensive use of the electronic format, regardless the faculty. However, some differences emerge when it comes to details. Such is the case of the usage rate of reference management software which is considerably more used by Engineering academics than Social Science ones. Generally, electronic journals meeting the information needs of its users and are increasingly used as a preferred means of research. Though, some particular differences in the use of them have emerged, when comparing academics from these two faculties.

99 GENERAL AND MISCELLANEOUS↗

PySAGES: flexible, advanced sampling methods accelerated with GPUs

Abstract Molecular simulations are an important tool for research in physics, chemistry, and biology. The capabilities of simulations can be greatly expanded by providing access to advanced sampling methods and techniques that permit calculation of the relevant underlying free energy landscapes. In this sense, software that can be seamlessly adapted to a broad range of complex systems is essential. Building on past efforts to provide open-source community-supported software for advanced sampling, we introduce PySAGES, a Python implementation of the Software Suite for Advanced General Ensemble Simulations (SSAGES) that provides full GPU support for massively parallel applications of enhanced sampling methods such as adaptive biasing forces, harmonic bias, or forward flux sampling in the context of molecular dynamics simulations. By providing an intuitive interface that facilitates the management of a system’s configuration, the inclusion of new collective variables, and the implementation of sophisticated free energy-based sampling methods, the PySAGES library serves as a general platform for the development and implementation of emerging simulation techniques. The capabilities, core features, and computational performance of this tool are demonstrated with clear and concise examples pertaining to different classes of molecular systems. We anticipate that PySAGES will provide the scientific community with a robust and easily accessible platform to accelerate simulations, improve sampling, and enable facile estimation of free energies for a wide range of materials and processes.

Chemistry↗

Graph Contractions for Calculating Correlation Functions in Lattice QCD

Computing correlation functions for many-particle systems in Lattice QCD is vital to extract nuclear physics observables like the energy spectrum of hadrons such as protons. However, this type of calculation has long been considered to be very challenging and computing-resource intensive because of the complex nature of a hadron composed of quarks with many degrees of freedom. In particular, a correlation function can be calculated through a sum of all possible pairs of quark contractions, each of which is a batched tensor contraction, dictated by Wick's theorem. Because the number of terms of this sum can be very large for any hadronic system of interest, fast evaluation of the sum faces several challenges: an extremely large number of contractions, a huge memory footprint at runtime, and the speed of tensor contractions. In this paper, we present a Lattice QCD analysis software suite, Redstar, which addresses these challenges by utilizing novel algorithmic and software engineering methods targeting modern computing platforms such as many-core CPUs and GPUs. In particular, Redstar represents every term in the sum of a correlation function by a graph, applies efficient graph algorithms to reduce the number of contractions to lower the cost of computations, and minimizes the total memory footprint. Moreover, Redstar carries out the contractions on either CPUs or GPUs utilizing an internal and highly efficient Hadron contraction library. Specifically, we illustrate some important algorithmic optimizations of Redstar, show various key design features of Hadron library, and present the speedup values due to the optimizations along with performance figures for calculating six correlations functions on four computing platforms.

Chen, Jie↗

Graph Contractions for Calculating Correlation Functions in Lattice QCD

Computing correlation functions for many-particle systems in Lattice QCD is vital to extract nuclear physics observables like the energy spectrum of hadrons such as protons. However, this type of calculation has long been considered to be very challenging and computing-resource intensive because of the complex nature of a hadron composed of quarks with many degrees of freedom. In particular, a correlation function can be calculated through a sum of all possible pairs of quark contractions, each of which is a batched tensor contraction, dictated by Wick's theorem. Because the number of terms of this sum can be very large for any hadronic system of interest, fast evaluation of the sum faces several challenges: an extremely large number of contractions, a huge memory footprint at runtime, and the speed of tensor contractions. In this paper, we present a Lattice QCD analysis software suite, Redstar, which addresses these challenges by utilizing novel algorithmic and software engineering methods targeting modern computing platforms such as many-core CPUs and GPUs. In particular, Redstar represents every term in the sum of a correlation function by a graph, applies efficient graph algorithms to reduce the number of contractions to lower the cost of computations, and minimizes the total memory footprint. Moreover, Redstar carries out the contractions on either CPUs or GPUs utilizing an internal and highly efficient Hadron contraction library. Specifically, we illustrate some important algorithmic optimizations of Redstar, show various key design features of Hadron library, and present the speedup values due to the optimizations along with performance figures for calculating six correlations functions on four computing platforms.

Chen, Jie↗

Automated Transformation of Lignin Topologies into Atomic Structures with LigninBuilder

Lignin is an abundant aromatic heteropolymer found in secondary plant cell walls and is a potential feedstock for conversion into bioderived fuels and chemicals. Lignin chemical diversity complicates traditional structural studies, and so, relatively little experimental evidence exists for how lignin structure exists in aqueous solution or how lignin polymers respond to changes in their chemical environment. Molecular modeling can address these concerns; however, prior computational structural lignin models typically did not capture lignin heterogeneity, as only a few polymers were considered. LigninBuilder creates a framework for building structural libraries for lignin from existing topological libraries, permitting significantly greater diversity of lignin structures to be sampled at atomic detail. As a demonstration of its capabilities, LigninBuilder was applied to three libraries of lignin from hardwood, softwood, and grass, and the resulting polymer structures were simulated in an aqueous environment. The lignins adopted compact globular structures, as would be expected for polymers in poor solvents. The differences between the libraries were largest when quantifying the packing of nonadjacent aromatic residues, with greater branching within the polymer resulting in poorer aromatic packing. Individual lignin polymers were also found to undergo rapid conformational changes, with the dwell time within a state growing as the square of the molecular weight. This first application of LigninBuilder demonstrates the potential for atomic-level insight into lignin interactions. LigninBuilder is distributed as a plugin to the visualization software VMD, lowering the barrier for modeling lignin structure in diverse environments.

09 BIOMASS FUELS↗

MBX: A many-body energy and force calculator for data-driven many-body simulations

Many-Body eXpansion (MBX) is a C++ library that implements many-body potential energy functions (PEFs) within the “many-body energy” (MB-nrg) formalism. MB-nrg PEFs integrate an underlying polarizable model with explicit machine-learned representations of many-body interactions to achieve chemical accuracy from the gas to the condensed phases. MBX can be employed either as a stand-alone package or as an energy/force engine that can be integrated with generic software for molecular dynamics and Monte Carlo simulations. MBX is parallelized internally using Open Multi-Processing and can utilize Message Passing Interface when available in interfaced molecular simulation software. In this study, MBX enables classical and quantum molecular simulations with MB-nrg PEFs, as well as hybrid simulations that combine conventional force fields and MB-nrg PEFs, for diverse systems ranging from small gas-phase clusters to aqueous solutions and molecular fluids to biomolecular systems and metal-organic frameworks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

TorchBraid: High-Performance Layer-Parallel Training of Deep Neural Networks with MPI and GPU Acceleration

TorchBraid is a high-performance implementation of layer-parallel training for deep neural networks (DNNs) supporting MPI-based parallelism and GPU acceleration. Layer-parallel training has been developed to overcome the serialization inherent in forward and backward propagation of DNNs that limits utilization of computational resources in the strong scaling limit. To achieve this, TorchBraid integrates the PyTorch neural network framework with the state-of-the-art XBraid time-parallel library. Furthermore, this article presents the use and performance of TorchBraid, in addition to solutions for overcoming the algorithmic challenges inherent in combining automatic differentiation with layer-parallel. Results are presented with and without GPU acceleration for the Tiny ImageNet and MNIST image classification data sets, as well as recurrent neural networks. Overall, TorchBraid enables fast training of DNNs, both in a strong and weak scaling context. In addition to the TorchBraid software, several new advances in applying layer-parallel algorithms are detailed. Integration of layer-parallel with data-parallel algorithms is presented for the first time, showing the computational advantages of the combination. Standard deep learning techniques, like batch-normalization, are developed for layer-parallel training. Finally, a new approach combining layer-parallel with spatial coarsening in order to accelerate training for 3D image classification shows roughly a 10× speedup over serial execution.

Layer-parallel↗

Scaling Irregular Applications through Data Aggregation and Software Multithreading

Bioinformatics, data analytics, semantic databases, knowledge discovery are emerging high performance application areas that exploit dynamic, linked data structures such as graphs, unbalanced trees or unstructured grids. These data structures usually are very large, requiring significantly more memory than available on single shared memory systems. Additionally, these data structures are difficult to partition on distributed memory systems. They also present poor spatial and temporal locality, thus generating unpredictable memory and network accesses. The Partitioned Global Address Space (PGAS) programming model seems suitable for these applications, because it allows using a shared memory abstraction across distributed-memory clusters. However, current PGAS languages and libraries are built to target regular remote data accesses and block transfers. Furthermore, they usually rely on the Single Program Multiple Data (SPMD) parallel control model, which is not well suited to the fine grained, dynamic and unbalanced parallelism of irregular applications. In this paper we present {\bf GMT} (Global Memory and Threading library), a custom runtime library that enables efficient execution of irregular applications on commodity clusters. GMT integrates a PGAS data substrate with simple fork/join parallelism and provides automatic load balancing on a per node basis. It implements multi-level aggregation and lightweight multithreading to maximize memory and network bandwidth with fine-grained data accesses and tolerate long data access latencies. A key innovation in the GMT runtime is its thread specialization (workers, helpers and communication threads) that realize the overall functionality. We compare our approach with other PGAS models, such as UPC running using GASNet, and hand-optimized MPI code on a set of typical large-scale irregular applications, demonstrating speedups of an order of magnitude.

Aggregation↗

A Formalization of Core Why3 in Coq

Intermediate verification languages like Why3 and Boogie have made it much easier to build program verifiers, transforming the process into a logic compilation problem rather than a proof automation one. Why3 in particular implements a rich logic for program specification with polymorphism, algebraic data types, recursive functions and predicates, and inductive predicates; it translates this logic to over a dozen solvers and proof assistants. Accordingly, it serves as a backend for many tools, including Frama-C, EasyCrypt, and GNATProve for Ada SPARK. But how can we be sure that these tools are correct? The alternate foundational approach, taken by tools like VST and CakeML, provides strong guarantees by implementing the entire toolchain in a proof assistant, but these tools are harder to build and cannot directly take advantage of SMT solver automation. As a first step toward enabling automated tools with similar foundational guarantees, we give a formal semantics in Coq for the logic fragment of Why3. We show that our semantics are useful by giving a correct-by-construction natural deduction proof system for this logic, using this proof system to verify parts of Why3's standard library, and proving sound two of Why3's transformations used to convert terms and formulas into the simpler logics supported by the backend solvers.

97 MATHEMATICS AND COMPUTING↗

The deal.II library, Version 9.6

Here, this paper provides an overview of the new features of the finite element library deal.II, version 9.6.

97 MATHEMATICS AND COMPUTING↗

SU-E-T-438: Commissioning of An In-Vivo Quality Assurance Method Using the Electronic Portal Imaging Device

Purpose: Patient specific pre-treatment quality assurance (QA) using arrays of detectors or film have been the standard approach to assure the correct treatment is delivered to the patient. This QA approach is expensive, labor intensive and does not guarantee or document that all remaining fractions were treated properly. The purpose of this abstract is to commission and evaluate the performance of a commercially available in-vivo QA software using the electronic portal imaging device (EPID) to record the daily treatments. Methods: The platform EPIgray V2.0.2 (Dosisoft), which machine model compares ratios of TMR with EPID signal to predict dose was commissioned for an Artiste (Siemens Oncology Care Systems) and a Truebeam (Varian medical systems) linear accelerator following the given instructions. The systems were then tested on three different phantoms (homogeneous stack of solid water, anthropomorphic head and pelvis) and on a library of patient cases. Simple and complex fields were delivered at different exposures and for different gantry angles. The effects of the table attenuation and the EPID sagging were evaluated. Gamma analysis of the measured dose was compared to the predicted dose for complex clinical IMRT cases. Results: Commissioning of the EPIgray system for two photon energies took 8 hours. The difference between the dose planned and the dose measured with EPIgray was better than 3% for all phantom scenarios tested. Preliminary results on patients demonstrate an accuracy of 5% is achievable in high dose regions for both 3DCRT and IMRT. Large discrepancies (>5%) were observed due to metallic structures or air cavities and in low dose areas. Flat panel sagging was visible and accounted for in the EPIgray model. Conclusion: The accuracy achieved by EPIgray is sufficient to document the safe delivery of complex IMRT treatments. Future work will evaluate EPIgray for VMAT and high dose rate deliveries. This work is supported by Dosisoft, Cachan, France.

60 APPLIED LIFE SCIENCES↗

Fully Implicit Conjugate Heat Transfer Analysis of the ARC-Class Vacuum Vessel

The coupled simulation of fusion reactor blankets including neutronics, thermal-hydraulics and thermo-mechanics is expected to speed up the design cycle of fusion reactor design concepts. In this work we demonstrate tight implicit coupling of conjugate heat transfer using the open-source Computational Fluid Dynamics software OpenFOAM for thermo-fluid mechanics and Diablo for thermo-solid mechanics. The heat transfer analysis is augmented by volumetric energy deposition from neutronic calculations using the Monte Carlo N-particle code on both solid and fluid parts of the vacuum vessel. An additional heat flux is imposed on the first wall estimated from the design power of the reactor. The tight coupling is realized through the open-source coupling library, preCICE, and tested on the vacuum vessel of the affordable, robust, compact reactor design by Commonwealth Fusion Systems. The features of the coupling and the influence of different coupling parameters such as coupling schemes, acceleration techniques and convergence criterion are discussed. The coupled simulation results are compared to a thermal-hydraulics simulation which includes only the fluid domains (the liquid immersion molten salt blanket and cooling channel) to demonstrate usefulness of a coupled simulation. Further analysis is performed to identify regions of hot spots for subsequent design improvement. This introduces the outline for integrating conjugate electromagnetics and fluid/solid mechanics (e.g., allow for deformation of the cooling channel walls) with our present approach for future analysis.

Sircar, Arpan↗

Real-Time Drilling Optimization System for Improved Overall Rate of Penetration and Reduced Cost Per Foot in Geothermal Drilling

The key to success in geothermal drilling is economic feasibility, and a major cost in the development of geothermal resources is the actual drilling of the wells. In this project, a real-time drilling optimization system for geothermal drilling was developed. The system couples three individual components while drilling. The first component is a drill stem vibration analysis model, the second is Mechanical Specific Energy (MSE) analyses, and the third is a detailed PDC Rate of Penetration (ROP) drill bit model for optimum RPM and WOB combinations. The benefit of the coupled system is that the range of WOB and RPM could be selected to avoid drill stem vibrations. Secondly, MSE is used as an efficiency measure and the detailed PDC drill bit model ensures the drill bit does not endure temperatures that exceed the temperature at which the PDC cutters experience accelerated wear. The new detailed PDC bit model is based on rock/bit interaction that physically tracks the PDC cutter wear flats as the bit drills ahead giving the capability to calculate the temperature being generated underneath the worn cutters to better advise on operational parameters to avoid accelerated cutter wear and failure and to ensure that operational parameters are applied so that overall ROP is maximized. By combining the drill stem vibrations and the detailed PDC bit cutter wear and “safe” non-accelerated cutter wear temperature and optimum ranges of operating parameters, it results in higher ROP and lower cost drilling. Single cutter PDC testing performed in different lithologies at Sandia was utilized to verify the PDC cutter forces and depth of cut for new and worn cutters. Based on single cutter PDC temperature modeling, verification using single cutter data from the testing done by National Oilwell Varco (NOV) was performed. Sandia’s Hard-Rock Drilling Facility (HRDF) was utilized to test different drill bit configurations with different cutter designs and wear status with different induced modes of vibration to obtain the critical bit RPM/WOB ranges resulting in ineffective drilling and low ROP. The collected test data were further used to verify and calibrate the full hole PDC ROP model that was developed based on single cutter interaction data. A full coupled drill stem vibration model was formulated and verified with geothermal field data from the Chocolate Mountain Aerial Gunnery Range (CMAGR). A graphical user interface (GUI) was developed using Tkinter library in the computer programming language Python, which integrates all the developed models in one system. The developed system consists mainly of the PDC ROP model, PDC bit wear model, PDC cutter temperature model, Mechanical Specific Energy (MSE) model, and drillstring vibration model integrated into one system. The developed system can be used for both, post well analysis and real-time optimization using different criteria such as ROP maximization or MSE minimization. The software uses Differential Evolution Algorithm (DEA) to find optimum values for operational parameters based on last foot drilled while avoiding the drillstring vibration and cutter temperature critical operating parameters.

15 GEOTHERMAL ENERGY↗