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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Investigation on Reaction Probabilities for Cyclic Voltammetry of Zirconium in LiCl-KCl Eutectic Molten Salt via Reverse-Engineering Method

This work focuses on a reverse-engineering program design for the cyclic voltammetry (CV) method to help illustrating, enhancing, and providing robustness in detection probability of the reactions during an electrorefining process. Analyzing the reactions from CV plots for 1.07, 2.49, and 4.98 wt% of zirconium chloride (ZrCl{sub 4}) in a LiCl-KCl molten salt at 773 K under different scan rates from Hoover's study has been accomplished. The results provide the current (Amp) versus potential (V) to investigate the optimum probability of the reactions at each peak. All simulations have been done through the commercial software MATLAB by using various experimental values existing in the literature as a first guess (e.g., the diffusion coefficients, and formal potentials). Preliminary results indicate that the zirconium reactions at each peak from cathodic to anodic side for 1.07 wt% are as followed: combination of Zr{sup +2}/Zr and Zr{sup +4}/ Zr{sup +2} at approximately -1.15 V, Zr{sup +2}/Zr at -1.5 V, Zr{sup +4}/Zr at -1.85 V, 70% Zr/Zr{sup +4} and 30% Zr/Zr{sup +2} at -0.5 V, and Zr{sup +2}/ Zr{sup +4} at - 0.25 V. These mentioned peaks are similar for 2.49 wt% at different voltages. In addition, the highest anodic peak for 4.98 wt% is found to be a combination of 30% Zr/Zr{sup +4} and 70% Zr/Zr{sup +2}. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

The Safe Guards Analysis Toolbox for Non-Proliferation Modeling

As nuclear technology becomes more widely available outside of nuclear weapon armed nations, the importance of preventing the illicit development of nuclear weapons has increased dramatically. A system of nuclear safeguards prevents proliferation; stopping the diversion of nuclear material from civilian facilities for the construction of nuclear weapons. These safeguards function through a combination of accountancy measurements that confirm the locations and amounts of nuclear material, and physical security at every nuclear facility. Safeguards Analysis (SGA) is a computational tool-kit developed using the MATLAB Simulink software suite that simulates a system of nuclear safeguards. The code allows users to examine current safeguards systems for weaknesses and to pinpoint the best locations for improvement. Developed as a modular framework, SGA does not directly model material flows, and relies on an outside code for this information. Instead, it represents a series of safeguards measurements that represent the values inspectors would see. Furthermore, SGA is designed to be modular in design, allowing the user to construct as complex or simple a model as needed: simulating an entire fuel-cycle comprising multiple facilities to a single measurement. SGA's use of outside codes for its material flow information makes it applicable to any model available to the user, as opposed to previous work that applied only to a single scenario. This paper discusses the initial development of SGA and the testing of these modules. The tests described in this paper demonstrate that the different pieces of SGA currently work as intended. Future work on the toolbox will develop and test SGA against a more realistic and complicated scenario using the various tools available in SGA. SGA currently has the pieces to perform analysis on safeguards systems, but needs to be tested on large, multi-MBA, systems to ensure the code is capable of handling such large models, to see the effect large models have on calculation speed, and if further work is needed to improve code efficiency. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

H-Canyon Tunnel Inspection System

Cable Analysis: The crawlers is controlled remotely via a long coax cable. In order to make sure we do not damage the cable which could cut our contact (data/power/control), we had to determine the bend radius of it at certain locations. This was done with dimensional analysis in Creo and inputting the information into a MATLAB program. Engineering Drawings: Since many of the crawlers parts, such as the main body, needed to be machined, part drawings were needed to communicate the details necessary for manufacture. Project Overview: The H-Canyon Ventilation Exhaust Tunnel due both to its long service history and exposure to aggressive process offgas chemistry has experienced degradation. In order to allow visual inspection of the condition of the radioactively-contaminated tunnel, a remote crawler was designed to deliver cameras as needed to return high definition video. My involvement included Vendor Searches, Cable Analysis, Thermal Analysis, and Engineering Drawings. Thermal Analysis: In order to see if internal active cooling is needed for the crawler, simplified thermal models were designed to approximate the internal temperatures: Conduction (assumes all heat sources are thermally-coupled to the crawler body); Convection (assumes no heat sources are thermally-bonded to crawler). Going Forward: Continue refining the cable analysis as the crawler model is updated and refined. Use my knowledge gained here to further my skills in engineering.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Radioactive Source Localization via Bayesian Particle Filter

In the event of a misplaced radioactive source or other emergency situation, measuring a radiation field, mapping its distribution, and determining a source location are essential tasks to ameliorating the situation. However, radiation fields may be extremely hazardous to human surveyors and minimizing received radiation doses is just as essential. Robots appear to be a potential solution to these problems. Beyond simply measuring radiation, the robot's computer processing capabilities offer a way to apply complex data analysis methods to radiation measurements in real-time. Methods which predict likely source locations can then feed this information into other processes, potentially improving path planning and enabling more efficient measurements. Given a robot mounted with a gamma-ray detector, can we: develop a methodology to account for detector performance across a wide range of source angles, distances, and photon energies? operate an autonomously navigating robot to effectively survey and characterize an area of interest? implement a data analysis method, conventionally used in measurements of motion, for source localization purposes? An open-source TurtleBot 3 robot, running Robot Operating System (ROS) on Ubuntu 16.04 LTS, was fitted with a Kromek GR1{sup R} Cadmium Zinc Telluride (CZT) solid-state gamma-ray detector. As a part of ROS, the packages OpenSlam, gmapping, and amcl were used to perform Simultaneous Localization and Mapping (SLAM), determining the robot's position and mapping the surrounding area. Data was acquired via Lidar mounted on top the TurtleBot 3. Detector Calibration Fit: The equation was fit to 365 counts of various energies, distances, and angles. A MATLAB{sup R} program was written to simulate measurements taken a robot on a random walk, with count data and positions discretized into finite element pixels. Using this program, a sample of 100 runs was performed on a map with a simulated source at the center, with a total of 200 of 2 pixels each. Similarly, multiple runs of the filter were performed on recorded robot measurement data. In both simulation and real tests, when corrected for errors (particles placed outside of bounds or on the robot, and simulation-specific errors), corresponding t-tests of predicted x and y-coordinates were within a 95% confidence interval of the actual position. For the real trial, these positions are slightly skewed right in the x-axis as the robot remained mainly to the left side of the source within the sample area. These simulations demonstrate potential validity for the usage of a particle filter as method of radioactive source localization. In the future, true real-time implementation and data fusion may further augment the performance of the robot to localize lost sources. Additionally, identification of multiple sources, determination of source types, and usage of a collimator are areas to potentially be explored.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

Development of a Smart Signal Detection Method for Cyclic Voltammetry via Artificial Neural Intelligence

The electrorefiner (ER) is the heart of pyroprocessing technology which is a high-temperature method for separating uranium from Used Nuclear Fuel (UNF). It is important to improve this technology with respect to nuclear materials detection and accountability. Artificial Neural Intelligence (ANI) is a novel data analysis and simulation method that can be applied to electrochemical data sets. This computational code, which has been performed using the commercial software MATLAB, can be trained to generalize adequate electrical current and potential simulated data sets for the unseen data with a high accuracy of prediction. For this purpose, a massive collection of cyclic voltammetry (CV) data sets by Hoover (2014), for 0.5, 1, 2.5, and 5 wt% of zirconium chloride in LiCl-KCl molten salt with different scan rates at 773 K has been considered. The computer is trained via ANI to predict the unseen data after providing suitable hidden layers and validation numbers. In addition, this work can trace the CV plot for a blind condition by interpolating between two simulated data set. The different hidden layers with various neurons (from 5 to 30) at several validation numbers (from 5 to 30) has been studied and the average percent error between experimental and theoretical data for 0.5 wt% with 200 and 450 mV/s has been calculated. Preliminary results demonstrate that if the number of hidden layers increases from one to three, the average error falls down from 44% to around 8%. The best condition which gives a minimum average percent error has been discovered and the simulated CV graph has been compared with the experimental data. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Neutronic-Thermal Hydraulic Benchmark using the SKINATH Point Kinetics Model

The SKINATH point kinetics program, whose name comes from: A Computer Program for Solving the Reactor Point Kinetics Equations with Simple Thermal Hydraulic-Feedback, had its origins at ORNL. The authors were interested in investigating low-level criticality events such as might occur in a waste repository of 55 gallon drums. We assume a relatively realistic thermal feedback for fissionable material in a container (called a reactor here), either critical or subcritical, with the potential to undergo a further excursion. The system is a homogeneous fissile fluid (mostly water), increasing in temperature as it absorbs all heat from the excursion. An ideal gas eternally cools the reactor at constant temperature T{sub c}. The SKINATH numerical concept in 1984 was to use an off-the-shelf Livermore ODE solver LSODE to solve the kinetic equations, much as one does today with MATLAB{sup TM}. Here, we develop a novel solver that uses the simplest of numerical methods including convergence acceleration. This work represents the initiation of a Multiphysics benchmark effort at INL. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Solution of the C5G7 Benchmark using the MOCUM Transport Code with ANSYS Unstructured Mesh

The objective of this research is to use finite element analysis (FEA) software to perform the unstructured meshing for reactor geometries and port the mesh information to MOCUM code. This workflow takes the advantage of the mature meshing module of the mainstream FEA; reduces the development effort for nuclear codes; and increases the nuclear code accuracy by using high quality meshes. In this work, ANSYS is used to build the reactor model and perform the unstructured meshing. The mesh information is then imported to the MOCUM code, a method of characteristic (MOC) transport program supporting unstructured meshes, using a developed MATLAB script. This procedure is verified by modeling the C5G7 benchmark problem. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

Performance of Energy Storage with Nuclear and Renewable Power Systems

Nuclear power represents a traditional baseload power component to the grid, while renewables output variable amounts of power at unpredictable times. Nuclear and renewable energy integration (NREI) systems have previously been explored as valuable additions to the existing power grid. Energy storage systems can allow power plants to store energy mechanically using pumped hydro or compressed air, thermally using liquid salt or steam accumulators, or chemically with batteries. The addition of energy storage to NREI systems can allow for the discharge of electricity to adjust to demands placed on the grid. The allocation of energy to adjust to grid demand is explored through a historical 8,760-hour analysis including nuclear, wind, and solar NREI systems with energy storage. The operation of storage systems is compared through simulations of different storage and discharge parameters based off the maximum daily grid demands, and hourly electricity prices. The overall profitability of storage systems is determined by comparing simulation results for the overall net revenue increases exceeding storage efficiency losses. Several NREI systems have increased the value of electricity for both solar PV and wind energy, in combination with nuclear power. Optimization of these systems is developed using MATLAB code and the simulations of the most profitable results are explored.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

Thermal Analysis of a Fuel Channel

To better study and understand the efficiency of fuel channels in reactors, a thermal analysis method was developed to conduct a study on any fuel channel to gain a better understanding of the temperature, heat rate, heat flux, and power profiles of such fuel channels. This study is conducted on a typical PWR (Pressurized Water Reactor) fuel assembly. The justification for specifically choosing a PWR fuel channel is because PWRs are the most widely used type of reactors around the world. This study was performed using NEWT code in SCALE 6, ANSYS, and SolidWorks. MATLAB and Excel were also utilized for some of the calculations. The power and heat rate of the fuel rod were calculated, alongside a 2D mapping of temperature, heat flux, directional heat flux, and thermal error on a 3D model of the fuel channel by analyzing the fuel rod and the surrounding water moderator independently.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Dynamic Analysis of the Next Generation Molten-Salt Breeder Reactor System

The conceptual design of a 1000 MWe two-fluid molten-salt breeder reactor (MSBR) station was performed during 1966-1970 at Oak Ridge National Laboratory (ORNL) with the goal of using the U-233 fuel cycle. This modular design with four reactor modules has an electrical output of 250 MW/module. The abundance of Th-232 was attractive in converting it to U-233 in a thermal breeder reactor. The two-fluid reactor is graphite-moderated 'with a {sup 7}LiF-BeF{sub 2}-UF{sub 4} fuel salt circulated through the core and a {sup 7}LiF-ThF{sub 4}-BeF{sub 2} blanket salt circulated through separate flow channels distributed throughout the core, as well as in a surrounding under-moderated region'. This reactor design has the advantages of not producing Pu isotopes and easy chemical separation of U-233 in the blanket material. The low fuel cost and the graphite replacement cost are attractive for such a design, with an estimated cost of $141/kW at 1968-dollar value. Potential non-proliferation features of MSBR, the low costs of raw materials, and established chemical separation of fission products are important factors in further development of MSBRs as an alternative to current light water reactors. These reactors are often referred to as liquid fluoride thorium reactor (LFTR). The current research on the modeling of the MSBR is based on an earlier work on the preliminary dynamic model of the MSBR concept. Revisions are made to the original model. The physics parameters have been modified for the current modeling effort based on additional information including more recent values of model parameters from and several typographical misprints in equations and values. Details of this effort will be presented in a full follow-up paper. This work presents the detailed nodal models of reactor fluid core, blanket fluid loop, graphite moderator, fuel fluid heat exchanger, and the blanket fluid heat exchanger. A set of simulations has been performed to validate the model using the commercial software MATLAB{sup TM}-Simulink. Simulations were performed for the cases of source perturbation and reactivity perturbation. The model is inherently nonlinear. A description of the MSBR conceptual design, along with a tabulation of relevant system parameters, is presented. The results of simulation demonstrate stable behavior of the MSBR dynamics, temperature feedback effects, and load-following capability of the MSBR system.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Integrating Sudo-Kaminaga Correlation to the Safety Analysis Code PARET-ANL

The accuracy and reliability of the reactor safety analysis code are critical because the specifications of the safety systems will depend on the analytical results for initiating events that could occur to the reactor. In safety analyses for power reactors, there exist many proficient and applicable computer codes such as RELAP5, CATHARE, RETRAN, or CATHENA. For non-power research reactors, an issue arises, however, due to lacking of such specialized codes. As a result, computer codes developed for the transient analysis of power reactors need to be applied carefully after proving their applicability to a specified research reactor. The safety analysis discussed in this paper concerns the prevention of fuel damage in hypothetical accidental scenarios for the new research reactor design at National Institute of Standards and Technology (NIST). One fuel integrity criterion is analyzed by investigating the Minimum Critical Heat Flux Ratio (MCHFR) in the core during the accidental transients. Some previous studies develop a relationship that correlates the nominal value of the minimum CHFR with the probability of CHF occurring and causing fuel damage. Using a statistical approach, the nominal value of the MCHFR for the low-enriched uranium (LEU) reactor needs to remain above the recommended limiting value to better ensure the safety of the research reactor. The probability of no fuel failure increases as the MCHFR increases. More importantly, as long as the MCHFR remains above the recommended limiting value of 1.301, the reactor will have a 90% probability that CHF is not reached. The safety analysis program PARET, developed by Argonne National Laboratory (ANL), is employed to obtain the CHFR information in this study. PARET/ANL is a digital computer programming code intended primarily for the analysis of test and research reactors that use plate-type (flat) fuel elements. This program has its own set of correlations that are used to calculate the CHFR. One of them is the Mirshak correlation, which was believed to provide the best estimation of the critical heat flux under the operating conditions of the new reactor at the time. Following the critical heat flux (CHF) experiments for vertical rectangular channels in the JRR-3 (Japan Research Reactor unit 3), the CHF calculations now use the Sudo- Kaminaga correlation. This correlation is of our interest because this method has an enhanced geometric similarity and an increased range of applicability that are more representative of the actual operating conditions of our current design, and has a more mechanistic approach [9]. Since the source code of PARET/ANL is inaccessible, a MATLAB-based utility has been written to integrate the Sudo-Kaminaga correlation as a substitute for the Mirshak correlation in the code. As a result, both the Mirshak and the Sudo-Kaminaga correlations are used to analyze the standard output from the PARET program to calculate the CHFR.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Equivalent-Circuit Models of Grid-Forming IBRs for Electromagnetic-Transient Simulations: Preprint

We derive equivalent-circuit models of control- and physical-layer subsystems of grid-forming (GFM) inverter-based resources (IBRs) for electromagnetic-transient (EMT) simulations. Three different primary controllers are considered: droop, virtual synchronous machine (VSM), and dispatchable virtual oscillator control (dVOC). In addition, the models include cascaded voltage- and current-control loops, and LCL output filters. Simulations for a single-inverter setup and for a network of five inverters in a modified IEEE 14-bus topology are presented. The equivalent-circuit models simulated with analog electronic circuit-simulator software (in our case, LTspice) offer the same accuracy with up to 150x lower computational burden compared to block-diagram-based implementations in commercial off-the-shelf EMT software (in our case, MATLAB-Simulink).

circuit-equivalent model↗
DOE OSTI2022

Development of a water source heat pump hardware-in-the-loop (HIL) testing facility for smart building applications

Over the last decade, the global fight against climate change through electrification has led to an increase in research on building heating, ventilation, and air conditioning (HVAC) systems that utilize intelligent control algorithms to provide demand-side grid service while maintaining the thermal comfort of building occupants. As the pivotalpoint between building electricity consumption and indoor thermal comfort, high-efficiency electrical heatpumps are at the center of these emerging studies, and various grid-interactive and occupant-comfort control algorithms have been developed for them. The impact of these algorithms on the heatpump operation andperformance under different weather, building load, and grid requests calls for investigation and verification via experimental tests with actual heat pumps integrated with real-time building and grid responses. This study presents a Water-Source Heat Pump Hardware-in-The-Loop (HIL) Test Facility developed with the capability to perform such tests. The hardware configuration for this testfacility introduces a hydronic system that emulates the conditions for the heat pump water-side, and a duct system that emulates conditions for the heat pump airside. Both data acquisition and emulator control are implemented through the National Instruments (NI) LabVIEW software running on an NI PXIplatform. The HIL mechanism based on the hardware-software integration that allows the testbed to communicate with a generic simulation environment is also discussed. Currently, the test facility setup includes a single heatpump and virtual building model in EnergyPlus coupled with an occupant behavioral model in MATLAB. Preliminary test results of the current setup demonstrate the building load emulator's ability to track the simulated gone temperature with a Root Mean Square Deviation (RSME) below 0.12°C (0.216°F). An uncertainty analysis based on sensor accuracies shows that the heat pump coefficient of performance (COP) can be measured with a relative uncertainty of 10.4% in cooling and 3.7% in heating. Apartfrom the current testing on a single heat pump, the test facility also provides the flexibility to include additional heat pumps to form a heat pump cluster, as well as coupling the heat pump with active thermal storage to provide enhanced demandflexibility.

Calfa, Caleb↗
DOE OSTI2024

Dynamic Modeling of Near Isothermal Compressor for Transcritical Carbon Dioxide Cycle

Compressors are the major energy consumption components in vapor compression systems, drawing much research effort in reducing carbon emissions. The isothermal compressor integrates the compressor chamber and gas cooler to achieve near isothermal compression, reaching up to 30% energy reduction compared to the traditional isentropic compression work. This paper presents a detailed isothermal compressor model combined with a generalized liquid piston model to account for the carbon dioxide (CO2) isothermal compression process. The model is established based on MATLAB environment. The model uses the real experimental data as boundary conditions and initial settings, which also considers the CO2 solubility in liquid piston (mineral oil) for designing, optimizing and customizing the compression chambers. The validation was carried out with experimental data using a prototype with 3.5 kW capacity. The results have demonstrated the accuracy of the dynamic model (6.2% relative error for chamber pressure and 0.5 K deviation for chamber temperature), which provide a guideline for designing and customizing the isothermal compression cycle.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Beam-water Interaction System for PFAS Treatment

Per- and polyfluoroalkyl substances (PFAS) are persistent environmental contaminants that resist conventional water treatment due to the strength of their carbon fluorine bonds. Electron-beam irradiation has shown potential for PFAS degradation; however, its effectiveness is limited by shallow energy penetration, restricting treatment to the upper region of the water surface. This project addresses this limitation through the design of a compact mixing chamber and recirculation system to enhance surface renewal and improve treatment efficiency. The system, constrained to approximately 10 cm 20 cm 20 cm, was developed using computer-aided design, computational fluid dynamics (CFD), and MATLAB-based hydraulic analysis to evaluate flow behavior, head loss, and pump compatibility. A perforated-baffle geometry and controlled flow range of 1 3 GPM were implemented to promote uniform mixing while maintaining a stable free surface for electron-beam exposure. Modeling results indicate improved flow distribution and reduced dead zones, supporting more consistent interaction between contaminated water and the beam. The final design enhances hydrodynamic performance within system constraints, contributing to scalable electron-beam water treatment technologies for improved PFAS degradation.

Ramos, Gerardo [Northern Illinois U.]↗
DOE OSTI2025

Work updates for Prof. Lee

I want to update my work which I have done so far in INL, including matlab GUI for neutron radiography and some fuel EBSD works.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Computational model for predicting the performance of fluid-powered longiudinal piston modules and motors.

A Matlab-based computational model line code has been developed for predicting the performance of fluid-powered longitudinal piston modules and motors. A system of equations has been developed based upon applying the conservation of mass to a control volume defined by a reciprocating piston. The method results in a differential equation that predicts the pressure gradient in the control volume as a function of rotor rotation angle. The model uses an integration scheme to predict the pressure state in the control volume. The orifice equation is used to predict the flowrate in and out of the control volume. From the resulting solution the torque and speed characteristics of the motor configuration can be predicted, optimized, or tailored for a specific requirement.SAND2019-13983 M Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Taylor, Michael↗
DOE OSTI2016

ALFA Station Keeping Results for Seabotix vLBV300 Underwater Vehicle near Newport, OR

This data set presents results testing the station keeping abilities of a tethered Seabotix vLBV300 underwater vehicle equipped with an inertial navigation system. These results are from an offshore deployment on April 20, 2016 off the coast of Newport, OR (44.678 degrees N, 124.109 degrees W). During the mission period, the sea state varied between 3 and 4, with an average significant wave height of 1.6 m. The vehicle utilizes an inertial navigation system based on a Gladiator Landmark 40 IMU coupled with a Teledyne Explorer Doppler Velocity Log to perform station keeping at a desired location and orientation.

16 TIDAL AND WAVE POWER↗