Search NASA⌕ Search

SEARCH · Search NASA

Results for “Matlab data”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Simulated impedance of diffusion in porous media

This paper describes the use of a frequency domain, finite-difference scheme to simulate the impedance spectra of diffusion in porous microstructures. We investigate both open and closed systems for a range of ideal geometries, as well as some randomly generated synthetic volumes and tomographically derived microstructural data. In many cases, the spectra deviate significantly from the conventional Warburg-type elements typically used to represent diffusion in equivalent circuit analysis. Furthermore, a key finding is that certain microstructures show multiple peaks in the complex plane, which may be misinterpreted as separate electrochemical processes in real impedance data. This is relevant to battery electrode design as the techniques for nano-scale fabrication become more widespread. This simulation tool is provided as an open-source MatLab application and is freely available online as part of the TauFactor platform.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Pushing the Limits of High-speed X-ray Tomography to See the Unknown

First-of-their kind datasets from a high-speed X-ray tomography system were collected, and a novel numerical effort utilizing temporal information to reduce measurement uncertainty was shown. The experimental campaign used three high-speed X-ray imaging systems to collect data at 100 kHz of a scene containing high-velocity objects. The scene was a group of known objects propelled by a 12-gauge shotgun shell reaching speeds of hundreds of meters per second. These data represent a known volume where the individual components are known, with experimental uncertainties that can be used for reconstruction algorithm validation. The numerical effort used synthetic volumes in MATLAB to produce projections along known lines of sight to perform tomographic reconstructions. These projections and reconstructions were performed on a single object at two orientations, representing two timesteps, to increase the reconstruction accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2020

SorpPropLib : An Open-Source Database for Sorption Equilibrium Properties

Hundreds of sorption working pair equilibrium measurements have been reported in the literature, and many experimental results have been fitted to equilibrium equations. However, an open database does not exist for the equilibrium equations. Therefore, this work compiled 438 published equilibrium equations of sorption working pairs into an open-source database called SorpPropLib. SorpPropLib is accessible from six programming environments commonly used in sorption cycle simulation (C++, Excel®, Python®, Matlab®, Modelica® and Lab-VIEW™). Plotting capabilities allow fast cycle visualization in each programming environment. Thereby, the open-source database SorpPropLib aims to enable and encourage easy data sharing to advance sorption research.

Engelpracht, Mirko↗
DOE OSTI2016

Fission gas bubble identification using MATLAB's image processing toolbox

Automated image processing routines have the potential to aid in the fuel performance evaluation process by eliminating bias in human judgment that may vary from person-to-person or sample-to-sample. This study presents several MATLAB based image analysis routines designed for fission gas void identification in post-irradiation examination of uranium molybdenum (U–Mo) monolithic-type plate fuels. Frequency domain filtration, enlisted as a pre-processing technique, can eliminate artifacts from the image without compromising the critical features of interest. This process is coupled with a bilateral filter, an edge-preserving noise removal technique aimed at preparing the image for optimal segmentation. Adaptive thresholding proved to be the most consistent gray-level feature segmentation technique for U–Mo fuel microstructures. The Sauvola adaptive threshold technique segments the image based on histogram weighting factors in stable contrast regions and local statistics in variable contrast regions. Once all processing is complete, the algorithm outputs the total fission gas void count, the mean void size, and the average porosity. The final results demonstrate an ability to extract fission gas void morphological data faster, more consistently, and at least as accurately as manual segmentation methods. - Highlights: •Automated image processing can aid in the fuel qualification process. •Routines are developed to characterize fission gas bubbles in irradiated U–Mo fuel. •Frequency domain filtration effectively eliminates FIB curtaining artifacts. •Adaptive thresholding proved to be the most accurate segmentation method. •The techniques established are ready to be applied to large scale data extraction testing.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The Polar System Analysis Package (Ver. 1 Specifications)

The Polar System Analysis Package (POLeSTAr) is a modern data fusion tool being designed to benchmark,evaluate and compare the polar physics and biogeochemistry of Earth System Models (ESMs). It is being written in Python to make use of libraries including Xarray, Matplotlib, SciPy, NumPy and PyTorch within algorithms developed specifically for polar system modeling. Initial POLeSTAr algorithms will be adapted from the MATLAB Ridgepack package that was designed to analyze the Model for Prediction Across Scales (MPAS) sea ice component of the Energy Exascale Earth System Model (E3SM) and sea ice model CICE within the Community Earth System Model (CESM).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

A Probability-based Model for Cost-effective Integration of Renewables into the Electricity Grid (Final Technical Report)

Funded by the U.S. Department of Energy during 2017-2019, the resource adequacy evaluation and production costing tool, called MATPLAN (MATlab and probability-based PLANning), was developed. MATPLAN can take into account the variable nature of renewable energy sources (both solar PV and wind farms), rather than considering deterministic data as in the current practice. It combines both probabilistic and optimization techniques to allow the determination of optimal system expansion policy and perform effective load carrying capacity (ELCC) analysis of a system with high renewable penetration. This is in contrast to the current practice that treats these variable sources as negative loads. This tool expects to benefit system planners by enabling them to consider renewable energy sources as options for their expansion plans. System operators can use MATPLAN to calculate the capacity credit for each solar/wind farm using probabilistic forecasts, thus enabling quantification of the operating reserve requirements to cope with uncertainty and variability of solar/wind output committed in the market.

14 SOLAR ENERGY↗
DOE OSTI2017

An Introduction to Goodness of Fit for PMU Parameter Estimation

New results of measurements of phasor-like signals are presented based on our previous work on the topic. In this document an improved estimation method is described. The algorithm (which is realized in MATLAB software) is discussed. We examine the effect of noisy and distorted signals on the Goodness of Fit metric. The estimation method is shown to be performing very well with clean data and with a measurement window as short as a half a cycle and as few as 5 samples per cycle. The Goodness of Fit decreases predictably with added phase noise, and seems to be acceptable even with visible distortion in the signal. While the exact results we obtain are specific to our method of estimation, the Goodness of Fit method could be implemented in any phasor measurement unit.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

SU-F-R-55: Early Detection of Treatment Induced Bone Marrow Injury During Chemoradiation Therapy Using Quantitative CT

Purpose: Acute hematologic toxicity associated with bone marrow injury is a common complication of chemoradiation therapy (CRT) for pelvic malignancies. In this work, we investigate the feasibility of using quantitative CT to detect bone marrow injury during CRT. Methods: Daily CTs were acquired during routine CT-guided radiation therapy using a CT-on-rails for 15 cervical cancer patients. All patients treated with a radiation dose of 45.0 to 50.4 Gy in 1.8 Gy/fraction along with chemotherapy. For each patient, the contours of bone marrow were generated in L4, L5 and sacrum on the first daily CT and then populated to other daily CTs by rigid registration using MIM (MIM Software Inc., Cleveland, OH) with manual editing if possible. A series of CT texture parameters, including Hunsfield Unit (HU) histogram, mean HU, entropy, energy, in bone marrow contours were calculated using MATLAB on each daily CT and were correlated with the completed blood counts (CBC) collected weekly for each patient. The correlations were analyzed with Pearson correlation tests. Results: For all patient data analyzed, mean HU in bone marrow decreased during CRT delivery. From the first to the last fraction the average mean HU reduction is 58.1 ± 13.6 HU (P<0.01). This decrease can be observed as early as after first 5 fractions and is strongly associated with the changes of most CBC quantities, such as the reductions of white and blood cell counts (r=0.97, P=0.001). The reduction of HU is spatially varied. Conclusion: Chemoradiation induced bone marrow injury can be detected during the delivery of CRT using quantitative CT. Chemoradiation results in reductions in mean HU, which are strongly associated with the change in the pretrial blood cell counts. Early detection of bone marrow injury with commonly available CT opens a door to improve bone marrow sparing, reducing risk of hematologic toxicity.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TU-D-209-03: Alignment of the Patient Graphic Model Using Fluoroscopic Images for Skin Dose Mapping

Purpose: The Dose Tracking System (DTS) was developed to provide realtime feedback of skin dose and dose rate during interventional fluoroscopic procedures. A color map on a 3D graphic of the patient represents the cumulative dose distribution on the skin. Automated image correlation algorithms are described which use the fluoroscopic procedure images to align and scale the patient graphic for more accurate dose mapping. Methods: Currently, the DTS employs manual patient graphic selection and alignment. To improve the accuracy of dose mapping and automate the software, various methods are explored to extract information about the beam location and patient morphology from the procedure images. To match patient anatomy with a reference projection image, preprocessing is first used, including edge enhancement, edge detection, and contour detection. Template matching algorithms from OpenCV are then employed to find the location of the beam. Once a match is found, the reference graphic is scaled and rotated to fit the patient, using image registration correlation functions in Matlab. The algorithm runs correlation functions for all points and maps all correlation confidences to a surface map. The highest point of correlation is used for alignment and scaling. The transformation data is saved for later model scaling. Results: Anatomic recognition is used to find matching features between model and image and image registration correlation provides for alignment and scaling at any rotation angle with less than onesecond runtime, and at noise levels in excess of 150% of those found in normal procedures. Conclusion: The algorithm provides the necessary scaling and alignment tools to improve the accuracy of dose distribution mapping on the patient graphic with the DTS. Partial support from NIH Grant R01-EB002873 and Toshiba Medical Systems Corp.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

JGCRI/hydro-emulator

The hydrological emulator (HE) was built to mimic complex global hydrological models (GHMs) at a range of spatial scales. The HE was written in MATLAB and utilizes a genetic algorithm to calibrate the a, b, c, d, m parameters from the ABCD runoff algorithm to observed runoff from GHMs when forced by the same climate data. The HE employs two methods to choose from when processing basin data: lumped (by basin) and distributed (gridded). Both methods can be evaluated using the built-in Kling-Gupta efficiency as a measure of model performance when compared to the outputs from the more complex GHM. Testing shows that the HE is 7 orders of magnitude faster then the widely used Variable Infiltration Capacity (VIC) GHM.

Vernon, Chris↗
DOE OSTI2024

2D Polyhedral Template Matching for Atomic Resolution Microscopy

SAND2024-13879O The 2D Polyhedral Template Matching for Atomic Resolution Microscopy is a suite of functions for analyzing atomic resolution electron microscopy images using the 2D polyhedral template matching (2D-PTM) method. This software analyzes atomic resolution microscopy data from electron microscopic imaging or scanning probe microscopies. The primary application is for identifying different crystal phases, crystal orientations, and defect structures obtained in such atomic resolution images. Written in MATLAB and starting from an atomic resolution image, the code identifies the positions of atomic intensity peaks. It then matches predefined structural templates to the local atomic environments. Outputs include the local structural identification, the template scaling factor and rotation angle, root-mean-squared deviation (RMSD), and centrosymmetry parameter.

Medlin, Douglas↗
DOE OSTI2022

HERO WEC V1.0 - WEC-Sim Model

This zip file contains the files that are needed to simulate NREL's HERO WEC (hydraulic and electric reverse osmosis wave energy converter). This requires the user to have already installed WEC-Sim. In addition to the standard toolboxes that are required to run WEC-Sim the user will also need the Simscape Fluids and Simscape Driveline packages. In the zip file you will find the following: - HEROV1_HPTO.slx: Simulink-based WEC Sim model of the first gen (V1.0) Hydraulic PTO (power take-off) that was designed for the HERO WEC - wecSimInputFile.m: Input file needed to run the model - userDefinedFunctionsMCR.m: MCR (multi condition run) script that is needed if a use wants to simulate multiple wave conditions. - geometry (folder): Includes the geometry file that is needed for visualization - hydroData (folder): Includes the required WAMIT data to run WEC-Sim

16 TIDAL AND WAVE POWER↗
DOE OSTI2016

iMatTOUGH: An open-source Matlab-based graphical user interface for pre- and post-processing of TOUGH2 and iTOUGH2 models

TOUGH2 and iTOUGH2 are powerful models that simulate the heat and fluid flows in porous and fracture media, and perform parameter estimation, sensitivity analysis and uncertainty propagation analysis. However, setting up the input files is not only tedious, but error prone, and processing output files is time consuming. Here, we present an open source Matlab-based tool (iMatTOUGH) that supports the generation of all necessary inputs for both TOUGH2 and iTOUGH2 and visualize their outputs. The tool links the inputs of TOUGH2 and iTOUGH2, making sure the two input files are consistent. It supports the generation of rectangular computational mesh, i.e., it automatically generates the elements and connections as well as their properties as required by TOUGH2. The tool also allows the specification of initial and time-dependent boundary conditions for better subsurface heat and water flow simulations. The effectiveness of the tool is illustrated by an example that uses TOUGH2 and iTOUGH2 to estimate soil hydrological and thermal properties from soil temperature data and simulate the heat and water flows at the Rifle site in Colorado.

58 GEOSCIENCES↗
DOE OSTI2016

SU-F-R-18: Updates to the Computational Environment for Radiological Research for Image Analysis

Purpose: To present new tools in CERR for Texture Analysis and Visualization. Method: (1) Quantitative Image Analysis: We added the ability to compute Haralick texture features based on local neighbourhood. The Texture features depend on many parameters used in their derivation. For example: (a) directionality, (b) quantization of image, (c) patch-size for the neighborhood, (d) handling of the edge voxels within the region of interest, (e) Averaging co-occurance matrix vs texture features for different directions etc. A graphical user interface was built to set these parameters and then visualize their impact on the resulting texture maps. The entire functionality was written in Matlab. Array indexing was used to speed up the texture calculation. The computation speed is very competitive with the ITK library. Moreover, our implementation works with multiple CPUs and the computation time can be further reduced by using multiple processor threads. In order to reduce the Haralick texture maps into scalar features, we propose the use of Texture Volume Histograms. This lets users make use of the entire distribution of texture values within the region of interest rather than using just the mean and the standard deviations. (2) Qualitative/Visualization tools: The derived texture maps are stored as a new scan (derived) within CERR’s planC data structure. A display that compares various scans was built to show the raw image and the derived texture maps side-by-side. These images are positionally linked and can be navigated together. CERR’s graphics handling was updated and sped-up to be compatible with the newer Matlab versions. As a result, the users can use (a) different window levels and colormaps for different viewports, (b) click-and-drag or use mouse scroll-wheel to navigate slices. Results: The new features and updates are available via https://www.github.com/adityaapte/cerr . Conclusion: Features added to CERR increase its utility in Radiomics and Outcomes modeling.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Dynamic Modeling for the Separation of Rare Earth Elements Using Solvent Extraction: Predicting Separation Performance Using Laboratory Equilibrium Data

Industrial rare earth element (REE) separation facilities utilize acidic cation exchange ligands such as 2-ethylhexylphosphonic acid mono-2-ethylhexyl ester (PC88A) for solvent extraction processes. REE separations are costly and difficult due to their chemical similarities and subsequent low separation factors. Several empirical correlations are available in the literature to predict steady state extraction equilibria for various solvent systems. However, complete solvent extraction flowsheet design for REE separations requires complex scrubbing and stripping circuits to separate and produce individual pure species. Furthermore, the sensitive process chemistry is susceptible to process upsets and may manifest as irreversible transients that complicate startup and continuous operation of a solvent extraction plant. Thus, dynamic modeling of REE separations will aid in process design, optimization, and management of process fluctuations. A time-dependent MATLAB/SIMULINK modeling framework has been developed for REE solvent extraction processes using laboratory equilibrium data. The model was used to develop a flowsheet that produced high purity neodymium from a 25 wt% praseodymium and 75 wt% neodymium feed. Laboratory mixer-settlers were used to verify and validate model performance. Results indicated that the model reasonably predicts the dynamic behavior of a counter-current REE separation process, and accurately predicts the steady state REE concentration profiles across the cascade. The modeling framework can be applied to any REE separation or solvent system provided adequate laboratory equilibrium data are available.

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

SU-D-BRD-01: An Automated Physics Weekly Chart Checking System Supporting ARIA

Purpose: A software tool was developed in this study to perform automatic weekly physics chart check on the patient data in ARIA. The tool accesses the electronic patient data directly from ARIA server and checks the accuracy of treatment deliveries, and generates reports which summarize the delivery history and highlight the errors. Methods: The tool has four modules. 1) The database interface is designed to directly access treatment delivery data from the ARIA database before reorganizing the data into the patient chart tree (PCT). 2) PCT is a core data structure designed to store and organize the data in logical hierarchies, and to be passed among functions. 3) The treatment data check module analyzes the organized data in PCT and stores the checking results into PCT. 4) Report generation module generates reports containing the treatment delivery summary, chart checking results and plots of daily treatment setup parameters (couch table positions, shifts of image guidance). The errors that are found by the tool are highlighted with colors. Results: The weekly check tool has been implemented in MATLAB and clinically tested at two major cancer centers. Javascript, cascading style sheets (CSS) and dynamic HTML were employed to create the user-interactive reports. It takes 0.06 second to search the delivery records of one beam with PCT and compare the delivery records with beam plan. The reports, saved in the HTML files on shared network folder, can be accessed by web browser on computers and mobile devices. Conclusion: The presented weekly check tool is useful to check the electronic patient treatment data in Varian ARIA system. It could be more efficient and reliable than the manually check by physicists. The work was partially supported by a research grant from Varian Medical System.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

The Binary Alkali Nitrate and Chloride Phase Diagrams: NaNO{sub 3}-KNO{sub 3}, LiNO{sub 3}-NaNO{sub 3}, LiNO{sub 3}-KNO{sub 3}, and NaCl-KCl

LiNO{sub 3}-NaNO{sub 3}, LiNO{sub 3}-KNO{sub 3,} NaNO{sub 3}-KNO{sub 3}, and NaCl-KCl phase relations were formulated using the regular solution model to obtain desired compositions for their potential application as a heat transfer fluid. The simultaneous nonlinear equations formed were solved numerically using a contour plot technique in Matlab. An algorithm was developed to determine the real roots of the coupled nonlinear equations and subsequent phase diagrams were constructed. The phase diagrams were compared with those predicted using FactSage thermochemical package and literature data from experimental studies. The calculated NaNO{sub 3}-KNO{sub 3} phase diagram shows that this system has a eutectic point at 0.5 mole fraction of KNO{sub 3} with a eutectic temperature of 220.85 °C, while the LiNO{sub 3}-KNO{sub 3} system was found to have a eutectic point at 0.44 mole fraction LiNO{sub 3} at 137 °C. The LiNO{sub 3}-NaNO{sub 3} system shows a eutectic at 0.524 mole fraction LiNO{sub 3} at 192.8 °C, and finally the NaCl-KCl system has a eutectic point at 0.48 mole fraction NaCl at 634 °C. The liquidus and solidus curves obtained fitted close against values obtained from the quasi-chemical model used in FactSage and experimental values from the literature. This study showed that the regular solution model combined with the contour plots approach can sufficiently describe the liquidus–solidus curves of the alkali binary nitrate and chloride systems and could be a useful method for predicting the phase diagrams for higher order alkali nitrates and chlorides.

36 MATERIALS SCIENCE↗