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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Advanced Small Modular Reactor Economics Model Development

The US Department of Energy Office of Nuclear Energy’s Advanced Small Modular Reactor (SMR) research and development activities focus on four key areas: Developing assessment methods for evaluating advanced SMR technologies and characteristics; and Developing and testing of materials, fuels and fabrication techniques; and Resolving key regulatory issues identified by US Nuclear Regulatory Commission and industry; and Developing advanced instrumentation and controls and human-machine interfaces. This report focuses on development of assessment methods to evaluate advanced SMR technologies and characteristics. Specifically, this report describes the expansion and application of the economic modeling effort at Oak Ridge National Laboratory. Analysis of the current modeling methods shows that one of the primary concerns for the modeling effort is the handling of uncertainty in cost estimates. Monte Carlo–based methods are commonly used to handle uncertainty, especially when implemented by a stand-alone script within a program such as Python or MATLAB. However, a script-based model requires each potential user to have access to a compiler and an executable capable of handling the script. Making the model accessible to multiple independent analysts is best accomplished by implementing the model in a common computing tool such as Microsoft Excel. Excel is readily available and accessible to most system analysts, but it is not designed for straightforward implementation of a Monte Carlo–based method. Using a Monte Carlo algorithm requires in-spreadsheet scripting and statistical analyses or the use of add-ons such as Crystal Ball. An alternative method uses propagation of error calculations in the existing Excel-based system to estimate system cost uncertainty. This method has the advantage of using Microsoft Excel as is, but it requires the use of simplifying assumptions. These assumptions do not necessarily bring into question the analytical results. In fact, the analysis shows that the propagation of error method introduces essentially negligible error, especially when compared to the uncertainty associated with some of the estimates themselves. The results of these uncertainty analyses generally quantify and identify the sources of uncertainty in the overall cost estimation. The obvious generalization—that capital cost uncertainty is the main driver—can be shown to be an accurate generalization for the current state of reactor cost analysis. However, the detailed analysis on a component-by-component basis helps to demonstrate which components would benefit most from research and development to decrease the uncertainty, as well as which components would benefit from research and development to decrease the absolute cost.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2013

Energy-Efficient and Comfortable Buildings through Multivariate Integrated Control (ECoMIC)

This project aims to develop an integrated control solution for enhanced energy efficiency and user comfort in commercial buildings. The developed technology is a zone-based control framework that minimizes energy usage while maintaining occupants’ visual and thermal comfort through control of electric lights, motorized venetian blinds and thermostats. The control framework is designed following a modular, scalable and flexible architecture to facilitate easy integration with exiting building management systems. The control framework contains two key algorithms: 1) the lighting load balancing algorithm and 2) the thermostat control algorithm. The lighting load balancing algorithm adopts a model-based closed-loop control approach to determine the optimal electric light and venetian blind settings. It is formulated into an optimization problem with minimizing lighting-related energy consumptions as the objective and delivering adequate task light and preventing daylight glare as the constraints. The thermostat control algorithm is based on a well-established thermal comfort model and formulated as a root-finding problem to dynamically determine the optimal thermostat setpoint for both energy savings and improved thermal comfort. To address building-wide scalability, a system architecture was developed for the zone-based control technology. Three levels of services are defined in the architecture: external services, facility level services and zone level services. The zone-level service includes the control algorithms described above as well as the corresponding interfaces, profiles, sensors and actuators to realize the zone controller. The facility level services connect to the zones through a backbone network, handle supervisory level information and controls, and thus facilitate building-wide scalability. The external services provide communication capability to entities outside of the building for grid interaction and remote access. Various aspects of the developed control technology were evaluated and verified through both simulations and testbed implementations. Simulations coupling a DOE medium office reference building in EnergyPlus building simulation software and a prototype controller in Matlab were performed. During summer time in a mixed-humid climate zone, the simulations revealed reductions of 27% and 42% in electric lighting load and cooling load, respectively, when compared to an advanced base case with daylight dimming and blinds automatically tilted to block direct sun. Two single-room testbeds were established. The testbed at Philips Lighting business building (Rosemont, IL) was designed for quantifying energy performance of integrated controls. This particular implementation achieved 40% and 79% savings on lighting and HVAC energy, respectively, compared to a relatively simple base case operated on predefined schedules. While the resulting energy savings was very encouraging, it should be noted that there may be several caveats associated with it. 1) The test was run during late spring and early summer, and the savings numbers might not be directly used to extrapolate the annual energy savings. 2) Due to the needs for separate control and metering of the small-scale demonstrator within a large building, the HVAC system, hence the corresponding savings, did not represent a typical energy code-compliant design. 3) The light level in the control case was regulated at a particular setpoint, which was lower than then the full-on light level in the base case, and the savings resulted from tuning down the light level to the setpoint was not attributable to the contribution of the developed technology. The testbed at the Lawrence Berkeley National Laboratory (Berkeley, CA) specifically focused on glare control integration, and has demonstrated the feasibility and capability of the glare detection and prevention technique. While the short one-month test in this testbed provided a functional indication of the developed technology, and it would require at least a full solstice-to-solstice cycle to ruinously quantify the performance, which was not possible within the project timeframe. There are certain limitations inherited from the operational assumptions, which could potentially affect the effectiveness and applicability of the developed control technologies. The system takes a typical ceiling-mounting approach for the photosensor locations, and therefore, the control performance relies on proper commissioning or the built-in intelligence of the photosensor for pertinent task light level estimations. For spaces where daylight penetration diminishes significantly deeper into the zone, certain modification to the control algorithms is required to accommodate multiple lighting control subzones and the corresponding sensors for providing a more uniform light level across the entire zone. Integrated control of visual and thermal comfort requires the lighting control zone and thermal control zone to coincide with each other. In other words, the area illuminated by a lighting circuit needs to be the same area served by the thermostat. Thus, the original zoning will potentially constrain the applicability of this technology in retrofitting projects. This project demonstrated the technical feasibility of a zone-based integrated control technology. From the simulation results and testbed implementations, up to 60% lighting energy savings in daylit areas relative to a “no-controls” case can easily be achieved. A 20% reduction of whole building energy consumption is also attainable. In the aspect of occupant comfort, the testbed demonstrated the ability to maintain specified light level on the workplane while promptly mitigate daylight glare 90% of the time. The control system also managed to maintain the thermal environment at a comfortable level 90% of the time. The aspect of system scalability was guaranteed by the system architecture design, based on which the testbeds were instantiated. Analysis on the aspect of economic benefit has yielded an about 6-year payback time for a medium-sized building, including the installation of all hardware and software, such as motorized blinds and LED luminaires. The payback time can be significantly reduced if part of the hardware is already in place for retrofitting projects. It needs to be noted that since the payback analysis was partly based on the testbed performance results, it is constrained by the caveats associated with the testbed implementations. The main uncertainty lies in the contribution from the space conditioning energy savings as it was non-trivial to realistically configure a room-size HVAC system for directly extrapolating whole-building HVAC energy savings. It is recommended to further evaluate the developed technology at a larger scale, where the lighting and HVAC energy consumption can be realistically measured at the building level, to more rigorously quantify the performance potentials.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2014

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↗
DOE OSTI2015

MueMex User's Guide

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

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

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↗
DOE OSTI2015

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↗
DOE OSTI2016

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↗
DOE OSTI2016

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↗
DOE OSTI2016

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↗
DOE OSTI2014

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↗
DOE OSTI2017

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↗
DOE OSTI2017

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↗
DOE OSTI2017

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↗
DOE OSTI2017

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↗
DOE OSTI2018

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↗
DOE OSTI2018

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↗
DOE OSTI2016

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↗
DOE OSTI2018

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↗