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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

SPND Sensitivity Calculations Using MCNP and Experimental Data from ACRR

The use of the Monte Carlo N-Particle Transport Code (MCNP) to calculate detector sensitivity for Self-Powered Neutron Detectors (SPNDs) in the Annular Core Research Reactor (ACRR) could be a vital tool in the effort to optimize the design of next-generation SPNDs. Next-generation SPND designs, which consider specific materials and geometry, may provide experimenters with capabilities for advanced mixed field dosimetry. These detectors will need to be optimized for configuration, materials, and geometries and the ability to model and iterate must be available in order to decide on the ideal. SPND design. SPNDs were modeled in MCNP which closely resembled the dimensions and location of actual detectors used in the ACRR. Tallies were used to calculate detector sensitivity. Using metrics from a previous report, oscilloscope data from pulses were manipulated in a Matrix Laboratory computing environment (MATLAB) script to calculate experimental detector sensitivity. This report outlines the process in which experimental data from ACRR pulses verified results from tallies in an MCNP ACRR model. The sensitivity values from experiments and MCNP calculations agreed within one standard deviation. Parametric studies were also performed with MCNP to investigate the effects of materials and dimensions of different SPNDs.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Simcenter TestLab Error in uploading reference profiles

Our team exports data from Matlab to a Simcenter TestLab formatted .mat file. The .mat file is then uploaded as a reference profile for Environmental Random Control vibration testing. Operations are applied to the data before export that introduce floating point errors. We are experiencing errors uploading the generated reference profiles if the frequency values have a certain level of floating point error (1 Hz vs 1.000000000000001 Hz). I generated some test cases of reference profiles that were giving us trouble, and outlined the errors we were seeing. Two different errors, depending on how the control frequency range is defined. Error is fixed by copy/pasting the frequency column in the Profile Editor.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Distributed Accelerometer IMU-2: An Interim Development Report

This report details the development of DAIMU-2, a MEMS accelerometer-based Gyro-Free Inertial Measurement Unit (GF-IMU). Previously, a GF-IMU called the Distributed Accelerometer IMU (DAIMU) was developed at LLNL using traditional analog accelerometers. This project leverages the experience gained from DAIMU and recent advances in sensor and embedded systems technology to develop DAIMU-2. The report introduces the theory and mathematics of a GF-IMU, followed by the design of an Unscented Kalman Filter (UKF) and simulation results. Finally, the DAIMU-2 prototypes developed to date are presented. After approximately one year of development work, the project was suspended, to be resumed in the future. Because of this, some efforts were partially completed, and this report attempts to indicate areas where further work is needed. MATLAB files, drawings, and other design documents have been archived for future project resumption.

42 ENGINEERING↗
DOE OSTI2016

A Study of Charge Propagation and Quantification of Noise Within the SuperCDMS Detector

This paper discuses two projects that were designed around the SuperCDMS experiment. The first is an attempt to understand the SuperCDMS simulation that was designed in the Matlab environment. The tests that were conducted were designed to study how electron and hole propagation is effected by the initial event location and the charge bias within the detector. The second experiment was an attempt to systematically quantify the level of noise produced within the SuperCDMS detectors by the dilution refrigerator. The results show only initial findings and no conclusive results were obtained.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

3D CFD Model Validation Using Benchmark Data of 1/16th Scaled VHTR Upper Plenum and Development of Wall Heat-Transfer Correlation For Laminar Flow

With support from the U.S. Department of Energy-Office of Nuclear Energy’s (DOE-NE’s) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, an effort has been pursued to support high-temperature gas-cooled reactor (HTGR) technology development and its modeling and simulation needs. There is a particular need for advanced modeling and simulation tools to predict thermal-fluid behavior in the nuclear reactor primary system, especially in the core and the lower and upper plena, during safety-related transients. In this report, two main such activities are presented relevant to the HTGRs: (1) three-dimensional (3D) computational fluid dynamics (CFD) validation using benchmark data from the upper plenum of Texas A&M University’s 1/16th scaled very-high-temperature gas-cooled reactor (VHTR), and (2) development of wall heat-transfer correlation for laminar flow in a wall-heated pipe. The CFD tool validation exercises can be helpful to choose the models and CFD tools to simulate and design specific components of the HTRGs such as upper plenum where jet mixing is a complex phenomenon. In a loss of forced circulation event, the laminar flow can be observed during the development of natural circulation flow. This work includes the development and validation of heat transfer correlations for laminar flow using the Nek5000 CFD code due to limited available experimental data for laminar flow conditions to guide low-order models (1D). In this report, the flow characteristics of a single isothermal jet discharging into the upper plenum was investigated using the Nek5000 Large-Eddy Simulation (LES) CFD tool. Several numerical simulations were performed for various jet-discharged Reynolds numbers ranging from 3,413 to 12,819. A grid-independent study was performed. The numerical results of mean velocity, root-mean-square fluctuating velocity, and Reynolds stress were compared against the benchmark data. Good agreement was obtained between simulated and measured data for axial mean velocities, except near the upper plenum hemisphere. The maximum predicted errors for axial mean velocities at various normalized coolant channel diameter heights of 1, 5, and 10 are 1.56%, 1.88%, and 3.82%, respectively. In addition, the predicted root-mean-square fluctuating velocity and Reynolds stress are qualitatively in agreement with the experimental data. The Nek5000 code was used to develop wall-heat transfer correlation for laminar flow in a cylindrical tube. Several simulations were performed for various Reynolds flow and wall-heat fluxes. A new heat transfer correlation was developed using data from Nek5000 simulation results and regression functions in Matlab. The developed heat transfer correlation is valid for various Reynolds flows from 200 to 2000. The predicted R² value for model fit was 0.875, which ensures that 87.5% of the model data lies on the Nek5000 data. Moreover, a machine learning (ML) tool was used to train and test the Nek5000 data. A good fit of the ML-based model was observed with the test data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Nested Comments in BNFC

Many languages such as Lua, Haskell and MATLAB support the ability to create nested comments. My proposal is to add a nested comment pragma to BNFC which would allow it to support the parsing of code. In this proposal I will provide sample implementations of nested comments for each of the lexer generators used by BNFC and outline the extension points by which I will add the functionality to the codebase.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Oak Ridge National Laboratory Current Magnitude and Frequency Response of SIPROTEC Relays

This project examined the current magnitude and frequency response of 7SJ61 and 7SJ62 SIPROTEC relays. The measured sample frequencies for the 7SJ61 and 7SJS62 relays were 800 and 960 Hz, respectively. The relays were tested at frequencies of up to 1 and 20.16 kHz. Two experiments were run in a lab environment. The first experiment was run at frequencies of up to 1 kHz every 10 Hz. The second one was run at frequencies of up to 20.16 kHz every 200 and 240 Hz (common factors of sample rate frequency) for the 7SJ61 and 7SJ62 SIPROTEC relays, respectively. Both experiments were set at a current magnitude of 0.1 and 0.5 Arms. A testbed with a Relay Test Equipment (RTE) was set to run the experiment at frequencies of up to 1 kHz, while a testbed with a RTE and High Frequency Power Source/ Amplifier (HFPS/A) was set to run at frequencies of up to 20.16 kHz. Currents of 0.1 and 0.5 Arms at different frequencies were connected at Phase A of the relays, and a current of 1 Arms (60 Hz) was connected at Phase B of relays. The 1 Arms current allowed the relays to trip and saved their record events by setting both relays at a definite-time overcurrent protection function. The record events were collected after the 7SJ61 and 7SJ62 relays tripped, and the current magnitude and frequency were plotted. The record events from the relays were collected using the DIGSI® 4 software. The SYNCHROWAVE® Event and MATLAB® software were then used to plot the current magnitude and frequency, respectively, from the relay record events. The current magnitude and frequency were plotted at time and frequency domain, respectively.The 7SJ61 and 7SJ62 relays measured source frequencies of up to 400 and 480 Hz, respectively. The percentage error of frequency for the relays increased after 400 and 480 Hz, respectively. The 7SJ62 relay with 0.1 and 0.5 Arms measured the alias frequency up to 1440 Hz. The 7SJ61 relay with 0.1 and 0.5 Arms measured the alias frequency up to 1000 and 1200 Hz, respectively. At frequencies higher than the alias frequency, the relays measured a DC amplitude current instead of an AC amplitude current. In the DC amplitude current zone, non-desired frequencies were measured by both relays. The current magnitude intensity (dB) for the 7SJ61 and 7SJ62 relays from 1 to 20.16 kHz were plotted. The current magnitude intensity for the 7SJ61 relay decreased up to 2600 and 1200 Hz for 0.5 and 0.1 Arms tests, respectively. However, the current magnitude intensity for the 7SJ62 relay decreased up to 4800 and 2640 Hz for 0.5 and 0.1 Arms tests, respectively. The current magnitude and frequency response study for the 7SJ61 and 7SJ62 relays was performed satisfactorily, observing the source, alias, and non-desired measured frequency zones for both relays.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Integrated Optimization and Control of a Hybrid Gas Turbine/sCO 2 Power System

During phase-I, the project team led by Echogen Power Systems (EPS) had two primary objectives based on investigating the application of gas turbines with supercritical carbon dioxide (sCO 2 ) power cycles. The first objective was to improve the overall efficiency and performance of a hybrid gas turbine/sCO 2 power system through a joint optimization of the two subsystems (gas turbine and sCO 2 power cycle) using non-linear optimization techniques that simultaneously evaluate thermal performance of the combined cycle. The hybrid power system included several points of interaction, including (but not limited to) gas turbine exhaust, fuel heating, inlet chilling and turbine cooling. The second objective was to establish a baseline transient response model of the hybrid power system and a notional microgrid and begin steps to integrate the control systems of the three major elements (gas turbine, sCO 2 cycle and grid controller). The project team established a baseline performance for a combined cycle power plant using a production gas turbine and scaled sCO 2 power cycle only utilizing exhaust heat recovery. Echogen’s non-linear techno-economic optimization code was extended by adding gas turbine component models derived from a in-house developed gas turbine design code. With the two cycles coupled by the gas turbine exhaust, design parameters of both cycles were allowed to vary simultaneously to determine performance opportunity versus isolated designs. Returning to the baseline gas turbine/sCO 2 power cycle transient models: Echogen had in-house developed sCO 2 cycle transient model in GT-Suite system simulation software, and had partnered with Siemens Finspång for gas turbine transient model, and Siemens PTI group to provide micro-grid load profile as well as hybrid power cycle generated load (power and frequency) analysis. The transient model for the SGT-750 Siemens gas turbine was a “black-box” functional mock-up interface (FMI) model developed by Siemens Industrial Turbomachinery in Finspång, Sweden. The SGT-750 is a twin-shaft gas turbine that produces 40 MW electricity with an efficiency of about 40% at ISO conditions. At 100% gas turbine throttle (load), the SGT-750 has average exhaust conditions of 114.6 kg/s and 469.8°C. The transient model for sCO 2 power cycle was developed by Echogen in GT-SUITE 1D system simulation software platform. The basic CO 2 flow circuit has single-shaft turbomachinery with net 11.5 MW electrical power output at design conditions. The power turbine has a double-ended shaft with one end connected to synchronous generator through a fixed-ratio gearbox. The other end of power turbine is connected to the compressor through a continuously variable transmission. The major components of the sCO 2 power cycle modeled include air cooled condenser/cooler, CO 2 compressor, recuperator, two waste heat exchanger coils, power turbine, continuous variable transmission, gearbox and generator. Integration of SGT-750 transient model and sCO 2 power cycle transient model was done in Matlab Simulink. In the integrated model, the gas turbine and sCO 2 power cycle interacted at two points, first one being the gas turbine exhaust gas flow rate and temperature, which were inputs to sCO 2 power cycle model. The second point was the distribution of micro-grid load demand signal between the SGT-750 generator and sCO 2 cycle generator. For a given combined-cycle load demand, the gas turbine load demand was equal to the total demand minus the sCO 2 cycle power generated. In the present study the integrated model was simulated for two cases of grid load demand: (i) for a step change, both positive-step and negative-step, in grid load demand (ii) for a micro-grid load demand curve provided by Siemens PTI group. Finally, the time series plots representing load demand versus integrated system response were presented including the sCO 2 power cycle control system performance plots. The actual generated power and frequency of both the generators, gas turbine and sCO 2 power cycle, was supplied to Siemens PTI group for dynamic grid assessment, results of which are provided in appendices.

03 NATURAL GAS↗
DOE OSTI2020

SRNL Analysis of ICCWR LCM and WAMS data for Corrosion and Cracking (FY2020 Progress Report)

The development of algorithms for machine learning and data analysis for the 3013 Surveillance Program is a collaborative effort by the Savannah River National Laboratory (SRNL) and the University of South Carolina (USC). For corrosion detection, Laser Confocal Microscope (LCM) or Wide Area 3D Measurement System (WAMS) data is extracted from large binary files, with software written to convert the data to physical attributes (e.g. height, color and grayscale values; all as functions of a location in a plane projection). It is the objective of this project to produce a user-friendly interface that incorporated all operations needed to perform surface examination. For this reason, a Matlab-based Graphical User Interface (GUI) was created to integrate data input with software developed for processing and evaluation. In summary, the GUI permits selective downloading of binary data, interrogation of attributes, data labeling, flagging of significant features, execution of Machine Learning (ML) algorithms, output of parameters for trained ML algorithms, reports of ML model accuracy with respect to labeled data, and generation of graphical representations of various analyses.

3013 corrosion↗
DOE OSTI2020

Parameter Sensitivity Analysis of the SparTen High Performance Sparse Tensor Decomposition Software (Extended Analysis)

Tensor decomposition models play an increasingly important role in modern data science applications. One problem of particular interest is fitting a low-rank Canonical Polyadic (CP) tensor decomposition model when the tensor has sparse structure and the tensor elements are nonnegative count data. SparTen is a high-performance C++ library which computes a low-rank decomposition using different solvers: a first-order quasi-Newton or a second-order damped Newton method, along with the appropriate choice of runtime parameters. Since default parameters in SparTen are tuned to experimental results in prior published work on a single real-world dataset conducted using MATLAB implementations of these methods, it remains unclear if the parameter defaults in SparTen are appropriate for general tensor data. Furthermore, it is unknown how sensitive algorithm convergence is to changes in the input parameter values. This report addresses these unresolved issues with large-scale experimentation on three benchmark tensor data sets. Experiments were conducted on several different CPU architectures and replicated with many initial states to establish generalized profiles of algorithm convergence behavior.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

The Polar System Analysis Package (Ver. 1 Specifications)

The Polar System Analysis Package (POLeSTAr) is a modern data fusion tool being designed to benchmark,evaluate and compare the polar physics and biogeochemistry of Earth System Models (ESMs). It is being written in Python to make use of libraries including Xarray, Matplotlib, SciPy, NumPy and PyTorch within algorithms developed specifically for polar system modeling. Initial POLeSTAr algorithms will be adapted from the MATLAB Ridgepack package that was designed to analyze the Model for Prediction Across Scales (MPAS) sea ice component of the Energy Exascale Earth System Model (E3SM) and sea ice model CICE within the Community Earth System Model (CESM).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Method for Calculating Delayed Gamma-Ray Response in the ACRR Central Cavity and FREC-II Cavity Using MCNP

This document presents the process for a new method developed for the characterization of the delayed gamma-ray radiation fields in pulse reactors like the Annular Core Research Reactor (ACRR) and the Fueled Ring External Cavity (FREC-II). The environments used to test this method in the ACRR were FF, LB44, PLG and CdPoly, and the environments used in the FREC-II were FF with rods-down, FF with rods-up, CdPoly with rods-down and CdPoly with rods-up. All environment configurations used the same fission product gamma-ray source energy spectrum. This method required the fission sites located in the MCNP KCODE source tapes. A FORTRAN script was written to translate and extract the coordinates for the fission sites. The 10K fission sites were then input it into an MCNP SOURCE mode script. Using a MATLAB script, a parametric analysis was done, and it helped determine that 10K fission sites are an appropriate number of coordinates to converge to the correct answer. The method gave excellent results and was tested in the ACRR, FREC-II and White Sands Missile Range (WSMR). This method can be applied to other pulse research reactors as well.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2019

Exactly and Easily Applying Experimental Boundary Conditions in Computational Structural Dynamics

Most experimental setups and environment specifications define acceleration loads on the component. However, Sierra Structural Dynamics cannot apply acceleration boundary conditions in modal transient analysis. Modal analysis of these systems and environments must be done through the application of a huge artificial force to a large fictitious point mass. Introducing a large mass into the analysis is a common source of numerical error. In this report we detail a mathematical procedure to directly apply acceleration boundary conditions in modal analyses without the requirement of adding a non-physical mass to the system. We prototype and demonstrate this procedure in Matlab and scope the work required to integrate this procedure into Sierra Structural Dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Non-Parametric Source Term Uncertainty Estimation

This paper describes the source term calculation documented in SAND2011-0128, ''Accident Source Terms for Light-Water Nuclear Power Plants Using High-Burnup or MOX Fuel'' and describes one method for implementing the calculation in MATLAB.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

NMDQi Thermal-hydraulic Optimization, Analysis, and Scoping Tool (TOAST)

The Thermal-hydraulic Optimization, Analysis, and Scoping Tool (TOAST) is developed to support the thermal-hydraulic aspects for the design of optimal irradiation vehicles for the Nuclear Materials Discovery and Qualification initiative (NMDQi). TOAST is currently a MATLAB-based tool that utilizes a simplified steady-state analytical model where a radial thermal resistance circuit is utilized to account for conduction through up to 2 capsules, a gas gap, a thermal bond, a sample cladding, and a sample, as well as the convection from a water coolant outside the capsules. The geometry is discretized axially along the user-defined height of a basket with a user-defined number of points/nodes, essentially turning TOAST into a semi-2D heat transfer analysis utility for cylindrical irradiation vehicles with practically any configuration, solving for temperatures radially at multiple axial locations. TOAST utilizes a Graphical User Interface (GUI) in-which a user can define the geometrical layout, the materials utilized for each component, the coolant characteristics, and the required solution. The user can define the geometry by inputting the diameters, thicknesses, and heights for each component, whereas the materials are defined via the user-provided constant or variable thermal conductivities. The user can select the coolant's inlet temperature, pressure, and the pressure drop across the height of the problem (which is used to calculate the velocity of the flow). Finally, the user can choose to solve for a sample heat generation rate limit by inputting temperature limits for the sample and capsules or can choose to purely solve for the axial temperature distributions of each component by inputting a pre-selected heat generation rate. Either way, TOAST provides the user with axial temperature distributions of all the components including the coolant, and the heat generation rate limit as well as the maximum outlet coolant temperature. The user can also choose to do one of 6 sensitivity analyses in TOAST. The sensitivity analyses yield plots of the sensitivity of the heat generation rate limit, the maximum coolant temperature, and the pressure limit for the annular components due to perturbations in 1-2 unknown variables based on the selected sensitivity analysis. Benchmarks between computations in TOAST and equivalent 2D axis symmetric ABAQUS models are presented, showing that TOAST results are within less than 3% of ABAQUS results in most cases, and a maximum of 8% difference in some cases. The benchmarks also revealed that this uncertainty is tied to the selection of an appropriate Nusselt number correlation and appropriate thermal conductivities, which the user can do from the GUI. Regardless, TOAST is demonstrated as a computationally efficient, highly accessible, and accurate utility for optimization and scoping s of studies different irradiation vehicles.

42 ENGINEERING↗
DOE OSTI2019

A Probability-based Model for Cost-effective Integration of Renewables into the Electricity Grid (Final Technical Report)

Funded by the U.S. Department of Energy during 2017-2019, the resource adequacy evaluation and production costing tool, called MATPLAN (MATlab and probability-based PLANning), was developed. MATPLAN can take into account the variable nature of renewable energy sources (both solar PV and wind farms), rather than considering deterministic data as in the current practice. It combines both probabilistic and optimization techniques to allow the determination of optimal system expansion policy and perform effective load carrying capacity (ELCC) analysis of a system with high renewable penetration. This is in contrast to the current practice that treats these variable sources as negative loads. This tool expects to benefit system planners by enabling them to consider renewable energy sources as options for their expansion plans. System operators can use MATPLAN to calculate the capacity credit for each solar/wind farm using probabilistic forecasts, thus enabling quantification of the operating reserve requirements to cope with uncertainty and variability of solar/wind output committed in the market.

14 SOLAR ENERGY↗
DOE OSTI2021

Load Deflection Data (LDDATA)

The current data set is load deflection data (Stress vs strain for compression stress and stress vs elongation for tension stress) from a material test. The goal of the test was to characterize this rubber-like material by looking at different variables: sample creation process, anisotropy, temperature, etc. The data will be used with the Texas A&M Online Research Experience for Undergraduates (O-REU) program. This program is joint with LANL. The student is not an official LANL student so they will not be able to access data and documents on the yellow. The load deflection data set is unclassified and will be sent to the student for analysis. Analysis will include curve fitting, scalar value analysis (stiffness, toughness, etc.) and analysis of variance to determine what variables statistically affect the response while correcting for covariates (temperature soak time and humidity). Data will be stored as Matlab structs and as HDF5 files. Same data, different formats.

36 MATERIALS SCIENCE↗