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

Versatile semi-batch apparatus for manometric measurement of gas–solid reaction rates

In this study, we describe the design, construction, operation, and performance of a simple and versatile semi-batch reactor that is especially useful for measurement of gas/solid reaction rates at pressures in the range of 1 mTorr to 1500 Torr. The reactor operates by repeatedly imposing small AC modulations of reactant gas pressure on top of a much larger DC pressure background. Based on the rate of pressure relaxation following each AC pulse, the reaction rate is determined. Our design is characterized by modular construction from off-the-shelf, ultra-high-vacuum-compatible components, which facilitate easy retrofitting and adaptation to a range of experimental conditions. Automated experiment control and data acquisition is accomplished via a custom National Instruments © LabView virtual instrument. Data analysis is automated using a custom series of Mathworks © Matlab scripts. We demonstrate reactor performance through measurements of hydrogenation kinetics for a composite H 2 getter material consisting of 1,4-bis(phenylethynyl)benzene mixed with a palladium/carbon catalyst.

47 OTHER INSTRUMENTATION↗
DOE OSTI2016

SU-F-T-563: Delivered Dose Reconstruction of Moving Targets for Gated Volumetric Modulated Arc Therapy (VMAT)

Purpose: Actual delivered dose of moving tumors treated with gated volumetric arc therapy (VMAT) may significantly differ from the planned dose assuming static target. In this study, we developed a method which reconstructs actual delivered dose distribution of moving target by taking into account both tumor motion and dynamic beam delivery of gated VMAT, and applied to abdominal tumors. Methods: Fifteen dual-arc VMAT plans (Eclipse, Varian Medical Systems) for 5 lung, 5 pancreatic, and 5 liver cancer patients treated with gated VMAT stereotactic body radiotherapy (SBRT) were studied. For reconstruction of the delivered dose distribution, we divided each original arc beam into control-point-wise sub-beams, and applied beam isocenter shifting to each sub-beam to reflect the tumor motion. The tumor positions as a function of beam delivery were estimated by synchronizing the beam delivery with the respiratory signal which acquired during treatment. For this purpose, an in-house program (MATLAB, Mathworks) was developed to convert the original DICOM plan data into motion-involved treatment plan. The motion-involved DICOM plan was imported into Eclipse for dose calculation. The reconstructed delivered dose was compared to the plan dose using the dose coverage of gross tumor volume (GTV) and dose distribution of organs at risk (OAR). Results: The mean GTV dose coverage difference between the reconstructed delivered dose and the plan dose was 0.2 % in lung and pancreas cases, and no difference in liver cases. Mean D1000cc of ipsilateral lungs was reduced (0.8 ± 1.4cGy). Conclusion: We successfully developed a method of delivered dose reconstruction taking into account both respiratory tumor motion and dynamic beam delivery, and applied it to abdominal tumors treated with gated VAMT. No significant deterioration of delivered dose distribution indicates that interplay effect would be minimal even in the case of gated SBRT. This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIP) (No. 2015038710)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Historical global Earth System Model simulations with E3SM's dynamic root module

This data package contains global-scale Earth System Model simulations run with DOE's E3SM Model. The package includes one set of simulations run in the default configuration and a second set of simulations run with the dynamic root module enabled. These two datasets can be used to assess the impact that dynamic roots have on Land Surface fluxes such as carbon and water. The outputs are provided here on a 0.5 x 0.5 global grid at monthly resolution. The simulations were run with atmospheric forcing from the Global Soil Wetness Project. We provide specifically here Gross Primary Production, Transpiration, relative root fraction per soil layer as well as Air Temperature, Precipitation and all the parameterizations used for the simulations. Our purpose for generating and analyzing these simulations was to assess the role that root dynamics and foraging for water has in the recovery timescale of ecosystems following climate perturbation. All data provided are in the "netcdf" format using standard CF-1 (Climate and Forecast) Convention. The data are broken up into 20 year chunks to facilitate easier read-in. The data are all machine readable using various software platforms including Matlab, NCO, Panoply and GrADS.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Earth System Model Simulations of Ameriflux sites using dynamic roots and dynamic leaf to root allocation

This data package contains point-scale Earth System Model simulations run with DOE's E3SM Model. The package includes four sets of simulations, one set of simulations run in the default configuration, a second set of simulations run with the dynamic root module enabled, the third set with dynamic leaf to fine root allocation, and the fourth with the combination of dynamic roots and dynamic allocation. These simulations were run in point mode a Ameriflux sites, US-MMS, US-Moz, US-Syv, and US-NR1. These datasets can be used to assess the impact that dynamic roots have on Land Surface fluxes such as carbon and water. The simulations were run with atmospheric forcing from the Global Soil Wetness Project. We provide specifically here Gross Primary Production, Transpiration, relative root fraction per soil layer as well as Air Temperature, Precipitation and all the parameterizations used for the simulations. Our purpose for generating and analyzing these simulations was to assess the role that root dynamics and foraging for water has in the recovery timescale of ecosystems following climate perturbation. All data provided are in the "netcdf" format using standard CF-1 (Climate and Forecast) Convention. The data are all machine readable using various software platforms including Matlab, NCO, Panoply and R.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2008

Automated Power Assessment for Helicopter Turboshaft Engines

An accurate indication of available power is required for helicopter mission planning purposes. Available power is currently estimated on U.S. Army Blackhawk helicopters by performing a Maximum Power Check (MPC), a manual procedure performed by maintenance pilots on a periodic basis. The MPC establishes Engine Torque Factor (ETF), an indication of available power. It is desirable to replace the current manual MPC procedure with an automated approach that will enable continuous real-time assessment of available power utilizing normal mission data. This report presents an automated power assessment approach which processes data currently collected within helicopter Health and Usage Monitoring System (HUMS) units. The overall approach consists of: 1) a steady-state data filter which identifies and extracts steady-state operating points within HUMS data sets; 2) engine performance curve trend monitoring and updating; and 3) automated ETF calculation. The algorithm is coded in MATLAB (The MathWorks, Inc.) and currently runs on a PC. Results from the application of this technique to HUMS mission data collected from UH-60L aircraft equipped with T700-GE-701C engines are presented and compared to manually calculated ETF values. Potential future enhancements are discussed.

Simon, Donald L.↗
NASA NTRS2012

Accelerated Creep Testing of High Strength Aramid Webbing

A series of preliminary accelerated creep tests were performed on four variants of 12K and 24K lbf rated Vectran webbing to help develop an accelerated creep test methodology and analysis capability for high strength aramid webbings. The variants included pristine, aged, folded and stitched samples. This class of webbings is used in the restraint layer of habitable, inflatable space structures, for which the lifetime properties are currently not well characterized. The Stepped Isothermal Method was used to accelerate the creep life of the webbings and a novel stereo photogrammetry system was used to measure the full-field strains. A custom MATLAB code is described, and used to reduce the strain data to produce master creep curves for the test samples. Initial results show good correlation between replicates; however, it is clear that a larger number of samples are needed to build confidence in the consistency of the results. It is noted that local fiber breaks affect the creep response in a similar manner to increasing the load, thus raising the creep rate and reducing the time to creep failure. The stitched webbings produced the highest variance between replicates, due to the combination of higher local stresses and thread-on-fiber damage. Large variability in the strength of the webbings is also shown to have an impact on the range of predicted creep life.

Jones, Thomas C.↗
DOE OSTI2017

Photovoltaic Hosting Capacity of Feeders with Reactive Power Control and Tap Changers

This paper proposes an algorithm to determine photovoltaic (PV) hosting capacity of power distribution networks as a function of number of PV injection nodes, reactive power support from the PVs, and the sub-station load tap changers (LTCs). In the proposed method, several minute by minute simulations are run based on randomly chosen PV injection nodes, daily PV output profiles, and daily load profiles from a pool of high-resolution realistic data set. The simulation setup is built using OpenDSS and MATLAB. The performance of the proposed method is investigated in the IEEE 123-node distribution feeder for multiple scenarios. The case studies are performed particularly for one, two, five and ten PV injection nodes, and looking at the maximum voltage deviations. Case studies show that the PV hosting capacity of the 123-node feeder greatly differs with the number of PV injection nodes. We have also observed that the PV hosting capacity increases with reactive power support and higher tap position of sub-station LTC.

14 SOLAR ENERGY↗
DOE OSTI2023

Simulink Model of a New Energy EVG-005 5kW Hydropower Turbine

Simulink model for a New Energy 5kW hydropower turbine. ADCP data ("ds_streamwise_7_13.nc") and DC voltage, DC current, and rotor rotation observed from the New Energy EVG-005 Current Energy Converter (CEC) ("electrical_7_13_10ohms.nc") were collected at the Tanana River Test Site (TRTS) near Nenana Alaska. - Data was collected on July 13th, 2023 with a constant 10 ohms resistance applied with a DC load bank. - Simulink model is meant to resemble the electrical setup at the TRTS. - Model is initialized by running the "NewEnergy_2023_10hz.m" Matlab script. Then the Simulink model ("New_Energy_Model_PMSM.slx") can be run. - Results are processed with the "NewEnergy_2023_processing.m" Matlab code. The Matlab results are also saved in the "New_Energy_7_13_model_validation_results.mat" Matlab file. This can be directly loaded into the Matlab Workspace using the "Load()" command. - The timetable variables "electrical_model_downsampled" is the model results and the "electrical_limited" is the experimental data from the TRTS.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

CROCUS Air Quality Dataset from the University of Illinois Chicago (UIC), July 2024

This dataset was collected by the measurement system in the Atmosphere, Climate, and Ecosystems (ACE) Lab at the University of Illinois Chicago (UIC) from July 12 to July 31, 2024, as part of the Community Research on Climate and Urban Science (CROCUS) Urban Integrated Field Laboratory (UIFL) project, led by Argonne National Laboratory.To enhance understanding of urban air quality dynamics in Chicago, and as part of the CROCUS 2024 Urban Canyon Intensive Observation Period (IOP), several instruments were set up to provide continuous measurements of air quality parameters in Chicago during July 2024. These measurements cover both aerosols and gas-phase species. It focuses on particle size distribution (2.5–478 nm) measured by two Scanning Mobility Particle Sizers (SMPS) at a 4-min resolution, total particle number concentrations at a 1-s resolution, and chemical composition from a High-Resolution Time-of-Flight Aerosol Mass Spectrometer (AMS) at a 1-min resolution. Key gas-phase species, including NO, NO₂, SO₂, and O₃, are measured at a 1-min resolution, along with high-resolution NO and dimethyl sulfide (DMS) data from a Chemical Ionization Mass Spectrometer (CIMS). Volatile organic compound (VOC) data for toluene, isoprene, and benzene are provided by a GC-PID with a time resolution of 25 minutes.The data are formatted as NetCDF (.nc) files, making them easily accessible using common software such as MATLAB, R, and Python. Each parameter is stored in an individual dataset, which includes detailed instrument information in the header, as well as the corresponding sample start time and concentration/distribution data for each sample.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2014

Acoustics Research of Propulsion Systems

The liftoff phase induces high acoustic loading over a broad frequency range for a launch vehicle. These external acoustic environments are used in the prediction of the internal vibration responses of the vehicle and components. Present liftoff vehicle acoustic environment prediction methods utilize stationary data from previously conducted hold-down tests to generate 1/3 octave band Sound Pressure Level (SPL) spectra. In an effort to update the accuracy and quality of liftoff acoustic loading predictions, non-stationary flight data from the Ares I-X were processed in PC-Signal in two flight phases: simulated hold-down and liftoff. In conjunction, the Prediction of Acoustic Vehicle Environments (PAVE) program was developed in MATLAB to allow for efficient predictions of sound pressure levels (SPLs) as a function of station number along the vehicle using semi-empirical methods. This consisted of generating the Dimensionless Spectrum Function (DSF) and Dimensionless Source Location (DSL) curves from the Ares I-X flight data. These are then used in the MATLAB program to generate the 1/3 octave band SPL spectra. Concluding results show major differences in SPLs between the hold-down test data and the processed Ares I-X flight data making the Ares I-X flight data more practical for future vehicle acoustic environment predictions.

Gao, Ximing↗
DOE OSTI2015

MO-F-CAMPUS-I-01: EIT Imaging to Monitor Human Salivary Gland Functionality: A Feasibility Study

Purpose: Clinically, there exists a need to develop a non-invasive technique for monitoring salivary activity. In this study, we investigate the feasibility of a using the electrical conductivity information from Electrical Impedance Tomography (EIT) to monitor salivary flow activity. Methods: To acquire EIT data, eight Ag/AgCl ECG electrodes were placed around the mandible of the subject. An EIT scan was obtained by injecting current at 50 KHz, 0.4 mA through each pair of electrodes and recording voltage across other electrode pairs. The functional conductivity image was obtained through reconstruction of the voltage data, using Electrical Impedance Tomography and Diffuse Optical Tomography Reconstruction Software (EIDORS) in Matlab. In using EIDORS, forward solution was obtained using a user-defined finite element model shape and inverse solution was obtained using one-step Gaussian solver. EIT scans of volunteer research team members were acquired for three different physiological states: pre-stimulation, stimulation and post-stimulation. For pre-stimulation phase, data were collected in intervals of 5 minutes for 15 minutes. The salivary glands were then stimulated in the subject using lemon and the data were collected immediately. Post-stimulation data were collected at 4 different timings after stimulation. Results: Variations were observed in the electrical conductivity patterns near parotid regions between the pre- and post-stimulation stages. The three images acquired during the 15 minute pre-stimulation phase showed no major changes in the conductivity. Immediately after stimulation, electrical conductivity increased near parotid regions and 15 minutes later slowly returned to pre-stimulation level. Conclusion: In the present study involving human subjects, the change in electrical conductivity pattern shown in the EIT images, acquired at different times with and without stimulation of salivary glands, appeared to be consistent with the change in salivary gland activity. The conductivity changes imaged through EIT are potentially useful for the purpose of salivary monitoring.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2014

Acoustics Research of Propulsion Systems

The liftoff phase induces some of the highest acoustic loading over a broad frequency for a launch vehicle. These external acoustic environments are used in the prediction of the internal vibration responses of the vehicle and components. Thus, predicting these liftoff acoustic environments is critical to the design requirements of any launch vehicle but there are challenges. Present liftoff vehicle acoustic environment prediction methods utilize stationary data from previously conducted hold-down tests; i.e. static firings conducted in the 1960's, to generate 1/3 octave band Sound Pressure Level (SPL) spectra. These data sets are used to predict the liftoff acoustic environments for launch vehicles. To facilitate the accuracy and quality of acoustic loading, predictions at liftoff for future launch vehicles such as the Space Launch System (SLS), non-stationary flight data from the Ares I-X were processed in PC-Signal in two forms which included a simulated hold-down phase and the entire launch phase. In conjunction, the Prediction of Acoustic Vehicle Environments (PAVE) program was developed in MATLAB to allow for efficient predictions of sound pressure levels (SPLs) as a function of station number along the vehicle using semiempirical methods. This consisted, initially, of generating the Dimensionless Spectrum Function (DSF) and Dimensionless Source Location (DSL) curves from the Ares I-X flight data. These are then used in the MATLAB program to generate the 1/3 octave band SPL spectra. Concluding results show major differences in SPLs between the hold-down test data and the processed Ares IX flight data making the Ares I-X flight data more practical for future vehicle acoustic environment predictions.

Gao, Ximing↗
DOE OSTI2019

Frio 2 : Raw CASSM Datasets

This submission contains the raw SEG-2 seismic datasets acquired in the CASSM experiment, part of the Frio-2 GCS pilot. The directory also includes the log files generated by the Geode seismograph (geometrics) used during acquisition. All SEG-2 files have a ".dat" suffix and log files have a ".log" suffix. Various supplementary files are also included (see below) including a set of flow and seismic models used in interpretation. The earlyGathers.tar.gz resource includes raw data acquired during system testing and pre-injection periods. All files correspond to the single source continuous active source seismic monitoring (CASSM) dataset described in the following papers. Daley, T.M., Solbau, R.D., Ajo-Franklin, J.B., and Benson, S.M. 2007, “Continuous Active-Source Seismic Monitoring of CO2 Injection in a Brine Aquifer”, Geophysics, Vol. 72, No. 5, p. A57-A61 Daley, T., Ajo-Franklin, J.B., and C. Doughty, 2011, “Constraining the reservoir model of an injected CO2 plume with crosswell CASSM at the Frio-II Brine Pilot”, International Journal of Greenhouse Gas Control (IJGGC), Vol. 5, No. 2, pp. 1022-1030. Zhu, T., Ajo-Franklin, J., and T. Daley, 2017, “Spatio-temporal changes in seismic attenuation caused by injected CO2 at the Frio-II pilot site, Dayton TX, USA”, Journal of Geophysical Research – Solid Earth, Vol. 122, No. 9, pp. 7156-7171. cassmDataRaw.tar : The main tar file containing all of the SEG-2 field records and the Geode log files. Stacking, time sampling (dt), and date information should be in the .log file. Each file corresponds to a single stacked gather for the single source deployed. Data was recorded on a 24 channel downhole hydrophone array. fieldNotes.tar : contains some of the scanned observer logs from the experiments as well as some planning images fieldPhotos.tar : contains a small set of photographs taken during the experiment frio_flow_seismic_models.mat : a matlab structure containing a set of flow and seismic models used in interpreting the CASSM data (see IJGGC paper).

GCS↗
DOE OSTI2020

Radioactive Source Localization via Bayesian Particle Filter

In the event of a misplaced radioactive source or other emergency situation, measuring a radiation field, mapping its distribution, and determining a source location are essential tasks to ameliorating the situation. However, radiation fields may be extremely hazardous to human surveyors and minimizing received radiation doses is just as essential. Robots appear to be a potential solution to these problems. Beyond simply measuring radiation, the robot's computer processing capabilities offer a way to apply complex data analysis methods to radiation measurements in real-time. Methods which predict likely source locations can then feed this information into other processes, potentially improving path planning and enabling more efficient measurements. Given a robot mounted with a gamma-ray detector, can we: develop a methodology to account for detector performance across a wide range of source angles, distances, and photon energies? operate an autonomously navigating robot to effectively survey and characterize an area of interest? implement a data analysis method, conventionally used in measurements of motion, for source localization purposes? An open-source TurtleBot 3 robot, running Robot Operating System (ROS) on Ubuntu 16.04 LTS, was fitted with a Kromek GR1{sup R} Cadmium Zinc Telluride (CZT) solid-state gamma-ray detector. As a part of ROS, the packages OpenSlam, gmapping, and amcl were used to perform Simultaneous Localization and Mapping (SLAM), determining the robot's position and mapping the surrounding area. Data was acquired via Lidar mounted on top the TurtleBot 3. Detector Calibration Fit: The equation was fit to 365 counts of various energies, distances, and angles. A MATLAB{sup R} program was written to simulate measurements taken a robot on a random walk, with count data and positions discretized into finite element pixels. Using this program, a sample of 100 runs was performed on a map with a simulated source at the center, with a total of 200 of 2 pixels each. Similarly, multiple runs of the filter were performed on recorded robot measurement data. In both simulation and real tests, when corrected for errors (particles placed outside of bounds or on the robot, and simulation-specific errors), corresponding t-tests of predicted x and y-coordinates were within a 95% confidence interval of the actual position. For the real trial, these positions are slightly skewed right in the x-axis as the robot remained mainly to the left side of the source within the sample area. These simulations demonstrate potential validity for the usage of a particle filter as method of radioactive source localization. In the future, true real-time implementation and data fusion may further augment the performance of the robot to localize lost sources. Additionally, identification of multiple sources, determination of source types, and usage of a collimator are areas to potentially be explored.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Evaluating the Feasibility of Transactive Approach for Voltage Management Using Inverters of a PV Plant

This article evaluates the feasibility of a double-auction-based Transactive Energy System (TES) for engaging the reactive power (Var) capability of inverters in a utility-scale pho¬tovoltaic (PV) plant for voltage management in distribution feed¬ers. In this approach, a PV plant owner (seller) provides Var supply curves in terms of price-quantity pairs consider¬ing the inverters’ P-Q capability, efficiency, opportunity cost for real power curtail¬ment, losses, etc. At the same time, the utility (buyer) proposes Var demand curves in price-quantity pairs exhibiting the marginal price of Var. Then, market clearing points are obtained from the intersection points of supply and demand curves which are then used to create Var dispatch sig-nals from the PV plant. To conduct a fea¬sibility study of this approach, an Australian distribution network with a utility-scale PV plant (capacity 3.275 MWP) is modeled in RSCAD and simu¬lated in real-time on MATLAB, RSCAD-RTDS co-simulation platform using real historical data. Consid¬ering a use case, namely, conservation voltage reduction (CVR), the operational feasibility along with cost-benefit analysis is demonstrated in software-in-the-loop (SIL) and hardware-in-the-loop (HIL) platform. Another use case, exploring the benefit from tap-change reduction of step-voltage-regulator, is investigated. These studies explore further use cases, e.g., PV hosting capacity, network reconfiguration, etc., to be incorporated in TES for experiencing substantial economic benefits.

Alam, Mollah Rezaul↗
DOE OSTI2016

Tutorial: Crystal orientations and EBSD — Or which way is up?

Electron backscatter diffraction (EBSD) is an automated technique that can measure the orientation of crystals in a sample very rapidly. There are many sophisticated software packages that present measured data. Unfortunately, due to crystal symmetry and differences in the set-up of microscope and EBSD software, there may be accuracy issues when linking the crystal orientation to a particular microstructural feature. In this paper we outline a series of conventions used to describe crystal orientations and coordinate systems. These conventions have been used to successfully demonstrate that a consistent frame of reference is used in the sample, unit cell, pole figure and diffraction pattern frames of reference. We establish a coordinate system rooted in measurement of the diffraction pattern and subsequently link this to all other coordinate systems. A fundamental outcome of this analysis is to note that the beamshift coordinate system needs to be precisely defined for consistent 3D microstructure analysis. This is supported through a series of case studies examining particular features of the microscope settings and/or unambiguous crystallographic features. These case studies can be generated easily in most laboratories and represent an opportunity to demonstrate confidence in use of recorded orientation data. Finally, we include a simple software tool, written in both MATLAB® and Python, which the reader can use to compare consistency with their own microscope set-up and which may act as a springboard for further offline analysis. - Highlights: • Presentation of conventions used to describe crystal orientations • Three case studies that outline how conventions are consistent • Demonstrates a pathway for calibration and validation of EBSD based orientation measurements • EBSD computer code supplied for validation by the reader.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dissolved gas recovery from water using a sidestream hollow-fiber membrane module: First principles model synthesis and steady-state validation

This paper presents a first-principles model for the recovery of dissolved gases from liquids using a sidestream hollow-fiber membrane module. The model avoids the use of new empirical coefficients, thus providing a parametric understanding of the process behavior for future design and optimization of membrane modules. This type of first-principles model could be particularly useful when gas recovery is beneficial to biological or chemical reactions of interest, such as the acetogenesis reactions in two-stage anaerobic digesters. The steady-state behavior of the model was validated against both new experimental data for the recovery of H 2 , CH 4 and H 2 –CH 4 mixtures from pure water, as well as existing published data. The modeled gas recovery predictions agreed with experimental data to an absolute average error of 13%, and an average R value of 0.98. Parametric analysis of mixed-gas recovery suggests possible key transition points in the composition of the recovered gases. For example, at 40 °C, increasing trans-membrane pressure while keeping hydraulic residence time (HRT) under 0.5 s will result in an increase in the ratio of H 2 to CH 4 recovered. Otherwise, increasing trans-membrane pressure will instead decrease the ratio of H 2 to CH 4 recovered. The model has potential to be extended to transient analysis, but has yet to be validated with transient experimental data. Furthermore, this model was successfully implemented in both Python and MATLAB, and provides valuable insights for future net-energy optimization for anaerobic digestion systems with in-situ gas recovery.

Anaerobic Digestion↗
DOE OSTI2023

Accurate real space iterative reconstruction (RESIRE) algorithm for tomography

Tomography has made a revolutionary impact on the physical, biological and medical sciences. The mathematical foundation of tomography is to reconstruct a three-dimensional (3D) object from a set of two-dimensional (2D) projections. As the number of projections that can be measured from a sample is usually limited by the tolerable radiation dose and/or the geometric constraint on the tilt range, a main challenge in tomography is to achieve the best possible 3D reconstruction from a limited number of projections with noise. Over the years, a number of tomographic reconstruction methods have been developed including direct inversion, real-space, and Fourier-based iterative algorithms. Here, we report the development of a real-space iterative reconstruction (RESIRE) algorithm for accurate tomographic reconstruction. RESIRE iterates between the update of a reconstructed 3D object and the measured projections using a forward and back projection step. The forward projection step is implemented by the Fourier slice theorem or the Radon transform, and the back projection step by a linear transformation. Our numerical and experimental results demonstrate that RESIRE performs more accurate 3D reconstructions than other existing tomographic algorithms, when there are a limited number of projections with noise. Furthermore, RESIRE can be used to reconstruct the 3D structure of extended objects as demonstrated by the determination of the 3D atomic structure of an amorphous Ta thin film. We expect that RESIRE can be widely employed in the tomography applications in different fields. Finally, to make the method accessible to the general user community, the MATLAB source code of RESIRE and all the simulated and experimental data are available at https://zenodo.org/record/7273314.

97 MATHEMATICS AND COMPUTING↗