Search NASA⌕ Search

SEARCH · Search NASA

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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Analysis of Inter Converter Interactions using Harmonic State Space Modeling

Power converters have become a key participant within the modern grid, facilitating the grid's interaction with renewable energy sources and various loads, among others. These converters can also interact with each other, creating harmonics that were previously absent in the grid. While some previous research studies have focused on inter-converter interactions, not much focus has been on the role of pulse width modulators in these inter-converter interactions. This paper focuses on the issue and builds upon harmonic state-space modeling (HSS) to observe inter-converter harmonic interactions. A new modulator model is developed, and its implementation within HSS modeling to observe harmonic coupling between multiple inverters is highlighted. A simulation testbench is developed, and to verify the advantages and limitations of the proposed modeling approach, it is compared with two other traditional modeling methods: namely, average and switching models. The proposed method is shown to be significantly more accurate and informative than conventional average modeling methods and has considerably lower computation costs and simulation time than detailed switching models. The results are verified using extensive simulations on MATLAB.

Xue, Yaosuo↗
DOE OSTI2022

DEFECT DETECTION USING DYNAMIC ANALYSIS FOR ADDITIVE MANUFACTURED METALS

Additive manufacturing (AM) has the ability to produce parts with complex geometries and internal features, however, for demanding applications such as the automotive and aerospace industries, it is crucial that the parts can meet the demanding functional and geometric requirements. Quality control for AM parts focuses on nondestructive methods of testing, but many of the current methods are expensive and time-consuming. The research presented in this report explores various methods of nondestructive evaluation (NDE) using dynamic analysis on stainless steel parts produced with selective laser melting (SLM). Methods include, but are not limited to, frequency response functions (FRF), impedance-based measurements, and scanning laser doppler vibrometry. Additionally, mode shape analysis was performed in MATLAB and FEA simulations were used for comparison with experimental results. The results indicate that dynamic analysis has the potential to be a feasible method of defect detection and NDE in AM parts and future work should focus on refining these methods, such as optimizing test parameters to improve sensitivity to defects.

Deonarain, Gita↗
DOE OSTI2022

Recent Developments in the WEC-Sim Open-Source Design Tool: Preprint

WEC-Sim (Wave Energy Converter SIMulator) is an open-source code for simulating wave energy converters, which has been actively developed and applied to simulate a wide variety of device archetypes and has become a popular tool since its initial release in 2014. WEC-Sim is developed jointly by the National Renewable Energy Laboratory (NREL) and Sandia National Laboratories (SNL) within the MATLAB/SIMULINK environment. Figure 1 illustrates a general wave-to-wire model which begins with a deployment site resource characterization, which is used to complete the hydrodynamic simulation of a single WEC (or array), with the power generation profile imported to a grid simulator to understand the influence on the local electrical network. While modelling the entire wave-to-wire is difficult and encompass multiple time scales and physics, WEC-Sim is focused on the hydrodynamics simulation to predict, analyze and optimize WEC dynamics and power performance. WEC-Sim simulations are performed in the time domain based on the radiation and diffraction method using hydrodynamics coefficients derived from boundary element method (BEM) based-frequency-domain potential flow solvers (e.g., WAMIT, NEMOH, Capytaine, or ANSYS-AQWA). Within this level of modeling fidelity, WEC-Sim can handle floating body hydrodynamics, mechanical and electrical power generation methods, advanced control implementation, mooring systems, and other unique applications such as desalination. Table 1 lists additional WEC-Sim functionalities, which are created using prebuilt Simulink blocks and MATLAB scripts that can simulate a wide range of floating systems and the corresponding auxiliary subsystems.

hydrodynamics modeling↗
DOE OSTI2023

MATBOX, an Open-Source Microstructure Analysis Toolbox for Meshing, Generation, Segmentation, and Characterization of 3D Heterogenous Volumes

Battery performance is strongly correlated with electrode microstructural properties. To account for its impact, lithium-ion battery (LIB) models either abstract the microstructural heterogeneity of composite electrodes using effective macroscopic properties (macro- or meso- scale models) or directly solve the system of equations on the microstructure geometry or mesh (microstructure-scale models). Therefore, to be adequate, both families of models require information from the microstructure geometry, which can be provided by the numerical tool presented in this work. MATBOX is a MATLAB open-source application [1] developed by NREL for performing various microstructure-related tasks including microstructure numerical generation, image filtering and microstructure segmentation, microstructure characterization and correlation, visualization, and microstructure meshing. MATBOX was originally developed for the analysis of LIB electrode microstructures; however, the algorithms provided by the toolbox are widely applicable to other heterogeneous materials. The toolbox provides a user-friendly experience thanks to a Graphical-User Interface, requires no coding by the user, and is well documented. This presentation will illustrate various MATBOX features for the characterization of a LIB electrode, including a fully automated Representative Volume Element (RVE) analysis, the numerical generation of complex 'virtual' microstructure, including dual-layer electrodes and carbon-binder additive phase, and the meshing of a complex NMC/graphite full cell microstructure suitable for 3D finite-element modeling. Other modules (segmentation, visualization, and correlation) will be briefly presented. Thanks to its modular, open-source approach, MATBOX can easily incorporate third-party algorithms to eventually build a standard in the field that will benefit the whole scientific community. Effective diffusion coefficient [2], additive phase numerical generation [3], and meshing [4] third-party algorithms have been already integrated in the toolbox with more to come.

DIRECT ENERGY CONVERSION,MATHEMATICS AND COMPUTING↗
DOE OSTI2023

New Developments and Capabilities Within WEC-Sim: Preprint

WEC-Sim is an open-source software for simulating wave energy converters, which has been actively developed and applied since its initial release in 2014 to simulate a wide variety of device archetypes. WEC-Sim is developed jointly by the National Renewable Energy Laboratory (NREL) and Sandia National Laboratories (Sandia) within the MATLAB/SIMULINK environment. A general wave-to-wire model begins with a deployment site resource characterization, which is used to complete the hydrodynamic simulation of wave energy converters (WEC), with the power generation profile imported to a grid simulator to understand the influence on the local electrical network. While modeling the entire wave-to-wire is difficult and encompasses multiple time scales and physics, WEC-Sim is focused on the hydrodynamics simulation to predict, analyze, and optimize WEC dynamics and power performance. WEC-Sim simulations are performed in the time domain based on the radiation and diffraction method using hydrodynamics coefficients derived from boundary element method (BEM)-based frequency-domain potential flow solvers (e.g., WAMIT, NEMOH, Capytaine, or ANSYS-AQWA). With this level of modeling fidelity, WEC-Sim can handle floating body hydrodynamics, mechanical and electrical power generation methods, advanced control implementation, mooring systems, and other unique applications such as desalination. Additional WEC-Sim functionalities include pre-built Simulink blocks and MATLAB scripts that can simulate a wide range of floating systems and the corresponding auxiliary subsystems. The developers of WEC-Sim continue to release new versions of the software, at least annually, with our latest release in September 2022. These releases include bug fixes, updates to software documentation, as well as new features to expand WEC-Sim's capabilities to model a wide range of WEC concepts. This publication will highlight the new features added to WEC-Sim between versions 4.1.0 to 5.0.1 which spans over a two year period from June 2020 to September 2022. New features to be described will include topics such as continuous integration checks, revised Morison Element and nonlinear hydro implementations, run directly from Simulink (required for hardware-in-the-loop execution), BEMIO updates to import Capytaine BEM hydrodynamics, addition of cable blocks, and new wave visualization features.

TIDAL AND WAVE POWER↗
DOE OSTI2023

Assessment of the Phase-to-Ground Fault Apparent Admittance Method with Phase/Ground Boundaries to Detect Types of Electrical Faults for Protective Relays Using Signature Library and Simulated Events

Protective relays in electric power grids recognize the types of electrical faults in a few seconds. The most common detection method to detect the types of electrical faults is based on measuring the angle between the zero and negative sequence currents. However, it is not completely accurate because the phase to phase ground and phase to ground electrical faults could have the same detection conditions. Therefore, engineers need to plot the events after an electrical fault to observe the nature of the incidents in detail. In this study, the phase to ground fault apparent (PGFA) admittance method with phase/ground boundaries identified the types of electrical faults located in distribution power lines and feeders. This method was based on measuring the PGFA admittance magnitudes for the faulted and non-faulted phases, resulting in greater than zero and near zero, respectively. The PGFA admittance algorithm was built with MATLAB/ Simulink software and tested with signature library and grid simulation events. The PGFA method with phase/ground boundaries was evaluated with the confusion matrix. The measured and predicted values matched in more than 90% of the tests, and the PGFA admittance method with phase/ground boundaries presented an accuracy of 94.3% and a precision of 100%.

Piesciorovsky, Emilio↗
DOE OSTI2023

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.

43 PARTICLE ACCELERATORS↗
DOE OSTI2012

Computation of neutron fluxes in clusters of fuel pins arranged in hexagonal assemblies (2D and 3D)

For computations of fluxes, we have used Carvik's method of collision probabilities. This method requires tracking algorithms. An algorithm to compute tracks (in 2D and 3D) has been developed for seven hexagonal geometries with cluster of fuel pins. This has been implemented in the NXT module of the code DRAGON. The flux distribution in cluster of pins has been computed by using this code. For testing the results, they are compared when possible with the EXCELT module of the code DRAGON. Tracks are plotted in the NXT module by using MATLAB, these plots are also presented here. Results are presented with increasing number of lines to show the convergence of these results. We have numerically computed volumes, surface areas and the percentage errors in these computations. These results show that 2D results converge faster than 3D results. The accuracy on the computation of fluxes up to second decimal is achieved with fewer lines. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2013

Spectroscopic and physicochemical measurements for on-line monitoring of used nuclear fuel separation processes

Separation processes for used nuclear fuel are often complicated and challenging due to the high constraints in purity of the products and safeguards of the process streams. In order to achieve a safe, secure and efficient separation process, the liquid streams in the separation process require close monitoring. Due to the high radiation environment, sampling of the materials is difficult. Availability of a detection technique that is remote, non-destructive and can avoid time-delay caused by retrieving samples would be beneficial and could minimize the exposure to personnel and provide material accountancy to avoid diversion (non-proliferation). For example, Ultra Violet (UV), Visible (Vis), Near-Infrared (NIR) and Raman spectroscopy that detect and quantify elements present in used nuclear fuel, e.g. lanthanides, actinides and molecules such as nitrate, can be used. In this work, we have carried out NIR and Raman spectroscopy to study aqueous solutions composed of different concentrations of nitric acid, sodium nitrate, and neodymium at varied temperatures. A chemometric model for online monitoring based on the PLS-Toolbox (MATLAB) software has been developed and validated to provide chemical composition of process streams based on spectroscopic data. In conclusion, both of our NIR and Raman spectra were useful for H{sup +} and NO{sub 3} prediction, and only NIR was helpful for the Nd{sup 3+} prediction.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2015

{sup 6}LiF:ZnS(Ag) Neutrons Scintillator Detector Configuration for Optimal Readout

A Chromatic Analysis Neutron Diffractometer Or Reflectometer (CANDOR) is under development at the NIST Center for Neutron Research (NCNR). The CANDOR neutron sensor will rely on scintillator material for detecting the neutrons scattered by the sample under test. It consists of {sup 6}LiF:ZnS(Ag) scintillator material into which wavelength shifting (WLS) fibers have been embedded. Solid state photo-sensors (silicon photomultipliers) coupled to the WLS fibers are used to detect the light produced by the neutron capture event ({sup 6}Li (n,α) {sup 3}H reaction) and ionization of the ZnS(Ag). This detector configuration has the potential to accomplish the CANDOR performance requirements for efficiency of 90% for 5 A (3.35 meV) neutrons with high gamma rejection (10{sup 7}) along with compact design, affordable cost and materials availability. However this novel design includes challenges for precise neutron detection. The recognizing of the neutron signature versus the noise event produce by gamma event cannot be easy overcome by pulse height discrimination obstacle as can be achieved with {sup 3}He gas tube. Furthermore the selection of silicon photomultipliers (SiPM) as the light sensor maintains the obstacle of dark noise that does not exist when a photomultiplier tube is coupled to the scintillator. A proper selection of SiPM should focus on increasing the output signal and reducing the dark noise in order to optimize the detection sensitivity and to provide a clean signal pulse shape discrimination. The main parameters for evaluation are: - Quantum Efficiency (QE) - matching the SiPM peak QE with the peak transmission wavelength emission of the WLS. - Recovery time - a short recovery time is preferred to minimize the pulse width beyond the intrinsic decay time of the scintillator crystal (improves the gamma rejection based output pulse shape (time)). - Diode dimensions -The dark noise is proportional to the diode active area while the signal is provided by the WLS fibers; therefore the diode area should ideally be only minimally larger than fiber bundle area. - Low dark noise - it is desirable to minimize the dark noise during the pulse integration period so as to minimize the background for pulse shape discrimination. - Photon Detection Efficiency - it is desirable to increase the SiPM PDE in order to enhance light collection. This will increase the likelihood of detecting neutron events with lower light production and will present a cleaner raw signal for pulse shape discrimination. We will present the SiPM optimization process and studies of dark noise and gamma and neutron sensitivity as a function of bias voltage and operating temperature that have enabled us to optimize the detector sensitivity and gamma rejection. The gamma rejection performance goal requires to overcome the challenge of discriminating between the light signature accepted by neutron event to the one received by the noise. In addition there is a huge variation between the number of light photons that reaching the WLS fibers for different neutron events caused by the heavy ions energy losses prior to ionizing the ZnS(Ag) and the high light attenuation of the scintillation mixture. This variation in the light signal along with the long decay time of the ZnS(Ag) (tens of microseconds) can cause double counting of the same neutron event in the case of high light output signature or preventing the detection of low sequential light output signature neutron event. We will presents the algorithm developed for {sup 6}LiF:ZnS(Ag) sensor readout and the results achieved by an off-line analysis by Matlab software code that successfully achieved both the high gamma rejection with a sensitive and accurate neutron event detection. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Advanced Features and Recent Developments in the WEC-Sim Open-Source Design Tool: Preprint

WEC-Sim is an open-source tool developed in Matlab/Simulink to simulate wave energy converters. WEC-Sim has been developed jointly by the National Renewable Energy Laboratory (NREL) and Sandia National Laboratories (SNL) since its initial release in 2014. WEC-Sim simulations are performed in the time domain based on the radiation and diffraction method using hydrodynamics coefficients derived from boundary element method (BEM)-based frequency-domain potential flow solvers (e.g., WAMIT, NEMOH, Capytaine, or ANSYS-AQWA). Within this level of modeling fidelity, WEC-Sim can handle floating body hydrodynamics, mechanical and electrical power generation methods, advanced control implementation, mooring systems, and other unique applications such as desalination. WEC-Sim functionalities are developed using Simulink blocks and Matlab functions. WEC-Sim is capabilities include simulation of different wave conditions, constraints and PTOs, body dynamics, mooring systems, among other advanced features such as nonlinear hydrodynamics, nonlinear PTOs, body-to-body interactions, etc.

hydrodynamics modeling↗
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↗