Search NASA⌕ Search

SEARCH · Search NASA

Results for “MatLab”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2024

Real space iterative reconstruction for vector tomography (RESIRE-V)

Tomography has had an important impact on the physical, biological, and medical sciences. To date, most tomographic applications have been focused on 3D scalar reconstructions. However, in some crucial applications, vector tomography is required to reconstruct 3D vector fields such as the electric and magnetic fields. Over the years, several vector tomography methods have been developed. Here, we present the mathematical foundation and algorithmic implementation of REal Space Iterative REconstruction for Vector tomography, termed RESIRE-V. RESIRE-V uses multiple tilt series of projections and iterates between the projections and a 3D reconstruction. Each iteration consists of a forward step using the Radon transform and a backward step using its transpose, then updates the object via gradient descent. Incorporating with a 3D support constraint, the algorithm iteratively minimizes an error metric, defined as the difference between the measured and calculated projections. The algorithm can also be used to refine the tilt angles and further improve the 3D reconstruction. To validate RESIRE-V, we first apply it to a simulated data set of the 3D magnetization vector field, consisting of two orthogonal tilt series, each with a missing wedge. Our quantitative analysis shows that the three components of the reconstructed magnetization vector field agree well with the ground-truth counterparts. We then use RESIRE-V to reconstruct the 3D magnetization vector field of a ferromagnetic meta-lattice consisting of three tilt series. Our 3D vector reconstruction reveals the existence of topological magnetic defects with positive and negative charges. We expect that RESIRE-V can be incorporated into different imaging modalities as a general vector tomography method. To make the algorithm accessible to a broad user community, we have made our RESIRE-V MATLAB source codes and the data freely available at https://github.com/minhpham0309/RESIRE-V.

47 OTHER INSTRUMENTATION↗
DOE OSTI2018

Auto-generated materials database of Curie and Néel temperatures via semi-supervised relationship extraction

article presents an auto-generated database of 39,822 records containing chemical compounds and their associated Curie and Néel magnetic phase transition temperatures. The database was produced using natural language processing and semi-supervised quaternary relationship extraction, applied to a corpus of 68,078 chemistry and physics articles. Evaluation of the database shows an estimated overall precision of 73%. Therein, records processed with the text-mining toolkit, ChemDataExtractor, were assisted by a modified Snowball algorithm, whose original binary relationship extraction capabilities were extended to quaternary relationship extraction. Consequently, its machine learning component can now train with ≤500 seeds, rather than the 4,000 originally used. Data processed with the modified Snowball algorithm affords 82% precision. Database records are available in MongoDB, CSV and JSON formats which can easily be read using Python, R, Java and MatLab. As a result, this makes the database easy to query for tackling big-data materials science initiatives and provides a basis for magnetic materials discovery.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Model systems for screening and investigation of lithium metal electrode chemistry and dendrite formation

The use of lithium metal as an electrode for electrochemical energy storage will provide a significant impact on practical energy storage technology. Unfortunately, the use of lithium metal is plagued with challenging chemical problems. Specifically, the formation of a solid electrolyte interphase layer and the nucleation and growth of lithium dendrites: both must be addressed and controlled in order to achieve a practically useable pure lithium metal electrode. Currently sophisticated experimental techniques and computationally expensive simulations are being used to probe these problems but these methods are arduous and time consuming which delays timely evaluation and insight into the rapidly changing field of advanced energy storage. Here, we report the use of DFT simulations of lithium nanoclusters to investigate and explore lithium metal chemistry with inexpensive computational methods to gain greater insight into electrochemical reductions and the nucleation and growth of dendrites. DME, LiTFSI, and LiFSI reduction energetics and structures with electrode effects from lithium metal are reported providing better physical description of the absolute reduction potential characterization. The electronic structure of the lithium nanoclusters were used to investigate the nucleation and growth of lithium dendrites from an ab initio perspective. The results demonstrate that kinetic processes have more control over non uniform growth than thermodynamic processes. Based on this information, a non ab initio model was created in Matlab that shows the initial stages of dendrite nucleation considering approximately 2000 atoms.

25 ENERGY STORAGE↗
DOE OSTI2022

Kinetic modelling of an environmentally friendly carbamazepine synthesis via urea and iminostilbene in batch and continuous processes

Accurate kinetic models for reaction systems allow for improved process understanding and greater quality control, which is particularly beneficial as the pharmaceutical industry shifts from batch to continuous manufacturing (CM). In this work, a first principles kinetic model has been developed for the synthesis of carbamazepine (CBZ) from iminostilbene and urea, starting in a batch reactor and subsequently in a continuous flow reactor. An eco-friendly reaction pathway using urea was selected to avoid the toxic reagents that are typically used for synthesis of CBZ. The kinetic parameters determined from batch reactions were utilized in a MATLAB based kinetic model to simulate the yield for the continuous process. Overall, good agreement between the model prediction and corresponding experimental values was observed for the batch and continuous reaction systems within the full factorial design space. However, the model slightly overpredicted the yield of the continuous reaction system for higher conversion values (>60%) since it did not account for the reverse reaction that can occur at the studied reaction conditions. Here, the use of broken order kinetics was compared with whole number orders, and it was determined that the whole number orders resulted in better agreement at all conversion values for the continuous system.

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

Chapter 6: Wind Turbine Control Design

In this chapter, the author introduced wind turbine control, discussing sensors and actuators, operating regions, and the operational controller loops. The author then described the different levels of models needed in the controller development process, emphasizing that the models needed for control design are a simplification of the detailed models used in control simulation testing. Scripts used in conjunction with MATLAB are used to synthesize the controller, and detailed models are used for controller simulations. We described the design of the basic operational controller based on simplified models, including the generator torque controller for optimizing energy capture, and the blade pitch controller for regulating turbine speed (or power). The author described how to account for the effects of turbine nonlinearities and actuator dynamics and saturation. We then gave an overview of modern state-space control design methods, which are useful when designing a controller to meet multiple control objectives. The author also described a variety of advanced multivariable control methods. The author gave an overview of the use of advanced sensors (such as lidar sensors for upwind wind-speed measurements) and actuators (outboard blade aerodynamic devices used in 'smart' blade technology). The use of individual blade pitch, although not considered an advanced actuator, was also described. The author concluded with some special control issues for offshore floating systems.

17 WIND ENERGY↗
DOE OSTI2021

Modular horizontal pendulum wave energy converter: Exploring feasibility to power ocean observation applications in the U.S. pacific northwest

Abstract Marine renewable energy as a power source for ocean observation applications has the potential to allow longer deployment operations due to the consistent, higher, and denser energy available from this resource. This additionally could encourage deployments in remote locations where maintenance is costly or resource availability is low if dependent on solar power. More importantly, gaps in spatial data could be filled. This paper examines the feasibility of a modular horizontal pendulum wave energy converter to power National Data Buoy Center's Self‐Contained Ocean Observations Payload (SCOOP) off the coast of Washington State, U.S. The effect on power output was studied when the pendulum's radius arm, mass, and power take‐off damping were varied. Results using Matlab toolbox WEC‐Sim revealed positive correlation between radius arm length and mass to power output, where power maximised for optimal damping values. Seasonal trends in power were not significant, where a 20 kg pendulum mass was needed to meet the SCOOP base power requirement of 5 W throughout the year.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

Evaluating Connected and Automated Vehicles in Co-Simulation Environment of Traffic Microsimulation and Vehicle Dynamics

Connected and automated vehicles (CAVs) have the potential to improve many aspects of the current transportation systems such as safety, mobility, and energy efficiency. In order to evaluate the benefits and impacts of a CAV, the CAV control algorithm is typically implemented on vehicles simulated in a traffic microsimulation environment. However, traffic microsimulation usually lacks detailed vehicle and powertrain dynamics, making it challenging to fully understand how a CAV control algorithm will perform and respond on an actual vehicle. Whether the same benefits measured in the simulation will also be observed in real-world remains an open question. One potential approach to fill in this gap is to conduct a co-simulation of traffic microsimulation with detailed vehicle and powertrain dynamics models, often developed in MATLAB Simulink. However, current microsimulation tools such as VISSIM and SUMO do not have a ready-to-use interface for co-simulation with vehicle dynamics and Simulink. Also, even if such an interface exists, it will be tool-specific, making it challenging to shift from one tool to another or test CAV controls in different tools. There are needs for tool-agnostic co-simulation as different microsimulation tools have their pros and cons, and researchers often need to use different tools based on the purposes of the simulation, project needs, and applications. In this work, Flexible Interface for X-in-the-loop Simulation (FIXS) is developed that can support the co-simulation of microsimulation, CAV control algorithm, and vehicle dynamics model in Simulink. Enabled by the FIXS, the benefit and performance of a CAV control algorithm can be better understood with the consideration of vehicle responses and dynamics. The connection to VISSIM and SUMO is handled internally by the interface, and users can easily switch tools by changing a configuration file. The co-simulation capability is demonstrated for a VISSIM eco-approach and departure CAV scenario and a SUMO cooperative merging scenario for both a passenger CAV and a class 8 heavy-duty connected and automated trucks.

Shao, Yunli↗
DOE OSTI2012

Measurement of parameters in Indus-2 synchrotron radiation source

The paper presents the measurement of optics parameters in Indus-2 synchrotron radiation source, which include betatron tune, beta function, dispersion function, natural chromaticity, corrected chromaticity, central RF frequency, momentum compaction factor, and linear betatron coupling. Two methods were used for beta function measurement; a conventional quadrupole scan method and a method using the fitting of the orbit response matrix. A robust Levenberg-Marquardt algorithm was used for nonlinear least square fitting of the orbit response matrix. In this paper, detailed methods for the parameter measurements are described. The measured results are discussed and compared with the theoretical values obtained using accelerator simulation code Accelerator Toolbox in MATLAB.

43 PARTICLE ACCELERATORS↗
DOE OSTI2013

Simplified modeling of liquid sodium medium with temperature and velocity gradient using real thermal-hydraulic data. Application to ultrasonic thermometry in sodium fast reactor

In the framework of the French R and D program for the Generation IV reactors and specifically for the sodium cooled fast reactors (SFR), studies are carried out on innovative instrumentation methods in order to improve safety and to simplify the monitoring of fundamental physical parameters during reactor operation. The aim of the present work is to develop an acoustic thermometry method to follow up the sodium temperature at the outlet of subassemblies. The medium is a turbulent flow of liquid sodium at 550 Degree-Sign C with temperature inhomogeneities. To understand the effect of disturbance created by this medium, numerical simulations are proposed. A ray tracing code has been developed with Matlab Copyright-Sign in order to predict acoustic paths in this medium. This complex medium is accurately described by thermal-hydraulic data which are issued from a simulation of a real experiment in Japan. The analysis of these results allows understanding the effects of medium inhomogeneities on the further thermometric acoustic measurement.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2013

Sensitivity of blackbody effective emissivity to wavelength and temperature: By genetic algorithm

A variable-temperature blackbody (VTBB) is used to calibrate an infrared radiation thermometer (pyrometer). The effective emissivity (ε{sub eff}) of a VTBB is dependent on temperature and wavelength other than the geometry of the VTBB. In the calibration process the effective emissivity is often assumed to be constant within the wavelength and temperature range. There are practical situations where the sensitivity of the effective emissivity needs to be known and correction has to be applied. We present a method using a genetic algorithm to investigate the sensitivity of the effective emissivity to wavelength and temperature variation. Two matlab® programs are generated: the first to model the radiance temperature calculation and the second to connect the model to the genetic algorithm optimization toolbox. The effective emissivity parameter is taken as a chromosome and optimized at each wavelength and temperature point. The difference between the contact temperature (reading from a platinum resistance thermometer or liquid in glass thermometer) and radiance temperature (calculated from the ε{sub eff} values) is used as an objective function where merit values are calculated and best fit ε{sub eff} values selected. The best fit ε{sub eff} values obtained as a solution show how sensitive they are to temperature and wavelength parameter variation. Uncertainty components that arise from wavelength and temperature variation are determined based on the sensitivity analysis. Numerical examples are considered for illustration.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2013

A reliable control system for measurement on film thickness in copper chemical mechanical planarization system

In recent years, a variety of film thickness measurement techniques for copper chemical mechanical planarization (CMP) are subsequently proposed. In this paper, the eddy-current technique is used. In the control system of the CMP tool developed in the State Key Laboratory of Tribology, there are in situ module and off-line module for measurement subsystem. The in situ module can get the thickness of copper film on wafer surface in real time, and accurately judge when the CMP process should stop. This is called end-point detection. The off-line module is used for multi-points measurement after CMP process, in order to know the thickness of remained copper film. The whole control system is structured with two levels, and the physical connection between the upper and the lower is achieved by the industrial Ethernet. The process flow includes calibration and measurement, and there are different algorithms for two modules. In the process of software development, C++ is chosen as the programming language, in combination with Qt OpenSource to design two modules’ GUI and OPC technology to implement the communication between the two levels. In addition, the drawing function is developed relying on Matlab, enriching the software functions of the off-line module. The result shows that the control system is running stably after repeated tests and practical operations for a long time.

42 ENGINEERING↗
DOE OSTI2013

Valuation of exotic options in the framework of Levy processes

In this paper we explore a straightforward procedure to price derivatives by using the Monte Carlo approach when the underlying process is a jump-diffusion. We have compared the Black-Scholes model with one of its extensions that is the Merton model. The latter model is better in capturing the market’s phenomena and is comparative to stochastic volatility models in terms of pricing accuracy. We have presented simulations of asset paths and pricing of barrier options for both Geometric Brownian motion and exponential Levy processes as it is the concrete case of the Merton model. A desired level of accuracy is obtained with simple computer operations in MATLAB for efficient computational time.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Improved algorithm for processing grating-based phase contrast interferometry image sets

Grating-based X-ray and neutron interferometry tomography using phase-stepping methods generates large data sets. An improved algorithm is presented for solving for the parameters to calculate transmissions, differential phase contrast, and dark-field images. The method takes advantage of the vectorization inherent in high-level languages such as Mathematica and MATLAB and can solve a 16 × 1k × 1k data set in less than a second. In addition, the algorithm can function with partial data sets. This is demonstrated with processing of a 16-step grating data set with partial use of the original data chosen without any restriction. Also, we have calculated the reduced chi-square for the fit and notice the effect of grating support structural elements upon the differential phase contrast image and have explored expanded basis set representations to mitigate the impact.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Consideration of a ultracold neutron source in two-dimensional cylindrical geometry by taking simulated boundaries

A new idea to calculate ultracold neutron (UCN) production by using Monte Carlo simulation method to calculate the cold neutron (CN) flux and an analytical approach to calculate the UCN production from the simulated CN flux was given. A super-thermal source (UCN source) was modeled based on an arrangement of D{sub 2}O and solid D{sub 2} (sD{sub 2}). The D{sub 2}O was investigated as the neutron moderator, and sD{sub 2} as the converter. In order to determine the required parameters, a two-dimensional (2D) neutron balance equation written in Matlab was combined with the MCNPX simulation code. The 2D neutron-transport equation in cylindrical (ρ − z) geometry was considered for 330 neutron energy groups in the sD{sub 2}. The 2D balance equation for UCN and CN was solved using simulated CN flux as boundary value. The UCN source dimensions were calculated for the development of the next UCN source. In the optimal condition, the UCN flux and the UCN production rate (averaged over the sD{sub 2} volume) equal to 6.79 × 10{sup 6} cm{sup −2}s{sup −1} and 2.20 ×10{sup 5} cm{sup −3}s{sup −1}, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Methods and preliminary measurement results of liquid Li wettability

A test of lithium wettability was performed in high vacuum (< 3 × 10{sup −4} Pa). High magnification images of Li droplets on stainless steel substrates were produced and processed using the MATLAB{sup ®} program to obtain clear image edge points. In contrast to the more standard “θ/2” or polynomial fitting methods, ellipse fitting of the complete Li droplet shape resulted in reliable contact angle measurements over a wide range of contact angles. Using the ellipse fitting method, it was observed that the contact angle of a liquid Li droplet on a stainless steel substrate gradually decreased with increasing substrate temperature. The critical wetting temperature of liquid Li on stainless steel was observed to be about 290 °C.

36 MATERIALS SCIENCE↗
DOE OSTI2014

New approach for calibration the efficiency of HpGe detectors

This work evaluates the efficiency calibrating of HpGe detector coupled with Canberra GC3018 with Genie 2000 software and Ortec GEM25-76-XLB-C with Gamma Vision software; available at Neutron activation analysis laboratory in Malaysian Nuclear Agency (NM). The efficiency calibration curve was constructed from measurement of an IAEA, standard gamma–point sources set composed by {sup 214}Am, {sup 57}Co, {sup 133}Ba, {sup 152}Eu, {sup 137}Cs and {sup 60}Co. The efficiency calibrations were performed for three different geometries: 5, 10 and 15 cm distances from the end cap detector. The polynomial parameters functions were simulated through a computer program, MATLAB in order to find an accurate fit to the experimental data points. The efficiency equation was established from the known fitted parameters which allow for the efficiency evaluation at particular energy of interest. The study shows that significant deviations in the efficiency, depending on the source-detector distance and photon energy.

42 ENGINEERING↗