Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

LineVISAR: A Fringe-Trace Data Analysis Program

The line-imaging ORVIS or VISAR provides velocity as a function of position and time for a line on an experimental setup via a streak camera record of interference fringes. This document describes a Matlab-based program which guides the user through the process of converting these fringe data to a velocity surface. The data reduction is of the "fringe trace" type, wherein the changes in velocity at a given position on the line are calculated based on fringe motion past that point. The analyst must establish the fringe behavior up front, aided by peak-finding routines in the program. However, the later work of using fringe jumps to compensate for phase problems in other analysis techniques is greatly reduced. This program is not a standard GUI construction, and is prescriptive. At various points it saves the progress, allowing later restarts from those points.

97 MATHEMATICS AND COMPUTING↗

MueMex User's Guide

This is the user guide for the Matlab interface in MUELU.

97 MATHEMATICS AND COMPUTING↗

Final Technical Report: Distributed Controls for High Penetrations of Renewables

The goal of this effort was to apply four potential control analysis/design approaches to the design of distributed grid control systems to address the impact of latency and communications uncertainty with high penetrations of photovoltaic (PV) generation. The four techniques considered were: optimal fixed structure control; Nyquist stability criterion; vector Lyapunov analysis; and Hamiltonian design methods. A reduced order model of the Western Electricity Coordinating Council (WECC) developed for the Matlab Power Systems Toolbox (PST) was employed for the study, as well as representative smaller systems (e.g., a two-area, three-area, and four-area power system). Excellent results were obtained with the optimal fixed structure approach, and the methodology we developed was published in a journal article. This approach is promising because it offers a method for designing optimal control systems with the feedback signals available from Phasor Measurement Unit (PMU) data as opposed to full state feedback or the design of an observer. The Nyquist approach inherently handles time delay and incorporates performance guarantees (e.g., gain and phase margin). We developed a technique that works for moderate sized systems, but the approach does not scale well to extremely large system because of computational complexity. The vector Lyapunov approach was applied to a two area model to demonstrate the utility for modeling communications uncertainty. Application to large power systems requires a method to automatically expand/contract the state space and partition the system so that communications uncertainty can be considered. The Hamiltonian Surface Shaping and Power Flow Control (HSSPFC) design methodology was selected to investigate grid systems for energy storage requirements to support high penetration of variable or stochastic generation (such as wind and PV) and loads. This method was applied to several small system models.

14 SOLAR ENERGY↗

High Penetration Solar PV Deployment Sunshine State Solar Grid Initiative (SUNGRIN)

The report provides results from the Sunshine State Solar Grid Initiative (SUNGRIN) high penetration solar PV deployment project led by Florida State University’s (FSU) Center for Advanced Power Systems (CAPS). FSU CAPS and industry and university partners have completed a five-year effort aimed at enabling effective integration of high penetration levels of grid-connected solar PV generation. SUNGRIN has made significant contributions in the development of simulation-assisted techniques, tools, insight and understanding associated with solar PV effects on electric power system (EPS) operation and the evaluation of mitigation options for maintaining reliable operation. An important element of the project was the partnership and participation of six major Florida utilities and the Florida Reliability Coordinating Council (FRCC). Utilities provided details and data associated with actual distribution circuits having high-penetration PV to use as case studies. The project also conducted foundational work supporting future investigations of effects at the transmission / bulk power system level. In the final phase of the project, four open-use models with built-in case studies were developed and released, along with synthetic solar PV data sets, and tools and techniques for model reduction and in-depth parametric studies of solar PV impact on distribution circuits. Along with models and data, at least 70 supporting MATLAB functions have been developed and made available, with complete documentation.

14 SOLAR ENERGY↗

Echo™ User Manual

Echo™ is a MATLAB-based software package designed for robust and scalable analysis of complex data workflows. An alternative to tedious, error-prone conventional processes, Echo is based on three transformative principles for data analysis: self-describing data, name-based indexing, and dynamic resource allocation. The software takes an object-oriented approach to data analysis, intimately connecting measurement data with associated metadata. Echo operations in an analysis workflow automatically track and merge metadata and computation parameters to provide a complete history of the process used to generate final results, while automated figure and report generation tools eliminate the potential to mislabel those results. History reporting and visualization methods provide straightforward auditability of analysis processes. Furthermore, name-based indexing on metadata greatly improves code readability for analyst collaboration and reduces opportunities for errors to occur. Echo efficiently manages large data sets using a framework that seamlessly allocates resources such that only the necessary computations to produce a given result are executed. Echo provides a versatile and extensible framework, allowing advanced users to add their own tools and data classes tailored to their own specific needs. Applying these transformative principles and powerful features, Echo greatly improves analyst efficiency and quality of results in many application areas.

97 MATHEMATICS AND COMPUTING↗

Analysis of the LSC microbunching instability in MaRIE linac reference design

In this report we estimate the effect of the microbunching instability in the MaRIE XFEL linac. The reference design for the linac is described in a separate report. The parameters of the L1, L2, and L3 linacs as well as BC1 and BC2 bunch compressors were the same as in the referenced report. The beam dynamics was assumed to be linear along the accelerator (which is a reasonable assumption for estimating the effect of the microbunching instability). The parameters of the bunch also match the parameters described in the referenced report. Additionally, it was assumed that the beam radius is equal to R = 100 m and does not change along linac. This assumption needs to be revisited at later studies. The beam dynamics during acceleration was accounted in the matrix formalism using a Matlab code. The input parameters for the linacs are: RF peak gradient, RF frequency, RF phase, linac length, and initial beam energy. The energy gain and the imposed chirp are calculated based on the RF parameters self-consistently. The bunch compressors are accounted in the matrix formalism as well. Each chicane is characterized by the beam energy and the R56 matrix element. It was confirmed that the linac and beam parameters described previously provide two-stage bunch compression with compression ratios of 10 and 20 resulting in the bunch of 3kA peak current.

43 PARTICLE ACCELERATORS↗

ROMUSE 2.0 User Manual

ROMUSE (Reduced Order Modeling Based Uncertainty/Sensitivity Estimator) is an effort within the Consortium for Advanced Simulation of Light water reactors (CASL) to provide an analysis tool to be used in conjunction with reactor core simulators, especially the Virtual Environment for Reactor Applications (VERA). ROMUSE is written in C++ and is currently capable of performing various types of parameters perturbations, uncertainty quantification, surrogate models construction and subspace analysis. Version 2.0 has the capability to interface with DAKOTA which gives ROMUSE access to the various algorithms implemented within DAKOTA. ROMUSE is mainly designed to interface with VERA and the Comprehensive Modeling and Simulation Suite for Nuclear Safety Analysis and Design (SCALE) [1,2,3], however, ROMUSE can interface with any general model (e.g. python and matlab) with Input/Output (I/O) format that follows the Hierarchical Data Format 5 (HDF5). In this brief user manual, the use of ROMUSE will be overviewed and example problems will be presented and briefly discussed. The algorithms provided here range from algorithms inspired by those discussed in Ref.[4] to nuclear-specific algorithms discussed in Ref. [3].

97 MATHEMATICS AND COMPUTING↗

SERIIUS-MAGEEP Visiting Scholars Program

Recent studies have assessed closed-loop supercritical carbon dioxide (s-CO 2 ) Brayton cycles to be a higher energy-density system in comparison to equivalent superheated steam Rankine systems. At turbine inlet conditions of 700°C and 20 MPa, a cycle thermal efficiency of ~50% can be achieved. Achieving these high efficiencies will help concentrating solar power (CSP) technologies to become a competitive alternative to current power generation methods. To incorporate an s-CO 2 Brayton power cycle in a solar power tower system, the development of a solar receiver capable of providing an outlet temperature of 700°C (at 20 MPa) is necessary. To satisfy the temperature requirements of an s-CO 2 Brayton cycle with recuperation and recompression, the s-CO 2 must undergo a temperature rise of ~200°C as it flows through the solar receiver. The main objective is to develop an optical-thermal-fluid and structural model to validate a tubular receiver that will receive a heat input ~0.33 MWth from the heliostat field at the National Solar Thermal Test Facility (NSTTF), Albuquerque, NM, USA. We also commenced the development of computational models and testing of air receivers being developed by the Indian Institute of Science (IISc) and the Indian Institute of Technology in Bombay (IIT-B). The helical tubular receiver is expected to counteract the effect of thermal expansion while using a cavity to reduce the radiative and convective losses. Initially, this receiver will be tested for a temperature range of 100-300°C under 1 MPa of pressurized air. The helical air receiver will be exposed to 10kWth to achieve a temperature rise of ~200°C. Preliminary tests to validate the modeling will be performed before the design and construction of a larger scale receiver. Lastly, I focused on the development of a new computational tool that would allow us to perform a nodal creep-fatigue analysis on the receivers and heat exchangers being developed. This tool was developed using MATLAB and is capable of processing the results obtained from ANSYS Fluent and Structural combined, which was limited when using commercial software. The main advantage of this code is that it can be modified to run in parallel making it more affordable and faster compared to commercial codes available. The code is in the process of validation and is currently being compared to nCode Design Life.

25 ENERGY STORAGE↗

Echo

This document is a white paper marketing proposal for Echo™ is a data analysis platform designed for efficient, robust, and scalable creation and execution of complex workflows. Echo’s analysis management system refers to the ability to track, understand, and reproduce workflows used for arriving at results and decisions. Echo improves on traditional scripted data analysis in MATLAB, Python, R, and other languages to allow analysts to make better use of their time. Additionally, the Echo platform provides a powerful data management and curation solution allowing analysts to quickly find, access, and consume datasets. After two years of development and a first release in early 2016, Echo is now available for use with many data types in a wide range of application domains. Echo provides tools that allow users to focus on data analysis and decisions with confidence that results are reported accurately.

97 MATHEMATICS AND COMPUTING↗

Multi-compartment Aerosol Transport Model

A simple aerosol transport model was developed for a multi-compartmented cleanroom. Each compartment was treated as a well-mixed volume with ventilating supply and return air. Gravitational settling, intercompartment transport, and leakage of exterior air into the system were included in the model. A set of first order, coupled, ordinary differential equations was derived from the conservation equations of aerosol mass and air mass. The system of ODEs was then solved in MATLAB using pre-existing numerical methods. The model was verified against cases of (1) constant inlet-duct concentration, and (2) exponentially decaying inlet-duct concentration. Numerical methods resulted in normalized error of less than 10 -9 when model solutions were compared to analytical solutions. The model was validated against experimental measurements from a single field test and showed good agreement in the shape and magnitude of the aerosol concentration profile with time.

54 ENVIRONMENTAL SCIENCES↗

Novel Infiltration Diagnostics based on Laser-line Scanning and Infrared Temperature Field Imaging

This project targets the building energy efficiency problems induced by building infiltration/leaks. The current infiltration inspection techniques often require extensive visual inspection and/or whole building pressure test. These current techniques cannot meet more than three of the below five criteria of ideal infiltration diagnostics: 1. location and extent diagnostics, 2. building-level application, 3. least surface preparation, 4. weather-proof, and 5. non-disruption to building occupants. These techniques are either too expensive or time consuming, and often lack accuracy and repeatability. They are hardly applicable to facades/facades section. The goal of the project was to develop a novel infiltration diagnostics technology based on laser line-scanning and simultaneous infrared temperature imaging. A laboratory scale experimental setup was designed to mimic a model house of well-defined pressure difference below or above the outside pressure. Algorithms and Matlab-based programs had been developed for recognition of the hole location in infrared images. Our experiment based on laser wavelengths of 450 and 1550 nm and laser beam diameters of 4-25 mm showed that the location of the holes could be identified using laser heating; the diagnostic approach however could not readily distinguish between infiltration and non-infiltration points. To significantly improve the scanning throughput and recognition accuracy, a second approach was explored, developed, and extensively tested. It incorporates a liquid spray on the surface to induce extra phase change cooling effect. In this spray method, we termed it as PECIT (Phase-change Enhanced Cooling Infrared Thermography), phase-change enhanced cooling was used, which significantly amplifies the effect of air flow (infiltration and exfiltration). This heat transfer method worked extremely well to identify infiltration and exfiltration locations with high accuracy and increased throughput. The PECIT technique was systematically developed and tested for through holes with diameters 1 mm to 2 mm, and diagonal lines of 0.5 mm width at different camera-wall distances of 46 cm to 200 cm, under different pressure differences from 5 Pa to 20 Pa, and under different wind conditions. The PECIT technique had either met or exceeded the goals proposed in the project. For exfiltration, we achieved 100% accuracy under a much lower pressure difference of 10 Pa (proposed one: 50 Pa with stretch goal of 15 Pa). For infiltration, we achieved >90% accuracy under a much lower pressure difference of 10 Pa (proposed one: 50 Pa with stretch goal of 15Pa). For exfiltration, we achieved 100% accuracy under a much lower pressure difference of 10 Pa. For infiltration, we achieved 100% accuracy under a much lower pressure difference of 10 Pa. The PECIT technique can reach a throughput of 120 m2/h, which is 4 times the proposed goal for the laser line-scanning and simultaneous infrared temperature imaging approach. For commercialization and market penetration, we had meetings with two companies for feedback collection and further improvement for practical use. Also, we have interacted with Office of Intellectual Property and Technology Transfer of Iowa State University for idea disclosure and patent application.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Analysis of BSRT Profiles in the LHC at Injection

The beam synchrotron radiation telescope (BSRT) at the LHC allows to take profiles of the transverse beam distribution, which can provide useful additional insight in the evolution of the transverse beam distribution. A python class has been developed, which allows to read in the BSRT profiles, usually stored in binary format, run different analysis tools and generate plots of the statistical parameters and profiles as well as videos of the the profiles. The detailed analysis will be described in this note. The analysis is based on the data obtained at injection energy (450 GeV) during MD1217 and MD1415 which will be also used as illustrative example. A similar approach is also taken with a MATLAB based analysis described in.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Transport in Random Media With Inhomogeneous Mixing Statistics

The Levermore-Pomraning (LP) model has been investigated for use in implicit Monte Carlo (IMC) algorithms to describe thermal photon transport computations in stochastically mixed cells with material energy coupling. For purely absorbing non-emitting materials, the LP-model is known to predict the exact material averaged fluxes when the mixing statistics are Markovian. In the presence of scattering or reemmission, or if the mixing is non-Markovian, the LP-model is approximate and has greatly diminished accuracy for optically thick, highly scattering materials. Nevertheless, the model enables a robust extension of the simple atomic mix approximation to accommodate spatial correlation effects through well characterized parameters such as mean chord lengths. In this work, we construct exact analytical solutions of the LP-equations for piecewise constant spatial variation of the material chord lengths. It is not possible in general to obtain closed form solutions when the chord lengths are continuous functions of space, although we briefly examine special conditions under which this is feasible. However, by allowing an arbitrarily fine spatial segmentation of the medium it is possible to construct highly refined solutions to fit desired spatial chord length profiles. Material volume fractions that are consistent with the chord lengths are first constructed by solving the associated Chapmann-Kolmogorov equations and applied to obtain the material averaged fluxes in two settings: when the materials are purely absorbing, and when scattering is present in rod model geometry. Illustrative numerical results for the volume fractions and material averaged fluxes are obtained by implementing the algorithms in MATLAB for spatially linear and quadratic variations of the chord lengths.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

System Advisor Model (SAM) General Description (Version 2017.9.5)

This document describes the capabilities of the System Advisor Model (SAM) developed and distributed by the U.S. Department of Energy's National Renewable Energy Laboratory. The document is for potential users and others wanting to learn about the model's capabilities. SAM is a techno-economic computer model that calculates performance and financial metrics of renewable energy projects. Project developers, policy makers, equipment manufacturers, and researchers use graphs and tables of SAM results in the process of evaluating financial, technology, and incentive options for renewable energy projects. SAM simulates the performance of photovoltaic, concentrating solar power, solar water heating, wind, geothermal, biomass, and conventional power systems. The financial models are for projects that either buy and sell electricity at retail rates (residential and commercial) or sell electricity at a price determined in a power purchase agreement (PPA). SAM's simulation tools facilitate parametric and sensitivity analyses, Monte Carlo simulation and weather variability (P50/P90) studies. SAM can also read input variables from Microsoft Excel worksheets. For software developers, the SAM software development kit (SDK) makes it possible to use SAM simulation modules in their applications written in C/C plus plus, C sharp, Java, Python, MATLAB, and other languages. NREL provides both SAM and the SDK as free downloads at https://sam.nrel.gov. SAM is an open source project, so its source code is available to the public. Researchers can study the code to understand the model algorithms, and software programmers can contribute their own models and enhancements to the project. Technical support and more information about the software are available on the website.

47 OTHER INSTRUMENTATION↗

Lagrangian particle tracking applied to high-speed tomographic particle imaging velocimetry

We are investigating the accuracy of Lagrangian analysis of time-resolved, 3D velocity measurements in turbulent flames obtained using high-speed tomographic particle image velocimetry (TPIV). A particle-tracking algorithm was implemented in MATLAB to extract Lagrangian statistics from the velocity fields. The accuracy of Lagrangian statistics depends on the accuracy of the estimated particle trajectories. Errors associated with the spatial and temporal resolutions of high-speed TPIV experiments were evaluated using a direct numerical simulation (DNS) of an isotropic turbulent flow. LPT was then implemented onto TPIV measurements of a laminar counterflow flame. Results calculated with LPT were compared to known results of quantities measured with an Eulerian approach.

47 OTHER INSTRUMENTATION↗

Boil-off losses along LH2 pathway

Losses along the LH 2 pathway are intrinsic to the utilization of a cryogenic fluid. They occur when the molecule is transferred between 2 vessels (liquefaction plant to trailer, trailer to station storage, station storage to pump or compressor, then fuel cell electric vehicles …) and when the fluid is warmed up due to heat transfer with the environment. Those losses can be estimated with good accuracy using thermodynamic models based on conservation of mass and energy, providing the thermodynamic states are correctly described. Indeed, the fluid undergoes various changes as it moves along the entire pathway (2 phase transition, super-heated warming, non-uniform temperature distributions across the saturation film) and accurate equations of state and 2 phase behavior implementations are essential. The balances of mass and energy during the various dynamics processes then enable to quantify the boil-off losses. In this work, a MATLAB code previously developed by NASA to simulate rocket loading is used as the basis for the LH 2 transfer model. This code implements complex physical phenomena such as the competition between condensation and evaporation and the convection vs. conduction heat transfer as a function of the relative temperatures on both sides of the saturated film. The original code was modified to consider real gas equations of state, and some semi-empirical relationships, such as between the heat of vaporization and the critical temperature, were also replaced by a REFPROP equivalent expression, assumed to be more accurate. Non-constant liquid temperature equations were added to simulate sub-cooled conditions. It is shown that although trailer depressurization and transfer losses from the receiving vessel during a LH 2 delivery may constitute the largest sources of losses, those can be greatly reduced and even eliminated if CFR/DOT regulations are properly applied and if no-vent fill methods are implemented, particularly for a stationary LH2 storage operating at reasonably low pressure (45 psia or so). As such, it is expected that the only remaining boil-off losses for a refueling station would come from the LH 2 pump (utilization and idling), pump vessel cool down/warm up, and environment heat transfer. Based on experimental data measured at LLNL on the Linde 875 bar LH 2 cryo-pump and results from the model, and extrapolating for refueling stations of various sizes, it can be shown that boil-off losses can vary from 15% of delivered LH 2 for a 100 kg/day, 5% at 400 kg/day, and down to less than 2% for stations above 1,800 kg/day (for a system similar to the one at LLNL). Less boil-off is to be expected for LH 2 pumps dispensing at 350 bar, so that less than 0.7% can be expected above 1,800 kg/day station capacities. At last, boil-off recovery solutions are briefly analyzed, including compressors, cryo-coolers and stationary fuel cell. Their economic value is strongly related to the costs of H2 and of the electricity needed to perform the recovery, on top of the initial capital expenditures. Effective boil-off could be reduced to 0.6-1% for a 2,000 kg/day station if adequate recovery solutions are used.

08 HYDROGEN↗

Distribution System V1G PEV Charging Impacts Report

Plug-in Electric Vehicles (PEVs) are a clean form of transportation that will require a new control paradigm between the energy supply and the transportation sectors as PEV adoption continues to grow. Meeting the projected PEV energy demands will require implementing PEV / grid integration methods beyond charging PEVs at incentivized times. MATLAB was used to simulate the daily vehicle trip demands and charging processes of hundreds of PEVs. The National Household Travel Survey (NHTS) data was used as the source of vehicle travel data, and EV Project provided measured PEV charging data to calibrate the PEV population models. GridLAB-D analyses used a prototypical feeder and PEV population models to evaluate residential PEV V1G (unidirectional) charging in multiple locations, V1G PEV distribution feeder impacts, Time-Of-Use (TOU) rate effects on PEV charging, and charging control methods. The key findings were: optimized charging using California TOU rates could result in savings of over $20 per month for non-TOU program participants; non-TOU participants consume 30-40% of PEV charging energy; about 50% of PEV drivers have their transportation fuel bill reduced by 50 dollars per month; 6.6kW charging will cause distribution transformer power limits to be exceeded on prototypical feeders if the household PEV adoption rate reaches 0.75 PEVs per home; distribution feeder loading is proportional to PEV charging rate; a method was devised to remotely identify when residential transformer overloads occur using changes in PEV charging rate; and the low-cost TOU period is composed of a distribution of PEVs needing small, medium and larger energy needs that can be more effectively shifted as groups to mitigate the TOU peak power and enable a higher PEV adoption rate per household without overloading distribution transformers.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Co-Decontamination Dynamic Modeling to Support the Experimental Campaign

The Co-Decontamination (CoDCon) Demonstration experiment at Pacific Northwest National Laboratory (PNNL) is designed to test the separation of a mixed U and Pu product from dissolved spent nuclear fuel. The primary purpose of the project is to demonstrate control of the Pu/U ratio throughout the entire process without producing a pure Pu stream. In addition, the project is quantifying the accuracy and precision to which a Pu/U mass ratio can be achieved. The system includes an on-line monitoring system using spectroscopy to monitor the ratios throughout the process. A dynamic model of the CoDCon flowsheet and the on-line monitoring system was developed to augment the experimental work. This model is based in MATLAB Simulink and provides the ability to expand the range of scenarios that can be examined for process control and determine overall measurement uncertainty. Experimental results have been used to inform and benchmark the model so that it can accurately simulate various transient scenarios. The results of the experimental benchmarking are presented here along with modeled scenarios to demonstrate the control and process monitoring of the system.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗