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 217 records · Page 12

In-Situ Blade Strain Measurements of a Crossflow Turbine Operating in a Tidal Flow

This data was collected between October 25 and December 12 of 2022 at the University of New Hampshire (UNH) and Atlantic Marine Energy Center (AMEC) turbine deployment platform (TDP). The goal was to collect blade strain data from a crossflow turbine operating in a tidal flow. A table in ('Deployment Schedule.PNG') outlines the various dates when each instrument was operational, and more details can be found via literature listed in 'Related Publications'.txt. This dataset includes zipped folders for each instrument containing data in .csv files for the relevant duration specific to each instrument, along with separate README file for each measurement. Some instrument files are quite large and can pose a challenge for a visual spreadsheet editor to open. A processing software like MATLAB or Python is recommended. All data contained in this submission is unfiltered/unprocessed data unless otherwise noted in the README file. Blade strain was measured using 8 foil-based strain gauges along the span of a single turbine blade. Water currents were measured using Acoustic Doppler Current Profilers (ADCP's) and Acoustic Doppler Velocimeters (ADV's) both upstream and downstream of the turbine for inflow, wake and turbulence measurements. Electrical power output was measured using the Voltsys rectifier. Shaft speed was calculated based on the Voltsys measurements of the permanent magnet three phase generator AC generation frequency, coupled directly to the cross flow turbine under test (i.e., no gear box). Platform motions were captured using a Yost IMU (inertial measurement unit). Turbine thrust loading was measured using a reaction arm about the turbine deployment platform spanning beam, where two bi-directional load cells were connected to the system via a pinned connection. The TDP is a floating structure moored on the Portsmouth facing side of Memorial Bridge pier #2, which spans the Piscataqua River between Portsmouth, NH and Kittery, ME. The Piscataqua River connects the Great Bay Estuary to the Gulf of Maine and the river currents are dominated by tidal flow with water velocities exceeding 2.5 m/s during spring ebb tides at this site which were previously characterized by Chancey 2019. The turbine under test was a modified New Energy Corporation (Calgary, CA) model EVG-025 4-blade H-Darrius type vertical axis cross flow turbine that rotates in the clockwise direction with a rotor diameter of 3.2m and blade length of 1.7m. The hydro-foil profile was a NACA 0021 with a 10 inch chord length and a blade preset pitch angle of +4deg with a positive angle corresponding with the toe in direction. The standard EVG-025 has a rotor diameter of 3.4m and its rated power output is 25kW at 3 m/s. The rotor diameter was reduced to accommodate the size of the existing TDP moon-pool. A single blade of this turbine was further modified to accommodate 8 full-bridge strain gauges (Bharath et al 2023, Bichanich et al 2024). For power performance and other relevant details on the turbine and its characteristics, see O'Byrne 2022.

16 TIDAL AND WAVE POWER↗

TEAMER: Triton Systems Oscillating Water Column Modeling Data and Report

This dataset provides the output of six Wave Energy Converter Simulator (WEC-Sim) simulations and accompanying documentation for the modeling of Triton Systems' oscillating water column (OWC) system at tank scale (validated using available data for tuning the model, Tests 1-2) and deployment scale (for which no validation data is available, Tests 4-6). Included are the output data in a MATLAB file structure, a comprehensive report on the modeling and design of the Triton OWC system, and a link to the WEC-Sim GitHub page. This work was supported by funding from TEAMER RFTS 5 (Request for Technical Support).

16 TIDAL AND WAVE POWER↗

DAISY Benchmark Performance Data

This repository contains the underlying data from benchmark experiments for Drifting Acoustic Instrumentation SYstems (DAISYs) in waves and currents described in "Performance of a Drifting Acoustic Instrumentation SYstem (DAISY) for Characterizing Radiated Noise from Marine Energy Converters" (https://link.springer.com/article/10.1007/s40722-024-00358-6). DAISYs consist of a surface expression connected to a hydrophone recording package by a tether. Both elements are instrumented to provide metadata (e.g., position, orientation, and depth). Information about how to build DAISYs is available at https://www.pmec.us/research-projects/daisy. The repository's primary content is three compressed archives (.zip format), each containing multiple MATLAB binary data files (.mat format). A table relating individual data files to figures in the paper, as well as the structure of each file, is included in the repository as a Word document (Data Description MHK-DR.docx). Most of the files contain time series information for a single DAISY deployment (file naming convention: [site]_DAISY_[Drift #].mat) consisting of processed hydrophone data and associated metadata. For a limited number of DAISY deployments, the hydrophone package was replaced with an acoustic Doppler velocimeter (file naming convention: [site]_DAISY_[Drift #]_ADV.mat). Data were collected over several years at three locations: (1) Sequim Bay at Pacific Northwest National Laboratory's Marine & Coastal Research Laboratory (MCRL) in Sequim, WA, the energetic tidal channel in Admiralty Inlet, WA (Admiralty Inlet), and the U.S. Navy's Wave Energy Test Site (WETS) in Kaneohe, HI. Brief descriptions of data files at each location follow. - MCRL - (1) Drift #4 and #16 contrast the performance of a DAISY and a reference hydrophone (icListen HF Reson), respectively, in the quiescent interior of Sequim Bay (September 2020). (2) Drift #152 and #153 are velocity measurements for a drifting acoustic Doppler velocimeter in in the tidally-energetic entrance channel inside a flow shield and exposed to the flow, respectively (January 2018). (3) Two non-standard files are also included: DAISY_data.mat corresponds to a subset of a DAISY drift over an Adaptable Monitoring Package (AMP) and AMP_data.mat corresponds to approximately co-temporal data for a stationary hydrophone on the AMP (February 2019). - Admiralty Inlet - (1) Drift #1-12 correspond to tests with flow shielded DAISYs, unshielded DAISYs, a reference hydrophone, and drifting acoustic Doppler velocimeter with 5, 10, and 15 m tether lengths between surface expression and hydrophone recording package (July 2022). (2) Drift #13-20 correspond to tests of flow shielded DAISYs with three different tether materials (rubber cord, nylon line, and faired nylon line) in lengths of 5, 10, and 15 m (July 2022). - WETS - (1) Drift #30-32 correspond to tests with a heave plate incorporated into the tether (standard configuration for wave sites), rubber cord only, and rubber cord, but with a flow shielded hydrophone (November 2022). (2) Drift #49-58 and Drift #65-68 correspond to measurements around mooring infrastructure at the 60 m berth where time-delay-of-arrival localization was demonstrated for different DAISY arrangements and hydrophone depths (November 2022).

16 TIDAL AND WAVE POWER↗

Drifting Acoustic Measurements around C-Power's SeaRay WEC

The repository contains underwater noise measurements and associated metadata collected around C-Power's SeaRay wave energy converter on July 15, 2024 and July 16, 2024 while it was deployed at the U.S. Navy's Wave Energy Test Site (WETS) in Kaneohe, HI. Measurements were obtained using Drifting Acoustic Instrumentation SYstems (DAISYs). DAISYs consist of a surface expression connected to a hydrophone recording package by a tether. Both elements are instrumented to provide metadata (e.g., position, orientation, and depth). Information about how to build DAISYs is available at https://www.pmec.us/research-projects/daisy. The repository's primary content is a compressed archive (.zip format), containing multiple MATLAB binary data files (.mat format). The structure of each file is included in the repository as a Word document (Data Description MHK-DR.docx). Each file contains time series information for a single DAISY deployment (file naming convention: WETS_DAISY_[Drift #].mat) consisting of processed hydrophone data and associated metadata. During these measurements, C-Power's SeaRay was located at approximately 21.48112 N, 157.74451 W.

16 TIDAL AND WAVE POWER↗

TEAMER: Twin Ocean Power Wave Energy Converter Comprehensive Overview

These files collectively provide a comprehensive overview of the testing process, data analysis, and validation for the Twin Ocean Power device tested at the O.H. Hinsdale Wave Research Laboratory, supported by TEAMER funding. This resource includes an overview of power results for a series of 7 trials. The files included in this comprehensive overview include a comprehensive log sheet for each trial, a summary of all trials, and processing scripts for the raw data. It includes all raw data in .tsv and MATLAB compatible formats, an average power chart, angular velocity charts for each trial, trial metrics, and power output files. This resource includes images of the Twin Ocean Power Wave Energy Converter device components and movement during testing and video recordings of each trial.

16 TIDAL AND WAVE POWER↗

1D Heat Loss Models Validation Experiment

Contains data from the model validation in the 1D Heat Loss Models to Predict the Aquifer Temperature Profile during Hot/Cold Water Injection Project. The data include two COMSOL models (2D axisymmetric benchmark model and 2D Vinsome model), one python code (1D Vinsome based FEM numerical simulation), one matlab main code (1D Newton analytical solution and all results comparison visualization), and output files generated from the above models.

1D↗

Simscape Modeling Verification in the Simulink Development Environment

The purpose of the Simulation Product Group of the Control and Data Systems division of the NASA Engineering branch at Kennedy Space Center is to provide a realtime model and simulation of the Ground Subsystems participating in vehicle launching activities. The simulation software is part of the Spaceport Command and Control System (SCCS) and is designed to support integrated launch operation software verification, and console operator training. Using Mathworks Simulink tools, modeling engineers currently build models from the custom-built blocks to accurately represent ground hardware. This is time consuming and costly due to required rigorous testing and peer reviews to be conducted for each custom-built block. Using Mathworks Simscape tools, modeling time can be reduced since there would be no custom-code developed. After careful research, the group came to the conclusion it is feasible to use Simscape's blocks in MatLab's Simulink. My project this fall was to verify the accuracy of the Crew Access Arm model developed using Simscape tools running in the Simulink development environment.

SIMSCAPE↗

The Development and Implementation of the Kennedy Space Center Umbilical Clearance Tool

In preparation for NASAs upcoming Space Launch System program, the Kennedy Space Center is currently developing subsystems to provide fuel, purges and communications to the flight vehicle, known as umbilicals. It is vital to the crew and mission that these umbilicals release at T-0 without re-contacting the vehicle as it is accelerating from the launch pad. To help ensure this requirement is met by the program, a methodology of evaluating the moving bodies was developed and implemented into a tool using MATLAB. The tool, known as the KSC Umbilical Clearance Tool, takes a given elevation of interest and an umbilical retract profile within the plane to evaluate the clearance between the umbilical arm and thousands of independent flight vehicle drift profiles from a Monte Carlo analysis. The presentation will delve into the challenges associated with developing and implementing the tool framed in the context of evaluating the clearance for one of the SLS umbilicals.

Clearance↗

Development and Application of a Tool for Optimizing Composite Matrix Viscoplastic Material Parameters

This document describes a recently developed analysis tool that enhances the resident capabilities of the Micromechanics Analysis Code with the Generalized Method of Cells (MAC/GMC) and its application. MAC/GMC is a composite material and laminate analysis software package developed at NASA Glenn Research Center. The primary focus of the current effort is to provide a graphical user interface (GUI) capability that helps users optimize highly nonlinear viscoplastic constitutive law parameters by fitting experimentally observed/measured stress-strain responses under various thermo-mechanical conditions for braided composites. The tool has been developed utilizing the MATrix LABoratory (MATLAB) (The Mathworks, Inc., Natick, MA) programming language. Illustrative examples shown are for a specific braided composite system wherein the matrix viscoplastic behavior is represented by a constitutive law described by seven parameters. The tool is general enough to fit any number of experimentally observed stress-strain responses of the material. The number of parameters to be optimized, as well as the importance given to each stress-strain response, are user choice. Three different optimization algorithms are included: (1) Optimization based on gradient method, (2) Genetic algorithm (GA) based optimization and (3) Particle Swarm Optimization (PSO). The user can mix and match the three algorithms. For example, one can start optimization with either 2 or 3 and then use the optimized solution to further fine tune with approach 1. The secondary focus of this paper is to demonstrate the application of this tool to optimize/calibrate parameters for a nonlinear viscoplastic matrix to predict stress-strain curves (for constituent and composite levels) at different rates, temperatures and/or loading conditions utilizing the Generalized Method of Cells. After preliminary validation of the tool through comparison with experimental results, a detailed virtual parametric study is presented wherein the combined effects of temperature and loading rate on the predicted response of a braided composite is investigated.

thermo-mechanical; MATLAB; GUI-Based Too↗

Swift Mission Gyro Patch and Re-Calibration Without A Calibration Campaign

The Swift spacecraft has three Two-Axis-Rate-Assemblies (TARAs) (designated as TARA 1, TARA 2, and TARA 3, and generically referred to as gyros) which until March of 2024 were all used in rate estimation as an operational workaround (dubbed “corrector-gyro”) of a software defect discovered in 2007. TARA 1 performance degraded throughout 2023 and 2024, making it an unreliable rate source. It was determined that shifting TARA 1 to be the corrector-gyro and using mainly TARAs 2 and 3, with the existing flight software (FSW) defect, would not meet performance requirements for science operations. Rather the FSW was patched to remove the defect found in 2007. This patch was delivered by the vendor in 2009, but was not fully tested or installed at that time. A decision was made in 2024 to install the patch despite limited simulation and testing capabilities. The problem facing the team was to come up with a set of alignment parameters that would work with the new patch, using data from the spacecraft operating normally in the old configuration. Calibration activities would have been possible in either the old or new FSW configuration, but were ultimately deemed unnecessary, saving staff time and pre-serving time on the spacecraft for science data collection. This paper describes the old and new FSW configurations related to gyro processing and alignment, and how the raw gyro and star tracker data from normal operation were used with a multidimensional unconstrained nonlinear minimization in MATLAB (Nelder-Mead1) to produce a fully calibrated set of gyro alignment parameters compatible with the patched FSW. Flight data is presented showing that not only did the process work, it improved performance significantly, rivaling the best slew performance of the mission to date.

MATLAB↗

Kraft black liquor concentration with graphene oxide membranes: Process simulations and technoeconomic analysis

Abstract Black liquor (BL) dewatering by multi‐effect evaporation in the kraft papermaking process is highly energy‐intensive. It was previously shown that graphene oxide (GO) nanofiltration membranes can remove lignin, other organics, and inorganic salts from BL while exhibiting stability in caustic BL conditions. Here, we design and simulate several candidate dewatering processes and evaluate their technoeconomic characteristics. All processes concentrate BL from 15 to 30 wt% solutes while producing aqueous permeate. Two process options were analyzed—option A including “last‐mile” permeate treatment to reduce solutes to 0.2 wt%, and option B excluding this treatment and producing a 3–4 wt% solutes stream. These processes were simulated in custom‐built ASPEN Plus flowsheets interfaced with Microsoft Excel and MATLAB. All processes deliver large (>40%) energy savings. Detailed technoeconomic analysis showed that option A processes are profitable in mills equipped with condensing turbines, but unprofitable with only purchased fuel savings. Option B processes are profitable in both situations, but require the caustic permeate to be utilized in other kraft process units. They are also profitable with electricity generation when operated at smaller scales matching the requirements of other process units. Monte‐Carlo sensitivity analysis shows that Option A can yield median 20‐year NPVs up to ~$10MM and Option B up to ~$25MM. Overall, GO membrane‐based BL dewatering is economically promising, assuming successful slipstream piloting and scale‐up campaigns. It would have immediate sustainability benefits from large energy savings, and broader implications for biorefinery processes due to the ability to fractionate biomass feedstock components under harsh conditions.

Wang, Zhongzhen↗

Open‐Source Anaerobic Digestion Modeling Platform, Anaerobic Digestion Model No. 1 Fast (ADM1F)

An open‐source modeling platform, called Anaerobic Digestion Model No. 1 Fast (ADM1F), is introduced to achieve fast and numerically stable simulations of anaerobic digestion processes. ADM1F is compatible with an iPython interface to facilitate model configuration, simulation, data analysis, and visualization. Faster simulations and more stable results are accomplished by implementing an advanced open‐source library of numerical methods called Portable Extensive Toolkit for Scientific Computation (PETSc) to solve the ADM1 system of equations. Leveraging PETSc, ADM1F can consistently complete a steady‐state simulation under 0.2 s, over 99% faster than a benchmark ADM1 model implemented with MATLAB while achieving agreement of model outputs within 1% of those obtained with the benchmark model. For dynamic simulations, however, ADM1F has a computational speed advantage only when the influent characteristics update more frequently than every 4 h. The ability of ADM1F to be useful as a tool to study anaerobic digestion systems is demonstrated through two example implementations of ADM1F: (1) a two‐phase co‐digestion scenario evaluating the impact of the organic loading rate and the substrate composition on reactor performance and stability, and (2) a conventional digester scenario assessing the effectiveness of recovery strategies after disruptions that led to instability. These examples demonstrate how the high simulation speed and the convenience of the iPython interface allow ADM1F to complete complex analyses within minutes, much faster than computational strategies currently reported in the literature.

anaerobic co-digestion↗

CP decomposition for tensors via alternating least squares with QR decomposition

The CP tensor decomposition is used in applications such as machine learning and signal processing to discover latent low-rank structure in multidimensional data. Computing a CP decomposition via an alternating least squares (ALS) method reduces the problem to several linear least squares problems. The standard way to solve these linear least squares subproblems is to use the normal equations, which inherit special tensor structure that can be exploited for computational efficiency. However, the normal equations are sensitive to numerical ill-conditioning, which can compromise the results of the decomposition. In this paper, we develop versions of the CP-ALS algorithm using the QR decomposition and the singular value decomposition, which are more numerically stable than the normal equations, to solve the linear least squares problems. Our algorithms utilize the tensor structure of the CP-ALS subproblems efficiently, have the same complexity as the standard CP-ALS algorithm when the input is dense and the rank is small, and are shown via examples to produce more stable results when ill-conditioning is present. Our MATLAB implementation achieves the same running time as the standard algorithm for small ranks, and we show that the new methods can obtain lower approximation error.

97 MATHEMATICS AND COMPUTING↗

Protection Scheme for Multiterminal HVDC System Based on Wavelet Transform Modulus Maxima

Fault protection is a major issue in the development of multiterminal high‐voltage direct current (HVDC) system based on voltage‐source converter. Owing to the rapid increase in fault current under DC faults, protection requires a fast fault detection method to isolate the fault current. However, the selectivity and reliability of the fast protection scheme should be primarily guaranteed. Therefore, a new protection strategy based on modulus maximum under wavelet transform is proposed in this paper. Fault detection and discrimination are achieved using single‐end signals of fault current. The external and internal faults can be discriminated via local information. The reliability and robustness of the proposed scheme against potential factors are verified. Numerous simulation cases are performed in PSCAD/EMTDC under different fault scenarios, and signal data are analyzed in MATLAB to demonstrate the effectiveness of the proposed fault protection scheme. © 2020 Institute of Electrical Engineers of Japan. Published by Wiley Periodicals LLC.

Huai, Qing↗

Deflection Shape Balancing: An Alternative to Modal Balancing

Rotor balancing is the final step in the rotor manufacturing process and may be necessary during rotor refurbishment. The process can be time consuming and costly, especially for supercritical monolithic rotors, impacting schedule and increasing costs. Modal balancing is a commonly used practice and is effective but requires multiple runs and at least one high-speed run. Deflection shape balancing is introduced as an alternative that may offer a more efficient method of balancing making the process faster with potentially fewer balancing runs. This paper explores the performance of the deflection shape fit balancing technique by utilizing MATLAB code simulations.

Morrell, William Jason↗

Response tailoring of elasto-plastic trusses

Abstract In this work we tailor the response of trusses loaded beyond their yield limit. The truss structures are modeled using finite strain theory and rate-independent elasto-plasticity. We design trusses with a tailored mechanical response that is between “elastic” and “elastic-ideal-plastic” subject to the volume constraint. The design updates are generated by the gradient-based Method of Moving Asymptotes (MMA) solver and the sensitivities of the response functions are computed using a path-dependent adjoint sensitivity analysis. The computations are performed in Matlab.

42 ENGINEERING↗

Initial conceptual demonstration of control co-design for WEC optimization

Abstract While some engineering fields have benefited from systematic design optimization studies, wave energy converters have yet to successfully incorporate such analyses into practical engineering workflows. The current iterative approach to wave energy converter design leads to sub-optimal solutions. This short paper presents an open-source MATLAB toolbox for performing design optimization studies on wave energy converters where power take-off behavior and realistic constraints can be easily included. This tool incorporates an adaptable control co-design approach, in that a constrained optimal controller is used to simulate device dynamics and populate an arbitrary objective function of the user’s choosing. A brief explanation of the tool’s structure and underlying theory is presented. To demonstrate the capabilities of the tool, verify its functionality, and begin to explore some basic wave energy converter design relationships, three conceptual case studies are presented. In particular, the importance of considering (and constraining) the magnitudes of device motion and forces in design optimization is shown.

16 TIDAL AND WAVE POWER↗