Results for “simulation output analysis” — Page 83

DOE OSTI2022

Estimating List-Mode Data Sensitivities to Nuclear Data with MCNP6

Nuclear data are a vital component of predictive simulations used in applications like experiment design, stockpile stewardship, nuclear nonproliferation/safeguards, health physics, and criticality safety. A singular simulation requires the coalescence of different areas of nuclear data such as cross sections, angular distributions, and energy distributions of emitted neutrons for different materials and energy ranges. Improving nuclear data and thus reducing the uncertainty in simulated parameters could enable smaller, better-informed safety factors and ultimately reduce operational and procedural costs. There is a constant effort to garner a better understanding of the physical quantities represented by nuclear data through experiments. Integral experiment benchmarks use simulated and measured results to validate current nuclear data values. In the past, benchmarks primarily focused on the effective multiplication factor (k eff ); however, this limited scope has caused compensating errors and areas of nuclear data that lack validation. Compensating errors are inaccuracies in nuclear data that are obfuscated by cancellation when observing integrated values such as k eff . Diverse integral benchmark experiments that look for quantities of interest other than k eff and include multiple responses minimize the possibility of compensating errors and provides validation to areas of nuclear data previously lacking experimental validation. Benchmark experiments can be optimized during the design process to be highly dependent on specific areas of nuclear data. The dependence of a response in an experiment to a specific area/type of nuclear data is defined as sensitivity. A larger sensitivity means that nuclear data uncertainties will play a larger role in the response(s) resulting in larger bias. Currently, the sensitivity capabilities of the Monte Carlo N-Particle (MCNP ®1 ) transport code are limited to responses of k eff and tallied values (e.g., flux, surface current). As a part of the EUCLID project, this work explores estimating list-mode nuclear data sensitivities that can be used to design experiments aimed to constrain and reduce compensating errors in nuclear data by focusing on responses other than k eff . Tallied values are ideal quantities that are estimated with detectors during experiments. List-mode data (a list of neutron collection times) are the direct output of detector systems in subcritical neutron noise experiments. Expanding MCNP sensitivity capabilities to include the sensitivity of responses estimated from list-mode data, such as the prompt neutron decay constant (α) and multiplicity estimates (S and D), enables more direct comparison of simulated and measured experimental quantities. Additionally, deterministic tools such as SENSMG are capable of obtaining sensitivities to a wide variety of responses; however, these tools cannot handle complex geometries due to the assumptions made in discretizing the phase-space variables of the Boltzman transport equation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2016

Global sensitivity analysis in stochastic simulators of uncertain reaction networks

Stochastic models of chemical systems are often subjected to uncertainties in kinetic parameters in addition to the inherent random nature of their dynamics. Uncertainty quantification in such systems is generally achieved by means of sensitivity analyses in which one characterizes the variability with the uncertain kinetic parameters of the first statistical moments of model predictions. In this work, we propose an original global sensitivity analysis method where the parametric and inherent variability sources are both treated through Sobol’s decomposition of the variance into contributions from arbitrary subset of uncertain parameters and stochastic reaction channels. The conceptual development only assumes that the inherent and parametric sources are independent, and considers the Poisson processes in the random-time-change representation of the state dynamics as the fundamental objects governing the inherent stochasticity. Here, a sampling algorithm is proposed to perform the global sensitivity analysis, and to estimate the partial variances and sensitivity indices characterizing the importance of the various sources of variability and their interactions. The birth-death and Schlögl models are used to illustrate both the implementation of the algorithm and the richness of the proposed analysis method. The output of the proposed sensitivity analysis is also contrasted with a local derivative-based sensitivity analysis method classically used for this type of systems.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2015

Experimental Investigation of Turbine Vane Heat Transfer for Alternative Fuels

The focus of this program was to experimentally investigate advanced gas turbine cooling schemes and the effects of and factors that contribute to surface deposition from particulate matter found in coal syngas exhaust flows on turbine airfoil heat transfer and film cooling, as well as to characterize surface roughness and determine the effects of surface deposition on turbine components. The program was a comprehensive, multi-disciplinary collaborative effort between aero-thermal and materials faculty researchers and the Department of Energy, National Energy Technology Laboratory (NETL). The primary technical objectives of the program were to evaluate the effects of combustion of syngas fuels on heat transfer to turbine vanes and blades in land-based power generation gas turbine engines. The primary questions to be answered by this investigation were; What are the factors that contribute to particulate deposition on film cooled gas turbine components? An experimental program was performed in a high-temperature and pressure combustion rig at the DOE NETL; What is the effect of coal syngas combustion and surface deposition on turbine airfoil film cooling? Deposition of particulate matter from the combustion gases can block film cooling holes, decreasing the flow of the film coolant and the film cooling effectiveness; How does surface deposition from coal syngas combustion affect turbine surface roughness? Increased surface roughness can increase aerodynamic losses and result in decreased turbine hot section efficiency, increasing engine fuel consumption to maintain desired power output. Convective heat transfer is also greatly affected by the surface roughness of the airfoil surface; Is there any significant effect of surface deposition or erosion on integrity of turbine airfoil thermal barrier coatings (TBC) and do surface deposits react with the TBC in any way to decrease its thermal insulating capability? Spallation and erosion of TBC is a persistent problem in modern turbine engines; and What advancements in film cooling hole geometry and design can increase effectiveness of film cooling in turbines burning high-hydrogen coal syngas due to the higher heat loads and mass flow rates of the core flow? Experimental and numerical investigations of advanced cooling geometries that can improve resistance to surface deposition were performed. The answers to these questions were investigated through experimental measurements of turbine blade surface temperature and coolant coverage (via infrared camera images and thermocouples) and time-varying surface roughness in the NETL high-pressure combustion rig with accelerated, simulated surface deposition and advanced cooling hole concepts, coupled with detailed materials analysis and characterization using conventional methods of Scanning Electron Microscopy (SEM), Transmission Electron Microscopy (TEM), X-Ray Diffraction (XRD), 3-D Surface Topography (using a 3-D stylus profilometer). Detailed surface temperatures and cooling effectiveness could not be measured due to issues with the NETL infrared camera system. In collaboration with faculty startup funding from the principal investigator, experimental and numerical investigations were performed of an advanced film cooling hole geometry, the anti-vortex hole (AVH), focusing on improving cooling effectiveness and decreasing the counter-rotating vortex of conventional cooling holes which can entrain mainstream particulate matter to the surface. The potential benefit of this program is in gaining a fundamental understanding of how the use of alternative fuels will effect the operation of modern gas turbine engines, providing valuable data for more effective cooling designs for future turbine systems utilizing alternative fuels.

10 SYNTHETIC FUELS↗

DOE OSTI2016

Transient dynamics of terrestrial carbon storage: Mathematical foundation and numeric examples

Terrestrial ecosystems absorb roughly 30% of anthropogenic CO 2 emissions since preindustrial era, but it is unclear whether this carbon (C) sink will endure into the future. Despite extensive modeling, experimental, and observational studies, what fundamentally determines transient dynamics of terrestrial C storage under climate change is still not very clear. Here we develop a new framework for understanding transient dynamics of terrestrial C storage through mathematical analysis and numerical experiments. Our analysis indicates that the ultimate force driving ecosystem C storage change is the C storage capacity, which is jointly determined by ecosystem C input (e.g., net primary production, NPP) and residence time. Since both C input and residence time vary with time, the C storage capacity is time-dependent and acts as a moving attractor that actual C storage chases. The rate of change in C storage is proportional to the C storage potential, the difference between the current storage and the storage capacity. The C storage capacity represents instantaneous responses of the land C cycle to external forcing, whereas the C storage potential represents the internal capability of the land C cycle to influence the C change trajectory in the next time step. The influence happens through redistribution of net C pool changes in a network of pools with different residence times. Furthermore, this and our other studies have demonstrated that one matrix equation can exactly replicate simulations of most land C cycle models (i.e., physical emulators). As a result, simulation outputs of those models can be placed into a three-dimensional (3D) parameter space to measure their differences. The latter can be decomposed into traceable components to track the origins of model uncertainty. Moreover, the emulators make data assimilation computationally feasible so that both C flux- and pool-related datasets can be used to better constrain model predictions of land C sequestration. We also propose that the C storage potential be the targeted variable for research, market trading, and government negotiation for C credits.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Recovery of Natural Gas Equipment Emissions into Gas Compression Engines for the Reduction of Potential Greenhouse Gas Emissions

Since the turn of the millennium, the United States (U.S.) oil and natural gas (ONG) industry has nearly doubled its natural gas production rate. As a result, the ONG industry has recently come under increasing scrutiny for its contributions to greenhouse gas (GHG) emissions. Consequently, various solutions to this problem have been proposed and formulated to reduce the impacts of GHG emissions on the environment. West Virginia University (WVU) have found it important to research the impacts of recovering vented gas streams into prime-mover engines. The U.S. Department of Energy (DOE) and National Energy Technology Laboratory (NETL) have granted WVU funding to research and develop a “Methane Mitigator” (M2) - a “Scalable Vent Mitigation Strategy to Simultaneously Reduce Methane Emissions and Fuel Consumption from the Compression Industry.” One of the main areas of interest for this research was the collection of emissions from natural gas equipment into a Caterpillar G3508J natural gas compression engine. The parameters being analyzed from the engine were brake-specific emissions and power output. The emissions sources considered for this research were pneumatic controllers (PCs), reciprocating compressor vents, and the engine’s open crankcase breather. The compressor vent and PC emissions were simulated using a mass flow controller (MFC) and flowed into the engine using two separate methods: (1) directly into the air intake, and (2) through a retrofitted closed crankcase ventilation system (CCV), serving as a buffer volume. The crankcase emissions were quantified without the CCV, and the impact on exhaust emissions from circulating the crankcase gases into the intake was measured. The simulated compressor vent and PC flows from the MFC had limited effect on the steady state operation of the engine and resulting performance. When the simulated flows were fed directly into the engine’s air intake, the changes within the engine’s continuous performance and emission parameters were larger but lasted for shorter durations. Conversely, when the simulated flows were fed into the CCV before entering the air intake, the changes in the engine’s performance and emission parameters were less pronounced for continuous analysis but lasted for longer durations. In either case, the continuous emission changes in both emissions and performance varied in size depending on the test scenario being run, but the cycle average changes in emissions and performance showed little impact overall compared to the engine’s baseline operation. As a result, the inclusion of a CCV shows a decrease in baseline carbon dioxide equivalent (CO2-eq.) engine emissions (from combined exhaust and open crankcase) of almost 4%. Likewise, the CCV inclusion reduced baseline total methane (CH4) from combined exhaust and open crankcase by upwards of 16%. These atmospheric emissions only decreased further with the inclusions of collected PC and compressor vent flows. The resulting changes in time-averaged rated exhaust behavior (or lack thereof) prove that the proposed M2 system could likely be deployed at sites with modern lean-burn natural gas engines as a viable option for reducing and eliminating potential GHG sources that would have otherwise been unutilized as energy sources.

03 NATURAL GAS↗

DOE OSTI2016

TH-B-BRC-00: How to Identify and Resolve Potential Clinical Errors Before They Impact Patients Treatment: Lessons Learned

Radiation treatment consists of a chain of events influenced by the quality of machine operation, beam data commissioning, machine calibration, patient specific data, simulation, treatment planning, imaging and treatment delivery. There is always a chance that the clinical medical physicist may make or fail to detect an error in one of the events that may impact on the patient’s treatment. In the clinical scenario, errors may be systematic and, without peer review, may have a low detectability because they are not part of routine QA procedures. During treatment, there might be errors on machine that needs attention. External reviews of some of the treatment delivery components by independent reviewers, like IROC, can detect errors, but may not be timely. The goal of this session is to help junior clinical physicists identify potential errors as well as the approach of quality assurance to perform a root cause analysis to find and eliminate an error and to continually monitor for errors. A compilation of potential errors will be presented by examples of the thought process required to spot the error and determine the root cause. Examples may include unusual machine operation, erratic electrometer reading, consistent lower electron output, variation in photon output, body parts inadvertently left in beam, unusual treatment plan, poor normalization, hot spots etc. Awareness of the possibility and detection of error in any link of the treatment process chain will help improve the safe and accurate delivery of radiation to patients. Four experts will discuss how to identify errors in four areas of clinical treatment. D. Followill, NIH grant CA 180803.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

TH-B-BRC-01: How to Identify and Resolve Potential Clinical Errors

Radiation treatment consists of a chain of events influenced by the quality of machine operation, beam data commissioning, machine calibration, patient specific data, simulation, treatment planning, imaging and treatment delivery. There is always a chance that the clinical medical physicist may make or fail to detect an error in one of the events that may impact on the patient’s treatment. In the clinical scenario, errors may be systematic and, without peer review, may have a low detectability because they are not part of routine QA procedures. During treatment, there might be errors on machine that needs attention. External reviews of some of the treatment delivery components by independent reviewers, like IROC, can detect errors, but may not be timely. The goal of this session is to help junior clinical physicists identify potential errors as well as the approach of quality assurance to perform a root cause analysis to find and eliminate an error and to continually monitor for errors. A compilation of potential errors will be presented by examples of the thought process required to spot the error and determine the root cause. Examples may include unusual machine operation, erratic electrometer reading, consistent lower electron output, variation in photon output, body parts inadvertently left in beam, unusual treatment plan, poor normalization, hot spots etc. Awareness of the possibility and detection of error in any link of the treatment process chain will help improve the safe and accurate delivery of radiation to patients. Four experts will discuss how to identify errors in four areas of clinical treatment. D. Followill, NIH grant CA 180803.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

TH-B-BRC-02: Troubleshooting Issues with Linear Accelerator

Radiation treatment consists of a chain of events influenced by the quality of machine operation, beam data commissioning, machine calibration, patient specific data, simulation, treatment planning, imaging and treatment delivery. There is always a chance that the clinical medical physicist may make or fail to detect an error in one of the events that may impact on the patient’s treatment. In the clinical scenario, errors may be systematic and, without peer review, may have a low detectability because they are not part of routine QA procedures. During treatment, there might be errors on machine that needs attention. External reviews of some of the treatment delivery components by independent reviewers, like IROC, can detect errors, but may not be timely. The goal of this session is to help junior clinical physicists identify potential errors as well as the approach of quality assurance to perform a root cause analysis to find and eliminate an error and to continually monitor for errors. A compilation of potential errors will be presented by examples of the thought process required to spot the error and determine the root cause. Examples may include unusual machine operation, erratic electrometer reading, consistent lower electron output, variation in photon output, body parts inadvertently left in beam, unusual treatment plan, poor normalization, hot spots etc. Awareness of the possibility and detection of error in any link of the treatment process chain will help improve the safe and accurate delivery of radiation to patients. Four experts will discuss how to identify errors in four areas of clinical treatment. D. Followill, NIH grant CA 180803.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

TH-B-BRC-03: Measurement and Calibration Discrepancy: Advanced Warning Guidance

Radiation treatment consists of a chain of events influenced by the quality of machine operation, beam data commissioning, machine calibration, patient specific data, simulation, treatment planning, imaging and treatment delivery. There is always a chance that the clinical medical physicist may make or fail to detect an error in one of the events that may impact on the patient’s treatment. In the clinical scenario, errors may be systematic and, without peer review, may have a low detectability because they are not part of routine QA procedures. During treatment, there might be errors on machine that needs attention. External reviews of some of the treatment delivery components by independent reviewers, like IROC, can detect errors, but may not be timely. The goal of this session is to help junior clinical physicists identify potential errors as well as the approach of quality assurance to perform a root cause analysis to find and eliminate an error and to continually monitor for errors. A compilation of potential errors will be presented by examples of the thought process required to spot the error and determine the root cause. Examples may include unusual machine operation, erratic electrometer reading, consistent lower electron output, variation in photon output, body parts inadvertently left in beam, unusual treatment plan, poor normalization, hot spots etc. Awareness of the possibility and detection of error in any link of the treatment process chain will help improve the safe and accurate delivery of radiation to patients. Four experts will discuss how to identify errors in four areas of clinical treatment. D. Followill, NIH grant CA 180803.

60 APPLIED LIFE SCIENCES↗

NASA NTRS2010

Utilizing In Situ Directional Hyperspectral Measurements to Validate Bio-Indicator Simulations for a Corn Crop Canopy

Two radiative transfer canopy models, SAIL and the two-layer Markov-Chain Canopy Reflectance Model (MCRM), were coupled with in situ leaf optical properties to simulate canopy-level spectral band ratio vegetation indices with the focus on the photochemical reflectance index in a cornfield. In situ hyperspectral measurements were made at both leaf and canopy levels. Leaf optical properties were obtained from both sunlit and shaded leaves. Canopy reflectance was acquired for eight different relative azimuth angles (psi) at three different view zenith angles (Theta (sub v)), and later used to validate model outputs. Field observations of photochemical reflectance index (PRI) for sunlit leaves exhibited lower values than shaded leaves, indicating higher light stress. Canopy PRI expressed obvious sensitivity to viewing geometry, as a function of both Theta (sub v) and psi . Overall, simulations from MCRM exhibited better agreements with in situ values than SAIL. When using only sunlit leaves as input, the MCRM-simulated PRI values showed satisfactory correlation and RMSE, as compared to in situ values. However, the performance of the MCRM model was significantly improved after defining a lower canopy layer comprised of shaded leaves beneath the upper sunlit leaf layer. Four other widely used band ratio vegetation indices were also studied and compared with the PRI results. MCRM simulations were able to generate satisfactory simulations for these other four indices when using only sunlit leaves as input; but unlike PRI, adding shaded leaves did not improve the performance of MCRM. These results support the hypothesis that the PRI is sensitive to physiological dynamics while the others detect static factors related to canopy structure. Sensitivity analysis was performed on MCRM in order to better understand the effects of structure related parameters on the PRI simulations. Leaf area index (LAI) showed the most significant impact on MCRM-simulated PRI among the parameters studied. This research shows the importance of hyperspectral and narrow band sensor studies, and especially the necessity of including the green wavelengths (e.g., 531 nm) on satellites proposing to monitor carbon dynamics of terrestrial ecosystems.

Cheng, Yen-Ben↗

DOE OSTI2022

Biorefinery upgrading of herbaceous biomass to renewable hydrocarbon fuels, Part 1: Process modeling and mass balance analysis

Process simulation has long been a well-established tool to track key operational, design, and mass and energy balance metrics for pre-commercial technologies such as advanced lignocellulosic biofuels. While tools such as this are well-documented in the public literature around 2nd-generation cellulosic ethanol technologies (which have been scaled up to commercial deployment to some degree over the past decade), such models and analysis information remain more sparse for more complex biorefinery pathways focused on producing drop-in hydrocarbon fuels and blend-stocks, particularly regarding information required to support air emissions or other environmental analysis. In this work, we summarize key details for an established "design case" modeling the conversion of herbaceous lignocellulosic biomass into a renewable diesel hydrocarbon blend-stock based on a representative lipid pathway from oleaginous yeast. The process is based on a biochemical deconstruction and upgrading approach utilizing deacetylation and dilute acid pretreatment, followed by enzymatic hydrolysis, fermentation, and catalytic upgrading of hydrolysate sugars to fuels. We provide key mass and energy balance outputs from the process models, with accompanying stream tables and component-level flowrates. A total of 12 model scenarios are presented spanning two feedstocks, three biorefinery scales, and two processing approaches for the lignin/residual solids waste streams. This "Part 1" manuscript presents the resulting impacts across the 12 cases on fuel yields and key output streams, focused here on direct biorefinery air emissions for selected components including CO 2 as well as sulfur (SOx) and nitrogen oxides (NOx). In the context of cleaner production, the latter focus on selected biorefinery air emission outputs establishes an initial baseline estimate and accompanying framework of the model cases, upon which an accompanying "Part 2" study will build to refine the values for these and other air pollutants across these scenarios, also considering mitigation options to comply with applicable regulatory standards. We also highlight further optimization opportunities based on potential tradeoffs identified here between air emissions versus life-cycle greenhouse gas profiles attributed to the disposition of lignin/residual solids.

09 BIOMASS FUELS↗

DOE OSTI2019

An integrated analytical and experimental study of contact acoustic nonlinearity at rough interfaces of fatigue cracks

This paper presents an analytical model that describes the relationship between interface stiffness and separation of rough interfaces such as fatigue cracks. The contact acoustic nonlinearity at the interface is simulated by a quasi-static model based on Hertzian contact theory. The model is validated using the results of dynamic acoustoelastic testing (DAET) with a Rayleigh wave probe on fatigue cracks in two aluminum alloy samples. One novel aspect of this work is that all the required geometrical parameters for the model are acquired directly from the aperture profile of real cracks extracted from their scanning electron microscopy (SEM) images. Also, the separation of the crack faces during dynamic perturbation is independently measured using a 3D laser Doppler vibrometer (LDV). In addition, using DAET allows for an unprecedented direct comparison between experimental and analytical results. Overall, this is because unlike conventional nonlinear ultrasonic tests that are based on measuring the amplitude of higher harmonics, DAET outputs directly the strain-dependency of transmission and time delay of ultrasonic waves propagating across the interface. The model is found in good qualitative agreement with the experimental results, although the predictions tend to underestimate the variation of transmission coefficient and time delay. We conduct a sensitivity analysis to investigate the influence of different assumptions and simplifications on model predictions.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Predicting defects in laser powder bed fusion using in-situ thermal imaging data and machine learning

Variation in the local thermal history during the Laser Powder Bed Fusion (LPBF) process in Additive Manufacturing (AM) can cause micropore defects, which add to the uncertainty of the mechanical properties (e.g., fatigue life, tensile strength) of the built materials. In-situ sensing has been proposed for monitoring the AM process to minimize defects, but successful minimization requires establishing a quantitative relationship between the sensing data and the porosity, which is particularly challenging with a large number of variables (e.g., laser speed, power, scan path, powder property). Physics-based modeling can simulate such an in-situ sensing-porosity relationship, but it is computationally costly. In this work, we develop Machine Learning (ML) models that can use in-situ thermographic data to predict the micropore of LPBF stainless steel materials. This work considers two identified key features from the thermal histories: the time above the apparent melting threshold ($\tau$) and the maximum radiance ( T max ). These features are computed, stored for each voxel in the built material, and then used as inputs. The binary state of each voxel, either defective or normal, is the output. Different ML models are trained and tested for the binary classification task. In addition to using the thermal features of each voxel to predict its own state, the thermal features of neighboring voxels are also included as inputs. This is shown to improve the prediction accuracy, which is consistent with thermal transport physics around each voxel contributing to its final state. Among the models trained, the F1 scores on test sets reach above 0.96 for Random Forests. Feature importance analysis based on the ML models shows that T max is more important to the voxel state than $\tau$ . The analysis also finds that the thermal history of the voxels above the present voxel is more influential than those beneath it. Our study significantly extends the capability of using in-situ thermographic data to predict porosity in LPBF materials. Finally, since ML models are fast, they may play integral roles in the optimization and control of such AM technologies.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Approach and Model Used to Represent a Timeline Analysis for Security Design Enhancements

Next-generation reactors will be able to use risk to inform and performance base the licensing of many aspects of the reactor, facility, and site design, including attributes of physical security. There are several factors related to security, including site topography, reactor design, and physical protection system components, to consider when designing the physical protection system into the overall facility design and plan of operation. With the versatility of advanced reactors, especially micro reactors, methods are needed to simplify and quickly evaluate potential timelines for designing a site configuration. This report describes an approach to generate qualitative and quantitative insights using a risk-informed simulation. The modeling process is described in detail, focusing on three aspects: (1) the facility mission time (the time required to control the plant until safe), (2) the attacker timeline (the time to potential sabotage), and (3) the response timeline (the time to counter the facility attack). The modeling capabilities also are extended to include facility physical phenomena such as thermal-hydraulics and heat transfer to capture realistic representation of dynamic changes to a facility. While the plant models and examples are hypothetical and do not represent a real facility, these modeling approaches could be used for future security-by-design engineering in advanced reactors. The outputs and insights from the modeling approach may be used to modify and optimize the security posture of a facility by efficiently making modifications to the model and seeing the overall impact from the modification. Lastly, use of the approach described in this report can also provide the technical basis for a physical protection program, describing how the facility and security strategy will cope with off-normal events.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2026

PowerAnalytics.jl: User-Centric Power Systems Analysis in Julia

The National Laboratory of the Rockies recently released version 1 of PowerAnalytics.jl, an analysis module for the outputs of its popular open-source electrical power systems modeling platform Sienna. It features an extensible framework - based on the flexible selecting of components, the execution of arbitrary metrics on them, and a familiar DataFrames-based output interface with embedded metadata - to process results in the Sienna style while keeping the interface as simple as possible for non-Julia experts. Here, I describe the package and where it fits into the Sienna ecosystem, how I harnessed user-centered design and Julia features to achieve beginner friendliness without sacrificing performance and expressibility, and what lessons might be drawn from the package's design and implementation.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1999

Automated Testing Experience of the Linear Aerospike SR-71 Experiment (LASRE) Controller

System controllers must be fail-safe, low cost, flexible to software changes, able to output health and status words, and permit rapid retest qualification. The system controller designed and tested for the aerospike engine program was an attempt to meet these requirements. This paper describes (1) the aerospike controller design, (2) the automated simulation testing techniques, and (3) the real time monitoring data visualization structure. Controller cost was minimized by design of a single-string system that used an off-the-shelf 486 central processing unit (CPU). A linked-list architecture, with states (nodes) defined in a user-friendly state table, accomplished software changes to the controller. Proven to be fail-safe, this system reported the abort cause and automatically reverted to a safe condition for any first failure. A real time simulation and test system automated the software checkout and retest requirements. A program requirement to decode all abort causes in real time during all ground and flight tests assured the safety of flight decisions and the proper execution of mission rules. The design also included health and status words, and provided a real time analysis interpretation for all health and status data.

Larson, Richard R.↗

DOE OSTI2022

Studsvik R2 materials test reactor Ad Hoc depletion strategy for the derivation of the fuel isotopic composition of the MPCMIV benchmark

The Ad Hoc Depletion Strategy elaborated by the NINE company, developed in support of the organization of the MPCMIV (Multi-physics Pellet Cladding Mechanical Interaction Validation) benchmark input and output specifications, is presented. This work aims at illustrating the strategy itself and then showing the results obtained with its application over the Studsvik R2 Testing Reactor, which is analyzed in the benchmark. The objective of the application of the strategy is to compute the fuel elements isotopic compositions at the beginning of some core loadings of interest for the benchmark. To this objective, it is necessary to implement first the simulation model of the three single assembly types and perform the infinite lattice depletions, then, to build the full core model and to perform the simulation of the core cycle. All the models and simulations were carried out with the use of the Monte Carlo particle transport code Serpent 2. Finally, the simulations results are assessed against Studsvik isotopic compositions of the fuel elements discharged from the R2 Reactor at the end of the core loading. Several assumptions were necessary during all the steps of the strategy, to overcome the lack of information regarding the core management. For this reason, the solution found at the end of the current analysis may not be completely optimized and further improvements regarding the model assumptions will be tested in a future work. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2013

Preliminary Scaling and controls Analysis of an FHR-HTSE System Idaho National Laboratory Summer 2013 Final Report

For new nuclear reactor system designs to be approved by regulatory agencies like the Nuclear Regulatory Commission (NRC), the details of system operation must be validated with respect to standards of safety, control, and output. A scaled experiment that replicates certain properties of the system can be used to validate compliance with regulatory standards, while avoiding the prohibitive cost and labor required to develop a fully functional prototype system; therefore, designing such an experiment is of special interest to current efforts to develop hybrid energy systems (HES) that integrate small modular reactors (SMRs), renewable energy systems, and industrial process applications such as hydrogen production and desalination. In addition, a scaled experiment can be an economical method of analyzing the interconnections between HES components and understanding the time constants associated between inter-component energy and information flows. This report discusses the results of a preliminary scaling analysis done for the primary loop of a 300 MWth Fluoride-Salt-Cooled High Temperature Reactor (FHR) that is coupled with a High-Temperature Steam Electrolysis system (HTSE), as well as the basic control logic that governs the primary components and the necessary hardware to achieve optimal functionality. The scaled facility will be a 1 MWth system that uses Dowtherm A as the simulant fluid for Flibe (the coolant of choice for the primary loop of molten salt reactors), and can validate the heat transfer and steady-state operational requirements of the 300 MWth prototype. The scaled facility matches the Prandtl and Reynolds numbers associated with steady-state operation of the FHR-HTSE’s primary loop without having to deal with very high temperatures, flow rates, or power inputs. This will allow the facility to run experiments that analyze various thermophysical and fluid-dynamic properties that characterize reactor operation, such as pressure drops, radial temperature distribution, heat exchanger conditions. The facility also has potential to integrate additional components of the prototype system, such as intermediate thermal-hydraulics loops, real-time grid-demand data, energy storage, and HTSE.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗