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 613 records · Page 34

Machine Learning Software for Cylindrical Battery Design and Performance Prediction

A software that delivers optimal design parameters and performance predictions for cylindrical cells, which can range in size from micro batteries to EV batteries, is developed. The Cylindrical battery design V1.0 is comprised of three types of cylindrical batteries, Micro battery (Primary), Micro battery (Secondary) and 18650/21700/xxxxx cylindrical battery. The software was developed in MATLAB. The software has the capability to output the cell design with the capacity ranges from several mAh to several million Ah. The software utilizes machine learning and includes a graphical user interface to enable rapid prototyping to accelerate energy storage research, development, and manufacturing.

25 ENERGY STORAGE↗

Pore-Scale Simulation of Spectral Induced Polarization

A new model was developed to simulate spectral induced polarization (SIP) in porous media. The model is based on a pore-network approximation of the soil or rock, in which bonds between pores are represented by equivalent circuits. In the circuit model, electrolytic and surface conduction occur in parallel, with the electrolyte behaving as a single resistor, and the mineral surface behaving as a resistor and capacitor in series. Based on the pore geometry, fluid conductivity, and mineral surface conductivity, an impedance network is generated. SIP spectra are calculated by solving the sinusoidal steady-state electrical problem for a range of frequencies. The new model provides a tool to bridge the gap between (1) our fundamental understanding of polarization processes at the scale of a single pore or grains and (2) existing simulation capabilities for macroscale systems that do not account for pore-scale properties. The model is implemented in Matlab and can solve two-dimensional and three-dimensional (2D and 3D) networks. Examples for regular cubic lattices are presented, but the method and code are flexible with respect to pore geometries and allow for future leveraging of existing open-source capabilities for generation of pore networks and extraction of network geometries from micro imaging of core samples.

58 GEOSCIENCES↗

Solute effects upon dislocation motion and recovery in Mg alloys (Final Report)

The objective the research was to develop a firmer understanding of the interactions between substitutional solute atoms and dislocations within Mg alloys. These interactions govern the absolute and relative mobilities of various types of dislocations in Mg (e.g., basal < a>, non-basal < a>, and pyramidal < c+a>). Furthermore, they have an impact on dislocation recovery processes (e.g., cross-glide, climb, rearrangement, and annihilation). Ultimately, solute-dislocation interactions strongly impact a) strain hardening, b) strain rate sensitivity, c) plastic anisotropy of textured polycrystals, d) texture evolution, e) dislocation substructures evolution, and even f) recrystallization behavior. It is broadly known that Mg alloys exhibit poor low temperature formability. Since formability is largely governed by properties (a) – (c) in this list, it is critical to develop a better understanding of solute-dislocation interactions, if one hopes to improve the situation. While the theory of static solute strengthening is well developed, especially for alloys with face centered cubic (FCC) crystal structures, there are outstanding questions related to applications to hexagonal close packed (HCP) crystal structures and dynamic strain aging (DSA) in materials of various structures. DSA has far-reaching implications for metal formability and, in the case of Mg alloys, appears to correlate with the so-called rare-earth (RE) texture which has been shown to benefit formability. The Portevin-Le Chatelier (PLC) effect, and associated negative strain rate sensitivity, occur at higher temperatures (>100℃) in Mg alloys as compared with similar Al alloys, even though they have similar melting points and solute diffusivities. Our preliminary research has shown that modern, physics-based models of DSA can be tuned to describe the behavior of Mg alloys if a rather higher activation enthalpy is assumed for cross-core diffusivity. While this partially explains the delay in DSA to higher temperatures, it is also hypothesized that this delay is due in part to the intrinsically more thermally activated (rate sensitive) nature of non-basal < a> dislocation motion which is required for macroscopic flow of Mg alloys, whereas octahedral slip in many FCC metals like aluminum is essentially athermal at room temperature. It was originally proposed to employ a combination of in-situ diffraction-based experimental characterization to validate existing theory. It was envisioned to perform in-situ transmission electron microscopy (TEM) to assess individual dislocation behavior and in-situ high-energy X-ray diffraction (HEXRD) techniques which were showing great promise for elucidating collective dislocation behavior, 2 including recovery, especially if the contributions to various forms of diffraction peak broadening (𝜂𝜂,𝜔𝜔,and 2𝜃𝜃) can be effectively integrated. Finally, it was envisioned to perform discrete dislocation dynamics (DDD) modeling approaches to aide in the interpretation of both TEM and HEXRD experiments. In the end, mechanical tests were performed on more complex Mg alloys which exhibited evidence of dynamic strain aging, and this led to the establishment of another project. Mechanical test data obtained at McMaster University served as the basis of an assessment of the applicability of Bazinski’s “stress equivalence” theory of solute strengthening to polycrystalline alloys of Mg. Although we did not succeed in applying the approach to Mg alloys, we did develop expertise with the HEXRD approach using a BCC, β-Ti alloy and demonstrated numerous new capabilities that may be applied to any polycrystalline material in collaboration with researchers at CHESS and around the world: (1) assessment of details of the elastoplastic transition (yielding) using a combination of HEXRD and full-field polycrystal plasticity modeling, (2) the first-ever experimental observation of strong stress rotation within the individual grains of a polycrystalline material, and (3) a comprehensive analysis of the grain-level dislocation density evolution based upon diffraction peak broadening along 𝜂𝜂,𝜔𝜔,and 2𝜃𝜃 directions. This final aspect allowed us to confirm that dislocations were gliding on multiple plane types and not restricted to {110} type planes, and it also provided clues as to why some grains were unloading during straining, with surprising implications for our understanding of the effects of geometrically necessary dislocations (GNDs). Finally, graduate student, Mohammed Shabana, developed a MATLAB code which confirmed the conclusions of Prof. Catalin Picu (Rensallear Polytechnic Institute, RPI) regarding the effect of solute-trapped, forest dislocations on the breaking stress of Lomer lock junctions in FCC metal alloys. He applied the same anisotropic line-tension model to a variety of dislocation junction configurations and found an inconsistency in the widely cited results of Dupuy and Fivel regarding the Hirth Lock, and he outlined an approach to extend these finding to HCP Mg alloys that we are still pursuing with discretionary fundings at UVA.

36 MATERIALS SCIENCE↗

Laser Diode Analysis and Verification (Professional Report)

The goal of this project at Lawrence Livermore National Laboratory (LLNL) was focused on analyzing and verifying laser diodes to ensure that the diodes meet requirements for LLNL mission applications. The different stages of this testing were documented in a process flow map which included location color coding so that each step was clearly defined in scope and location of appropriate testing facility. Data were collected and analyzed and compared to minimum acceptable values to see if requirements were met. The process map documents the initial receipt, inspection, and testing of the laser diodes. Initial inspections started with Keyence Microscope imaging and then moved on to High Potential, Ramp, and Burst Testing. Data from the diode testing were processed through MATLAB and Python codes to verify various metrics such as slope efficiency, threshold current, back irradiance, and beam divergence met requirements. These metrics were then recorded in Excel summary reports. Approximately 95% of the laser diodes passed all tests. Presentations were given to Lawrence Livermore’s internal leadership team, an external partner, and to a Lab-wide audience. The data released for this report was constrained by information protection considerations of LLNL’s national security missions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Preliminary Enhanced LAW Glass Formulation Algorithm

This report summarizes the Preliminary Enhanced Low-Activity Waste Glass Formulation Algorithm (GFA), its background information, and the calculations it performs. The Preliminary Enhanced Low-Activity Waste GFA is a tool developed in MATLAB to formulate glass for a given waste composition while attempting to maximize waste loading. It is intended for use at the Waste Treatment and Immobilization Plant, where nuclear waste will be vitrified into glass. The formulated glass is required to satisfy several processing and product quality constraints. In addition, calculations must account for associated uncertainties in constraint prediction and measurement. The GFA adheres to Pacific Northwest National Laboratory nuclear quality assurance procedures and has been validated and verified. In this second revision to the report, two constraints: the lithium and zirconia target constraints were amended to Table 3.1 and programmed into the algorithm. Additionally, three sections were added. Section 3.2.1 includes the Enhanced LAW Correlation Rules, an addition to the glass optimality criteria. Section 4.0 discusses a stepped process approach taken to ensure the algorithm can adequately react to differences between its recommended output and the actual measurements and transfers made during plant operation, and calculation of specific radionuclide concentrations. Section 5.2 describes changes to the user interaction with the algorithm program including the addition of a graphical user interface and compilation of the software.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Verification and Validation of the PLTEMP/ANL Code for Thermal-Hydraulic Analysis of Experimental and Test Reactors, Volume 1

This is Volume 1 of a two-volume document that collects the verification and validation (V&V) works done for the PLTEMP/ANL code during the years of its development and improvement. Volume 1 describes the V&V of sixteen capabilities of the PLTEMP/ANL code that were identified by research reactor analysts as frequently used in their thermal-hydraulic analysis. Volume 2 describes the V&V of developments and improvements since the release of PLTEMP/ANL Version 4.3. Each chapter of the document focuses on verifying or validating a specific part of the software that calculates a specific phenomenon, e.g., channel flow calculation, coolant property calculation, heat transfer calculation, and flow instability calculation. Software verification is performed by comparing the code with hand calculation, Microsoft spreadsheet calculation, Mathematica calculation, or MATLAB calculation. The software validation is done by comparing the code with experimental data or a widely tested code like the RELAP5 code. In addition, some PLTEMP/ANL V&V works that are available in the open literature are simply cited in Volume 1 of the document. PLTEMP/ANL has been used in the safety analysis reports of several US and foreign research reactors licensed and converted from highly enriched uranium fuel to low-enriched uranium fuel. A list of such reactors is given in Volume 1 of the document.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Verification and Validation of the PLTEMP/ANL Code for Thermal-Hydraulic Analysis of Experimental and Test Reactors, Volume 2

This is Volume 2 of a two-volume document that collects the verification and validation (V&V) works done for the PLTEMP/ANL code during the years of its development and improvement. Volume 1 describes the V&V of sixteen capabilities of the PLTEMP/ANL code that were identified by research reactor analysts as frequently used in their thermal-hydraulic analysis. Volume 2 describes the V&V of developments and improvements since the release of PLTEMP/ANL Version 4.3. Each chapter of the document focuses on verifying or validating a specific part of the software that calculates a particular phenomenon, e.g., channel flow calculation, coolant property calculation, heat transfer calculation, and flow instability calculation. Software verification is performed by comparing the code with a hand calculation, Microsoft spreadsheet calculation, Mathematica calculation, or MATLAB calculation. The software validation is achieved by comparing the code with experimental data or a widely tested code like the RELAP5 code. In addition, some PLTEMP/ANL V&V works that are available in the open literature are simply referenced in Volume 1 of the document. PLTEMP/ANL has been used in safety analysis reports of several US and foreign research reactors licensed and converted from highly enriched uranium fuel to low-enriched uranium fuel. A list of such reactors is given in Volume 1 of the document.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Simulation of the Photoionization of H- Together With the Subsequent Tracking of the Liberated Electrons

The Proton Improvement Plan - II (PIP-II) is a new linear accelerator (LINAC) complex being built at Fermilab. It is based on superconducting radiofrequency cavities and will accelerate H- ions to 800 MeV kinetic energy before injection into the existing Booster ring. Measurements of the profile of the beam along the LINAC must be done by non-intercepting methods due to the superconducting cavities. The method chosen is photoionization of a small number of H- by a focused infrared laser, aka laserwire. The number of ionized electrons is measured as a function of laser position within the H- beam. To aid in the design of the collection mechanism, a simulation was written in MATLAB with input from the commercial electromagnetic simulation, CST. This simulation calculates the number and positions of the liberated electrons and tracks them through the magnetic collection and H- beam fields to the collection point. Results from this simulation for various points along the LINAC will be shown.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A Tutorial for Generating Correlated Random Samples in the Context of Replica Cross Section Data Used in the Propagation of Uncertainty (Second Edition)

The following sample problem write-ups are designed as a tutorial for generating correlated random samples (e.g., replica multi-group cross section data) for use in propagation of uncertainty problems. These problems were set up and solved in MATLAB, but any programming environment with basic statistical functions can be used to generate similar results. Since linear sample problems were chosen, helpful comparisons to the “sandwich” rule propagation of uncertainty are available and were used.

97 MATHEMATICS AND COMPUTING↗

Augmenting the Structural Dynamics Workflow to Reduce Development Times [Slides]

The structural dynamics workflow at LANL stands to be improved by implementing methods and techniques related to reduced-order finite-element modeling (FEM), pretest predictions, and FEM-test correlation and updating. The benefits of these industry-standard practices are demonstrated in this work. Software used for this project included Ansys 2023 R2 (Mechanical and Spaceclaim), Creo Parametric 8.0, MathCAD Prime 8.0, Matlab R2023b.

42 ENGINEERING↗

Automated Calculation of U-10Mo Fuel and Zr Cladding Thickness

Hand calculation and verification of cladding and fuel thickness in a hot-isostatically pressed uranium-molybdenum (U-10Mo) alloy can incur human errors and longer image processing time when analyzing a large set of cross-sectional images. To help alleviate both time cost and errors made in hand measurements, an automated image processing procedure was developed using Octave, an open-source MATLAB alternative, to repeatedly determine the thickness of Zr and U-10Mo layers. U-10Mo specimens were imaged using standard secondary electron and backscattered-electron imaging at 250× magnification to capture the various layers present in U-10Mo. Further image processing used pixel-by-pixel calculation of the Zr and U-10Mo layers to gather quantitative statistics on the thickness variations associated with each of the layers.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Progress Report on Model Development for Aerosol Transport through Divergent Cracks Paths

This report summarizes the progress in developing a phenomenological model of aerosol transport, deposition, and plugging through microchannels. The purpose of this effort is to introduce to a user community—involving researchers, regulators, and industry—a generic, reliable numerical model for the prediction of aerosol transport while accounting for potential deposition and plugging of the leak paths to model spent nuclear fuel (SNF) release from postulated Stress Corrosion Cracks in canisters under storage or transportation. Current work focuses on expanding the model to predict aerosol (and gas) flow through complex microchannel (nozzle) geometries in addition to the rectangular and cylindrical geometries, as we approach more realistic Stress Corrosion Crack conditions. In this regard, a divergent nozzle geometry, used at Sandia National Laboratories (SNL) for testing aerosol release and retention (Durbin et al. 2021), was added to the model. The model was then validated with blowdown data for the particular microchannel, from experiments conducted at SNL. The report also presents preliminary non-benchmarked aerosol penetration fraction and mass flow rate for a monodisperse 10-micron (AED) particle concentration of 1.7e-08 kg/m 3 being released through the divergent microchannel geometry from the SNL’s aerosol experimental tank setup.Future work will involve validating the aerosol model and updating the model’s Graphical User Interface (GUI) to include the divergent nozzle geometry. This is expected to help stakeholders perform quick and easy first principles calculations without needing to understand the underlying MATLAB script.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Parametric and Sensitivity Analysis of a Steam Generator Model Using Python and Machine-Learning Tools

For this study, we used Python and machine-learning tools to perform a comprehensive parametric and sensitivity analysis on a steam generator (SG) model. (The Python model was based on a previously completed MATLAB framework for the Holtec SMR-160 SG.) We investigated the influence of various input parameters (e.g., heat transfer coefficient [HTC], Nusselt number, and heat exchanger effectiveness) on the system’s output. With machine-learning tools such as the Risk Analysis Virtual Environment (RAVEN), which was developed at Idaho National Laboratory, we were then able to perform an automated analysis of the SG inputs’ effect on the HTC. The analysis results give valuable insights into the performance and optimization of SG systems. We found the inlet mass flow rate (MFR) to have the greatest impact on the HTC, followed closely by the inlet temperature, and then pressure. Shifting of the input parameters causes the location of the maximum HTC along the SG length to change incrementally. The cold leg (CL) MFR was also found to impact the HTC magnitude as well as the location of the maximum HTC. At between 0.4–0.9 of the total SG length, the input parameters experience maximum impact on the HTC, leading us to suggest that sensors be efficiently placed on the SG so as to closely and effectively monitor thermal-hydraulic properties during reactor operation. We also found that the sensitivity data calculated manually agrees with the RAVEN – based data, confirming the same range of maximum sensitivity. However, the RAVEN-based analysis showed that cold leg pressure and hot leg temperature have a greater impact on the heat transfer coefficient than the mass flow rate, implying that a manual sensitivity study taking only two samples is not accurate.

20 FOSSIL-FUELED POWER PLANTS↗

The Digital Engineering Vision for DOME: Facilitating Design, Deployment, and Operations [Poster]

DOME is a planned microreactor test facility at INL’s Materials and Fuels Complex. It is a complex system with several interdependent sub-systems such as the reactor (up to 20 MWth), radioactive confinement, temperature and pressure regulation system, ventilation system, etc. The engineering design process for such a system traditionally involves several documents from various sources and the system information is scattered across these documents. Digital engineering represents a paradigm shift through which systems are designed using digital models and integrated data. The digital engineering vision for DOME utilizes a model-based systems engineering (MBSE) approach. The system architecture, physical components, control logic, and verification experiments are all designed using MathWorks MATLAB and Simulink. This hierarchical model can combine data from multiple sources at various levels of abstraction. It can be used to simulate the facility’s operations and to test the system using different sets of parameters. Its capabilities can be expanded by interfacing it with high-fidelity multi-physics models, risk analysis tools, etc. The same model can evolve into a digital twin that can monitor operations and conduct predictive analysis using real-time sensor data from the facility. The eventual goal of this effort is to transform the end-to-end engineering of nuclear facilities in every phase of their lifecycle, including design, deployment, and operations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Finite Element Analysis System Workflow Tools

A collection of MATLAB functions and class definitions called System Workflow Tools (SWFT) are available to semi-automate steps in the simulation process. Some of these steps are often simple and routine for smaller finite element models, but if done directly by an analyst can quickly become labor intensive, cumbersome, and error prone for larger, system level models. Some of SWFT’s capabilities demonstrated in this report includes writing Sierra input decks and processing Quantities of Interest (QOI) from results files. SWFT also writes scripts in order to utilize other software programs such as Cubit (separating system level CAD into subassemblies and components, creating nodesets and sidesets), DAKOTA (ensemble management), and ParaView (contour plots and animations). Detailed commands and workflows from mesh generation to report generation are provided as examples for analysts to utilize SWFT capabilities.

97 MATHEMATICS AND COMPUTING↗

Characterization of Lab-Grown Cracks for Aerosol Transmission Testing

Motivation: Determine the length and opening of two lab-grown cracks, designated as LT-14 and LT-28, representative of stress corrosion cracks in spent nuclear fuel dry storage casks to supplement future testing of gas and aerosol transport. Problem: The extreme aspect ratio of the crack length to opening requires that imaging occurs in stages with the results merged before final analysis. Method: High magnification (1500x) optical images of both sides of the two plates were acquired. 20x stitched images with LSCM were acquired, fully stitched along the length, and leveled with newly developed PLATES Method in MATLAB®. Conclusion for LT-14: Side 1 is 47.25 mm long and has 366 separate crack features with an average length of 23.50 µm and an average opening of 8.27 µm. Side 2 is 69.44 mm long and has 550 separate crack features with an average length of 81.63 µm and an average opening of 67.70 µm. Conclusion for LT-28: Side 1 is 71.95 mm long and has 1,127 separate crack features with an average length of 42.27 µm and an average opening of 10.31 µm. Side 2 is 74.88 mm long and has 520 separate crack features with an average length of 98.13 µm and an average opening of 14.99 µm. The adjacent crack on side 1 is 18.95 mm long and has 37 separate crack features with an average length of 17.46 µm and an average opening of 10.42 µm. The adjacent crack on side 2 is 26.40 mm long and has 55 separate crack features with an average length of 87.26 µm and an average opening of 48.29 µm. Each adjacent crack is approximately 26 mm from the main crack.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

SinhaRoy_TechPresentation_2024 [Slides]

DOME is a planned microreactor test facility at INL’s Materials and Fuels Complex. It is a complex system with several interdependent sub-systems such as the reactor (up to 20 MWth), radioactive confinement, temperature and pressure regulation system, ventilation system, etc. The engineering design process for such a system traditionally involves several documents from various sources and the system information is scattered across these documents. Digital engineering represents a paradigm shift through which systems are designed using digital models and integrated data. The digital engineering vision for DOME utilizes a model-based systems engineering (MBSE) approach. The system architecture, physical components, control logic, and verification experiments are all designed using MathWorks MATLAB and Simulink. This hierarchical model can combine data from multiple sources at various levels of abstraction. It can be used to simulate the facility’s operations and to test the system using different sets of parameters. Its capabilities can be expanded by interfacing it with high-fidelity multi-physics models, risk analysis tools, etc. The same model can evolve into a digital twin that can monitor operations and conduct predictive analysis using real-time sensor data from the facility. The eventual goal of this effort is to transform the end-to-end engineering of nuclear facilities in every phase of their lifecycle, including design, deployment, and operations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗