Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “library design”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

kynema-fmb [SWR-23-07]

Kynema-FMB (FKA: Kynema) is an open-source performance portable flexible multibody (FMB) dynamics solver designed for time-domain simulations. While originally tailored for wind turbine structural dynamics, the formulation and implementation are those of a general flexible-multidbody dynamics solver that can readily be applied to a wide range of systems. Kynema was designed with a narrow focus, namely to provide a lightweight, fast, accurate FMD solver for coupling to computational-fluid-dynamics (CFD) codes, especially the CFD codes in the Kynema suite, for fluid-structure-interaction (FSI) simulations. Kynema-FMB is equipped to model systems that can be represented as a collection of beams and rigid bodies that are connected through constraints. Degrees of freedom are defined in the inertial/global frame of reference and include displacements and rotations (formally as rotation matrices, but stored as quaternions). The underlying formulation is built on a Lie-group time integrator designed for index-3 differential-algebraic equations, which is second-order accurate in time (Bruls et al., 2012). Beam models are based on geometrically exact beam theory and are discretized as high-order spectral finite elements similar to those in BeamDyn (Wang et al., 2017). The governing equations for a FMD system like a wind turbine constitute a highly nonlinear system of constrained partial-differential equations. Kynema-FMB uses analytical Jacobians in the nonlinear-system solves in each time step. Linear systems use sparse storage and several third-party sparse-linear-system solvers are enabled. Ill conditioning of linear systems is mitigated with preconditioning described in Bottasso et al, 2008. Kynema-FMB is integrated with a simple open-source controller (ROSCO). There is an application programming interface (API) for coupling to geometry-resolved CFD (like that in Sharma et al., 2023) and actuator-force CFD (like that in Kuhn et al., 2025). In the latter, for actuator-line models, Kynema-FMB includes an internal blade-element solver that depends on user-provided lookup tables for coefficients of lift and drag, i.e., aerodynamic polars. Kynema-FMB is written in C++ and leverages Kokkos and Kokkos-Kernels (KokkosEcosystem) as its performance portability layer enabling simulations on both CPU and GPU systems. The repository is equipped with extensive automated testing at the unit and regression/system levels. The following describes the high-level development objectives conceived for Kynema: *Kynema will follow modern software development best practices, including test-driven development (TDD), version control, hierarchical automated testing, and continuous integration (CI) for a robust development environment. *The core data structures are memory efficient and enable vectorization and parallelization at multiple levels. *Data structures are data-oriented to exploit methods for accelerated computing including high utilization of chip resources (e.g., single instruction multiple data (SIMD) instruction sets) and parallelization using GP-GPUs. *The computational algorithms incorporate robust open-source libraries for mathematical operations, resource allocation, and data management. *The API design considers multiple stakeholder needs and ensure integration with existing and future ecosystems for data science, machine learning, and AI. *Kynema-FMB is written in modern C++ and leverages Kokkos as its performance-portability library with inspiration from the kynema stack.

Sprague, MichaelA.↗
DOE OSTI2024

Conceptual Design of Neutron Sensor Qualification Device

The Advanced Sensors and Instrumentation Program at Idaho National Laboratory has been formulating strategies to qualify sensors for use in nuclear environments, particularly in irradiation experiments and advanced reactors. When qualifying neutron sensors for use in high-temperature environments, the wide range of neutron flux levels and representative energy spectra presents significant challenges. This paper discusses the development of the Neutron Sensor Qualification Device (NQD), which is designed to test neutron sensors in high temperature controlled environments with known neutron spectra, addressing the spatial and spectral complexities of neutron fluxes in reactor cores. The proposed NQD will be situated in the exposure room at the Armed Forces Radiobiology Research Institute, thus affording a unique capability to expose sensors to high neutron and gamma fluxes. To achieve thermal control, the device will utilize a radiation-hardened tube furnace, accommodating multiple sensors and neutron activation dosimetry wires. Titanium, iron, and cobalt dosimeter wires are chosen from the American Society for Testing and Materials and International Reactor Dosimetry and Fusion File libraries as references for providing energy-dependent fluence measurements. The design ensures precise sensor positioning to minimize mutual shielding and flux perturbation, which are evaluated via Monte Carlo N particle Transport Code (MCNP) simulations. These simulations have informed the development of guidelines on sensor placement within the NQD. The NQD is essential to the qualification of neutron sensors for advanced reactor technologies. It enables controlled testing of a statistically significant number of sensors, thereby supporting assessments of sensor performance across various neutron flux levels and temperatures. This paper highlights the detailed planning for the NQD prototype, along with its inaugural irradiation (scheduled for fiscal year [FY] 2025). The results from this initial testing will be fundamental in evaluating the device’s performance and establishing measurement uncertainty for in-pile neutron sensor measurements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

The Kokkos OpenMPTarget Backend: Implementation and Lessons Learned

As the supercomputing landscape diversifies, solutions such as Kokkos to write vendor agnostic applications and libraries have risen in popularity. Kokkos provides a programming model designed for performance portability, which allows developers to write a single source implementation that can run efficiently on various architectures. At its heart, Kokkos maps parallel algorithms to architecture and vendor specific backends written in lower level programming models such as CUDA and HIP. Another approach to writing vendor agnostic parallel code is using OpenMP’s directives based approach, which lets developers annotate code to express parallelism. It is implemented at the compiler level and is supported by all major high performance computing vendors, as well as the primary Open Source toolchains GNU and LLVM. Since its inception, Kokkos has used OpenMP to parallelize on CPU architectures. In this paper, we explore leveraging OpenMP for a GPU backend and discuss the challenges we encountered when mapping the Kokkos APIs and semantics to OpenMP target constructs. As an exemplar workload we chose a simple conjugate gradient solver for sparse matrices. We find that performance on NVIDIA and AMD GPUs varies widely based on details of the implementation strategy and the chosen compiler. Furthermore, the performance of the OpenMP implementations decreases with increasing complexity of the investigated algorithms.

Gayatri, Rahulkumar↗
DOE OSTI2023

Using heterogeneous GPU nodes with a Cabana-based implementation of MPCD

In this study, the Kokkos based library Cabana, which has been developed in the Co-design Center for Particle Applications (CoPA), is used for the implementation of Multi-Particle Collision Dynamics (MPCD), a particle-based description of hydrodynamic interactions. Cabana allows for a function portable implementation, which has been used to study the interplay between CPU and GPU usage on a multi-node system as well as analysis of said interplay with performance analysis tools. As a result, we see most advantages in a homogeneous GPU usage, but we also discuss the extent to which heterogeneous applications might be more performant, using both CPU and GPU concurrently.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Horizontal Resolution Sensitivity of the Simple Convection-Permitting E3SM Atmosphere Model in a Doubly-Periodic Configuration

We develop a doubly periodic version of the Simple Convection-Permitting E3SM Atmosphere Model (SCREAM) to provide an efficient configuration for this global convection permitting model (GCPM), akin to a single column model often found in conventional general circulation models. The design details are explained, in addition to the extensive case library associated with the doubly periodic SCREAM (DP-SCREAM) configuration. We demonstrate that doubly periodic cloud resolving models are useful tools to explore the horizontal resolution sensitivity of GCPMs, in addition to replicating biases seen in the global models. Using DP-SCREAM, we show that SCREAM exhibits behaviors of a scale aware model as it is able to naturally partition between sub-grid scale (SGS) and resolved vertical transport across the gray zone of turbulence. We show that SCREAM is reasonably scale insensitive when run at resolutions from 1 to 5 km, but can exhibit sensitivity, particularly for the shallow convective regime, when run at resolutions approaching that of large eddy simulations. We conclude that SGS parameterization improvements are likely needed to reduce this scale sensitivity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

How R Developers explain their Package Choice: A Survey

Background: Contemporary software development relies heavily on reusing already implemented functionality, usually in the form of packages. Aims: We aim to shed light on developers' preferences when selecting packages in R language. Method: To do that, we create and administer a survey to over 1000 developers who have added one of two common dataframe enhancement libraries in R to their projects: data.table or tidyr. We design a questionnaire using the Social Contagion Theory (SCT) following prior work on technology adoption and ensure that key dimensions affecting developer choice are considered. Results: Of the 1085 developers we contacted, 803 completed the survey asking them to prioritize various factors known to affect developer perceptions of package quality and to provide their background. Most developers self-identified as data scientists with two to five years of work experience. We found significant differences between the preferences of developers who chose data.table and tidyr. Surprisingly, package reputation based on easy-to-see measures, such as the number of stars on GitHub, was not an important factor for either group. Conclusions: Our findings demonstrate the inherently social nature of package adoption. They can help design future studies on how different populations of developers make decisions on which software packages to use in their projects. Finally, package developers and maintainers can benefit by better understanding the prime concerns of the users of their packages.

Malviya Thakur, Addi↗
DOE OSTI2024

hypredrive: high-level interface for solving linear systems with hypre

This software introduces a high-level interface designed to simplify solving linear systems using hypre, a renowned library for such computational challenges. It is crafted to be accessible and user-friendly, making the powerful capabilities of hypre available to a broader audience without requiring in-depth technical knowledge. The interface is characterized by its use of YAML for input, a format celebrated for its structured yet straightforward readability. This choice ensures that users can easily configure the software to meet their specific needs. Additionally, the software boasts an intuitive API that encapsulates hypre's functionalities, making it easier for users to interact with the process of solving linear systems. It is particularly beneficial for prototyping, offering a quick and efficient means to test various solver and preconditioner configurations. Furthermore, the software allows for the creation of an offline testing framework in which predefined linear systems are read from files and benchmarked with user-defined solution strategies. This makes it an invaluable tool for developers and researchers exploring and validating their computational models. Overall, the software serves as a bridge, bringing the advanced computational capabilities of hypre closer to users who may need more specialized technical expertise, thereby facilitating innovation and exploration in the field of numerical linear algebra.

Paludetto Magri, Victor↗
DOE OSTI2022

National Security Research Center (NSRC) Overview

The National Security Research Center (NSRC) is Los Alamos' classified library. It houses 75+ years of scientific and engineering research, designs, procedures, videos, photos, and other reports.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2023

High-throughput printing of combinatorial materials from aerosols

The development of new materials and their compositional and microstructural optimization are essential in regard to next-generation technologies such as clean energy and environmental sustainability. However, materials discovery and optimization have been a frustratingly slow process. The Edisonian trial-and-error process is time consuming and resource inefficient, particularly when contrasted with vast materials design spaces 1 . Whereas traditional combinatorial deposition methods can generate material libraries 2,3 , these suffer from limited material options and inability to leverage major breakthroughs in nanomaterial synthesis. Here we report a high-throughput combinatorial printing method capable of fabricating materials with compositional gradients at microscale spatial resolution. In situ mixing and printing in the aerosol phase allows instantaneous tuning of the mixing ratio of a broad range of materials on the fly, which is an important feature unobtainable in conventional multimaterials printing using feedstocks in liquid–liquid or solid–solid phases 4–6 . We demonstrate a variety of high-throughput printing strategies and applications in combinatorial doping, functional grading and chemical reaction, enabling materials exploration of doped chalcogenides and compositionally graded materials with gradient properties. The ability to combine the top-down design freedom of additive manufacturing with bottom-up control over local material compositions promises the development of compositionally complex materials inaccessible via conventional manufacturing approaches.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Enabling particle applications for exascale computing platforms

The Exascale Computing Project (ECP) is invested in co-design to assure that key applications are ready for exascale computing. Within ECP, the Co-design Center for Particle Applications (CoPA) is addressing challenges faced by particle-based applications across four “sub-motifs”: short-range particle–particle interactions (e.g., those which often dominate molecular dynamics (MD) and smoothed particle hydrodynamics (SPH) methods), long-range particle–particle interactions (e.g., electrostatic MD and gravitational N-body), particle-in-cell (PIC) methods, and linear-scaling electronic structure and quantum molecular dynamics (QMD) algorithms. Our crosscutting co-designed technologies fall into two categories: proxy applications (or “apps”) and libraries. Proxy apps are vehicles used to evaluate the viability of incorporating various types of algorithms, data structures, and architecture-specific optimizations and the associated trade-offs; examples include ExaMiniMD, CabanaMD, CabanaPIC, and ExaSP2. Libraries are modular instantiations that multiple applications can utilize or be built upon; CoPA has developed the Cabana particle library, PROGRESS/BML libraries for QMD, and the SWFFT and fftMPI parallel FFT libraries. Success is measured by identifiable “lessons learned” that are translated either directly into parent production application codes or into libraries, with demonstrated performance and/or productivity improvement. The libraries and their use in CoPA’s ECP application partner codes are also addressed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

GP Cosmology Surrogate v1.0

GP Cosmology Surrogate is a Python library for building and training a generalized multi-output Gaussian process (GP) framework of @takhtaganov2021cosmic. In this approach, the surrogate is constructed sequentially, guided by a Bayesian optimization acquisition function that targets reduction of emulation error in the regions most consistent with the observational data. This adaptive design concentrates computational resources where they have the greatest impact on inference accuracy. The library supports efficient training for separable GP kernels, which allows the use of Kronecker algebra to handle high-dimensional input spaces and large numbers of correlated outputs. This makes it well suited for applications such as modeling cosmological power spectra, large-scale physical simulations, and multi-output hyperparameter tuning. By combining scalable multi-output GP modeling with data-driven adaptive sampling, GPsurrogate enables parameter inference and optimization with substantially fewer simulations than conventional space-filling designs.

Lukic, Zarija [Lawrence Berkeley National Laborato↗
DOE OSTI2020

Safety Risk Reliability Model Library

SR2ML is a software package which contains a set of safety and reliability models designed to be interfaced with the INL developed RAVEN code. These models can be employed to perform both static and dynamic system risk analysis and determine risk importance of specific elements of the considered system. Two classes of reliability models have been developed; the first class includes all classical reliability models (Fault-Trees, Event-Trees, Markov models and Reliability Block Diagrams) which have been extended to deal not only with Boolean logic values but also time dependent values. The second class includes several components aging models. Models of these two classes are designed to be included in a RAVEN ensemble model to perform time dependent system reliability analysis (dynamic analysis). Similarly, these models can be interfaced with system analysis codes to determine failure time of systems and evaluate accident progression (static analysis).

Wang, Congjian↗
DOE OSTI2020

NRC Multiphysics Analysis Capability Deployment FY2020 - Part 4

This report details progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Sup- port for Advanced Non-Light Water Reactors.” The tasks completed for this report are: • Task 3a: This task demonstrates the effectiveness of Griffin’s methods for interpolating cross sections and SPH correction factors as a function of control drum rotation. The Monte Carlo code Serpent is used to generate reference results. The accuracy of Griffin-calculated eigenvalues and reaction rates are studied as a function of control drum rotation for cross section and SPH factor libraries of varying fidelity. Higher-fidelity libraries require less interpolation and are therefore inherently more accurate; however, they require more computational resources to create. Overall, Griffin is able to accurately model eigenvalue and reaction rates over the full range of drum rotation using a library generated from three to five discrete drum rotations for a particular microreactor design. • Task 3c: This preliminary task successfully demonstrates the performance of a coupled neutron-photon transport calculation with Griffin using the spherical harmonics approximation. The Griffin solutions are compared to reference solutions from the continuous energy Monte Carlo codes MCNP and Serpent. Furthermore, this work also demonstrates how to prepare the neutron and photon libraries using NJOY, MCNP, and Serpent with the latest ENDF/B-VIII data for 235 U and graphite. The Griffin total energy deposition in the active core region is calculated within 0.3% of the reference, but it degrades in the reflector regions with a maximum difference of -1.5%. Overall, the results look promising and future analysis should use an improved data stream and focus on determining the potential cancellation of error within each energy group.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Impact of nuclear data covariance libraries on uncertainty quantification of sodium cooled fast reactor simulation

Uncertainty quantification in core modelling to assess reactor performance is a regular practice to identify design safety limits. Best-estimate plus uncertainty (BEPU) methods are commonly used to propagate uncertainties to core responses from input parameters, characterized by variance-covariance libraries. This study evaluates the impact of different covariance libraries on the uncertainty associated with core safety parameters and their propagation through reactor transient modelling for a metallic core (ABR-1000). The COMAC-1.0, COMMARA-2.0, and COMMARA-2.1 covariance libraries are used to propagate nuclear data uncertainties to neutron feedback coefficients using the Generalized Perturbation Theory, then through transient safety analyses to establish the confidence interval for safety performance of the reactor using stochastic sampling. The considered transients include Unprotected Transient Over Power and Unprotected Loss Of Flow. Most influential nuclide and reaction pairs in different variance-covariance libraries are identified and their impact on the uncertainty of safety parameters are evaluated. Significant differences were observed between uncertainties from COMAC compared to COMMARA libraries. Uncertainties from COMMARA libraries were consistently larger than COMAC for all cases. For structure, coolant density feedback coefficients and control rod worth, total nuclear data uncertainties were approximately 50% larger for COMMARA than COMAC. Differences were also observed in identification of top contributors of uncertainty in terms of reaction-cross sections. For instance, largest contribution to total uncertainty of radial expansion coefficient from COMMARA-2.0 and 2.1 originates from Na{sup 23} P1 Elastic cross section (0.81%) but from Pu{sup 239} Fission (0.53%) with COM C-1.0. Such differences are evaluated in this research for all quantities of interest. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

CINDER90 and CINDER2008 Comparison for Second Target Station Analysis

A preliminary target activation analysis for the Second Target Station (STS) at the Spallation Neutron Source (SNS) located at Oak Ridge National Laboratory (ORNL) was preformed using Monte Carlo N-Particle X (MCNPX), CINDER90, and the corresponding multi-group cross section libraries distributed with CINDER90. This preliminary analysis has been used to guide some design decisions of the STS target design. Using the MCNPX output files from this preliminary analysis, a comparison of the activation analysis using CINDER90 and CINDER2008 was performed. While the CINDER90 code has been validated for use in spallation source activation and transmutation calculations, and compares well with measurements, the transition to the modernized version of the code (CINDER2008) and the updated cross section libraries is desirable. The goal of this comparison is to evaluate the differences and implications for the STS activation analysis. Overall, the total activity differs between CINDER2008 and CINDER90 by no more than 20% for the materials evaluated in this report. The largest discrepancies were observed in the decay gamma intensity densities with CINDER90 calculating a factor of 2 higher than CINDER2008.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Signal Decomposition for Intrusion Detection in Reliability Assessment in Cyber Resilience (Summary Report)

The complexity of assuring cyber resilience for physical process interactions in connected systems such as energy grids increases dramatically as the coupling between processes becomes more direct and responsive. An example of this growing complexity is provided by Integrated Energy Systems (IES), in which various processes such as nuclear heat generation and commodity production are being directly coupled for increased responsiveness to highly variable signals such as market pricing or electricity demand. As such, the potential attack surface of the coupled processes is larger than the two processes independently. Securing these complex systems requires two-fold monitoring: cybersecure monitoring for potential malicious incursion, and physics monitoring for system tampering. Physics monitoring includes analyzing the behavior of the signals within the system for anomalous behavior. This analysis has been shown to be insufficient if approached by only data-driven machine learning and artificial intelligence (MLAI) techniques or only low-level model comparison. Previous efforts at Purdue University suggested combining high-fidelity models with MLAI algorithms as a basis for a software tool for detecting anomalies in physical processes. This work built on that suggestion, developing an advanced library for signal decomposition and analysis using both MLAI and high-fidelity physics algorithms for greatly improved anomaly detection, especially false data injection. This software can be used as part of a secure imbedded intelligence (SEI) system designed under Consequence-driven Cyber-informed Engineering (CCE) for complex coupled systems. This library established a foundation for online and posteriori analysis of digital signals for the purpose of detecting potential malicious tampering in digital signals representing physical processes. Demonstrations carried out throughout the development highlight the effective use of characterization algorithms to detect signal perturbations, particularly triangle attack-style perturbations, in three wide-ranging applications: seismic monitoring, nuclear thermal hydraulics system simulation, and custom manufacturing.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Bridging semantics, control specifications and assessment: A library for scalable demand flexibility controls

There is growing recognition that Demand Flexibility (DF) can play a major role in enhancing grid reliability, with building control applications emerging as key enablers for DF. However, the traditional approach to deploying new control applications in buildings, including those for DF, remains largely manual and tailored to individual buildings, making it difficult to scale. While research efforts have explored semantics-driven portability, DF controls specification, and assessment approaches, these initiatives are fragmented and limited in scope. This paper proposes a novel methodology, grounded in design science research, to integrate these elements and create a comprehensive DF controls library for both industry and academia. This approach is applied to develop the Demand FLEXibility controls LIBrary using Semantics (DFLEXLIBS), an extensible open-source library that provides DF controls for HVAC systems in Python. DFLEXLIBS enables portable, easy-to-deploy controls that abstract building-specific data points, facilitating assessment across diverse buildings. DFLEXLIBS features nine different control applications, and it is successfully implemented and tested across four virtual and two real buildings, bridging the gap between semantics-driven portability, DF controls specification, and rigorous performance assessment. Its benefits are measured by a reusability ratio greater than 90% and a functional overlap ratio of around 70% for the most common functions used in the library, significantly reducing time for deploying new controls.

Controls library↗