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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

SISSO (Sure-independence-screening sparsifying-operator regressor)

Implementation of the SISSO regression algorithm in MATLAB. The SISSO regression algorithm iteratively selects model features from candidates, converging even when the number of possible features is much greater than the number of available data points. Includes a test script that validates the algorithm by replicating the results using data and procedure from https://analytics-toolkit.nomad-coe.eu/hub/user-redirect/notebooks/tutorials/compressed_sensing.ipynb. Python code from the SISSO regressor in 'sisso.py', from the link above, was used as the basis for developing the MATLAB implementation. The SISSO regression algorithm is detailed by the original authors in R. Ouyang, S. Curtarolo, E. Ahmetcik et al., Phys. Rev. Mater. 2, 083802 (2018), R. Ouyang, E. Ahmetcik, C. Carbogno, M. Scheffler, and L. M. Ghiringhelli, J. Phys.: Mater. 2, 024002 (2019).

Gasper, Paul↗
DOE OSTI2019

Calculation of the cross and auto power spectral densities for low neutron counting from pulse mode detectors

Here, this report introduces a methodology to calculate the prompt neutron decay constant of a subcritical assembly from the cross and auto power spectral densities of detectors operating in pulse mode when the neutron counting is low. The pulse operation mode allows the easy discrimination of the gamma contribution to the detector signal by the pulse height or shape techniques. The methodology is based on a MATLAB script that allows to change arbitrary parameters at any time after the experiment is performed. This technique requires list mode data and is much more flexible than programming the digital acquisition system because in the latter case the arbitrary parameters are set prior to the experiment and cannot be changed any longer after the experiment is performed. The MATLAB script processes all the detectors timestamps from several thousand pulses of the particle accelerator and redistributes them into a time interval that spreads over only 10 particle accelerator pulses. This pulses collapsing procedure has been performed because of the low neutron counting.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

SU-E-T-385: Evaluation of DVH Change for PTV Due to Patient Weight Loss in Prostate VMAT Using Gaussian Error Function

Purpose: We proposed a method to predict the change of dose-volume histogram (DVH) for PTV due to patient weight loss in prostate volumetric modulated arc therapy (VMAT). This method is based on a pre-calculated patient dataset and DVH curve fitting using the Gaussian error function (GEF). Methods: Pre-calculated dose-volume data from patients having weight loss in prostate VMAT was employed to predict the change of PTV coverage due to reduced depth in external contour. The effect of patient weight loss in treatment was described by a prostate dose-volume factor (PDVF), which was evaluated by the prostate PTV. Along with the PDVF, the GEF was used to fit into the DVH curve for the PTV. To predict a new DVH due to weight loss, parameters from the GEF describing the shape of DVH curve were determined. Since the parameters were related to the PDVF as per the specific reduced depth, we could first predict the PDVF at a reduced depth based on the prostate size from the pre-calculated dataset. Then parameters of the GEF could be determined from the PDVF to plot the new DVH for the PTV corresponding to the reduced depth. Results: A MATLAB program was built basing on the patient dataset with different prostate sizes. We input data of the prostate size and reduced depth of the patient into the program. The program then calculated the PDVF and DVH for the PTV considering the patient weight loss. The program was verified by different patient cases with various reduced depths. Conclusion: Our method can estimate the change of DVH for the PTV due to patient weight loss quickly without CT rescan and replan. This would help the radiation staff to predict the change of PTV coverage, when patient’s external contour reduced in prostate VMAT.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Comparison of plug flow and multi-node stratified tank modeling approaches regarding computational efficiency and accuracy

Residential water heaters contain water stratified by temperature-driven density differences. This implies that a water tank can reach a state in which the top and bottom sections have different temperatures, unless mixing happens. A high degree of thermal stratification can improve the efficiency of some water heaters, by saving the amount of energy required for the heat-up process. Studies of stratification became popular in the 1970s and it remains an active research topic today. The research has led to the development of different models and techniques to better predict and define a stratified tanks behavior. By comparing these models and techniques used previously to describe thermal stratification, the phenomenon could be better understood, exploited, and used to increase efficiency and thermal energy capacity in modern water tanks. From the existing models, we found the one-dimensional standard plug-flow and a multi node model to be appropriate for analyzing the processes of the heat up and cool-down in a water tank. These two models are based on energy balances. This work involved comparing the accuracy and computational effort needed to implement these models. To assess accuracy, we compared both types of existing models to experimental data (also collected in this work) which included a heat up process using an external heat pump. This external process included a layering process that has an eddy diffusivity at five times the rate of thermal diffusion. For this project, we implemented the models in MATLAB, the multi-paradigm numerical computing environment. We quantified model accuracy using the root mean squared error between modeled data and experimental data for six measured tank temperatures. Comparing the accuracy and the computational time taken to run the simulation provides a method to contrast the performance of each model and a way to rate it. The multi node model was run using from 6 to 96 spatial nodes; the plug flow model was run using 1 to 0.001 º C temperature bin sizes. Additionally, timesteps were varied from 4 to 236 s. The results quantify the tradeoff between accuracy and computational time, providing guidance for simulations to intelligently select the best model type and simulation parameters. This research can be used to validate the pre-existing models and possibly improve the modern water tank.

Bulnes, Fernando Karg↗
DOE OSTI2019

Assessment of Calibration Approaches for the Stress Relaxation Test

Advanced manufacturing technologies are enabling the next generation of superalloys for extreme temperatures, pressures, and longevity. There is a need for Accelerated Creep Test (ACT) methods to reduce the time needed to implement new creep resistant superalloys. Without the development of ACTs, the qualification of new superalloys using conventional creep tests can take 11+ years corresponding to the 100,000-hour service life of long-lived turbomachinery components. A refined stress relaxation test (SRT) is introduced to succeed real-time conventional creep tests (CCTs). The stress relaxation test generates extensive stress versus creep-strain-rate data within a single specimen. Since the specimen is held below the elastic limit and relatively little creep deformation has occurred, several SRTs at different temperatures can be performed using a single specimen. To obtain data, specimens are rapidly loaded in force-control to just below the elastic limit. The specimen is switched to displacement-control and held at a fixed displacement while stress relaxation occurs. This generates a stress versus time curve that must be converted into stress versus creep-strain-rate. The stress versus creep-strain-rate curves are utilized to calculate creep activation energy and produce a creep-strain-rate master curve. The objective of this study is to determine the best calibration approach for SRT: regression of stress, regression of creep strain, and the finite difference of creep strain. It is recommended that the differentiation be applied; however, finite difference applied to the raw data produces inconsistent calculations of the creep-strain-rate. In this study, a MATLAB algorithm is written to evaluate all three calibration options where, stress and creep deformation are regressed into natural logarithm functions. The resulting creep-strain-rate curves are compared to determine the most reliable and consistent method. Analyzing and evaluating the results the regression of stress and regression or creep strain calibration methods prove to provide accurate results. These methods will help replace CCTs by conducting less expensive ACTs, while increasing the overall efficiency of the development of new materials.

36 MATERIALS SCIENCE↗
DOE OSTI2016

TU-FG-201-07: Development of SRS Conical Collimator Collision Prediction Software for Radiation Treatment Safety

Purpose: Stereotactic radiosurgery (SRS) treatments using conical collimators can potentially result in gantry collision with treatment table due to limited collision-clear spaces. An in-house software was developed to help the SRS treatment planner mitigate potential SRS conical collimator (Varian Medical System, Palo Alto, CA) collisions with the treatment table. This software was designed to remove treatment re-planning secondary to unexpected collisions. Methods: A BrainLAB SRS ICT Frameless Extension used for SRS treatments in our clinic was mathematically modelled using surface points registered to the 3D co-ordinate space of the couch extension. The surface points are transformed based on the treatment isocenter point and potential collisions are determined in 3D space for couch and gantry angle combinations. The distance between the SRS conical collimators and LINAC isocenter is known. The collision detection model was programmed in MATLAB (Mathwork, Natick, MA) to display graphical plots of the calculations, and the plotted data is used to avoid the gantry and couch angle combinations that would likely result in a collision. We have utilized the cone collision tool for 23 SRS cone treatment plans (8 retrospective and 15 prospective for 10 patients). Results: Twenty one plans strongly agreed with the software tool prediction for collision. However, in two plans, a collision was observed with a 0.5 cm margin when the software predicted no collision. Therefore, additional margins were added to the clearance criteria in the program to achieve a lower risk of actual collisions. Conclusion: Our in-house developed collision check software successfully avoided SRS cone re-planning by 91.3% due to a reduction in cone collisions with the treatment table. Future developments to our software will include a CT image data set based collision prediction model as well as a beam angle optimization tool to avoid normal critical tissues as well as previously treated lesions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Fault detection in PV systems

A practical fault detection approach for PV systems intended for online implementation is developed. The fault detection model here is built using artificial neural network. initially the photovoltaic system is simulated using MATLAB software and output power is collected for various combinations of irradiance and temperature. Data is first collected for normal operating condition and then four types of faults are simulated and data are collected for faulty conditions. Four faults are considered here and they are: Line to Line faults with a small voltage difference, Line to line faults with a large voltage difference, degradation fault and open-circuit fault. This data is then used to train the neural network and to develop the fault detection model. The fault detection model takes irradiance, temperature and power as the input and accurately gives the type of fault in the PV system as the output. This system is a generalized one as any PV module datasheet can be used to simulate the Photovoltaic system and also this fault detection system can be implemented online with the use of data acquisition system.

14 SOLAR ENERGY↗
DOE OSTI2018

DASH Slow Strain Rates from Brady Hot Springs Geothermal Field during PoroTomo Deployment Period

This submission contains slow strain rates summed to radians over 30 second intervals [rad/s] derived from horizontal distributed acoustic sensing measurements (DASH) of Brady geothermal field during PoroTomo deployment (2016-Mar-14 to 2016-Mar-26). There is one file corresponding to each day written in *.mat format for use with Matlab. The format for the binary Matlab .mat files are defined at: https://www.mathworks.com/help/pdf_doc/matlab/matfile_format.pdf. One such file includes the following variables: 'flist': list of raw DASH files used in the summation 'time_tag_mdt': sample time tag in datetime format with hours given in 24-hr format (yyyy/MM/dd HH:mm:ss.SSSSSSS) 'time_tag_uts': sample time tag in Unix time 'strain_rate_summed_over30s_in_radians_per_second': slow strain rates summed over 30 second intervals in units rad/s 'sample_standard_deviation_in_radians_per_second': corresponding sample standard deviation of slow strain rates in units rad/s The PoroTomo final technical report, raw DASH data, and software repository are also available through the links below.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Parameter Sensitivity Analysis of the SparTen High Performance Sparse Tensor Decomposition Software (Extended Analysis)

Tensor decomposition models play an increasingly important role in modern data science applications. One problem of particular interest is fitting a low-rank Canonical Polyadic (CP) tensor decomposition model when the tensor has sparse structure and the tensor elements are nonnegative count data. SparTen is a high-performance C++ library which computes a low-rank decomposition using different solvers: a first-order quasi-Newton or a second-order damped Newton method, along with the appropriate choice of runtime parameters. Since default parameters in SparTen are tuned to experimental results in prior published work on a single real-world dataset conducted using MATLAB implementations of these methods, it remains unclear if the parameter defaults in SparTen are appropriate for general tensor data. Furthermore, it is unknown how sensitive algorithm convergence is to changes in the input parameter values. This report addresses these unresolved issues with large-scale experimentation on three benchmark tensor data sets. Experiments were conducted on several different CPU architectures and replicated with many initial states to establish generalized profiles of algorithm convergence behavior.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Analytic reconstruction algorithms for triple-source CT with horizontal data truncation

Purpose: This paper explores a triple-source imaging method with horizontal data truncation to enlarge the field of view (FOV) for big objects. Methods: The study is conducted by using theoretical analysis, mathematical deduction, and numerical simulations. The proposed algorithms are implemented in c + + and MATLAB. While the basic platform is constructed in MATLAB, the computationally intensive segments are coded in c + +, which are linked via a MEX interface. Results: A triple-source circular scanning configuration with horizontal data truncation is developed, where three pairs of x-ray sources and detectors are unevenly distributed on the same circle to cover the whole imaging object. For this triple-source configuration, a fan-beam filtered backprojection-type algorithm is derived for truncated full-scan projections without data rebinning. The algorithm is also extended for horizontally truncated half-scan projections and cone-beam projections in a Feldkamp-type framework. Using their method, the FOV is enlarged twofold to threefold to scan bigger objects with high speed and quality. The numerical simulation results confirm the correctness and effectiveness of the developed algorithms. Conclusions: The triple-source scanning configuration with horizontal data truncation cannot only keep most of the advantages of a traditional multisource system but also cover a larger FOV for big imaging objects. In addition, because the filtering is shift-invariant, the proposed algorithms are very fast and easily parallelized on graphic processing units.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Mass separation of deuterium and helium with conventional quadrupole mass spectrometer by using varied ionization energy

Deuterium pressure in deuterium-helium mixture gas is successfully measured by a common quadrupole mass spectrometer (model: RGA200) with a resolution of ∼0.5 atomic mass unit (AMU), by using varied ionization energy together with new developed software and dedicated calibration for RGA200. The new software is developed by using MATLAB with the new functions: electron energy (EE) scanning, deuterium partial pressure measurement, and automatic data saving. RGA200 with new software is calibrated in pure deuterium and pure helium 1.0 × 10{sup −6}–5.0 × 10{sup −2} Pa, and the relation between pressure and ion current of AMU4 under EE = 25 eV and EE = 70 eV is obtained. From the calibration result and RGA200 scanning with varied ionization energy in deuterium and helium mixture gas, both deuterium partial pressures (P{sub D{sub 2}}) and helium partial pressure (P{sub He}) could be obtained. The result shows that deuterium partial pressure could be measured if P{sub D{sub 2}} > 10{sup −6} Pa (limited by ultimate pressure of calibration vessel), and helium pressure could be measured only if P{sub He}/P{sub D{sub 2}} > 0.45, and the measurement error is evaluated as 15%. This method is successfully employed in EAST 2015 summer campaign to monitor deuterium outgassing/desorption during helium discharge cleaning.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Unsupervised Machine Learning for Evaluation of Aging in Explosive Pressed Pellets

This project is an evaluation of the changes in the dimensional geometry of pressed explosive powder pellets when stored for prolonged periods under unmonitored storage conditions. Through the application of Structural Health Monitoring (SHM) techniques, questions to be answered by this investigation are what, if any, effect has the long-term storage of the pellets had on their physical attributes and, if changes have occurred, are the resulting differences of significant magnitude to affect the performance of the pellets. This resulting report includes a summary of the relevant background information, explanation of the SHM approach used to analyze the available data, presentation and discussion of the results, concluding remarks, and an accompanying appendix containing the MATLAB program developed for the project.

36 MATERIALS SCIENCE↗
DOE OSTI2016

SU-F-T-266: Dynalogs Based Evaluation of Different Dose Rate IMRT Using DVH and Gamma Index

Purpose: This work investigates the impact of low and high dose rate on IMRT through Dynalogs by evaluating Gamma Index and Dose Volume Histogram. Methods: The Eclipse™ treatment planning software was used to generate plans on prostate and head and neck sites. A range of dose rates 300 MU/min and 600 MU/min were applied to each plan in order to investigate their effect on the beam ON time, efficiency and accuracy. Each plan had distinct monitor units per fraction, delivery time, mean dose rate and leaf speed. The DVH data was used in the assessment of the conformity and plan quality.The treatments were delivered on Varian™ Clinac 2100C accelerator equipped with 120 leaf millennium MLC. Dynalogs of each plan were analyzed by MATLAB™ program. Fluence measurements were performed using the Sun Nuclear™ 2D diode array and results were assessed, based on Gamma analysis of dose fluence maps, beam delivery statistics and Dynalogs data. Results: Minor differences found by adjusted R-squared analysis of DVH’s for all the plans with different dose rates. It has been also found that more and larger fields have greater time reduction at high dose rate and there was a sharp decrease in number of control points observed in dynalog files by switching dose rate from 300 MU/min to 600 MU/min. Gamma Analysis of all plans passes the confidence limit of ≥95% with greater number of passing points in 300 MU/min dose rate plans. Conclusion: The dynalog files are compatible tool for software based IMRT QA. It can work perfectly parallel to measurement based QA setup and stand-by procedure for pre and post delivery of treatment plan.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Downscaled precipitation and mean air temperature datasets; East-Taylor subbasin; 2008-2019; daily temporal resolution; 400 m spatial resolution

This dataset provides gridded meteorological forcing data (specifically, daily precipitation and daily mean air temperature). The dataset has been generated by downscaling the Parameter-elevation Regressions on Independent Slopes Model (PRISM) dataset from a spatial resolution of 800 m to 400 m. The time period is 2008-2019, and the mapped area is East Taylor subbasin in Upper Colorado. The data are in the form of NetCDF files, arranged by year. Compared to the PRISM dataset, the dates of the downscaled dataset are shifted backwards by one (e.g., downscaled data for May 25 corresponds to PRISM data for May 26). This temporal shifting makes the meteorological forcing correspond more closely to the prescribed date. The NetCDF format is a standard raster format that can be read using any Geographic Information System (GIS) software; plenty of modules exist in popular scripting languages (such as Python, R and Matlab) that can also be used to read NetCDF files.East_Taylor_Tavg_PRISM400 corresponds to mean air temperature.East_Taylor_Precip_PRISM400.zip corresponds to precipitation.The proprietary PRISM data (800 m resolution) were purchased with funding from the Watershed Function Scientific Focus Area supported by U.S. Department of Energy, Office of Science, Office of Biological and Environmental Research under award no. DE-AC02-05CH11231.Dataset update on March 9, 2022:Uploaded new data files that are identical to the earlier files but are now in the NetCDF format, arranged by year. Earlier, the files were in the GeoTiff format, arranged by date.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

2017 Western Passage Tidal Energy Resource Characterization Measurements

These data are from tidal resource characterization measurements collected between April and July 2017 in Western Passage near Eastport, Maine, USA. The dataset contains the following four sub-datasets, each of which is described in more detail in the README.pdf. 1. A bottom-mounted Teledyne RDI Workhorse 600 kHz acoustic Doppler current profiler (ADCP) was deployed at 44.92015 N, 66.98915 W in ~50 m of water from 3 April to 18 July (106 days). Data were recorded in 6-minute increments in the ENU (East, magnetic North, Up) coordinate system with bin-mapping enabled. 2. A bottom-mounted Nortek Signature 500 kHz ADCP was deployed at 44.92192 N, 66.98913 W in ~50 m of water from 4 April to 18 July (105 days). Data were sampled and recorded at 2 Hz and recorded in the ENU (East, magnetic North, Up) coordinate system. 3. Between those stations along a cross-channel transect, a Stable Tidal Turbulence Mooring (STTM) positioned ~10 m above the seabed was deployed for one week during a spring tide. The STTM was outfitted with two Nortek Vector acoustic Doppler velocimeters equipped with inertial motion units (ADVs), a bottom-tracking downward-looking Teledyne RDI Workhorse 600 kHz ADCP to provide motion-corrected flow and turbulence characteristics at high temporal resolution, and an upward-looking Teledyne RDI Sentinel V20 ADCP. The STTM was deployed at 44.92098 N, 66.98922 W from 24-31 May. 4. A vessel-mounted Teledyne RDI Workhorse 300 kHz ADCP collected current data along three transects over two days, 4-5 April. The data processing used DOLfYN version 0.11.2. All hdf5 files (i.e., files ending in `.h5`) contained here can be opened using that version of DOLfYN (e.g., `dat = dolfyn.load('')`). All distances are in meters (e.g., depth, range, MLLW, hab, eta, z_), and all velocities in m/s. See the DOLfYN documentation https://lkilcher.github.io/dolfyn/), and/or the Nortek and Teledyne RDI documentation for additional details. Additional details on the dataset can be found in the README.pdf, including: - Format details of each data file. - How to regenerate the data-processing (using the files in the `wp2017_processing.zip` archive).

16 TIDAL AND WAVE POWER↗
DOE OSTI2018

Computational examination of two-phase microchannel heat transfer correlations with conjugate heat spreading

Microchannel flow boiling is an attractive thermal management strategy for ever-growing volumetric heat dissipation demands associated with electronic systems. Due to difficulties related to measurement at the microscales, the majority of researchers have chosen to study relatively simple situations with uniform heat flux. However, many applications involve local hotspots which give rise to highly non-uniform heat flux and temperature gradients due to heat spreading. This necessitates the consideration of conjugate heat transfer for accurate analysis. The current work is aimed at investigating conjugate heat transfer in a two-phase microchannel array. Experimental data was collected on R134a flow boiling heat transfer for very small hydraulic diameter (<100 µm) silicon microchannels with very high heat fluxes (>1 kW cm -2 ) applied via platinum strip heaters with a footprint size of 1 cm × 1 mm. The collected experimental data was then combined with detailed computational modeling utilizing finite element modeling with COMSOL Multiphysics and MATLAB to examine the applicability of five published heat transfer correlations for use in determining local heat transfer coefficients. A two-phase correlation from Agostini and Bontemps, developed for markedly different test parameters, provided the best computational agreement with an RMS temperature difference from experiment of 3.3 °C and a predicted peak perimeter-averaged heat transfer coefficient of 116 MW m -2 K -1 . Modeling confirms the presence of highly non-uniform local heat flux and correspondingly non-uniform local heat transfer coefficient. Finally, the results of this study make clear the need for better micro-scale two-phase correlations developed to predict local heat transfer coefficients at these small scales and high local heat fluxes.

42 ENGINEERING↗
DOE OSTI2017

Quantitative analysis of total reflection X-ray fluorescence from finely layered structures using XeRay

Total reflection x-ray fluorescence (TXRF) is a widely applicable experimental technique for studying chemical element distributions across finely layered structures at extremely high sensitivity. To promote and facilitate scientific discovery using TXRF, we developed a MATLAB-based software package with a graphical user interface, named XeRay, for quick, accurate, and intuitive data analysis. XeRay lets the user model any layered system, each layer with its independent chemical composition and thickness, and enables fine-tuned data fitting. The accuracy of XeRay has been tested in the analysis of TXRF data from both air/liquid interface and liquid/liquid interfacial studies and has been compared to literature results. In an air/liquid interface study, Ca 2+ sequestration was measured at a Langmuir monolayer of 1-stearoyl-2-oleoyl-sn-glycero-3-phosphatidic acid (SOPA) on a buffer solution of 1 mM CaCl2 at pH 7.5. Data analysis with XeRay reveals that each 1 nm2 of interfacial area contains 2.38 ± 0.06 Ca 2+ ions, which corresponds to a 1:1 ratio between SOPA headgroups and Ca 2+ ions, consistent with several earlier reports. For the liquid/liquid interface study of Sr 2+ enrichment at the dodecane/surfactant/water interface, analysis using XeRay gives a surface enrichment of Sr 2+ at $68^{+6}_{-5}$ Å 2 per ion, consistent with the result published for the same dataset.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗