Search NASA⌕ Search

SEARCH · Search NASA

Results for “Matlab data”

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 199 records · Page 11

Data from: "Towards CONUS-Wide ML-Augmented Conceptually-Interpretable Modeling of Catchment-Scale Precipitation-Storage-Runoff Dynamics"

This data package was generated to support the manuscript “Towards CONUS-Wide Machine Learning-Augmented Conceptually Interpretable Modeling of Catchment-Scale Precipitation-Storage-Runoff Dynamics.” It provides input files, model outputs, plotting data, scripts, notebooks, and documentation used to develop, evaluate, and reproduce Mass-Conserving Perceptron (MCP)-based hydrologic modeling experiments across 513 selected Catchment Attributes and Meteorology for Large-sample Studies in the United States (CAMELS-US) basins. The files are organized by modeling component and analysis purpose, including rainfall–runoff experiments, snow module experiments, coupled hydrologic-snow experiments, Long Short-Term Memory (LSTM) benchmark results, model skill metrics, initialization and epoch records, cell-state normalization files, Akaike Information Criterion (AIC)-based model comparison files, and data used to generate manuscript figures. Tabular files can be opened using standard spreadsheet software or Python/R data-analysis tools. Python scripts, Jupyter notebooks, and selected MATLAB scripts are included for model execution, postprocessing, plotting, and statistical analysis. Quality assurance and quality control were conducted through the source-data selection and modeling workflow. Meteorological forcing, streamflow, and static catchment attributes were derived from the CAMELS-US dataset, and snow water equivalent data were derived from the University of Arizona (UA) Snow Water Equivalent dataset. Selected basins and time periods were screened during the associated research workflow to avoid missing observations or poor-quality cases. Static geospatial features were processed primarily using Quantum Geographic Information System (QGIS) and Geospatial Data Abstraction Library (GDAL) workflows. Additional details are provided in the associated manuscript and documentation.

ESS-DIVE CSV File Formatting Guidelines Reporting ↗

CalWave - TEAMER UMaine 2021 Tank Testing

This is the data submission to the teamer report for RFTS1 during which tank testing at UMaine W2 tank occurred. This submission includes the relevant information and recorded data of CalWave's UMaine wave tank testing conducted under TEAMER Request for Support Round 1. During the testing CalWave assessed an autonomous control approach for changing PTO gain parameters based on "on board" sea state detection. The data submitted includes an overview of normalized PTO gains based on sea state values and for a set of 8 relevant sea states the WEC motion. All data was collected on CalWave's SCADA system and UMaine's motion tracking system. The data submitted suffices the creation of figures as shown in the PDF report submitted alongside the data and the submission includes excel sheets (parameter gains), MATLAB structures (absorber motion), and the mentioned PDF report.

16 TIDAL AND WAVE POWER↗

JHTDB-wind: a web-accessible large-eddy simulation database of a wind farm with virtual sensor querying

This paper introduces JHTDB-wind (https://turbulence.idies.jhu.edu/datasets/windfarms, last access: 11 November 2025), a publicly accessible database containing large-eddy simulation (LES) data from wind farms. Building on the framework of the Johns Hopkins Turbulence Database (JHTDB), which hosts direct numerical simulation (DNS) and some LES datasets of canonical turbulent flows, JHTDB-wind stores the 4D space–time history of the flow and provides users the ability to access and query the data via a web-based virtual sensor interface. The initial dataset comprises LES results from a large wind farm with 10×6 turbines, modeled using a filtered actuator line method, under conventionally neutral atmospheric conditions. These data comprise 1 h (hour) of flow field data (velocity, pressure, potential temperature deviation, subgrid-scale (SGS) eddy viscosity, and turbine forces, approximately 15 TB (terabytes) and wind turbine data – including both turbine-level operational quantities and blade-level aerodynamic quantities (approximately 1.3 TB) – stored in Zarr and Parquet formats, respectively. Data retrieval is facilitated by the giverny Python package, allowing remote users to query the database in Python or MATLAB (C and Fortran support are available for flow field data). This paper details the simulation setup and demonstrates data access through examples that analyze wind farm flow structures and turbine performance. The framework is extensible to future datasets, including the JHTDB-wind diurnal cycle simulation analyzed in Xiao et al. (2025).

17 WIND ENERGY↗

ELM-Wet: Inclusion of a wetland landunit with sub-grid representation of eco-hydrological patches and hydrological forcing in E3SM Land Model (ELM)

Wetlands emit the most biogenic methane (CH4) and present the greatest uncertainty in the global CH4 budget. Modeling these emissions is challenging due to the temporal and spatial variability in wetland structure and CH4 flux rates, along with complex interactions among hydrological, ecological, meteorological, and microbial processes that govern CH4 dynamics. To address these issues, we aim to enhance the accuracy of wetland representation in the U.S. Department of Energy’s Exascale Earth System Model (E3SM) Land Model, ELM. This effort led to the development of ELM-Wet, which incorporates a dedicated wetland landunit with subgrid representation of eco-hydrological patch types. We implement wetland-specific hydrology by imposing site-specific constraints on surface water levels, thereby allowing different patches to sustain varying depths of inundation. Additionally, we refined the calculation of aerenchyma transport diffusivity based on observed conductance across different vegetation types. We validated these enhancements through site-specific simulations of a coastal freshwater wetland, the Salvador WMA Freshwater Marsh (Ameriflux, site ID US-LA2), located in the coast of Louisiana (29.85N,90.29W) at an elevation of 0 m. The site was simulated with ELM-Wet and the default version ELMv1. We use Bayesian Optimization to parameterize CO2 and CH4 fluxes. Eddy covariance observations of CO2 and CH4 fluxes from 2012-2013 were used to train the model and data from 2021 were used for validation. In this repository, we include all the output of all the simulations performed using different versions of ELMv1 and ELM-Wet, all input required to the model, and the Matlab scripts we used to processed the output data. ELM-Wet_flmd.csv includes a detailed description of the datasets files.

54 ENVIRONMENTAL SCIENCES↗

TEAMER: Water Tunnel Data from Testing the Pterofin Skimmer Concept

Pterofin's Skimmer concept relies on a flapping and pitching hydrofoil to extract hydrokinetic energy from water flows. The concept aims to utilize unsteady fluid dynamics phenomena (added mass, shed vorticity, and unsteady boundary layer development) to achieve higher lift coefficients, enabling increased power density of the hydrokinetic device and a fundamental shift in the rpm/torque scaling of the power take off compared with turbines. The Applied Research Laboratory at Penn State, in collaboration with Pterofin, designed and built a proof-of-concept flapping/pitching mechanism which was subsequently tested in ARL's 12-inch water tunnel facility. The mechanical power supplied to or extracted from the mechanism was measured for a range of hydrofoils provided by Pterofin over operating conditions including reduced frequency, Reynolds number, and the ratio between pitching and flapping amplitudes. The power lost to friction in the mechanism was removed from the net power measurement by means of a bare hub tare, with the resultant hydrodynamic power being used to calculate a mechanism-independent and non-dimensional power coefficient. The product of this effort is a dataset describing the power coefficient of a hydrofoil having simultaneous pitching and flapping motions, both of which are approximately sinusoidal. Power coefficients were collected for a range of primary design variables including: - Reduced frequency: 0.01 to 0.95 - Pitching/flapping peak angle ratio: 1.5 to 3.0 - Chord-based Reynolds number: 60,000 to 560,000 Secondary design variables relating to the hydrofoil geometry were explored including: - Aspect ratio - Planform shape - Section thickness distribution - Hydrofoil position relative to the pitching axis - Hydrofoil sweep angle relative to the pitching axis Measured data are provided in mean and time series formats. MATLAB scripts are provided which can be used to generate figures of time-averaged and phase-averaged hydrodynamic power coefficients calculated from the measured data. A complete description of the experiment and data reduction can be found in the Post Access Report for the Pterofin Skimmer test effort which will be available on the TEAMER website. This work was supported by the Pacific Energy Ocean Trust via a TEAMER award.

16 TIDAL AND WAVE POWER↗

Trick Simulation Environment 07

The Trick Simulation Environment is a generic simulation toolkit used for constructing and running simulations. This release includes a Monte Carlo analysis simulation framework and a data analysis package. It produces all auto documentation in XML. Also, the software is capable of inserting a malfunction at any point during the simulation. Trick 07 adds variable server output options and error messaging and is capable of using and manipulating wide characters for international support. Wide character strings are available as a fundamental type for variables processed by Trick. A Trick Monte Carlo simulation uses a statistically generated, or predetermined, set of inputs to iteratively drive the simulation. Also, there is a framework in place for optimization and solution finding where developers may iteratively modify the inputs per run based on some analysis of the outputs. The data analysis package is capable of reading data from external simulation packages such as MATLAB and Octave, as well as the common comma-separated values (CSV) format used by Excel, without the use of external converters. The file formats for MATLAB and Octave were obtained from their documentation sets, and Trick maintains generic file readers for each format. XML tags store the fields in the Trick header comments. For header files, XML tags for structures and enumerations, and the members within are stored in the auto documentation. For source code files, XML tags for each function and the calling arguments are stored in the auto documentation. When a simulation is built, a top level XML file, which includes all of the header and source code XML auto documentation files, is created in the simulation directory. Trick 07 provides an XML to TeX converter. The converter reads in header and source code XML documentation files and converts the data to TeX labels and tables suitable for inclusion in TeX documents. A malfunction insertion capability allows users to override the value of any simulation variable, or call a malfunction job, at any time during the simulation. Users may specify conditions, use the return value of a malfunction trigger job, or manually activate a malfunction. The malfunction action may consist of executing a block of input file statements in an action block, setting simulation variable values, call a malfunction job, or turn on/off simulation jobs.

Lin, Alexander S.↗

Connecting Users and Applications with Po.daac Hosted GHRSST Data

The 80+ GHRSST public datasets represent a rich resource for sea surface temperature research and applications given their time series length, resolution, spatial coverage, varying measurement types and processing levels, and availability in the full spectrum of PO.DAAC tools and services ecosystem. The PO.DAAC has created a publicly accessible recipe suite for the user community to perform straightforward yet powerful computations on GHRSST data using python recipes, Jupyter notebooks, R, Matlab, and the NCO programming language. These recipes include numerical computations for regional and global SST trends, anomaly derivations, EOF analysis, climate signal reproduction, and ocean phenology. For example, one recipe reproduces a famous SST based warming figure from the Fourth National Climate Assessment (USA) while another focuses on quantifying the regional changes in ocean SST phenology. Most are python-based while some contain hybrid calls and leverage the NCO programming interface too. All are available on the PO.DAAC user forum (https://podaac.jpl.nasa.gov/forum/) and/or via the open source NASA GitHub repository (https://github.com/nasa/podaac_tools_and_services). Several are available in the Jupyter notebook framework including podaacypy (https://github.com/nasa/podaacpy), a recipe for GHRSST granule metadata discovery and application, and more recently a Jupyter notebook developed to support data analysis and visualization of a cloud-based Zarr formatted Level 4 MUR dataset in the AWS Open Data Registry. Throughout the summer of 2020, the PO.DAAC intends to add and migrate more of its numerical recipes to the Jupyter notebook framework and publish them on its open source GitHub repository.

Gentemann, Chelle↗

TEAMER - AquaHarmonics High Fidelity WEC Sim PTO and Control Model Validation, Sim Model

Collaborative effort between AquaHarmonics, Sandia National Laboratories (SNL), and the National Renewable Energy Laboratory (NREL) to revise and validate Aquaharmonics' full wave to wire model, allowing for reduced uncertainty and increased understanding of design requirements of a utility scale wave energy converter (WEC). SNL and NREL in collaboration with AquaHarmonics, will set up and run WEC Simulator (WEC-Sim) models of the AquaHarmonics WEC, building off past model developments for inclusion of custom PTO (power take-off) dynamics. The intent is to review, update, and verify or validate a new WEC-Sim model against wave tank experimental data. Furthermore, the WEC-Sim model will be coupled to an energy storage system model to better understand the wave-to-wire functionality. Project is part of the TEAMER (Testing & Expertise for Marine Energy) RFTS (request for technical support) system of WEC research projects. Testing data can be found in the associated MHKDR link below.

16 TIDAL AND WAVE POWER↗

Zircon, Baddeleyite, and Reidite Found in Ries Crater Suevite

Introduction: Impact events can generate superheated impact melts and even vapor [1]. In the past decade, evidence of the high temperatures and high pressures of the impact process has been found in the impact melt from terrestrial craters (e.g., Mistastin Lake [2], Meteor Crater [3]) and even on the Moon [4–5]. Often, these studies involve zircon or zirconium-bearing phases. Zircon is a particularly useful mineral, due to its robustness and durability against weathering. It is used principally for chronology, but it has a multitude of geologic applications, including geothermometry and fingerprinting magma sources. Here we focus on its capability of recording impact conditions. Grains of zircon (ZrSiO4) are converted at high temperature and/or pressure during an impact event to reidite (a high-pressure polymorph of ZrSiO4) or to tetragonal-ZrO2 + SiO2 [6]. Reidite and tertragonal-ZrO2 leave identifying relicts in the rocks, markers of the high pressures and/or temperatures these rocks underwent. Ries Crater is a 26-km-diameter peak ring crater in southern Germany that formed approximately 15 Ma [7–9]. Ries is the type locality of the polymict impact breccia known as suevite. Here, we report a microanalytical study of suevite from the Ries impact structure. These analyses are used to inform our understanding of the pressure and temperature conditions involved in creating a polymict impact breccia such as suevite in a peak ring crater. Sample Description: The sample, denoted 16RS08, originates from Otting Quarry, at 48.8777° N, 10.7921° E, approximately 17 km from the center of the Ries Crater in Germany (approximately 4 km outside the crater rim). The material recorded in sample 16RS08 is therefore considered an outer suevite. The outer suevite is a discontinuous layer of polymict impact breccia that occurs outside the central ring of Ries, up to 22 km from the center of the crater [9]. Lithic clasts in the outer suevite consist primarily of crystalline basement rocks (gneiss, granite, amphibolite), with less than 5% of lithic clasts being overlying sedimentary rocks (limestone, sandstone, shale) [9]. The thin section studied contains variably shocked lithic and mineral clasts, impact glass, and interstitial minerals that make up the matrix of the breccia (Fig. 1) [9]. Methods: We used a Cameca SX100 electron probe microanalyzer (EPMA) located in the Kuiper Materials Imaging and Characterization Facility (KMICF) at the University of Arizona to obtain 15 elemental X-ray maps of 16RS08. We next used the JEOL 7900F SEM at the Astromaterials Research & Exploration Science (ARES) at NASA Johnson Space Center (JSC) to obtain electron backscatter diffraction (EBSD) maps and energy dispersive X-ray spectroscopy (EDS) maps of select portions of the section. The EBSD data were collected under beam conditions of 20 kV, and ~9 μA, with step sizes varying from 0.05 to 2 μm. Following EBSD data collection, we processed the data using AZtecCrystal and MTEX, a free MATLAB toolbox. Results: We used the elemental X-ray maps to identify the areas of interest in the section, particularly phosphates and Zr-bearing grains. These areas of interest were then targeted for follow up EBSD and EDS analyses. We have identified, through combined EBSD and EDS analysis, the presence of zircon, reidite, and monoclinic-ZrO2 (baddeleyite) in 16RS08. In 16RS08, we have found singular grains of zircon, zircon with a vermicular baddeleyite halo (Fig. 2a), and granular zircon with reidite (Fig. 2b). The different Zr-rich phases and their corresponding textures signifies that this sample underwent a broad spectrum of pressure and temperature conditions during the impact event. For example, zircons surrounded by a vermicular baddeleyite and SiO2 intergrowth (i.e., Fig. 2a) have been shown to preserve evidence of the extremely high temperatures of impact melt, upwards of 2370 °C [2, 6]. Similarly, reidite and granular zircon (i.e., Fig. 2b) have been shown to preserve evidence of high pressure, as the transition to reidite occurs >30 GPa [3, 6]. Future Work: Next, we will process the EBSD and EDS data for these Zr-rich grains, specifically looking for indicators of cubic- or tetragonal-ZrO2 in the baddeleyite remnants and of shock-precursors to the reidite. To further inform our work, we will also obtain BSE images of these grains using a Hitachi S-4800 SEM in KMICF at the University of Arizona, as well as geochemical spot analyses via EPMA. The data collected will be used to constrain the formation conditions of the Ries Crater outer suevite. Acknowledgments: We thank Ken Domanik and Jerry Chang for their support with data collection. This work was supported by a University of Arizona RII Core Facilities Pilot Program grant and start-up funds to JJB. TME thanks A. Cavosie and N. Timms for assistance during field sampling. We acknowledge support from NASA’s Planetary Science Research program for analysis performed at JSC. References: [1] Melosh H. J. (1989) Oxf. U. Press. [2] Timms et al. (2017) EPSL 477, 52–58. [3] Cavosie et al. (2016) Geology 44:9, 703–706. [4] White et al. (2020) Nature Astr. 4, 974–978. [5] Crow C. A. et al. (2017) GCA 202, 264–284. [6] Timms et al. (2017) Earth-Sci. Rev. 165, 185–202. [7] Schmieder M. et al. (2018) GCA 220, 146–157. [8] Schwarz W. H. et al. (2020) M&PS 55:2, 312–325. [9] Stöffler et al. (2013) M&PS 43:4, 515–589.

Zircon↗

Extraction of Qualitative Features from Sensor Data Using Windowed Fourier Transform

In this paper, we use Matlab to model the health monitoring of a system through the information gathered from sensors. This implies assessment of the condition of the system components. Once a normal mode of operation is established any deviation from the normal behavior indicates a change. This change may be due to a malfunction of an element, a qualitative change, or a change due to a problem with another element in the network. For example, if one sensor indicates that the temperature in the tank has experienced a step change then a pressure sensor associated with the process in the tank should also experience a step change. The step up and step down as well as sensor disturbances are assumed to be exponential. An RC network is used to model the main process, which is step-up (charging), drift, and step-down (discharging). The sensor disturbances and spike are added while the system is in drift. The system is allowed to run for a period equal to three time constant of the main process before changes occur. Then each point of the signal is selected with a trailing data collected previously. Two trailing lengths of data are selected, one equal to two time constants of the main process and the other equal to two time constants of the sensor disturbance. Next, the DC is removed from each set of data and then the data are passed through a window followed by calculation of spectra for each set. In order to extract features the signal power, peak, and spectrum are plotted vs time. The results indicate distinct shapes corresponding to each process. The study is also carried out for a number of Gaussian distributed noisy cases.

Amini, Abolfazl M.↗

Pink Spinel in Apollo Impact Melt Rock 68815: Implications for Mg-Suite Magmatism

Introduction: Magnesian rocks from the lunar highlands are collectively termed the Mg-suite. Characterized by high (>60) Mg# (molar 100×Mg/[Mg+Fe]) in mafic phases and calcic plagioclase, these rocks are plutonic to hypabyssal in origin, and include a range of bulk mineralogies such as troctolites, dunites, norites, gabbronorites, and spinel troctolites [1]. These Mg-suite lithologies have distinct trace element concentrations and ratios that differentiate them from other lunar rock types. These rocks are ancient, generally dated to between 4.5 and 4.1 Ga, although it is unknown if this represents the full range of Mg-suite ages [1–2]. Among the Mg-suite lithologies, the spinel troctolites are relatively rare, to date only found in polymict breccias [3]. Spinel troctolites, as their name suggests, consist of calcic plagioclase and forsteritic olivine, with minor amounts of spinel (MgAl2O4), ± pyroxene and cordierite [1,4]. This form of spinel is often called ‘pink’ spinel because of its appearance in thin section under plane polarized light (PPL; Fig. 1), due to minor amounts of Cr. Spinel troctolites are generally plutonic or hypabyssal in origin (subsequently exhumed and incorporated into polymict breccias), or formed through impact processes (e.g., crystalline impact melt) [5]. A spinel-rich lithology has also been found in the Moscoviense region of the Moon via the Moon Mineralogy Mapper (M3) and lacks other mafic phases [6]. Finally, while the Mg-suite sampled thus far consists of plutonic (or hypabyssal) rocks, the question remains if such magmas could have erupted on the surface of the Moon [7]. These magmas have much lower density than mare basalts, but little sample or remote sensing evidence has been found to support the idea that extrusive Mg-suite volcanism occurred [7]. Here, we present a coordinated microanalytical study of spinel-bearing lithic and mineral clasts found in Apollo sample 68815. These data will be used to understand their petrogenesis (magmatic or impact) and modification histories, and to shed light on the existence of volcanic Mg-suite rocks. Sample Description: Apollo sample 68815 is a polymict impact melt breccia containing a variety of lithic and mineral fragments embedded in devitrified impact melt. This sample was chipped off the top of a boulder at Station 8 during the Apollo 16 mission and had an original weight of nearly 1.8 kg. In this study, we investigated two polished thin sections of 68815: 68815,17 and 68815,148, both containing spinel. Methods: The thin sections of 68815 were studied using optical light microscopy (PPL, cross-polarized light, and reflected light) with a Keyence VHX-7100 Digital Microscope. Each section was then X-ray mapped for 13–14 elements using a Cameca SX100 electron probe microanalyzer (EPMA) located in the Kuiper Materials Imaging and Characterization Facility (KMICF) at the University of Arizona. We have obtained geochemical information about the phases (olivine, plagioclase, spinel, pyroxene) in the thin sections also using the EPMA. In addition, we have used ThermoScientific Helios NanoLab 660 Focused-Ion-Beam Scanning-Electron Microscope (FIB-SEM) and a Hitachi S-4800 SEM (both in KMICF) to obtain backscattered electron (BSE) images and energy dispersive Xray spectrometry (EDS) maps of areas of interest. Using a JEOL 7900F SEM at the Astromaterials Research & Exploration Science (ARES) at NASA Johnson Space Center (JSC), we have obtained electron backscatter diffraction (EBSD) maps of the spinel-bearing portions of the thin sections. The EBSD data were collected under beam conditions of 20 kV, and ~90 μA, with step sizes varying from 0.05 to 2 μm. Following EBSD data collection, we processed the data using AZtecCrystal and MTEX, a free MATLAB toolbox. Results: We have found clasts with subophitic textures, that consist of primarily olivine and plagioclase, with minor amounts of pink Mg-Al spinel and pyroxene (Fig. 1). These clasts are up to ~1 mm in length and contain spinels up to 50 μm across. We have additionally identified pink Mg-Al spinels within the impact melt (i.e., not contained in lithic clasts) in both thin sections. In one instance, a single spinel grain is approximately 300 μm across (Fig. 1b, 2). The spinel fragments embedded in impact melt have varying compositions, typically distinct from the compositions of spinels in the lithic clasts. Spinel-Bearing Clasts: Ten lithic clasts with similar textures and mineral compositions were identified between 68815,17 (two clasts) and ,148 (eight clasts). These clasts fall into two groups. The first has skeletal olivine with intergranular plagioclase, with minor amounts of pyroxene and spinel (Fig. 1a, 2c, 2d). The spinel in these clasts are found amid the plagioclase. The second group have an intergranular texture of olivine and plagioclase, again with minor spinel and pyroxene. The second group may contain spinels surrounded by plagioclase, and spinels enclosed in olivine. Spinels located within both clast types range from no apparent Cr zoning, to reverse zoning (Cr-enrichment inward; Fig. 2d), to normal zoning (Cr-enrichment outward). In the clasts thus far investigated with EPMA, plagioclase compositions range from An# (molar 100×Ca/[Ca+Na+K]) 92–96. Olivine Mg# ranged from 77 to 94, while pyroxene had Mg# from 54–84. Spinel in the clasts have Cr# (molar 100×Cr/[Cr+Al]) 2–4 and Mg# 88–91, which is within the range of pristine and plutonic spinel troctolites [8]. Isolated Spinels: These crystals are generally euhedral to subhedral, and can exhibit reverse Cr zoning (Cr enrichment inward) or no apparent Cr zoning. The spinels thus far investigated via EPMA have Cr# 9–14 and Mg# 65–82. The Cr# for these spinels is within the range reported by [8], but have lower Mg#. Future Work: We will continue to process the EBSD data for these lithic and mineral clasts. We will also continue to characterize these clasts using EPMA and SEM. By thoroughly characterizing the various spinels and spinel-bearing clasts, we aim to constrain the petrogenesis of these minerals and rock fragments. Acknowledgments: We thank NASA for the loan of these thin sections. We thank Ken Domanik and Jerry Chang for their support with data collection. Work was supported by a University of Arizona RII Core Facilities Pilot Program grant and start-up funds to JJB. We acknowledge support from NASA’s Planetary Science Research program for analysis performed at JSC. References: [1] Shearer C. K. et al. (2015) Am. Min. 100, 294–325. [2] Borg L. E. et al. (2020) GCA 290, 312–332. [3] Warren P. H. (1993) Am. Min. 78, 360–376. [4] Dymek R. F. et al. (1976) LPS VII, 2335–2378. [5] Treiman et al. (2019) Am. Min. 104, 370–384. [6] Pieters et al. (2011) JGR: Plan. 116:E00G08. [7] Prissel et al. (2016) Icarus 277, 319–329. [8] Prissel et al. (2016) Am. Min. 101, 1624–1635.

spinel↗

Wave Measurements taken NW of Culebra Is., PR, 2023

Wave and sea surface temperature measurements collected by a Sofar Spotter buoy in 2023. The buoy was deployed on July 27, 2023 at 11:30 UTC northwest of Culebra Island, Puerto Rico, (18.3878 N, 65.3899 W) and recovered on Nov 5, 2023 at 12:45 UTC. Data are saved here in netCDF format, organized by month, and include directional wave statistics, GPS, and SST measurements at 30-minute intervals. Figures produced from these data are provided here as well. They include timeseries of wave height/period/direction and SST, GPS location, wave roses, and directional spectra. Additionally, raw CSV files from the Spotter's memory card can also be found below. NetCDF files can be read in python using the netCDF4 or Xarray packages, or through MATLAB using the "ncread()" command.

16 TIDAL AND WAVE POWER↗

Data for Post-Access Report: Tidal Energy Resource Characterization and Model Validation via the Assessment of Cost-effective Spatially Distributed Current Measurements from microFloat Swarms

This project evaluated how high-resolution, spatially distributed field data can be used to refine and validate site-scale hydrodynamic simulations of tidal channels. Use of such spatially-distributed field observations or site-scale hydrodynamic simulations will be needed for producing accurate predictions of tidal energy production over larger arrays of tidal turbines. The data contained in this archive includes field observations of water velocity from a series of short surveys in Agate Pass taken on August 20th, 2020 using station-keeping ADCP (acoustic doppler current profiler) measurements, drifting downward-looking ADCPs, and microFloats, gathered by researchers at APL-UW under a previous project. It also includes high-resolution hydrodynamic model runs for the site during that time generated under TEAMER support. An explanatory file (.pdf format) describes the survey and field data structure. Processed data is stored in .mat format. Survey summary: start time, number of floats, water level. microFloat data: - Cleaned data for each float in survey, - Time series of position (lat., lon.), depth (relative to surface), and horizontal velocity (u, v) Drifting ADCP data: - Cleaned data for each drifting ADCP in survey - Time series of position (lat., lon.), sample depth (relative to surface), and observed water velocity (u,v,w) Stationary ADCP data: - Cleaned data for each station-keeping data set in survey - Time series of position (lat., lon.), sample depth (relative to surface), and observed water velocity (u,v,w) FVCOM (Finite Volume Coastal Ocean Model) model outputs: Water level and velocity, as well as derived variables (e.g., power density) from the final refined model simulations for the Agate Pass subdomain. Contained in NetCDF and ASCII files for the duration of August 18th and 20th 2020 field surveys. Additional scripts for processing the data and generating the figures in the report are included.

16 TIDAL AND WAVE POWER↗

Predictive model using artificial neural network to design phase change material-based ocean thermal energy harvesting systems for powering uncrewed underwater vehicles

Uncrewed Underwater Vehicles (UUVs) are a major beneficiary of the phase change material (PCM)-based ocean thermal energy harvesting technology for their mission needs. However, this technology relies on different parameters and energy conversion steps that could be critical to the general energy generation efficiency. Sea trials showed that the design performed lower than their laboratory design specifications. This underperformance results from different factors, mainly the UUV’s trajectory, travel time, underwater ocean currents, temperature fluctuations, and biofouling on the heat exchanger due to long term underwater operations. Therefore, there exists a need to continuously monitor the ambient energy harvesting system and predict system performance, for mission planning purposes. Two major parameters influencing the energy harvesting system include the final pressure inside the hydraulic energy storage vessel or accumulator, and the electrical load value. Here, this work focuses on the hydraulic to electric energy conversion system. Therefore, a combination of numerical model and experimental testing is used to develop a predictive model using artificial neural network using MATLAB. After validation with experimental testing, 1000 data samples obtained from the numerical model are used to train the ANN. Compared to the experimental results, the developed ANN model can predict in less than a second the designed benchtop system’s total efficiency with less than 15 percent maximum error range. This predictive model development represents a cost-effective way for optimization and a computational energy efficient mode aboard UUVs for mission planning for deployed UUVs using PCM-based ocean thermal energy harvesting technology.

30 DIRECT ENERGY CONVERSION↗

SU-D-204-07: Comparison of AAPM TG150 Draft Image Receptor Tests with Vendor Automated QC Tests for Five Mobile DR Units

Purpose: To determine the consistency of digital detector (DR) tests recommended by AAPM TG150 and tests provided by commercially available DirectView Total Quality Tool (TQT). Methods: The DR tests recommended by the TG150 Detector Subgroup[1] were performed on 4 new Carestream DRX-Revolution and one Carestream DRX1C retrofit of a GE AMX-4 that had been in service for three years. After detector calibration, flat-field images plus images of two bar patterns oriented parallel and perpendicular to the A-C axis, were acquired at conditions recommended by TG150. Raw images were harvested and then analyzed using a MATLAB software previously validated[2,3,4]. Data were analyzed using ROIs of two different dimensions: 1) 128 x 128 ROIs matching the detector electronics; and 2) 256 x 256 ROIs, each including 4 adjacent smaller ROIs. TG150 metrics from 128 x 128 ROIs were compared to TQT metrics, which are also obtained from 128 x 128 ROIs[5]. Results: The results show that both TG150 and TQT measurements were consistent among these detectors. Differences between TG150 and TQT values appear systematic. Compared with 128 x 128 ROIs, noise and SNR non-uniformity were lower with 256 x 256 ROIs, although signal non-uniformity was similar, indicating detectors were appropriately calibrated for gain and offset. MTF of the retrofit unit remained essentially the same between 2012 and 2015, but was inferior to the new units. The older generator focal spot is smaller (0.75mm vs. 1.2mm), and the SID for acquisition is 182cm as well, so focal spot dimensions cannot explain the difference. The difference in MTF may be secondary to differences in generator X-ray spectrum or by unannounced changes in detector architecture. Further investigation is needed. Conclusion: The study shows that both TG150 and TQT tests are consistent. The numerical value of some metrics are dependent on ROI size.

60 APPLIED LIFE SCIENCES↗

Rapid Process to Generate Beam Envelopes for Optical System Analysis

The task of evaluating obstructions in the optical throughput of an optical system requires the use of two disciplines, and hence, two models: optical models for the details of optical propagation, and mechanical models for determining the actual structure that exists in the optical system. Previous analysis methods for creating beam envelopes (or cones of light) for use in this obstruction analysis were found to be cumbersome to calculate and take significant time and resources to complete. A new process was developed that takes less time to complete beam envelope analysis, is more accurate and less dependent upon manual node tracking to create the beam envelopes, and eases the burden on the mechanical CAD (computer-aided design) designers to form the beam solids. This algorithm allows rapid generation of beam envelopes for optical system obstruction analysis. Ray trace information is taken from optical design software and used to generate CAD objects that represent the boundary of the beam envelopes for detailed analysis in mechanical CAD software. Matlab is used to call ray trace data from the optical model for all fields and entrance pupil points of interest. These are chosen to be the edge of each space, so that these rays produce the bounding volume for the beam. The x and y global coordinate data is collected on the surface planes of interest, typically an image of the field and entrance pupil internal of the optical system. This x and y coordinate data is then evaluated using a convex hull algorithm, which removes any internal points, which are unnecessary to produce the bounding volume of interest. At this point, tolerances can be applied to expand the size of either the field or aperture, depending on the allocations. Once this minimum set of coordinates on the pupil and field is obtained, a new set of rays is generated between the field plane and aperture plane (or vice-versa). These rays are then evaluated at planes between the aperture and field, at a desired number of steps perceived necessary to build up the bounding volume or cone shape. At each plane, the ray coordinates are again evaluated using the convex hull algorithm to reduce the data to a minimal set. When all of the coordinates of interest are obtained for every plane of the propagation, the data is formatted into an xyz file suitable for FRED optical analysis software to import and create a STEP file of the data. This results in a spiral-like structure that is easily imported by mechanical CAD users who can then use an automated algorithm to wrap a skin around it and create a solid that represents the beam.

Howard, Joseph↗

Thermochemical Ablation Analysis of the Orion Heatshield

The Orion Multi-Purpose Crew Vehicle will one day carry astronauts to the Moon and beyond, and Orion's heatshield is a critical component in ensuring their safe return to Earth. The Orion heatshield is the structural component responsible for absorbing the intense heating environment caused by re-entry to Earth's atmosphere. The heatshield is primarily composed of Avcoat, an ablative material that is consumed during the re-entry process. Ablation is primarily characterized by two processes: pyrolysis and recession. The decomposition of in-depth virgin material is known as pyrolysis. Recession occurs when the exposed surface of the heatshield reacts with the surrounding flow. The Orion heatshield design was changed from an individually filled Avcoat honeycomb to a molded block Avcoat design. The molded block Avcoat heatshield relies on an adhesive bond to keep it attached to the capsule. In some locations on the heatshield, the integrity of the adhesive bond cannot be verified. For these locations, a mechanical retention device was proposed. Avcoat ablation was modelled in CHAR and the in-depth virgin material temperatures were used in a Thermal Desktop model of the mechanical retention device. The retention device was analyzed and shown to cause a large increase in the maximum bondline temperature. In order to study the impact of individual ablation modelling parameters on the heatshield sizing process, a Monte Carlo simulation of the sizing process was proposed. The simulation will give the sensitivity of the ablation model to each of its input parameters. As part of the Monte Carlo simulation, statistical uncertainties on material properties were required for Avcoat. Several properties were difficult to acquire uncertainties for: the pyrolysis gas enthalpy, non-dimensional mass loss rate (B´c), and Arrhenius equation parameters. Variability in the elemental composition of Avcoat was used as the basis for determining the statistical uncertainty in pyrolysis gas enthalpy and B´c. A MATLAB program was developed to allow for faster, more accurate and automated computation of Arrhenius reaction parameters. These parameters are required for a material model to be used in the CHAR ablation analysis program. This MATLAB program, along with thermogravimetric analysis (TGA) data, was used to generate uncertainties on the Arrhenius parameters for Avcoat. In addition, the TGA fitting program was developed to provide Arrhenius parameters for the ablation model of the gap filler material, RTV silicone.

Sixel, William↗

Testing and characterization of an IR CBIRD FPA

High performance IR cameras are in great demand for a variety of applications, including defense (e.g. night vision, missile detection) and space (e.g. imaging at IR wavelengths). IR focal plane arrays (FPAs) are an essential component of IR cameras. An FPA is a collection of detectors at the focal plane of an imaging device, where each detector can be thought of a pixel in the image it is detecting. The FPA in this study is based on Complimentary Barrier Infrared Detectors (CBIRDs). The CBIRD FPA was tested by imaging flat black body targets at three different temperatures, 20ºC, 25ºC, and 30ºC, with images captured using SEIR software. The important figures of merit for an IR FPA, Noise Equivalent Difference Temperature (NEDT), Noise Equivalent Irradiance (NEI), detectivity (D*), responsivity, and quantum Efficiency (QE), were extracted from the image data for each individual detector using a MATLAB program, as was the array Modulation Transfer Function MTF. The test results show a mean NEDT of 16.98mK and an MTF of about 0.34 at the Nyquist frequency. Overall, the array performs well in the MWIR range and takes quality IR images.

Das, Tanya↗