Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “open source software”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Evaluation of diffusion kurtosis and diffusivity from baseline staging MRI as predictive biomarkers for response to neoadjuvant chemoradiation in locally advanced rectal cancer

Purpose: To evaluate the role of diffusion kurtosis and diffusivity as potential imaging biomarkers to predict response to neoadjuvant chemoradiation therapy (CRT) from baseline staging magnetic resonance imaging (MRI) in locally advanced rectal cancer (LARC). Materials and methods: This retrospective study included 45 consecutive patients (31 male/14 female) who underwent baseline MRI with high b-value sequences (up to 1500 mm/s{sup 2}) for LARC followed by neoadjuvant chemoradiation and surgical resection. The mean age was 57.4 years (range 34.2–72.9). An abdominal radiologist using open source software manually segmented T2-weighted images. Segmentations were used to derive diffusion kurtosis and diffusivity from diffusion-weighted images as well as volumetric data. These data were analyzed with regard to tumor regression grade (TRG) using the four-tier American Joint Committee on Cancer (AJCC) classification, TRG 0–3. Proportional odds regression was used to analyze the four-level ordinal outcome. A sensitivity analysis was performed using univariable logistic regression for binary TRG groups, TRG 0/1 (> 90% response), or TRG 2/3 (< 90% response). p < 0.05 was considered significant throughout. Results: In the univariable proportional odds regression analysis, higher diffusivity summary (D{sub sum}) values were observed to be significantly associated with higher odds of being in one or more favorable TRG group (TRG 0 or 1). In other words, on average, patients with higher D{sub sum} values were more likely to be in a more favorable TRG group. These results are mostly consistent with the sensitivity analysis, in which higher values for most D{sub sum} values [all but region of interest (ROI)-max D median (p = 0.08)] were observed to be significantly associated with higher odds of being TRG 0 or 1. Tumor volume of interest (VOI) and ROI volume, ROI kurtosis mean and median, and VOI kurtosis mean and median were not significantly associated with TRG. Conclusion: Diffusivity derived from the baseline staging MRI, but not diffusion kurtosis or volumetric data, is associated with TRG and therefore shows promise as a potential imaging biomarker to predict the response to neoadjuvant chemotherapy in LARC. Clinical relevance statement: Diffusivity shows promise as a potential imaging biomarker to predict AJCC TRG following neoadjuvant CRT, which has implications for risk stratification. Patients with TRG 0/1 have 5-year disease-free survival (DFS) of 90–98%, as opposed to those who are TRG 2/3 with 5-year DFS of 68–73%.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2017

3D Histopathology—a Lung Tissue Segmentation Workflow for Microfocus X-ray-Computed Tomography Scans

Lung histopathology is currently based on the analysis of 2D sections of tissue samples. The use of microfocus X-ray-computed tomography imaging of unstained soft tissue can provide high-resolution 3D image datasets in the range of 2–10 μm without affecting the current diagnostic workflow. Important details of structural features such as the tubular networks of airways and blood vessels are contained in these datasets but are difficult and time-consuming to identify by manual image segmentation. Providing 3D structures permits a better understanding of tissue functions and structural interrelationships. It also provides a more complete picture of heterogeneous samples. In addition, 3D analysis of tissue structure provides the potential for an entirely new level of quantitative measurements of this structure that have previously been based only on extrapolation from 2D sections. In this paper, a workflow for segmenting such 3D images semi-automatically has been created using and extending the ImageJ open-source software and key steps of the workflow have been integrated into a new ImageJ plug-in called LungJ. Results indicate an improved workflow with a modular organization of steps facilitating the optimization for different sample and scan properties with expert input as required. This allows for incremental and independent optimization of algorithms leading to faster segmentation. Representation of the tubular networks in samples of human lung, building on those segmentations, has been demonstrated using this approach.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2018

Retrospective Analysis of Clinical Performance of an Estonian Speech Recognition System for Radiology: Effects of Different Acoustic and Language Models

The aim of this study was to analyze retrospectively the influence of different acoustic and language models in order to determine the most important effects to the clinical performance of an Estonian language-based non-commercial radiology-oriented automatic speech recognition (ASR) system. An ASR system was developed for Estonian language in radiology domain by utilizing open-source software components (Kaldi toolkit, Thrax). The ASR system was trained with the real radiology text reports and dictations collected during development phases. The final version of the ASR system was tested by 11 radiologists who dictated 219 reports in total, in spontaneous manner in a real clinical environment. The audio files collected in the final phase were used to measure the performance of different versions of the ASR system retrospectively. ASR system versions were evaluated by word error rate (WER) for each speaker and modality and by WER difference for the first and the last version of the ASR system. Total average WER for the final version throughout all material was improved from 18.4% of the first version (v1) to 5.8% of the last (v8) version which corresponds to relative improvement of 68.5%. WER improvement was strongly related to modality and radiologist. In summary, the performance of the final ASR system version was close to optimal, delivering similar results to all modalities and being independent on user, the complexity of the radiology reports, user experience, and speech characteristics.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2024

pyEGAF: Modernization of the EGAF database

One of the most comprehensive resources for thermal neutron-capture data is the Evaluated Gamma-ray Activation File (EGAF), containing data from prompt gamma activation analysis measurements carried out at the Budapest Research Reactor for 245 isotopes. Although these valuable datasets have been freely available for many years, the outdated and cryptic adopted format makes it difficult to utilize the data and it is not generally suitable for modern computational technologies. Furthermore, to help overcome these challenges, we have converted the datasets into an open standard JSON format. Additionally, we have developed a Python implementation of an open-source software package, pyEGAF, designed to interact with the JSON data structures for general purpose access, manipulation, and rapid assessment of the capture-gamma data in EGAF.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

An efficient three-dimensional rhizosphere modeling capability to study the effect of root system architecture on soil water and reactive transport

This research was to develop a rhizosphere modeling capability for plant-soil interactions by integrating water and ion uptake and release from three dimensional (3D) root systems and 3D variably saturated flow and multicomponent reactive transport in soil. We combined open source software for simulating plant and soil interactions with parallel processing computational technologies to address highly-resolved root system architecture and soil hydrobiogeochemical processes. Simulation capability was demonstrated on Brachypodium distachyon . Our demonstration problem showed that availability of water and plant nutrients is controlled by the interplay between 1) transpiration-driven cycles of water uptake, root zone saturation and desaturation; 2) sorption and desorption via competitive ion exchange; 3) buildup of ions not taken up during kinetic uptake of nutrients/ions from pore water; and 4) advection, dispersion, and diffusion of ions in the soil. The uptake of water and ions by individual roots leads to dynamic, local gradients in ion concentrations. The resulting higher solution ionic strength can lead to increasing sorption of K+ and NH4+ on the soil exchange sites, despite their preferential absorption by roots. The modeling framework makes it possible to explore alternative conceptual models and improve the understanding of the biogeochemical and physical environment within the rhizosphere.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Scalable branching on dual decomposition of stochastic mixed-integer programming problems

In this work, we present a scalable branching method for the dual decomposition of stochastic mixed-integer programming. Our new branching method is based on the branching method proposed by Caroe and Schultz that creates branching disjunctions on first-stage variables only. We propose improvements to the process for creating branching disjunctions, including (1) branching on the optimal solutions of the Dantzig-Wolfe reformulation of the restricted master problem and (2) using a more comprehensive (yet simple) measure for the dispersions associated with subproblem solution infeasibility. We prove that the proposed branching process leads to an algorithm that terminates finitely, and we provide conditions under which globally optimal solutions can be identified after termination. We have implemented our new branching method, as well as the Caroe-Schultz method and a branch-and-price method, in the open-source software package DSP. Using SIPLIB test instances, we present extensive numerical results to demonstrate that the proposed branching method significantly reduces the number of node subproblems and solution times.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Rucio: Scientific Data Management

Rucio is an open-source software framework that provides scientific collaborations with the functionality to organize, manage, and access their data at scale. The data can be distributed across heterogeneous data centers at widely distributed locations. Rucio was originally developed to meet the requirements of the high-energy physics experiment ATLAS, and now is continuously extended to support the LHC experiments and other diverse scientific communities. In this article, we detail the fundamental concepts of Rucio, describe the architecture along with implementation details, and report operational experience from production usage.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2018

Role of scan strategies on thermal gradient and solidification rate in electron beam powder bed fusion

Local microstructure control in electron beam powder bed fusion (EB-PBF) is of great interest to the additive manufacturing community to realize complex part geometry with targeted performance. The local microstructure control relies on having a detailed understanding of local melt pool physics (e.g., 3-D melt pool shape as well as spatial and temporal variations of thermal gradient (G) and solidification rate (R)). In this research, a new scan strategy referred to as ghost beam is numerically evaluated as a candidate to achieve the targeted G and R of IN718 alloy. The boundary conditions for simulations, including the speed (490 mm/s) and spatial locations of the beam within a given layer, are obtained by using series of snapshot images, recorded at 12,000 frames per second, using a high-speed camera. The heat transfer simulations were performed using TRUCHAS an open-source software deployed within a high-performance computational infrastructure. The simulation results showed that reheating at short beam on-time and time delay decreases both G and R. Local variation of R at the center of the melt pool trailing edge showed periodic temporal fluctuations. Finally, the ghost beam scan strategy was compared to other existing raster and spot scan strategies.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Westcott g factors extended to arbitrary neutron energy spectra

Westcott 𝑔 factors are used in Neutron Activation Analysis (NAA) and Prompt Gamma-ray Activation Analysis (PGAA) to evaluate the impact of non-1∕𝑣 behavior in the neutron-capture cross sections of certain nuclei on activation product yields. This non-1∕𝑣 behavior arises from the presence of neutron resonances in the neutron- capture cross sections that overlap with the source neutron spectrum at low (< 5 eV) energies. Historically, Westcott 𝑔 factors that have been cataloged for NAA and PGAA applications are the result of calculations that assume a Maxwellian neutron flux distribution with a given temperature. In this work, we use this approach with updated neutron-capture cross sections from the Evaluated Nuclear Data File, version VIII.1 (ENDF/B-VIII.1) to tabulate Westcott 𝑔 factor values for a broad range of Maxwellian distribution temperatures, comparing the results against currently-available 𝑔 factors from International Atomic Energy Agency tables and other sources. Here, it was discovered during this analysis that the use of guided thermal and cold-neutron beams at certain facilities necessitates an approach for evaluating Westcott 𝑔 factors based on arbitrary non-Maxwellian spectra. In this paper, we present an approach for calculating 𝑔 factors with user-specified neutron spectra, and we demonstrate these methods to obtain Westcott 𝑔-factors for guided- and cold-neutron beams at the Budapest Research Reactor and the Forschungsreaktor München II reactor. As part of this work, open-source software has been developed that can be used to perform these calculations for applications in PGAA and NAA experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Definition and solution of a stochastic inverse problem for the Manning’s n parameter field in hydrodynamic models

The uncertainty in spatially heterogeneous Manning’s n fields is quantified using a novel formulation and numerical solution of stochastic inverse problems for physics-based models. The uncertainty is quantified in terms of a probability measure and the physics-based model considered here is the state-of-the-art ADCIRC model although the presented methodology applies to other hydrodynamic models. An accessible overview of the formulation and solution of the stochastic inverse problem in a mathematically rigorous framework based on measure theory is presented in this paper. Technical details that arise in practice by applying the framework to determine the Manning’s n parameter field in a shallow water equation model used for coastal hydrodynamics are presented and an efficient computational algorithm and open source software package are developed. A new notion of “condition” for the stochastic inverse problem is defined and analyzed as it relates to the computation of probabilities. Finally, this notion of condition is investigated to determine effective output quantities of interest of maximum water elevations to use for the inverse problem for the Manning’s n parameter and the effect on model predictions is analyzed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Generation of random geological models using multi-randomization for machine learning

Generating high-fidelity geological models is essential for advancing machine learning (ML) methods in automated seismic interpretation. For instance, seismic images paired with corresponding fault labels are foundational for ML-based fault detection from seismic migration sections. While several open-access datasets of random geological models exist, open-source tools specifically designed to produce large volumes of such models for ML applications remain scarce. To address this gap, we present RGM (Random Geological Model), an open-source software package for efficiently generating 2D and 3D synthetic geological models tailored for ML workflows. RGM supports the creation of diverse model components, including medium property distributions (P-/S-wave velocities and density), seismic reflectivity images (i.e., synthetic migration sections), relative geological time, and discrete fault attributes such as probability, dip, strike, rake, and displacement. It also accommodates the creation of complex geological features such as salt bodies and unconformities. The model generation algorithm employs a multi-randomization strategy, yielding an effectively infinite-dimensional model space that encompasses a wide range of geological scenarios and associated seismic features. Furthermore, RGM incorporates a method to generate synthetic elastic migration images using analytical elastic reflection coefficients combined with frequency-dependent scaling. This functionality enables the creation of training datasets for ML models that leverage elastic seismic images. RGM is implemented in modern object-oriented Fortran, allowing users to flexibly control statistical parameters governing model variability. We demonstrate the capability, performance, and geological realism of the package through comprehensive 2D and 3D examples.

58 GEOSCIENCES↗
DOE OSTI2025

Thermodynamic modeling of countercurrent chemical looping reverse water gas shift process for redox material screening

The reverse water gas shift (RWGS) reaction is a key pathway for CO 2 utilization, particularly within Power-to-X process chains aimed at sustainable fuel and chemical production. Countercurrent chemical looping (CL-RWGS) using non-stoichiometric oxides can overcome equilibrium limitations of conventional RWGS reactors, enabling significantly higher CO 2 conversions. However, modeling the limiting performance of such systems is challenging due to their multiphase nature and coupled spatial and temporal variation in chemical composition. In this work, we present a discretized batch equilibrium model that simulates CL-RWGS reactors as a series of localized equilibrium exchanges between gas and solid elements. The model is numerically stable, computationally efficient, and free of kinetic source terms, making it well-suited for parametric studies and system-level integration. It is validated against established convection–diffusion models and shown to predict reasonable upper bounds on experimental results. Application of the model to a range of oxygen carrier materials identifies cerium–zirconium solid solutions, particularly Ce 0.80 Zr 0.20 O 2 , as a promising class offering superior oxygen storage characteristics compared to state-of-the-art La 0.6 Sr 0.4 FeO 3 . This framework provides a robust platform for materials screening, reactor sizing, and performance optimization in chemical looping systems. The model implementation is available as open-source software to support further research and development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Computational studies of impurity migration during induction stirring of molten uranium

Understanding and controlling impurity behavior in actinide metal casting processes are foundational for efficient part production yet remain major challenges for researchers and industry. To help provide insight regarding impurity distribution during actinide metal casting, we have developed computational fluid dynamics (CFD) models for a laboratory-scale system using commercial and open-source code. Multiple simulation frameworks allow for improved confidence in the resulting outputs, while taking advantage of the maturity and convenience offered by commercial platforms and simultaneously maintaining the transparency and flexibility often provided by open-source software. Here, we describe multiple experiment-informed models designed to simulate a specific laboratory system in which uranium melt, containing a known starting concentration of carbon impurity, is electromagnetically stirred in an induction furnace. Here, the goal of the simulations is to predict the motion of uranium carbide microparticles in the velocity field of the melt. Prior to simulating the uranium-carbon system, numerical models were validated using a previously published nonradioactive experimental system. Effects of the size and shape of impurity particles in the models were investigated and agree with experimental findings. Simulation of smaller particles (< 50 µm) shows more homogenous distribution throughout the stirred melt. With increased particle size (100 µm), the body forces, which include the buoyancy force, dominate over the drag force, causing larger particles to move toward the crucible walls and upward in the system.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Highly-scalable GPU-accelerated compressible reacting flow solver for modeling high-speed flows

Emerging supercomputing systems utilize a combination of central processing units (CPUs) and graphics processing units (GPUs) in an effort to reach exascale capabilities while minimizing the energy footprint of operating such systems. Such heterogeneous machines introduce new challenges for fluids solvers because the hardware architecture and operation of a GPU are fundamentally different from conventional CPUs. In this work, a general approach for efficient implementation of finite-volume based reacting flow solvers on such heterogeneous systems is presented. Three main challenges, namely, data access pattern, thread divergence, and thread safety, are addressed. Since compressible reacting flows require special methods to deal with chemical reactions, hyperbolic and nonlinear convection terms, and the presence of turbulence, specific algorithms that ensure GPU-based efficiency are developed. The approach is demonstrated on the widely available OpenFOAM open source software by modifying core algorithms for GPU accessibility. The scalability of the resulting solver, is demonstrated using practical test cases, including flow through a scramjet engine and the dynamics of a rotating detonation engine. Here, the solver provides near-ideal scaleup on a large number of GPUs (>3000), and extremely efficient use of the GPUs, with throughput nearly a constant even when processing a large number of control volumes.

42 ENGINEERING↗
DOE OSTI2017

Improved treatment of exact exchange in Quantum ESPRESSO

Here, we present an algorithm and implementation for the parallel computation of exact exchange in Quantum ESPRESSO (QE) that exhibits greatly improved strong scaling. QE is an open-source software package for electronic structure calculations using plane wave density functional theory, and supports the use of local, semi-local, and hybrid DFT functionals. Wider application of hybrid functionals is desirable for the improved simulation of electronic band energy alignments and thermodynamic properties, but the computational complexity of evaluating the exact exchange potential limits the practical application of hybrid functionals to large systems and requires efficient implementations. We demonstrate that existing implementations of hybrid DFT that utilize a single data structure for both the local and exact exchange regions of the code are significantly limited in the degree of parallelization achievable. We present a band-pair parallelization approach, in which the calculation of exact exchange is parallelized and evaluated independently from the parallelization of the remainder of the calculation, with the wavefunction data being efficiently transformed on-the-fly into a form that is optimal for each part of the calculation. For a 64 water molecule supercell, our new algorithm reduces the overall time to solution by nearly an order of magnitude.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

ELSI — An open infrastructure for electronic structure solvers

Routine applications of electronic structure theory to molecules and periodic systems need to compute the electron density from given Hamiltonian and, in case of non-orthogonal basis sets, overlap matrices. System sizes can range from few to thousands or, in some examples, millions of atoms. Different discretization schemes (basis sets) and different system geometries (finite non-periodic vs. infinite periodic boundary conditions) yield matrices with different structures. The ELectronic Structure Infrastructure (ELSI) project provides an open-source software interface to facilitate the implementation and optimal use of high-performance solver libraries covering cubic scaling eigensolvers, linear scaling density-matrix-based algorithms, and other reduced scaling methods in between. In this paper, we present recent improvements and developments inside ELSI, mainly covering (1) new solvers connected to the interface, (2) matrix layout and communication adapted for parallel calculations of periodic and/or spin-polarized systems, (3) routines for density matrix extrapolation in geometry optimization and molecular dynamics calculations, and (4) general utilities such as parallel matrix I/O and JSON output. The ELSI interface has been integrated into four electronic structure code projects (DFTB+, DGDFT, FHI-aims, SIESTA), allowing us to rigorously benchmark the performance of the solvers on an equal footing. Based on results of a systematic set of large-scale benchmarks performed with Kohn–Sham density-functional theory and density-functional tight-binding theory, we identify factors that strongly affect the efficiency of the solvers, and propose a decision layer that assists with the solver selection process. As a result, we describe a reverse communication interface encoding matrix-free iterative solver strategies that are amenable, e.g., for use with planewave basis sets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

NESSi : The N on- E quilibrium S ystems S imulation package

The nonequilibrium dynamics of correlated many-particle systems is of interest in connection with pump–probe experiments on molecular systems and solids, as well as theoretical investigations of transport properties and relaxation processes. Nonequilibrium Green’s functions are a powerful tool to study interaction effects in quantum many-particle systems out of equilibrium, and to extract physically relevant information for the interpretation of experiments. Here, we present the open-source software package NESSi (The Non-Equilibrium Systems Simulation package) which allows to perform many-body dynamics simulations based on Green’s functions on the L-shaped Kadanoff–Baym contour. NESSi contains the library libcntr which implements tools for basic operations on these nonequilibrium Green’s functions, for constructing Feynman diagrams, and for the solution of integral and integro-differential equations involving contour Green’s functions. The library employs a discretization of the Kadanoff–Baym contour into time points and a high-order implementation of integration routines. The total integrated error scales up to $\mathcal{O}(N^{-7})$, which is important since the numerical effort increases at least cubically with the simulation time. A distributed-memory parallelization over reciprocal space allows large-scale simulations of lattice systems. We provide a collection of example programs ranging from dynamics in simple two-level systems to problems relevant in contemporary condensed matter physics, including Hubbard clusters and Hubbard or Holstein lattice models. The libcntr library is the basis of a follow-up software package for nonequilibrium dynamical mean-field theory calculations based on strong-coupling perturbative impurity solvers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Methodology for physics-informed generation of synthetic neutron time-of-flight measurement data

Accurate neutron cross section data are a vital input to the simulation of nuclear systems for a wide range of applications from energy production to national security. The evaluation of experimental data is a key step in producing accurate cross sections. There is a widely recognized lack of reproducibility in the evaluation process due to its artisanal nature and therefore there is a call for improvement within the nuclear data community. This can be realized by automating/standardizing viable parts of the process, namely, parameter estimation by fitting theoretical models to experimental data. This automation effort could greatly benefit from a synthetic data resource. This work leverages problem-specific physics, Monte Carlo sampling, and a general methodology for data synthesis to generate unlimited, labelled experimental cross-section data that is statistically indistinguishable to the observed data. Heuristic and, where applicable, rigorous statistical comparisons to observed data support this claim. The demonstration is based on/limited to transmission measurements at Rensselaer Polytechnic Institute (RPI) and energy-differential cross sections in the resolved resonance region (RRR). An open-source software is published alongside this article that executes the complete methodology to produce high-utility synthetic datasets. The goal of this work is to provide an approach and corresponding tool that will allow the evaluation community to begin exploring more data-driven, ML-based solutions to long-standing challenges in the field.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗