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 685 records · Page 38

Thermal Storage System for Electric Vehicle Cabin Heating Component and System Analysis

Cabin heating of current electric vehicle (EV) designs is typically provided using electrical energy from the traction battery, since waste heat is not available from an engine as in the case of a conventional automobile. In very cold climatic conditions, the power required for space heating of an EV can be of a similar magnitude to that required for propulsion of the vehicle. As a result, its driving range can be reduced very significantly during the winter season, which limits consumer acceptance of EVs and results in increased battery costs to achieve a minimum range while ensuring comfort to the EV driver. To minimize the range penalty associated with EV cabin heating, a novel climate control system that includes thermal energy storage from an advanced phase change material (PCM) has been designed for use in EVs and plug-in hybrid electric vehicles (PHEVs). The present paper focuses on the modeling and analysis of this electrical PCM-Assisted Thermal Heating System (ePATHS) and is a companion to the paper Design and Testing of a Thermal Storage System for Electric Vehicle Cabin Heating. A detailed heat transfer model was developed to simulate the PCM heat exchanger that is at the heart of the ePATHS and was subsequently used to analyze and optimize its design. The results from this analysis were integrated into a MATLAB Simulink system model to simulate the fluid flow, pressure drop and heat transfer in all components of the ePATHS. The system model was then used to predict the performance of the climate control system in the vehicle and to evaluate control strategies needed to achieve the desired temperature control in the cabin. The analysis performed to design the ePATHS is described in detail and the system s predicted performance in a vehicle HVAC system is presented.

cabin heating↗

Cortexsys v. 3.0

This code implements a variety of "deep learning" algorithms from openly published academic journals. These deep learning algorithms allow the user to train neural networks to do various tasks. Such tasks include predicting future values in a time sequence, categorizing images or compressing information. The code is mostly written in the easy to understand Matlab / GNU Octave language, which enables rapid prototyping and understanding for research and educational purposes.

Cox, Jonathan A.↗

SWinzip v. 1.0

SWinzip is a Matlab and C++ library for scientific lossy data compression and reconstruction using compressed sensing and tree-wavelets transforms. These methods are known for their large compression and usefulness in data analytics such as features extraction. Compressed sensing and wavelets methods rely heavily on sparse and dense linear algebra operations implemented through the Boost codes in our library. SWinzip accommodates data represented on both regular grids (e.g. image data) and point-clouds (e.g. unstructured meshes).

Salloum, Maher↗

Super Strypi HWIL 6DOF (Hardware-In-Loop six-degree-of-freedom) Rev. 2175

The Super Strypi HWIL is a six degree-of-freedom (6DOF) simulation for the Super Strypi Launch Vehicle. The simulation is used to test the NGC flight software including the navigation software. Aerodynamic and propulsive forces, mass properties, ACS (attitude control system) parameters are defined in input files. Output parameters are saved to a Matlab mat file.

Gilkey, Jeff C.↗

First indirect x-ray imaging tests with an 88-mm diameter single crystal

Using the 1-BM-C beamline at the Advanced Photon Source (APS), we have performed the initial indirect x - ray imaging point-spread-function (PSF) test of a unique 88-mm diameter YAG:Ce single crystal of only 100 - micron thickness. The crystal was bonded to a fiber optic plat e (FOP) for mechanical support and to allow the option for FO coupling to a large format camera. This configuration resolution was compared to that of self - supported 25-mm diameter crystals, with and without an Al reflective coating. An upstream monochromator was used to select 17-keV x-rays from the broadband APS bending magnet source of synchrotron radiation. The upstream , adjustable Mo collimators were then used to provide a series of x-ray source transverse sizes from 200 microns down to about 15-20 microns (FWHM) at the crystal surface. The emitted scintillator radiation was in this case lens coupled to the ANDOR Neo sCMOS camera, and the indirect x-ray images were processed offline by a MATLAB - based image processing program. Based on single Gaussian peak fits to the x-ray image projected profiles, we observed a 10.5 micron PSF. This sample thus exhibited superior spatial resolution to standard P43 polycrystalline phosphors of the same thickness which would have about a 100-micron PSF. Lastly, this single crystal resolution combined with the 88-mm diameter makes it a candidate to support future x-ray diffraction or wafer topography experiments.

43 PARTICLE ACCELERATORS↗

General Mode Scanning Probe Microscopy

A critical part of SPM measurements is the information transfer from the probe-sample junction to the measurement system. Current information transfer methods heavily compress the information-rich data stream by averaging the data over a time interval, or via heterodyne detection approaches such as lock-in amplifiers and phase-locked loops. As a consequence, highly valuable information at the sub-microsecond time scales or information from frequencies outside the measurement band is lost. We have developed a fundamentally new approach called General Mode (G-mode), where we can capture the complete information stream from the detectors in the microscope. The availability of the complete information allows the microscope operator to analyze the data via information-theory analysis or comprehensive physical models. Furthermore, the complete data stream enables advanced data-driven filtering algorithms, multi-resolution imaging, ultrafast spectroscropic imaging, spatial mapping of multidimensional variability in material properties, etc. Though we applied this approach to scanning probe microscopy, the general philosophy of G-mode can be applied to many other modes of microscopy. G-mode data is captured by completely custom software written in LabVIEW and Matlab. The software generates the waveforms to electrically, thermally, or mechanically excite the SPM probe. It handles real-time communications with the microscope software for operations such as moving the SPM probe position and also controls other instrumentation hardware. The software also controls multiple variants of high-speed data acquisition cards to excite the SPM probe with the excitation waveform and simultaneously measure multiple channels of information from the microscope detectors at sampling rates of 1-100 MHz. The software also saves the raw data to the computer and allows the microscope operator to visualize processed or filtered data during the experiment. The software performs all these features while offering a user-friendly interface.

Somnath, Suhas↗

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

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

14 SOLAR ENERGY↗

Study of a Novel Oscillating Surge Wave Energy Converter: Preprint

This study investigates the performance of an oscillating surge wave energy converter (OSWEC) that utilizes adjustable geometry as a means of controlling the hydrodynamic coefficients, a concept originally proposed by [1]. The body of the device consists of a bottom-hinged solid rectangular frame with five horizontal flaps spanning the interior of the frame. The flaps can rotate independently about their center of rotation within the frame like a large window shutter. Changing the orientation of the flaps alters the hydrodynamic coefficients and natural frequency of the device as well as the ability to shed or absorb structural loads accordingly. This ability may allow the device to operate in a wider range of sea states than other current wave energy converter designs. This paper presents and compares the results of numerical simulations and experimental testing of the OSWEC's response to regular waves with all five of the horizontal fin configurations sharing the same orientation of 0 degrees (fully closed interior) and 90 degrees (fully open). The numerical simulations were performed using WAMIT, which calculates hydrodynamic coefficients using a boundary element method code to solve the linear potential flow problem, and WEC-Sim, a MATLAB-based tool that simulates multibody devices in the time domain by solving the governing equations of motion. A 1:14 scale model of the device was built for experimental evaluation in an 8-m-long, 1-m wide wave tank, which supports a water depth of 0.7 m. The OSWEC motion in different wave conditions was measured with displacement sensors while nonlinear wave-structure interaction effects like slamming and overtopping were captured using a high-speed camera and used to understand differences between the simulation and experiments.

13 HYDRO ENERGY↗

right-sized dimple evaluator

The code (aka computer program written as a Matlab script) uses a unique set of n independent equations to solve for n turbulence variables. The code requires the input of a characteristic dimension, a characteristic fluid velocity, the fluid dynamic viscosity, and the fluid density. Most importantly, the code estimates the size of three key turbulent eddies: Kolmogorov, Taylor, and integral. Based on the eddy sizes, dimples dimensions are prescribed such that the key eddies (principally Taylor, and sometimes Kolmogorov), can be generated by the dimple rim and flow unimpeded through the dimple’s concave cavity. It is hypothesized that turbulent eddies are generated by the dimple rim at the dimple-surface interface. The newly-generated eddies in turn entrain the movement of surrounding regions of fluid, creating more mixing. The eddies also generate lift near the wall surrounding the dimple, as they accelerate and reduce pressure in the regions near and at the dimple cavity, thereby minimizing the fluid drag.

Rodriguez, Sal↗

Beam Based RF Voltage Measurements and Longitudinal Beam Tomography at the Fermilab Booster

Increasing proton beam power on neutrino production targets is one of the major goals of the Fermilab long term accelerator programs. In this effort, the Fermilab 8 GeV Booster synchrotron plays a critical role for at least the next two decades. Therefore, understanding the Booster in great detail is important as we continue to improve its performance. For example, it is important to know accurately the available RF power in the Booster by carrying out beam-based measurements in order to specify the needed upgrades to the Booster RF system. Since the Booster magnetic field is changing continuously measuring/calibrating the RF voltage is not a trivial task. Here, we present a beam based method for the RF voltage measurements. Data analysis is carried out using computer programs developed in Python and MATLAB. The method presented here is applicable to any RCS which do not have flat-bottom and flat-top in the acceleration magnetic ramps. We have also carried out longitudinal beam tomography at injection and extraction energies with the data used for RF voltage measurements. Beam based RF voltage measurements and beam tomography were never done before for the Fermilab Booster. The results from these investigations will be very useful in future intensity upgrades.

43 PARTICLE ACCELERATORS↗

Stereo-Optic High Definition Imaging: A New Technology to Understand Bird and Bat Avoidance of Wind Turbines

There is a critical need to develop monitoring tools to track aerofauna (birds and bats) in three dimensions around wind turbines. New monitoring systems will reduce permitting uncertainty by increasing the understanding of how birds and bats are interacting with wind turbines, which will improve the accuracy of impact predictions. Biodiversity Research Institute (BRI), The University of Maine Orono School of Computing and Information Science (UMaine SCIS), HiDef Aerial Surveying Limited (HiDef), and SunEdison, Inc. (formerly First Wind) responded to this need by using stereo-optic cameras with near-infrared (nIR) technology to investigate new methods for documenting aerofauna behavior around wind turbines. The stereo-optic camera system used two synchronized high-definition video cameras with fisheye lenses and processing software that detected moving objects, which could be identified in post-processing. The stereo- optic imaging system offered the ability to extract 3-D position information from pairs of images captured from different viewpoints. Fisheye lenses allowed for a greater field of view, but required more complex image rectification to contend with fisheye distortion. The ability to obtain 3-D positions provided crucial data on the trajectory (speed and direction) of a target, which, when the technology is fully developed, will provide data on how animals are responding to and interacting with wind turbines. This project was focused on testing the performance of the camera system, improving video review processing time, advancing the 3-D tracking technology, and moving the system from Technology Readiness Level 4 to 5. To achieve these objectives, we determined the size and distance at which aerofauna (particularly eagles) could be detected and identified, created efficient data management systems, improved the video post-processing viewer, and attempted refinement of 3-D modeling with respect to fisheye lenses. The 29-megapixel camera system successfully captured 16,173 five-minute video segments in the field. During nighttime field trials using nIR, we found that bat-sized objects could not be detected more than 60 m from the camera system. This led to a decision to focus research efforts exclusively on daytime monitoring and to redirect resources towards improving the video post- processing viewer. We redesigned the bird event post-processing viewer, which substantially decreased the review time necessary to detect and identify flying objects. During daytime field trials, we determine that eagles could be detected up to 500 m away using the fisheye wide-angle lenses, and eagle-sized targets could be identified to species within 350 m of the camera system. We used distance sampling survey methods to describe the probability of detecting and identifying eagles and other aerofauna as a function of distance from the system. The previously developed 3-D algorithm for object isolation and tracking was tested, but the image rectification (flattening) required to obtain accurate distance measurements with fish-eye lenses was determined to be insufficient for distant eagles. We used MATLAB and OpenCV to improve fisheye lens rectification towards the center of the image, but accurate measurements towards the image corners could not be achieved. We believe that changing the fisheye lens to rectilinear lens would greatly improve position estimation, but doing so would result in a decrease in viewing angle and depth of field. Finally, we generated simplified shape profiles of birds to look for similarities between unknown animals and known species. With further development, this method could provide a mechanism for filtering large numbers of shapes to reduce data storage and processing. These advancements further refined the camera system and brought this new technology closer to market. Once commercialized, the stereo-optic camera system technology could be used to: a) research how different species interact with wind turbines in order to refine collision risk models and inform mitigation solutions; and b) monitor aerofauna interactions with terrestrial and offshore wind farms replacing costly human observers and allowing for long-term monitoring in the offshore environment. The camera system will provide developers and regulators with data on the risk that wind turbines present to aerofauna, which will reduce uncertainty in the environmental permitting process.

Adams, Evan↗

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↗