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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Surface Water Quality Data from Beaver-Impacted Streams; Trail Creek and East River, Colorado 2025

This data package contains surface water chemistry measurements collected in 2025 to evaluate how beaver damming and low-tech process-based stream restoration influence water quality and metal mobility in mountainous headwater systems of the Upper Colorado River Basin. Sampling was conducted at Trail Creek (Taylor Park watershed, Colorado), a tributary undergoing restoration through installation of low-tech process-based structures (i.e., beaver dam analogs), and at off-channel beaver ponds within the East River floodplain (East River watershed, Colorado). Samples were collected along longitudinal transects spanning upstream control reaches, beaver-influenced ponded reaches, and downstream segments. Additional samples were collected from near-surface pore waters within a beaver dam seepage face. The dataset includes concentrations of major and trace elements measured by inductively coupled plasma–mass spectrometry (ICP-MS) and inductively coupled plasma–optical emission spectrometry (ICP-OES), major anions measured by ion chromatography (IC), and dissolved organic carbon (DOC; reported as non-purgeable organic carbon, NPOC). Samples were size-fractionated at 0.45 micrometers (µm), 0.22 µm, and 0.02 µm to distinguish particulate (>0.45 µm), colloidal (0.22–0.02 µm), and dissolved (<0.02 µm) fractions. The data package consists of comma-separated value (.csv) files containing tabulated chemical concentration data, sample metadata (site identifiers, geographic coordinates, sampling dates, fraction type), and quality control flags. All files are provided in open, non-proprietary formats that can be accessed using standard data analysis software such as Microsoft Excel, R, Python, MATLAB, or other programs capable of reading .csv files. Units, detection limits, and analytical methods are documented in accompanying metadata files. The dataset is designed to support analyses of (1) how beaver impoundment and restoration structures alter elemental partitioning and transport, (2) the role of iron and organic carbon in mediating trace metal mobility, and (3) reach-scale changes in water quality across restoration gradients. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231. Part of this work was performed at SLAC Accelerator Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-76SF00515.

Anions↗
DOE OSTI2026

Chamber flux measurements at US-ORv located in Columbus, OH from 2010 to 2014

Closed and flow through chamber methane and carbon dioxide flux measurements were taken at different water depths, vegetation cover types, and times of day at the Wilma H. Schiermeier Olentangy River Wetland Research Park (ORWRP) from 2010 to 2014 to investigate the effects of these variables on gas exchange. The ORWRP consists of two urban experimental wetland basins constructed in 1994 within the footprint of the AmeriFlux tower US-ORv that was functional from 2011 to 2016. Chamber measurements are separated into flow through and closed chamber files (flow_through.csv, closed_chamber.csv) with location labels corresponding with coordinates listed in the metadata file (Location_ID.csv). All files are .csv and could be viewed or analyzed with Excel, MATLAB, Python, or R.

EARTH SCIENCE > ATMOSPHERE > ATMOSPHERIC CHEMISTRY↗
DOE OSTI2026

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

Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 2. Evaluating Controls on Flow Persistence in an Urbanized Catchment”

This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in an urbanized catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, distributed temperature sensing (DTS), continuous self-potential (SP) monitoring, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Field_Application subfolder contains the ATS XML input scripts, data files, output data for the SP site. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. The flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.m can only be used with COMSOL with MATLAB) is executed using the ATS output data to simulate the potential field. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) DTS Contains collated DTS data including raw Stokes and anti-Stokes measurement (provided as .h5 file). It also includes DTS processing.ipynb, a Jupyter notebook for calibrating the DTS data using dts_calibration Python package. cooler_calibration.csv is the DTS calibration CSV used in the calibration sequence. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion. 6) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 7) SP Contains the SP data collected in field at the SP sites (provided as CSV files). 8) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). Note: Code files (.ipynb, .py, .xml) can be opened in any standard code editor, .exo file can be viewed using Paraview, .h5 files can be opened using HDFView software and h5py Python package, and .resipy file can be opened with the open-source ResIPy software.

ATS↗
DOE OSTI2019

Frio 2 : Raw CASSM Datasets

This submission contains the raw SEG-2 seismic datasets acquired in the CASSM experiment, part of the Frio-2 GCS pilot. The directory also includes the log files generated by the Geode seismograph (geometrics) used during acquisition. All SEG-2 files have a ".dat" suffix and log files have a ".log" suffix. Various supplementary files are also included (see below) including a set of flow and seismic models used in interpretation. The earlyGathers.tar.gz resource includes raw data acquired during system testing and pre-injection periods. All files correspond to the single source continuous active source seismic monitoring (CASSM) dataset described in the following papers. Daley, T.M., Solbau, R.D., Ajo-Franklin, J.B., and Benson, S.M. 2007, “Continuous Active-Source Seismic Monitoring of CO2 Injection in a Brine Aquifer”, Geophysics, Vol. 72, No. 5, p. A57-A61 Daley, T., Ajo-Franklin, J.B., and C. Doughty, 2011, “Constraining the reservoir model of an injected CO2 plume with crosswell CASSM at the Frio-II Brine Pilot”, International Journal of Greenhouse Gas Control (IJGGC), Vol. 5, No. 2, pp. 1022-1030. Zhu, T., Ajo-Franklin, J., and T. Daley, 2017, “Spatio-temporal changes in seismic attenuation caused by injected CO2 at the Frio-II pilot site, Dayton TX, USA”, Journal of Geophysical Research – Solid Earth, Vol. 122, No. 9, pp. 7156-7171. cassmDataRaw.tar : The main tar file containing all of the SEG-2 field records and the Geode log files. Stacking, time sampling (dt), and date information should be in the .log file. Each file corresponds to a single stacked gather for the single source deployed. Data was recorded on a 24 channel downhole hydrophone array. fieldNotes.tar : contains some of the scanned observer logs from the experiments as well as some planning images fieldPhotos.tar : contains a small set of photographs taken during the experiment frio_flow_seismic_models.mat : a matlab structure containing a set of flow and seismic models used in interpreting the CASSM data (see IJGGC paper).

GCS↗
DOE OSTI2019

CIAM Climatological Isolation and Attraction Model–Climatological Lagrangian Coherent Structures

Matlab code to compute climatological Lagrangian Coherent Structures. Climatological Lagrangian Coherent Structures have been shown to be effective in identifying regions that persistently attract or isolate trajectories at the sea surface, while also identifying recurrent transport patterns and transport barriers. It is therefore known as the Climatological Isolation and Attraction Model within NETL DoE. This repository also includes the information necessary to replicate the cCLS published in the paper which is freely available at: https://www.nature.com/articles/s41598-018-23121-y The method is developed in: Duran, R., F. J. Beron-Vera, M. J. Olascoaga (2018). Extracting quasi-steady Lagrangian transport patterns from the ocean circulation: An application to the Gulf of Mexico. Scientific Reports, 8(1), 5218. https://www.nature.com/articles/s41598-018-23121-y The method is successfully evaluated with a large dataset of satellite-tracked drifters, and with a different ocean model in: Gough M.K., F. J. Beron-Vera, M. J. Olascoaga, J. Sheinbaum, J. Jouenno, R. Duran (2019). Persistent Lagrangian transport patterns in the northwestern Gulf of Mexico. J. Phys. Oceanogr., 49, 353–367, https://doi.org/10.1175/JPO-D-17-0207.1 and with floats and yet a different model in: Maslo, A., Azevedo Correia de Souza, J. M., Andrade-Canto, F., & Rodríguez Outerelo, J. (2020). Connectivity of deep waters in the Gulf of Mexico. Journal of Marine Systems, 203, 103267. https://doi.org/10.1016/j.jmarsys.2019.103267

CIAM↗
DOE OSTI2022

Task 5 Field Scale Stress Modeling Scripts and Data

Project data and Matlab scripts for Task 5 of the project “A Non-Invasive Approach for Elucidating the Spatial Distribution of In Situ Stress in Deep Subsurface Geologic Formations Considered for CO2 Storage”. Project includes the static earth models, grids and other data used for simulation and calibration of the two case studies addressed in the project, referred to herein as the “FutureGen2” model (Illinois) and the “Perch” model (Otsego county, Michigan). Large Dataset, please contact EDXSupport@netl.doe.gov

FutureGen2↗
DOE OSTI2025

AP4

The NETL-CMU Air Pollution Emissions Experiments and Policy version 4 (AP4) MATLAB to Octave model.

AS↗
DOE OSTI2025

Proton pulse charge calculation algorithm in Beam Power Limiting System at Spallation Neutron Source

A proton pulse charge calculation algorithm in the Beam Power Limiting System (BPLS) at the Spallation Neutron Source (SNS) was developed and implemented in an FPGA. The algorithm calculates one-minute running average of the pulse charges and issues a fault to the Personal Protection System (PPS) and the Machine Protection System (MPS) when a limit is reached. A bit-accurate model of the algorithm was first developed and tested in Matlab® and then implemented and simulated in VHDL using Vivado® design environment. Finally, the algorithm was verified on a µTCA-based hardware platform.

Bobrek, Miljko [ORNL] (ORCID:0000000332763451)↗
DOE OSTI2023

Prototype HB650 Cryomodule Heat Loads Simulations

During the design stages of the PIP-II cryomodules, many analytical calculations and FEA have been performed on simpler geometry in order to estimate the heat loads and also to optimize the design. To better analyze the cryomodule cold tests, simulations have been performed with MATLAB to determine the temperature of the main components during cool down and to determine the heat loads of the cryomodule. These simulations have been applied to the High Beta 650 MHz prototype cryomodule design and compared to the cold tests performed on it.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Automatic Calibration of a Geomechanical Model from Sparse Data for Estimating Stress in Deep Geological Formations

Summary In this study, we demonstrate geomechanical modeling with fully automatic parameter calibration to estimate the full geomechanical stress fields of a prospective US carbon dioxide (CO2) storage site, based on sparse measurement data. The goal is to compute full stress tensor field estimates (principal stresses and orientations) that are maximally compatible with observations within the constraints of the model assumptions, thereby extending pointwise, incomplete partial stress measurement to a simulated full formation stress field, as well as a rough assessment of the associated error. We use the Perch site, located in Otsego County, Michigan, USA, as our case study. The input data consist of partial stress tensor information inferred from in-situ borehole tests, geophysical well logs, and processing of seismic data. A static earth model (SEM) of the site was developed, and geomechanical simulation functionality of the open-source MATLAB Reservoir Simulation Toolbox (MRST) was used to model the stress field. Adjoint-based nonlinear optimization was used to adjust boundary conditions and material properties to calibrate simulated results of observations. Results were interpreted through a Bayesian framework. The focus of this paper is to demonstrate how the fully automatic calibration procedure works and discuss the results obtained; it does not attempt a detailed analysis of the stress field in the context of the proposed CO2 storage initiatives. Our work is part of a larger effort to noninvasively determine in-situ stresses in deep formations considered for CO2 storage. Guided by previously published research on geomechanical model calibration, our work presents a novel calibration approach supporting a potentially large number of linear or nonlinear calibration parameters to produce results optimally agreeing with available measurements and thus extend partial pointwise estimates to full tensor fields compatible with the physics of the site.

Engineering↗
DOE OSTI2024

QSpace - An open-source tensor library for Abelian and non-Abelian symmetries

This is the documentation for the tensor library QSpace (v4.0), a toolbox to exploit ‘quan tum symmetry spaces’ in tensor network states in the quantum many-body context. QSpace permits arbitrary combinations of symmetries including the abelian symmetries $\mathbb{Z}_n$ and U(1), as well as all non-abelian symmetries based on the semisimple classical Lie algebras: A n , B n , C n , and D n , or respectively, the special unitary group SU(n), the odd orthogonal group SO(2n+1), the symplectic group Sp(2n), and the even orthogonal group SO(2n). The code (C++ embedded via the MEX interface into Matlab) is available open source as of QSpace v4.0 on bitbucket under the Apache 2.0 license. QSpace is designed as a bottom-up approach for non-abelian symmetries. It starts from the defining representation and the respective Lie algebra. By explicitly comput ing and tabulating generalized Clebsch-Gordan coefficient tensors, QSpace is versatile in the type of operations that it can perform across all symmetries. At the level of an ap plication, much of the symmetry-related details are hidden within the QSpace C++ core libraries. Hence when developing tensor network algorithms with QSpace, these can be coded (nearly) as if there are no symmetries at all, despite being able to fully exploit general non-abelian symmetries.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Final Technical Report: Distributed Controls for High Penetrations of Renewables

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

14 SOLAR ENERGY↗
DOE OSTI2016

Echo™ User Manual

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

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Analysis of the LSC microbunching instability in MaRIE linac reference design

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

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

ROMUSE 2.0 User Manual

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

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Echo

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

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Multi-compartment Aerosol Transport Model

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

54 ENVIRONMENTAL SCIENCES↗