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 1,027 records · Page 57

Radioisotope Power System Dose Estimation Tool

The Radioisotope Power System Dose Estimation Tool (RPS-DET) is used to create, perform and analyze particle-transport simulations of radioisotope power systems (RPS). RPS-DET consists of a Matlab graphical user interface (GUI) to generate and analyze Monte Carlo N-Particle (MCNP) 6 input and output files. Simulation geometries are created using a library-based method of combining multiple RPS units and environmental scenarios while source terms are produced from databases of time-dependent photon and neutron spectra generated by the software packages IsoChain and ORIGEN. Simulation results are saved in a three-dimensional format for neutron and photon flux and dose in tissue or silicon and may be interactively inspected via the GUI. RPS-DET allows flexibility in the simulation-space by providing the ability to mix-and-match RPS units with relevant environments along with capturing the effects of fuel-aging for the particle-energy-spectra and intensity. This work has been benchmarked against measurements of existing RPS units and found to comparably predict neutron and gamma flux, absorbed dose, and dose equivalent in humans and electronics (silicon).

Smith, Michael↗

Aeroelastic oscillations of a cantilever with structural nonlinearities: theory and numerical simulation.

Our study details the derivation of the nonlinear equations of motion for the axial, biaxial bending and torsional vibrations of an aeroelastic cantilever undergoing rigid body (pitch) rotation at the base. The primary attenstion is focussed on the geometric nonlinearities of the system, whereby the aeroelastic load is modeled by the theory of linear quasisteady aerodynamics. This modelling effort is intended to mimic the wind-tunnel experimental setup at the Royal Military College of Canada. While the derivation closely follows the work of Hodges and Dowell [1] for rotor blades, this aeroelastic system contains new inertial terms which stem from the fundamentally different kinematics than those exhibited by helicopter or wind turbine blades. Using the Hamilton’s principle, a set of coupled nonlinear partial differential equations (PDEs) and an ordinary differential equation (ODE) are derived which describes the coupled axial-bending-bending-torsion-pitch motion of the aeroelastic cantilever with the pitch rotation. The finite dimensional approximation of the coupled system of PDEs are obtained using the Galerkin projection, leading to a coupled system of ODEs. Subsequently, these nonlinear ODEs are solved numerically using the built-in MATLAB implicit ODE solver and the associated numerical results are compared with those obtained using Houbolt’s method. It is demonstrated that the system undergoes coalescence flutter, leading to a limit cycle oscillation (LCO) due to coupling between the rigid body pitching mode and teh flexible mode arising from the flapwise bending motion.

42 ENGINEERING↗

SFDBSI_GLS v. 1.0

Matlab code for inversion of frequency domain, electrostatic geophysical data in terms of scalar scattering amplitudes in the subsurface. The data is assumed to be the difference between two measurements: electric field measurements prior to the injection of an electrically conductive proppant, and the electric field measurements after proppant injection. The proppant is injected into the subsurface via a well, and its purpose is to prop open fractures created by hydraulic fracturing. In both cases the illuminating electric field is assumed to be a vertically incident plane wave. The inversion strategy is to solve a set of linear system of equations, where each equation defines the amplitude of a candidate scattering volume. The model space is defined by M potential scattering locations and the frequency domain (of which there are k frequencies) data are recorded on N receivers. The solution thus solves a kN x M system of linear equations for M scalar amplitudes within the user-defined solution space. Practical Application: Oilfield environments where observed electrostatic geophysical data can reasonably be assumed to be scattered by subsurface proppant volumes. No field validation examples have so far been provided.

Poppeliers, Christian↗

MUFFSgenMC: An Open Source Muon Flexible Framework for Spectral Generation for Monte Carlo Applications

A cosmic ray Muon Flexible Framework for Spectral Generation for Monte Carlo Applications (MUFFSgenMC) has been developed to support state-of-the-art cosmic ray muon detection and tomographic applications. The flexible framework allows for easy and fast creation of source terms for popular Monte Carlo applications such as GEANT4 and MCNP, which simplifies the process of simulations. The flexible framework contains a variety of analytical, numerical, and parametric models that capture the main characteristics of the muon energy and angular distributions in the range 0 to 90 and at arbitrary energies. The primary model for muon energy distribution is an implementation of the Smith and Duller phenomenological model. The predictions for muon energies in the range from 1 GeV to 1 TeV and zenith angles 0 to 90 are validated against an extensive series of experimental spectrum measurements. Finally, the open source code available on Matlab’s Mathworks File Exchange (https://www.mathworks.com/matlabcentral/fileexchange/65585-muffsgenmc) was written to allow users to easily modify and expand the code. MUFFSgenMC can facilitate the study of cosmic ray muons for nuclear nonproliferation applications and will enable the development of new muon-based detection and imaging techniques by providing an easy source term generator.

Greulich, Christopher↗

Operation Data from the UW Natural Circulation FLiBe Flow Loop

Operating measurements from three multi-week test campaigns of the natural circulation FLiBe loop. This system is investigating thermal hydraulic behavior of the molten salt BeF_2 -LiF (33 - 67 mol %). The system behaves in a stable fashion, but shows unusual local transients, such as flow break-down in the riser and thermal jumps at the cooler exit. Heat transfer shows promising trends that FLiBe will behave as a normal heat transfer fluid, as long as salt purity can be maintained. The most recent test shows heat transfer degradation that is likely a result of oxides and impurities. The data archive contains as-built dimensions, Matlab analysis codes, and the raw data files.

42 ENGINEERING↗

Development of an Automated J-R Curve Analysis Software Based on the ASTM E1820-18 Normalization Method

The J-integral versus crack-growth resistance curve (J-R curve) is a useful tool for evaluating material fracture toughness in the ductile region. A widely accepted practice for conducting J-R curve testing is ASTM standard E1820, which introduces the normalization method. This method does not require compliance measurements during testing, significantly simplifying J-R curve testing compared with the conventional elastic unloading compliance (EUC) method. Despite the advantage in testing, analysis procedures for the normalization method are very complicated. Also, there is no open source software currently available to perform such analysis.The authors developed an open source automated J-R curve-analysis software based on the normalization method in ASTM E1820-18. Source codes were written in MATLAB and the compiled software, in the form of a stand-alone executable, has a user-friendly graphic interface and is readily compatible with modern Windows operating systems. Both source codes and the compiled executable file are available for download at https://code.ornl.gov/xc8/ANJR.git. Comparisons between resulting J-R curves from the software and manual analysis procedures show essentially identical results, indicating the suitability of applying the automated software for J-R curve analysis.

36 MATERIALS SCIENCE↗

Blackstart of Power Grids with Inverter-Based Resources: Preprint

This paper presents the findings of our investigation into inverter-based resource- (IBR-) driven blackstart of electric grids. Four potential black-start configurations with different setups are presented. To evaluate the technical feasibility of IBR-driven black start in the four configurations, a behavioral model of inverters that mimics current-limited inverter operation is developed using variable resistors in the MATLAB Simulink/Simscape environment. The inverter model is connected to an induction motor through transformers and a transmission line to simulate its startup. Simulation results show that even with the limited current supply capability of inverters because of their physical constraints, IBRs can black-start a motor under certain conditions. Results also show that by using soft-start techniques, such as ramped supply voltage, inrush currents can be reduced, thereby expanding the conditions under which IBRs can provide black-start support. Simulation results with different scenarios lead to discussions and key takeaways that could be valuable for further IBR-driven black-start research.

14 SOLAR ENERGY↗

Dynamic Equivalence of Large-Scale Power Systems Based on Boundary Measurements

Parallel computing helped speed up many tasks that can be done independently. The advance in gaming industry did not produce a technology that helps the power industry to reduce the computation time for dynamic simulation in individual cases. To perform faster than real-time simulation for the purpose of predicting power system dynamic trajectory, power industry continues to struggle to reduce the system size and simulation time of large-scale power systems while keeping its dynamic behavior under various disturbances. There is also an acute need for real-time dynamic model reduction as more renewables enter the generation mix with dramatic changes in the generation outputs. All existing model reduction solutions are based on having access to a detailed dynamic model of the system. The system changes by the minutes, dynamic models are only updated annually. To solve this problem, wide-area measurements obtained by the phasor measurement unit (PMU) at the boundaries between the reduced system and the study system is used to represent the external area. An artificial neuro-fuzzy inference system (ANFIS) is established to perform the mapping of the measurement to the external equivalent model. Here the external area is regarded as a black-box. Model reduction studies are conducted on the Northeast Power Coordinating Council (NPCC) under various types of contingencies, by using a co-simulation approach between PSS/E and MATLAB. The results look very promising and will be discussed in this paper.

Tong, Ning↗

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↗

Actinide Mixing Optimization Solver

The Actinide Mixing Optimization Solver (AMOS) is a set of programs written in MATLAB that computes theoretical mixes of different uranium, plutonium, thorium, and other actinide materials given a list of isotopic reference materials. Constraints are applied depending on the desired product and a mixture is computed to meet the specified requirements. The optimal mixture consists of a product that uses the lowest possible number of initial actinide materials and has the lowest cost. AMOS is currently applied to particular reference material requests from the sponsoring agency (NNSA) on behalf of the IAEA.

Baldwin, Aaron Taylor↗

Comparison of Simulated and Experimentally Measured Detector Impulse Responses to Pulsed X-rays

The Gamma Array Simulation Toolkit (GAST) is a multi-physics software stack comprised of Geant4, MATLAB, and experimentally measured photomultiplier tube (PMT) impulse response waveforms that is used to model and inform the end-to-end performance of radiation detector systems. In order to validate the ability of GAST to predict the impulse response of a modeled detector configuration operated in current mode, experimental measurements using the Nevada National Security Site Transformational Diagnostics and Imaging 2 MeV endpoint energy electron linear accelerator (linac) with a tungsten target were performed. Seven different detector configurations consisting of unique combinations of scintillator material, size, reflectivity conditions, presence or absence of light guide, and PMT were evaluated. The resultant impulse response of each detector configuration was recorded and compared to a GAST simulation that emulated the respective linac experiment conditions. The comparison of the experimental and simulated detector impulse responses indicated that GAST can reliably predict detector impulse responses.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Development and Validation of a Simulation Testbed for the Intelligent Building Agents Laboratory (IBAL) Using TRNSYS

This paper documents the development and validation of a dynamic primary cooling and thermal storage system simulation testbed. The system simulation testbed, sIBAL, is based on the Intelligent Building Agents Laboratory (IBAL) at the National Institute of Standards and Technology (NIST), which is a research infrastructure and testbed for the development, evaluation, and demonstration of intelligent control algorithms. The sIBAL testbed developed in this project will serve as a virtual twin of the real facility for future control algorithm development. The details of the methodologies used to develop and validate the simulation testbed, which replicates the dynamic behaviors of the primary cooling and ice storage system in the IBAL facility, are presented. The simulation testbed was developed in TRNSYS using built-in component models and MATLAB functions to replicate the two water-cooled chillers, a thermal storage tank, pumps, valves and other components for four different operation modes. Experiments on IBAL components were designed and executed to generate experimental data for model development and verification of the simulation platform. The validation of the simulation results was carried out in two phases: 1) independent component simulations for the chillers and thermal storage tank, and 2) a combined testbed simulation of the entire hydronic system. Comparison of simulation results to the experimental data obtained from the IBAL facility showed errors within 1 °C for the temperatures outputs of both the chiller and the thermal storage model. The error is within an acceptable range for further intelligent control algorithms development. The findings from the study are summarized and presented along with areas where additional research is needed. In addition, data filtering procedures and model refinement measures utilized to improve the accuracy and accelerate the computation time of the simulation are presented.

cooling systems↗

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↗

AXEAP (ARGONNE X-RAY EMISSION PACKAGE)

Argonne X-ray Emission Package (AXEAP), a singular purpose software package for processing X-ray emission (XES) images collected with a 2-dimensional position sensitive pixel array detector, has been developed. AXEAP can rapidly convert XES image files into a spectral form by applying parallel computation and unsupervised machine learning to compute vast amount of image data. Special focus has been placed on designing user-friendly-interface for processing multiple edges, non-resonant and resonant x-ray emission image analysis, in order to make data processing quick and easy. AXEAP is free software and is written in MATLAB, armed with powerful libraries and toolboxes. The software runs on all common operating systems such as Linux, Window, and Mac.

SUN, CHENGJUN↗

The Marine and Hydrokinetic ToolKit for Data Quality Control and Analysis: Preprint

The ability to handle data is critical at all stages of marine energy (ME) development. The marine hydrokinetic toolkit (MHKiT) is an open-source marine energy software, which includes modules for ingesting, applying quality control, processing, visualizing, and managing data. MHKiT-Python and MHKiT-MATLAB provide robust and verified functions that are needed by the ME community to standardize data processing. Calculations and visualizations adhere to International Electrotechnical Commission (IEC) technical specifications and other guidelines. A resource assessment of NDBC buoy 46050 near PACWAVE is performed using MHKiT and discusses comparisons to the resource assessment provided performed by Dunkel et al.

marine energy↗

Blockchain-Enabled Cyber-Secure Microgrid Control Using Consensus Algorithm: Preprint

This paper proposes a distributed control method based on a consensus algorithm for distributed energy resources (DERs) using blockchain as a secure communication medium for cyber resilience. Each DER communicates with a local blockchain server that is maintained by the distributed ledger technology to allow for securely sharing local measurements among neighboring assets to achieve the global control objectives, i.e., voltage and frequency regulation as well as accurate power sharing among the DERs, including collective grid-forming capability. To prove that the distributed control can retain system stability under a blockchain-induced variable communication delay, Lyapunov function-based stability analysis is carried out. This paper demonstrates the concept on an 11-bus test case developed in MATLAB Simulink, which has been modified from the IEEE 9-bus test case, to study dynamic operations of the five inverter-based DERs working with a blockchain-induced variable delay. The results validate the superior performance of the proposed control method, compared to heavily compromised operations of the other test cases without the security measure, affected by heavy communication delays and communication interruptions.

blockchain↗

The Marine and Hydrokinetic Toolkit (Mhkit) for Data Quality Control and Analysis

The ability to handle data is critical at all stages of marine energy (ME) development. The marine hydrokinetic toolkit (MHKiT) is an open-source marine energy software, which includes modules for ingesting, applying quality control, processing, visualizing, and managing data. MHKiT-Python and MHKiT-MATLAB provide robust and verified functions that are needed by the ME community to standardize data processing. Calculations and visualizations adhere to International Electrotechnical Commission (IEC) technical specifications and other guidelines. A resource assessment of NDBC buoy 46050 near PACWAVE is performed using MHKiT and discusses comparisons to the resource assessment provided performed by Dunkel et al.

marine energy↗

Geothermal play-fairway analysis of Washington State prospects

Herein, analysis of existing geologic, geophysical, and geochemical data revealed areas with elevated heat and permeability, defining three promising plays along the central axis of the magmatic arc of Washington State: Mount St. Helens, Wind River valley, and the southeast flank of Mount Baker. These areas are geothermal ‘fairways’, or locations with high geothermal resource potential based on modeling. This project applies innovative data analyses (Poly3D) in conjunction with methods proven in previous geothermal favorability studies (ArcGIS and MATLAB) to extract new value from existing public data. Heat and permeability potential models for each play are weighted using the Analytical Hierarchy Process, and combined to map geothermal resource potential and locations for further data collection at each site. Heat potential is based on: (1) temperature gradients, (2) volcanic vents, (3) Quaternary intrusive rocks, (4) spring temperature, and (5) spring geothermometry. Permeability potential is based on: (1) slip and (2) dilation tendency on mapped and seismic faults (3) maximum shear strain rate, and (4) dilatational strain rate at the surface, (5) modeled fault displacement distribution, and (6) displacement gradient, (7) shear, and (8) tensile fracture density, and (9) local geology and geophysical data. This information is vital for revealing the small scale heat and permeability potential of each study area, locating areas of undiscovered or untapped resources, and reducing the risk and cost involved in greenfield exploration and development. Preliminary analysis of the three play areas shows promise for geothermal development. Uncertainty modeling is underway and the results will guide future exploration plans.

15 GEOTHERMAL ENERGY↗