Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reference systems”

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

Principles of Corrosion in Nuclear Systems: Theory and Analytical Methods

This article is meant to provide the reader with the fundamental electrochemistry necessary to critically evaluate the data and discussions in the articles on corrosion that follow. As such it has two main sections, theory and analytical methods. In the section on theory, a review of half cell reactions, cell potentials and the Nernst equation, and reference electrodes in nuclear systems are discussed. In addition, a review of the thermodynamics from room temperature to the PWR and kinetics of dissolution and passive film formation is also presented. In the section on analytical methods, we review the principles of potentiodynamic polarization, electrochemical impedance spectroscopy, and we introduce a more seldom used technique, Mott-Schottky analysis.

Lillard, R Scott↗
DOE OSTI2016

SU-F-T-640: Feasibility of Using a Commercially Available Surface Guided Radiotherapy System with An Open-Face SRS Immobilization System

Purpose: To evaluate the inherent accuracy of using a surface guided radiotherapy system (SGRT) in the setup and monitoring of patients receiving stereotactic radiosurgery with an open-face SRS immobilization system. Methods: An anthropomorphic head phantom was set up using the Qfix Encompass SRS Immobilization System on a Varian Edge with OSMS and Varian TrueBeam with AlignRT. The phantom was positioned at 0° gantry and couch. A reference image was acquired using the SGRT system and an ROI was created over the mask opening. The couch and gantry were rotated to different combinations focusing on clinically used SRS gantry/couch combinations and those blocking the SGRT cameras. Perceived surface deviation by the SGRT system from the reference image was recorded. A Winston-Lutz test was performed on couch angles tested and used to exclude couch walkout. The deviation magnitude was calculated using translational values and rotational raw values were recorded. Results: The maximum couch walkouts were: 0.4mm (Edge) and 0.5mm (TB). Solely rotating the gantry resulted in a median couch deviation of 0.2mm and range of 0.1–0.3mm for both linacs. Only rotating the couch (0° gantry) resulted in median deviations of 0.6mm and 0.5mm with ranges of 0.3–1.0mm and 0.3–0.7mm for the Edge and TB, respectively. Combining gantry and couch rotations, the median deviations were 0.7mm and 0.9mm with ranges of 0.3–1.1mm and 0.2–1.9mm for the Edge and TB, respectively. Including all combinations, rotation, roll, and pitch median deviations ranged from 0.1–0.3° with pitch demonstrating consistently higher values and a maximum deviation of 1.0° (both linacs). Conclusion: SGRT is a reliable monitoring tool, though taking into account system fluctuations, 1mm is too restrictive a site tolerance to use with the Qfix Encompass mask. Gantry rotation has little effect on system fluctuation even with camera blockage, whereas couch rotation has a larger effect.

43 PARTICLE ACCELERATORS↗
DOE OSTI2013

Time-dependent resonant UHF CI approach for the photo-induced dynamics of the multi-electron system confined in 2D QD

We extend the static multi-reference description (resonant UHF) to the dynamic system in order to include the correlation effect over time, and simplify the TD Schrödinger equation (TD-CI) into a time-developed rate equation where the TD external field Ĥ′(t) is then incorporated directly in the Hamiltonian without any approximations. We apply this TD-CI method to the two-electron ground state of a 2D quantum dot (QD) under photon injection and study the resulting two-electron Rabi oscillation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Grid-Supportive Loads - A New Approach to Increasing Renewable Energy in Power Systems

This paper demonstrates the potential of inverter-based loads to support grid reliability during power system transients thereby enabling reliable integration of renewable energy in power systems. Such loads are referred to in this paper as grid-supportive loads (GSLs). A new GSL model is developed that simulates the transient response capabilities that can be programmed in electronic loads. The model’s design enables it to be easily integrated in widely used commercial power system transient analysis software. Theoretical expressions are derived that explain the workings of the GSL model. The performance, numerical stability, and impact of the GSL model is validated on 9-bus and 2000-bus synthetic power system models using generator tripping and bus fault disturbances. Results on the 2000 bus system show that in the absence of frequency support from wind/solar generation resources, just 20% of loads with grid-supportive capabilities can improve frequency response by up to 2000 MW/0.1 Hz and reduce deviation in frequency at nadir by up to 60% compared to the situation when GSLs are absent. Power system reliability also improves under fault events. Here, it is further shown that GSLs can aid in integrating more renewable generation without degrading the overall transient response of the power system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2016

Tutorial: Crystal orientations and EBSD — Or which way is up?

Electron backscatter diffraction (EBSD) is an automated technique that can measure the orientation of crystals in a sample very rapidly. There are many sophisticated software packages that present measured data. Unfortunately, due to crystal symmetry and differences in the set-up of microscope and EBSD software, there may be accuracy issues when linking the crystal orientation to a particular microstructural feature. In this paper we outline a series of conventions used to describe crystal orientations and coordinate systems. These conventions have been used to successfully demonstrate that a consistent frame of reference is used in the sample, unit cell, pole figure and diffraction pattern frames of reference. We establish a coordinate system rooted in measurement of the diffraction pattern and subsequently link this to all other coordinate systems. A fundamental outcome of this analysis is to note that the beamshift coordinate system needs to be precisely defined for consistent 3D microstructure analysis. This is supported through a series of case studies examining particular features of the microscope settings and/or unambiguous crystallographic features. These case studies can be generated easily in most laboratories and represent an opportunity to demonstrate confidence in use of recorded orientation data. Finally, we include a simple software tool, written in both MATLAB® and Python, which the reader can use to compare consistency with their own microscope set-up and which may act as a springboard for further offline analysis. - Highlights: • Presentation of conventions used to describe crystal orientations • Three case studies that outline how conventions are consistent • Demonstrates a pathway for calibration and validation of EBSD based orientation measurements • EBSD computer code supplied for validation by the reader.

36 MATERIALS SCIENCE↗
DOE OSTI2019

New Strategies for Time Delay Estimation during System Calibration for UAV-Based GNSS/INS-Assisted Imaging Systems

The need for accurate 3D spatial information is growing rapidly in many of today’s key industries, such as precision agriculture, emergency management, infrastructure monitoring, and defense. Unmanned aerial vehicles (UAVs) equipped with global navigation satellite systems/inertial navigation systems (GNSS/INS) and consumer-grade digital imaging sensors are capable of providing accurate 3D spatial information at a relatively low cost. However, with the use of consumer-grade sensors, system calibration is critical for accurate 3D reconstruction. In this study, ‘consumer-grade’ refers to cameras that require system calibration by the user instead of by the manufacturer or other high-end laboratory settings, as well as relatively low-cost GNSS/INS units. In addition to classical spatial system calibration, many consumer-grade sensors also need temporal calibration for accurate 3D reconstruction. This study examines the accuracy impact of time delay in the synchronization between the GNSS/INS unit and cameras on-board UAV-based mapping systems. After reviewing existing strategies, this study presents two approaches (direct and indirect) to correct for time delay between GNSS/INS recorded event markers and actual time of image exposure. Our results show that both approaches are capable of handling and correcting this time delay, with the direct approach being more rigorous. When a time delay exists and the direct or indirect approach is applied, horizontal accuracy of 1–3 times the ground sampling distance (GSD) can be achieved without either the use of any ground control points (GCPs) or adjusting the original GNSS/INS trajectory information.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Spectral Effects in Albedo and Rearside Irradiance Measurement for Bifacial Performance Estimation

We investigate the impact of spectral dependence of ground surface reflectivity on albedo and rearside irradiance measurements necessary for bifacial photovoltaic (PV) module performance estimation and monitoring. Because PV modules are spectrally selective, albedo and irradiance measurements performed with common irradiance sensors may require spectral mismatch corrections when used for performance prediction. We investigate via simulation the differences in spectrally responsive albedo measured with thermopile pyranometers and crystalline silicon PV reference cells in comparison to a typical crystalline-silicon bifacial PV module. Simulations are performed for nine different representative ground surface materials using simulated solar spectra together with spectral reflectivity data distributed with the SMARTS simulation software. For the materials considered, the results show that albedo spectral mismatch relative to the bifacial module is distributed over a range of ±9.2% for thermopile pyranometers versus only ±3.7% for a typical PV reference cell. We consider the impact of this spectrally-responsive albedo mismatch on bifacial PV module rearside irradiance measurements. Using synthesized rearside spectral irradiance distributions, we find that for the nine different ground surface materials the predicted rearside irradiance measurement deviates from the effective irradiance observed by the PV module by on the order of 16.5 W/m 2 for the pyranometer and 3.6 W/m 2 for the PV reference cell. We discuss the implications for bifacial albedo and irradiance measurement.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2025

Co-Located Wave Energy Converter (WEC) and Aquaculture System Annotated Bibliography

This annotated bibliography includes references that could aid in the design of a co-located WEC and aquaculture system off the coast of Guam. The breadth of this work covers multiple co-location archetypes such as: 1. WEC seawater desalination system a) Nearshore and deepwater WEC deployment b) Onshore and offshore aquaculture 2. WEC powering an offshore aquaculture platform a) Nearshore or deepwater WEC deployment 3. WEC powering an onshore aquaculture system a) Nearshore WEC deployment 4. Wave powered seawater pump a) Nearshore WEC deployment b) Onshore aquaculture There are two archetypes that may be of immediate interest to the community in Guam are to service the existing Fadian Hatchery (Mangilao) and to support freshwater aquaculture activities. First, the seawater pump at the hatchery that fills the facility’s seawater storage unit is broken. A nearshore seawater pumping WEC could be a solution to this issue. In addition, due to the high frequency of typhoons/extreme conditions and the island’s bathymetry, the likelihood of community support for an offshore aquaculture platform or WEC deployed in deepwater is low. Proactive and resilient solutions not just for power, but for freshwater are of interest as well to support any freshwater aquaculture activities. Therefore, a nearshore WEC desalination system is another archetype to consider.

16 TIDAL AND WAVE POWER↗
DOE OSTI2017

Mapping of orthogonal 2D flux coordinates for two nearby magnetic X-points to logically rectangular domains

The behavior of plasmas in configurations with magnetic null points is of potential interest for a number of systems. Such nulls, sometimes referred to as X-points, a name strictly valid only for wellseparated nulls where magnetic separatrices cross at a right angle, may occur in only one component of the magnetic field, such as the poloidal field in a tokamak. Numerical simulations are often used to study and predict turbulence and transport in such systems. For cases with strong magnetic field such that the ion and electron gyro-frequencies are much larger that particle collision frequencies, there is a strong anisotropy in the plasma transport along the magnetic field compared to across the field. Because of the strong anisotropy, it is often advantageous for simulations to use a coordinate system based on magnetic flux surfaces together with a coordinate orthogonal to the flux surfaces. This report presents the topological considerations needed to construct the mapping on an orthogonal flux-surface mesh onto a logically rectangular domain often used in simulation codes of the divertor region of tokamaks having two nearby X-points, thus extending the mapping now common single-null or well-separated double nullpoint cases with one null located at the top and bottom of the device. This extension is especially relevant to study the capability of snowflake divertors in tokamaks for reducing the peak heat flux of exhaust plasma on the target plates [1-3]. The analysis presented here follows the analytic description of the poloidal magnetic field, B pol , near two closely-spaced X-points given by Ryutov [1-3]. This analysis is applicable to experimental devices where the two X-points are generated by distributed poloidal field coils and the plasma core current. This analytic description is used to identify topological features of an orthogonal mesh. A key parameter that specifies the various topology changes is the poloidal angle between the primary X-point and the second X-point, and the changes are associated with a set of critical angles. Understanding these changes leads to the final step in the mesh construction by identifying the locations of various regions or patches in configuration space and in a corresponding logically rectangular index space. UEDGE [4] and other 2-D fluid plasma simulation codes [5-7] describe the edge plasma of toroidally symmetric tokamak configurations with one or more X-points in a divertor region. The 2-D spatial mesh is usually based on a magnetic configuration produced by an MHD equilibrium code such as EFIT and written out as one or more “eqdsk” files. As illustrated in Fig. 1, the domain of the simulated plasma consists of several flux tubes bounded by poloidal magnetic flux surfaces. The number and spacing of the flux surfaces defines the “radial” resolution of the mesh. The “radial” direction is locally orthogonal to flux surfaces; the “poloidal” direction is locally parallel to flux surfaces. The “poloidal” resolution of the mesh is defined by the number of cells in the poloidal direction between adjacent flux surfaces. For boundary plasma transport codes used for tokamak configurations such as UEDGE it is convenient, though not necessary, to assume that all flux tubes have the same number of poloidal cells so the mesh is logically rectangular in index space, but does possess internal boundaries linking different sub-rectangles in the computational domain. The following sections of this report describe indexing schemes that map tokamak magnetic configurations with one or more x-points in the divertor region(s) onto a rectangular index space. Simulation results from the UEDGE code illustrate the schemes for a few configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

SU-F-J-21: Clinical Evaluation of Surface Scanning Systems in Different Treatment Locations

Purpose: To reduce imaging dose in fractionated IGRT, the ability of optical surface imaging systems (OSIS) to detect setup errors was tested. Therefore, clinical studies to evaluate for different treatment locations setup corrections derived by OSIS in comparison with x-ray image guidance in fractionated radiation therapy was performed. Methods: The setup correction accuracy of an OSIS system (AlignRT, VisionRT, London, UK) will be analysed for the 4 tumour locations Pelvis, Upper Abdomen, Thorax and Breast, 20 patients for each location in comparison to a different system (Sentinel, C-RAD, SE). For each patient, the setup corrections of the cone-beam computed tomography (CBCT) of an Elekta Versa HD linear accelerator (Elekta, Crawley, UK) is considered as gold-standard and then compared with those of the OSIS for the first ten fractions retrospectively. There were no clinical decisions made based on the surrogate system. For the OSIS, the reference surface is highly important as it represents the actual ground truth. It can be obtained either with the system itself or the surface structure delineated in the planning CT can be imported via DICOM interface. In this paper, the first results for the treatment region thorax are presented. The reference image modalities were compared. Results: Table 1 displays the difference between the setup corrections obtained with OSIS and CBCT in lateral (LAT), longitudinal (LNG) and vertical (VRT) direction for the DICOM reference image. While the median deviations are within a few millimeters, some outliers showed large deviations. Generally, the mean deviation as well as the spread was smallest in lateral and largest in vertical direction. Conclusion: Although the system allows fast, simple and non-invasive determination of setup corrections, it should be evaluated treatment region dependant. Therefore, the study is ongoing. The application of OSIS may help to reduce the imaging dose for the patient. We gratefully acknowledge the support by our colleagues from the Radiological University Clinic Heidelberg, where the study was performed. This work was funded by the Federal Ministry of Education and Research (BMBF) Germany, grant number 01IB13001B.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

A Local Macroscopic Conservative (LoMaC) Low Rank Tensor Method for the Vlasov Dynamics

Abstract In this paper, we propose a novel Local Macroscopic Conservative (LoMaC) low rank tensor method for simulating the Vlasov-Poisson (VP) system. The LoMaC property refers to the exact local conservation of macroscopic mass, momentum and energy at the discrete level. This is a follow-up work of our previous development of a conservative low rank tensor approach for Vlasov dynamics ( arXiv:2201.10397 ). In that work, we applied a low rank tensor method with a conservative singular value decomposition to the high dimensional VP system to mitigate the curse of dimensionality, while maintaining the local conservation of mass and momentum. However, energy conservation is not guaranteed, which is a critical property to avoid unphysical plasma self-heating or cooling. The new ingredient in the LoMaC low rank tensor algorithm is that we simultaneously evolve the macroscopic conservation laws of mass, momentum and energy using a flux-difference form with kinetic flux vector splitting; then the LoMaC property is realized by projecting the low rank kinetic solution onto a subspace that shares the same macroscopic observables by a conservative orthogonal projection. The algorithm is extended to the high dimensional problems by hierarchical Tuck decomposition of solution tensors and a corresponding conservative projection algorithm. Extensive numerical tests on the VP system are showcased for the algorithm’s efficacy.

Guo, Wei↗
DOE OSTI2019

Moore vs. Murphy: Tradeoffs between complexity and reliability in distributed energy system scheduling using software-as-a-service

Software-based optimization of building control strategies, including scheduling, has the potential to improve the performance of existing complex heating, ventilation, and air conditioning (HVAC), storage, and other systems—especially if temporally variable energy production, such as solar thermal or photovoltaics, is included. If reductions in energy bills can be achieved using optimized control strategies that take advantage of cost-saving opportunities, such as time-of-use pricing, the additional bill savings can cover further efficiency investment costs. As computer processing becomes cheaper over time (Moore's Law), opportunities to perform complex control optimization become more abundant, and these can be performed remotely as software-as-a-service (SaaS). However, by “perfecting” our control strategies, we run an increased risk that when something unexpected happens (Murphy's Law), the consequences of failure are greater. This study used simulation to explore the potential benefits of HVAC schedule optimization, delivery, and implementation using a SaaS paradigm, at various levels of complexity. Implementing optimal schedules in a model of an efficient building's HVAC system, the study predicts energy cost savings of up to 10% compared to the naïve reference control strategy. Optimizing more system control variables increases the potential energy cost savings; however, these savings could be compromised by failures in communication inherent in delivering schedules via SaaS. The additional cost of energy resulting from the risk of increased demand charges generally increased with increased communication failure to a much larger extent than the risk of increased energy use charges. This work suggests that moderate improvements in performance, achieved at low cost by simple means, may be more effective than highly optimized schemes, which are more susceptible to failure due to their dependence on complex interactions between systems.

42 ENGINEERING↗
DOE OSTI2024

Corrigendum: Simulations of stand-off runaway electron beam termination by tungsten particulates for tokamak disruption mitigation (2024 Nucl. Fusion 64 056019)

In the original article, we inadvertently omitted a reference which previously introduced the concept of tungsten injection for 'Disruption Mitigation in Tokamak Reactor via Reducing the Seed Electrons of Avalanche.' In that reference, the authors have proposed injection of a tungsten cylinder of 80 mm by 8 mm dimensions, by rail guns, to absorb the runaway seed population. This can be contrasted with another interesting idea of using small tungsten pellets coated with a low Z material for depleting the runaway seeds as an option for upgrade of the ITER Disruption Mitigation System. Interested readers are referred to those two papers on the tungsten injection for runaway seed removal/reduction, and the related issue of not shortening the current quench excessively. The application focus of was not on reducing the runaway seeds in the initial Ohmic-to-runaway current conversion phase, but on safe termination of a fully formed and likely decaying runaway beam that is about to scrape off against the first wall, which can be accelerated by vertical displacement events. This scheme thus serves as a last line of defense against potential wall damage. The choice here is locally released tungsten particulates, similar to the previous idea of using tungsten particulates as a dust shield for the divertor. The most significant finding of is that strong pitch angle scattering, in addition to energy attenuation and absorption, of high-energy runaways by tungsten particulates, provides a transport mechanism by which runaways damage of the first wall can be mitigated by reducing and spreading the runaway wall load. We regret the omission in the original article, and thank the authors of the cited article for bringing this matter to our attention. We also find the possibility of using the injected tungsten rod from to terminate the runaways, as opposed to the cloud of tungsten particulates in intriguing. In that context, it is of interest to mention another possibility of a retractable tungsten metal arm that can be swung out for deployment in runaway termination, which would remove the need for post-mitigation recovery of tungsten debris in injection schemes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

SU-G-JeP1-01: A Combination of Real Time Electromagnetic Localization and Tracking with Cone Beam Computed Tomography in Stereotactic Radiosurgery for Brain Tumors

Purpose: 1. online verification of patient position during treatment using calypso electromagnetic localization and tracking system. 2. Verification and comparison of positional accuracy between cone beam computed tomography and calypso system. 3. Presenting the advantage of continuation localization in Stereotactic radiosurgery treatments. Methods: Ten brain tumor cases were taken for this study. Patients with head mask were under gone Computed Tomography (CT). Before scanning, mask was cut on the fore head area to keep surface beacons on the skin. Slice thickness of 0.65 mm were taken for this study. x, y, z coordinates of these beacons in TPS were entered into tracking station. Varian True Beam accelerator, equipped with On Board Imager was used to take Cone beam Computed Tomography (CBCT) to localize the patient. Simultaneously Surface beacons were used to localize and track the patient throughout the treatment. The localization values were compared in both systems. For localization CBCT considered as reference. Tracking was done throughout the treatment using Calypso tracking system using electromagnetic array. This array was in tracking position during imaging and treatment. Flattening Filter free beams of 6MV photons along with Volumetric Modulated Arc Therapy was used for the treatment. The patient movement was observed throughout the treatment ranging from 2 min to 4 min. Results: The average variation observed between calypso system and CBCT localization was less than 0.5 mm. These variations were due to manual errors while keeping beacon on the patient. Less than 0.05 cm intra-fraction motion was observed throughout the treatment with the help of continuous tracking. Conclusion: Calypso target localization system is one of the finest tools to perform radiosurgery in combination with CBCT. This non radiographic method of tracking is a real beneficial method to treat patients confidently while observing real-time motion information of the patient.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Time Dependent Analysis of Laser Heating of Metal Parts in Contact with HE

We consider the time dependent heat transfer of laser energy impinging on metal parts that are in contact with high explosives (HE). The predecessor to this calculation was TN-J4-17-2109, “Heat Transfer Analysis from Laser Energy on Metal Parts in Contact with HE”, which considered the time independent heat transfer and compared to experimental results. The intent here is to extend those results to the time dependent case. The steady-state results were developed because they were the “worst case” from a safety perspective: the steady state represents the highest possible temperature excursion. However, we found that the temperatures achieved were well above the safety limits established by the DOE High Explosives Safety Manual. Moreover, they were well above the safety limit published by the LANL / Weapons Facilities Operations (WFO) in WFO-OS-ES-48, which establishes an upper temperature limit of 70°C for metal in contact with High Explosive (HE). The J-4 response to the original results was to design, build, and deploy an interlock system for DARHT that would ensure no laser energy on the firing point when both personnel and HE were present. This system was tested and proved effective, refer TP-DARHT-1150, (12/19/2016). It was reviewed and approved by the LANL Explosives Instrumentation Review Committee, refer PRES-J4-16-2140, “DARHT Laser Diagnostics Interlock System: Presentation to the EIRC”, (2016). Although the interlock system prevented personnel from being on the firing point under conditions where the diagnostic lasers would heat metal above the WFO limit, the WFO limit was still exceeded during the Multiplexed Photonic Doppler Velocimetry (MPDV) tuning process. The only practical solution is to limit the amount of time that the MPDV lasers are engaged at full power. The diagnostics interlock system was implemented with laser shutters that can fully interrupt the beams, and these shutters can be used to limit the amount of time that the lasers will impinge on the hydrotest assembly. It is the purpose of this paper to provide quantitative guidance on what length of time can be permitted during this tuning process. Therefore, we will solve the time-dependent heat transfer equation to determine the temperature increase as a function of time. This function can be inverted to determine the maximum allowable time for energizing the lasers based on the maximum allowable temperature at the metal to HE interface.

42 ENGINEERING↗
DOE OSTI2017

Systems Security Architecture Research and Development

The general technical goals of the project were to better understand the assumptions of security guarantees using secure commercial off-the-shelf (COTS) technologies, to identify possible security weaknesses introduced during system integration, and to create improved security hardware and software which result in well understood reference architectures for critical computing systems and applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Implementation of a pressure correction term in the CTF sub-channel code

CTF [1] is a modernized and improved version of the legacy subchannel code, COBRA-TF [2], that is being jointly developed and maintained by Oak Ridge National Laboratory (ORNL) and North Carolina State University. The code was adopted by ORNL for use in the Consortium for Advanced Simulation of Light Water Reactors (CASL) in 2012 for aiding in addressing CASL challenge problems. Since that time, activities related to CTF have included implementing software quality assurance measures, implementing new features and models, performing validation and verification testing, establishing and supporting a CTF User Group, and developing the code for use in coupled applications [3, 4, 5]. CTF is applicable to single- and two-phase flows in light water reactor geometries at normal and accident operating conditions. The two-phase flow model consists of the six-equation two-phase flow model [6] augmented by equations for the droplet field and also includes appropriate source terms. This work focuses on the numerical stability of the six-equation two-phase flow model and is part of a recent eort to improve CTF capabilities for simulating flows in boiling nuclear reactors (BWRs). As it is of common knowledge among the thermodynamic community, the six-equation two-phase flow model is not strongly hyperbolic, meaning some of its eigenvalues are either not real or not distinct [7]. In a Cauchy sense and for an initial-value problem, strong hyperbolicity was found to be equivalent to well-posedness which ensures uniqueness of the solution in the space-time used. This concept is also of importance when designing high-order numerical methods and also when implementing boundary conditions that rely on real eigenvalues such as Riemann solvers. A tremendous amount of research is available in the literature on how to hyperbolize the six-equation two-phase flow model by adding ad hoc terms to the system of equations often referred to as virtual mass or interfacial pressure correction term. From a mathematical prospective, the main idea driving this approach is to add terms to the phasic momentum equation containing first-order partial derivatives of the unknowns to modify the Jacobian matrix and thus canceling the imaginary part of the eigenvalues.For this work it was chosen to implement a interfacial pressure correction (IPC) term on the same model as what was done for the CATHARE system code [8], that allows to recover a two-pressure formulation. In the following, the six-equation two-phase flow model is recalled along with the IPC term and its influence on the eigenvalues. Then, numerical results and convergence studies are presented for a water-faucet problem [9] at a pressure of 10 bar. Finally, conclusions and future work are addressed.

Delchini, Marco↗
DOE OSTI2023

Distributed Energy Resource Visual Emulator: Phase 1

To help federal energy managers assess, monitor, and manage cybersecurity while achieving decarbonization, the National Renewable Energy Laboratory's Distributed Energy Resource Cybersecurity Framework (DER-CF) offers a comprehensive, web-based assessment tool focusing on cyber governance or policies, technical management, and physical security. The DER-CF currently presents users with a series of pertinent cybersecurity questions, which are used to generate a site-specific report and recommendations. This paper outlines a technical approach to integrate the DER-CF with another key asset—NREL's Advanced Research on Integrated Energy Systems (ARIES) Cyber Range—to visualize cybersecurity resilience and compliance and to enhance the usability and accessibility of the DER-CF. The result is a new tool called the Distributed Energy Resource Visual Emulator (DER-VE). Its development will include regular conversations with stakeholders to assess the effectiveness of these efforts, refine the visualization capability, and ensure its value to our partners. Phase 0 of the integration project was concluded in 2021. Phase 1, completed in 2022, has two components: The first is developing a working visualization of system compliance using the DER-CF, and the second is planning the design of a server application that takes input data from the DER-CF and creates a personal emulated environment of the user's system or a selected reference architect. Major components that were addressed in this phase are the DER-CF output, compliance visualization, data model, and compliance server design.

24 POWER TRANSMISSION AND DISTRIBUTION↗