SEARCH · Search NASA
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.
WNTR is a python package designed to simulate and analyze resilience of water distribution networks. The software includes: - Pressure driven and demand driven hydraulic simulation - Water quality simulation to track concentration, trace, and water age - Conditional controls to simulate power outages - Models to simulate pipe breaks - A wide range of resilience metrics - Analysis and visualization tools
This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in a forested catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, ground penetrating radar (GPR), continuous self-potential (SP) monitoring, electromagnetic (EM) imaging, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Contains two subfolders: Synthetic and Field_Application subfolder. Synthetic subfolder contains the ATS XML input script (can be opened using any code editor) for the four synthetic hydrological cases tested (Connected and gaining, Connected and losing, Disconnected and losing, and dry stream). It also includes other experimental cases to test the influence of precipitation and concentration gradient. For each synthetic case, the flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.mph can be opened with the commercial software COMSOL and requires a license) is executed using the ATS output data to simulate the potential field. It also includes the Synthetic_model_plot.ipynb (can be opened using any code editor) to visualize the SP result and generate manuscript figures. The data subfolder contains mesh files to run both the ATS (.exo and .stl files can be viewed using Paraview; .h5 files can be opened using HDFView software and h5py Python package) and COMSOL models. Field_Application subfolder contains two subfolders: ES_MDA_inversion and Final_Model. ES_MDA_inversion contains the Python script (.py can be opened using any code editor) and SP observation data used to run the Ensemble Smoother with Multiple Data Assimilation (ES-MDA) inversion sequence to get the optimal model parameters. The Final_model subfolder contains the ATS XML input scripts, data files, output data for the two SP sites. The same workflow steps outlined for the Synthetic subfolder apply here. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) EM Contains the CSV file of the EM data from the DUALEM-42, including spatial coordinates (x, y, z), apparent conductivity, and in-phase measurements at 2 m coil separations for horizontal coplanar (HCP) and perpendicular (PRP) geometries. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion (.resipy can be opened with the open-source ResIPy software). 6) GPR Includes GPR field datasets collected at 100 MHz and 250 MHz antenna frequencies, along with the processing/interpretation project file (GPR_process.gpz can be viewed using EKKO_Project 6, a commercial software by Sensors & Software that requires a license). 7) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 8) SP Contains the SP data collected in field at the two SP sites (one in the perennial reach and the other in the intermittent reach), provided as DAT files. 9) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). It also includes DTW.ipynb, a Jupyter notebook containing the code for the dynamic time warping (DTW) with sliding window to evaluate SP signal synchronicity.
High-resolution simulations of a gas-solids jet in a 0.3 m diameter and 15.9 m tall circulating fluidized bed (CFB) riser were conducted with the open source software-MFIX. In the numerical simulations, both gas and solids injected through a 1.6 cm diameter radial-directed tube 4.3 m above the bottom distributor were tracked as tracers, which enable the analysis of the characteristics of a two-phase jet. Two jetting gas velocities of 16.6 and 37.2 m/s were studied with the other operating conditions fixed. Reasonable flow hydrodynamics with respect to overall pressure drop, voidage, and solids velocity distributions were predicted. Due to the different dynamic responses of gas and particles to the crossflow, a significant separation of gas and solids within the jet region was predicted for both cases. In addition, the jet characteristics based on tracer concentration and tracer mass fraction profiles at different downstream levels are discussed. Overall, the numerical predictions compare favorably to the experimental measurements made at NETL.
Detector simulation is a key component for studies on prospective future high-energy colliders, the design, optimization, testing and operation of particle physics experiments, and the analysis of the data collected to perform physics measurements. This review starts from the current state of the art technology applied to detector simulation in high-energy physics and elaborates on the evolution of software tools developed to address the challenges posed by future accelerator programs beyond the HL-LHC era, into the 2030–2050 period. New accelerator, detector, and computing technologies set the stage for an exercise in how detector simulation will serve the needs of the high-energy physics programs of the mid 21st century, and its potential impact on other research domains.
Sandia National Laboratories develops and maintains several open-source software packages to support material accountancy analyses. This includes the Material Accountancy Performance Indicator Toolkit (MAPIT), the Fissile Facility Flow Modeler (F3M) and the Separation and Safeguards Performance Model Library (SSPM-L). MAPIT is responsible for performing statistical safeguards analyses on bulk and itemized data from nuclear fuel cycle facilities and can operate on real or synthetic data. MAPIT is the only open-source software for such analyses. F3M is a library of modules, built in MATLAB Simulink, that contain pre made blocks to represent different generic fuel cycle processes. These blocks can be used together in a modular fashion to represent and simulate nuclear fuel cycle processes with the goal of improving facility-level accountancy during the design phase. F3M is also an open-source library. Finally, the SSPM-L library is a series of completed models built from F3M. The library includes facility models such as a generic PUREX facility and a fuel fabrication facility. The SSPM-L library is not open source, but is available to collaborators with a relevant use case. These tools include modeling and simulation pipelines to simulate nuclear fuel cycle facilities and the underlying software needed to simulate measurement uncertainty and perform statistical analyses. Together, these tools can perform end-to-end nuclear material accountancy analyses. This report documents the various improvements made to these tools in FY25. Specifically, we added new statistical test, new statistical modeling capabilities, new fuel cycle facility models, and launched a new open-source model component library.
Purpose: Substantial target motion during the delivery of radiation therapy (RT) for pancreatic cancer is well recognized as a major limiting factor on RT effectiveness. The aim of this work is to monitor intra-fractional motion of the pancreas using ultrasound during RT delivery. Methods: Transabdominal Ultrasound B-mode images were collected from 5 volunteers using a research version of the Clarity Autoscan System (Elekta). The autoscan transducer with center frequency of 5 MHz was utilized for the scans. Imaging parameters were adjusted to acquire images at the desired depth with good contrast and a wide sweep angle. Since well-defined boundaries of the pancreas can be difficult to find on ultrasound B-mode images, the portal vein was selected as a surrogate for motion estimation of the head of the pancreas. The selection was due to its anatomical location posterior to the neck of the pancreas and close proximity to the pancreas head. The portal vein was contoured on the ultrasound images acquired during simulation using the Clarity Research AFC Workstation software. Volunteers were set up in a similar manner to the simulation for their monitoring session and the ultrasound transducer was mounted on an arm fixed to the couch. A video segment of the portal vein motion was captured. Results: The portal vein was visualized and segmented. Successful monitoring sessions of the portal vein were observed. In addition, our results showed that the ultrasound transducer itself reduces breathing related motion. This is analogous to the use of a compression plate to suppress respiration motion during thorax or abdominal irradiation. Conclusion: We demonstrate the feasibility of tracking the pancreas through the localization of the portal vein using abdominal ultrasound. This will allow for real-time tracking of the intra-fractional motion to justify PTV-margin and to account for unusual motions, thus, improving normal tissue sparing. This research was funding in part by Elekta Inc.
This report describes improvements made to the solid mechanics formulation in the MOOSE open source finite element simulation environment and the open source Nuclear Engineering Material model Library (NEML) for mechanical constitutive models. The focus of these improvements is to improve the usability and performance of simulations involving one or both pieces of software. Specifically, this work completes a new system for solid mechanics simulations in the MOOSE ecosystem providing exact linearizations and optimal (quadratic) convergence, for a variety of coordinate systems and material types, including large deformation simulations. This work then provides users a framework to build highly efficient mechanical simulations of structures or materials or to couple in additional MOOSE physics modules to build complex, scalable multiphysics simulations.
This is a MATLAB-based tool to post process simulation results in the EPRI software TRANSCARE, for massive cascading failure analysis following severe disturbances. There are a few key modules available in this tool, including: 1. automatically creating a contingency list to run TRANSCARE simulations, including substation outages above a certain kV threshold, N-k (1, 2 or 3) generator outages and branch outages; 2. read in and analyze a CKO file of PCG definition, an initiating event list, and a CDN file; 3. post process all the simulation results saved in a CDN file and perform critical event corridor analysis; 4. provide a summary of TRANSCARE simulations; 5. Identify the most frequently occurring event corridors in the system; and 6. Rank the contingencies using a user defined security index to quantify consequences in terms of total load loss, total number of cascades, etc.
The Particle/Ray Interaction Simulation Manager (PRISM) is a software tool developed by Triple Ring Technologies (TRT), http://www.tripleringtech.com/, under a Small Business Innovation Research Award from the Department of Homeland Security (DHS) to investigate the utility of simulations in X-ray-based explosive detection systems. The PRISM tool is intended to provide a standardized and intuitive user interface to control and configure realistic simulations of X-ray imaging systems. It includes a command-line interface that enables the user to configure simulated X-ray radiography and computed tomography (CT) systems. It allows the import of computer-aided design (CAD) models of objects and it supports Monte Carlo simulation engines such as GEANT4. PRISM is designed to simulate realistic X-ray projection images including scatter, noise, and various detector effects such as quantum efficiency, lag, and veiling glare. At the request of DHS, the authors assessed PRISM’s capabilities to simulate data generated by a Micro-CT test-bed system (MCT-TB) housed at the Lawrence Livermore National Laboratory (LLNL). This document describes the results of that assessment.
Purpose: The low contrast detectability (LCD) of a CT scanner is its ability to detect and display faint lesions. The current approach to quantify LCD is achieved using vendor-specific methods and phantoms, typically by subjectively observing the smallest size object at a contrast level above phantom background. However, this approach does not yield clinically applicable values for LCD. The current study proposes a statistical LCD metric using software tools to not only to assess scanner performance, but also to quantify the key factors affecting LCD. This approach was developed using uniform QC phantoms, and its applicability was then extended under simulated clinical conditions. Methods: MATLAB software was developed to compute LCD using a uniform image of a QC phantom. For a given virtual object size, the software randomly samples the image within a selected area, and uses statistical analysis based on Student’s t-distribution to compute the LCD as the minimal Hounsfield Unit’s that can be distinguished from the background at the 95% confidence level. Its validity was assessed by comparison with the behavior of a known QC phantom under various scan protocols and a tissue-mimicking phantom. The contributions of beam quality and scattered radiation upon the computed LCD were quantified by using various external beam-hardening filters and phantom lengths. Results: As expected, the LCD was inversely related to object size under all scan conditions. The type of image reconstruction kernel filter and tissue/organ type strongly influenced the background noise characteristics and therefore, the computed LCD for the associated image. Conclusion: The proposed metric and its associated software tools are vendor-independent and can be used to analyze any LCD scanner performance. Furthermore, the method employed can be used in conjunction with the relationships established in this study between LCD and tissue type to extend these concepts to patients’ clinical CT images.
Ontario Power Generation's (OPG) Darlington Nuclear Generating Station (DNGS) is currently executing a major refurbishment of the first of their four Candu reactors. As expected, the volumes of Low and Intermediate Level Waste (LILW) generated from this large-scale project have been significant. As a result, unique solutions were needed for the safe and efficient transportation of radioactive waste from the Darlington Nuclear Generating Station to OPG's Western Waste Management Facility for volume reduction and interim storage. OPG's Nuclear Waste Management Division and the Darlington Refurbishment Project worked together to identify unique radioactive waste transportation solutions. These solutions included staging and temporary storage facilities, new waste containers and specialized transportation packages, additional human resources and enhanced planning and logistics required to transport large volumes of routine and non-routine radioactive waste in a manner that minimized the impact on employee and public safety, project schedule and project costs. From a facilities perspective a shielded 'Marshalling Area' that allowed for loaded transportation packages to be staged was built on the Darlington site and a temporary licensed off-site radioactive waste storage facility close to the DNGS was contracted through a third party vendor. These two facilities allowed for more consistent and efficient schedules that limited radioactive waste shipments during off-hours and avoided peak traffic periods thereby limiting the risk to Refurbishment critical path activities. Special purpose-built transportation packages and waste containers were designed and fabricated to handle unique waste streams such as feeder pipes and other miscellaneous intermediate level waste items associated with reactor dismantling. Twenty-four Type A qualified transportation packages mounted to six flat deck trailers were designed and fabricated to support the shipment of contaminated feeder pipes. An Intermediate Level Waste Utility Container (ILWUC) to be used in conjunction with the Type B qualified Trillium Transportation Package was designed and fabricated to handle unique intermediate level waste streams arising from the reactor dismantling process. Logistics software was used to plan the shipments to ensure waste transportation packages could be loaded, transported, unloaded and returned safely to the project site in a timely manner without impacting project schedules. Several shipping scenarios were simulated and adjusted using the software to ensure that the most efficient use of equipment and resources for the radioactive waste transportation aspects of the project. (authors)
TRISQD software is developed for the computer simulation of processes in which radiation defects are formed under the corpuscular irradiation of semiconductor heterostructures with lenticular nanoinclusions of various shapes. The computer program is used to study defect-formation processes in p-i-n diodes with the i region having a built-in 20-period lattice of self-assembled Ge(Si) nanoislands formed under irradiation with high-energy neutrons. It is found that the fraction of Ge(Si) nanoislands in which point radiation defects are formed under the impact of atomic-displacement cascades is ≤3% of their total number in the lattice. More than 94% of the defects are localized in the bulk of the p, n, and i regions of the diode and in silicon layers that separate sheets of Ge(Si) nanoislands.
The appearance of single photon avalanche diodes (SPADs) in the field of PET detector modules made it necessary to apply more complex optical design methods to refine the performance of such assemblies. We developed a combined simulation tool that is capable to model complex detector structures including scintillation material, light guide, light collection optics and sensor, correctly taking into account the statistical behavior of emission of scintillation light and its absorbance in SPADs. As a validation we compared simulation results obtained by our software and another optical design program. Calculations were performed for a simple PET detector arrangement used for testing purposes. According to the results, deviation of center of gravity coordinates between the two simulations is 0.0195 mm, the average ratio of total counts 1.0052. We investigated the error resulting from finite sampling in wavelength space and we found that 20 nm pitch is sufficient for the simulation in case of the given spectral dependencies. (authors)
A Functional Mock-up Interface (FMI) defines a standardized interface to be used in computer simulations to develop complex cyber-physical systems. FMI implementation by a software modeling tool enables the creation of a simulation model that can be interconnected, or the creation of a software library called a Functional Mock-up Unit (FMU). This report describes an FMU wrapper implementation that imports FMUs into a C++ environment and uses an Euler solver that executes FMUs in parallel using Open Multi-Processing (OpenMP). The purpose of this report is to elucidate the runtime performance of the solver when a multi-component system is imported as a single FMU (for the whole system) or as multiple FMUs (for different groups of components as sub-systems). This performance comparison is conducted using two test cases: (1) a simple, multi-tank problem; and (2) a more realistic use case based on the Modelica Buildings Library. In both test cases, the performance gains are promising when each FMU consists of a large number of states and state events that are wrapped in a single FMU. Load balancing is demonstrated to be a critical factor in speeding up parallel execution of multiple FMUs.
We present an overview of a software system, Bellerophon, built to support a production-level HPC application called CHIMERA, which simulates core-collapse supernova events at the petascale. Developed over the last four years, Bellerophon enables CHIMERA's geographically dispersed team of collaborators to perform data analysis in near real-time. Its n-tier architecture provides an encapsulated, end-to-end software solution that enables the CHIMERA team to quickly and easily access highly customizable animated and static views of results from anywhere in the world via a web- deliverable, cross-platform desktop application. In addition, Bellerophon addresses software engineering tasks for the CHIMERA team by providing an automated mechanism for performing regression testing on a variety of supercomputing platforms. Elements of the team's workflow management needs are met with software tools that dynamically generate code repository statistics, access important online resources, and monitor the current status of several supercomputing resources.
We present an overview of a software system, Bellerophon, built to support a production-level HPC application called CHIMERA, which simulates core-collapse supernova events at the petascale. Developed over the last four years, Bellerophon enables CHIMERA s geographically dispersed team of collaborators to perform data analysis in near real-time. Its n-tier architecture provides an encapsulated, end-to-end software solution that enables the CHIMERA team to quickly and easily access highly customizable animated and static views of results from anywhere in the world via a web-deliverable, cross-platform desktop application. In addition, Bellerophon addresses software engineering tasks for the CHIMERA team by providing an automated mechanism for performing regression testing on a variety of supercomputing platforms. Elements of the team s workflow management needs are met with software tools that dynamically generate code repository statistics, access important online resources, and monitor the current status of several supercomputing resources.
The AirflowNetwork software provides the ability to simulate air movement and the resulting transport of contaminants and other materials within a building. The performance of air distribution systems, including supply and return leaks, and multizone airflows driven by outdoor wind and forced air during HVAC system operation are simulated using a nodal network representation of the building. The pressure network model implemented here is based on George Walton's AIRNET and was modified for inclusion in the EnergyPlus building energy modeling software. The algorithm were then extracted from EnergyPlus into a stand-alone package to facility use without the building energy model and more extensive testing and validation of the underlying operations.
This report documents the results of Task 3 of Project Action Sheet PP05 between the United States Department of Energy (DOE) and the Republic of Korea (ROK) Ministry of Education, Science, and Technology (MEST) for Support with Review of an ROK Risk Evaluation Process. This task was to have Sandia National Laboratories collaborate with the Korea Institute of Nuclear Nonproliferation and Control (KINAC) on several activities concerning how to determine the Probability of Neutralization, P N , and the Probability of System Effectiveness, P E , to include: providing descriptions on how combat simulations are used to determine P N and P E ; comparisons of the strengths and weaknesses of two neutralization models (the Neutralization.xls spreadsheet model versus the Brief Adversary Threat-Loss Estimator (BATLE) software); and demonstrating how computer simulations can be used to determine P N . Note that the computer simulation used for the demonstration was the Scenario Toolkit And Generation Environment (STAGE) simulation, which is a stand-alone synthetic tactical simulation sold by Presagis Canada Incorporated. The demonstration is provided in a separate Audio Video Interleave (.AVI) file.