Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Cross Flow”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

High spatial resolution and contrast radiography of hydrodynamic instabilities at the National Ignition Facility

In this work, we are developing techniques for studying the Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities in a planar geometry at high-energy-densities at the National Ignition Facility (NIF). In particular, through the improvement of experimental imaging quality, we are progressing toward the study of the turbulent regime of the mixing regions in capsule implosion experiments for inertial confinement fusion, which requires few micrometers resolution. Using 60 NIF beams, a solid shock tube is driven launching a shock wave that crosses the interface between a dense and a light material pre-machined in the target to obtain sinusoidal ripples, which results in RM and RT instabilities that are imaged using the NIF Crystal Backlighter Imager. High-quality images were obtained with a mean resolution of 7 μm and improved contrast. While the obtained resolution does not allow the observation of the smallest scale of the “turbulent” energy spectrum, the generated image encompasses 63% of the total flow energy, a 50% improvement over previous studies, which is observed for the first time a roll-up feature in a high energy density-type RT experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Surface Properties of Double-Fin Generated Shock-Wave/Boundary-Layer Interactions

Shock-Wave/Boundary-Layer Interaction (SBLI) is ubiquitous in high-speed air vehicles. SBLI surface flowfield properties of crossing shocks induced by symmetric Double-Fins (DFs) at fin angles of 8° and 10° at Mach 2 are experimentally investigated; a 10° Single-Fin (SF) is also explored as a comparison. Surface visualization is used to capture distinct flowfield structures such as separation and upstream influence. Steady and unsteady Pressure-Sensitive Paints (PSPs) are used to obtain global surface pressure fields. Characteristics of DF-SBLI surface topology are discussed by comparing the corresponding SF-SBLI. Here, the effect of fin angle on DF-SBLI is characterized, and unsteady dynamics of the pressure field are examined. Spectral Proper Orthogonal Decomposition (SPOD) is performed on unsteady PSP, and traveling surface pressure waves and their frequency-dependent behavior are discussed. Dispersion relation of surface pressure waves is revealed by performing two-dimensional space–time Fourier transform on unsteady PSP, and the results are compared with the phase velocity distributions obtained from SPOD. Accompanying group velocity behavior is examined for distinct SBLI regions wherein wave propagation is found dependent on its frequency. Steady surface pressures along the centerline are compared with previous experimental data. Fin-tip spacing is found to be an essential parameter for scaling of upstream influence length.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2026

Numerical simulation of vortex-induced vibration response of a single IEA 10-MW wind turbine blade

Three-dimensional simulation of vortex-induced vibration (VIV) of a single International Energy Agency (IEA) 10-MW reference wind turbine blade with a length of 97.325 m is performed using the ExaWind stack, an open-source suite of codes. This study aims to illustrate the spanwise VIV response characteristics and cross-validate the results with an existing commercial framework. Five near-body meshes and three time steps are selected for the convergence study. To improve computational efficiency, several VIV triggering methods are also compared to shorten the VIV development period. The ExaWind-based VIV simulation strategy for a single IEA 10-MW blade is determined. First, the modal shape is validated against published results. Then, spanwise VIV responses of four blade configurations under a fixed and varied incoming flow velocity are analyzed. Results show that the VIV response is dominated by the first edgewise (second overall) mode. Little first-mode contributions appear near the second-mode node, producing a pi phase jump, and a higher harmonics response occurs near the blade root. Rotational degrees of freedom are minor compared with translational motion. The response versus reduced velocity is analyzed, showing a two-branch behavior similar to that of VIV for a bluff cylinder. Across all tested cases, the dominant frequency remains locked to the natural frequency of the second mode with no observed desynchronization. A mild deviation is observed for the case of 90-degree pitch and 310-degree azimuth rotation near a reduced velocity of 6, which will be examined with additional cases in future work. These findings indicate that severe VIV responses can arise under specific configurations and flow conditions, thereby increasing the potential for VIV fatigue damage and requiring greater attention during operation.

17 WIND ENERGY↗
DOE OSTI2023

Less Is More: Oligomer Extraction and Hydrothermal Annealing Increase PDMS Adhesion Forces for Materials Studies and for Biology-Focused Microfluidic Applications

Cues in the micro-environment are key determinants in the emergence of complex cellular morphologies and functions. Primary among these is the presence of neighboring cells that form networks. For high-resolution analysis, it is crucial to develop micro-environments that permit exquisite control of network formation. This is especially true in cell science, tissue engineering, and clinical biology. We introduce a new approach for assembling polydimethylsiloxane (PDMS)-based microfluidic environments that enhances cell network formation and analyses. We report that the combined processes of PDMS solvent-extraction and hydrothermal annealing create unique conditions that produce high-strength bonds between solvent-extracted PDMS (E-PDMS) and glass—properties not associated with conventional PDMS. Extraction followed by hydrothermal annealing removes unbound oligomers, promotes polymer cross-linking, facilitates covalent bond formation with glass, and retains the highest biocompatibility. Herein, our extraction protocol accelerates oligomer removal from 5 to 2 days. Resulting microfluidic platforms are uniquely suited for cell-network studies owing to high adhesion forces, effectively corralling cellular extensions and eliminating harmful oligomers. We demonstrate the simple, simultaneous actuation of multiple microfluidic domains for invoking ATP- and glutamate-induced Ca 2+ signaling in glial-cell networks. These E-PDMS modifications and flow manipulations further enable microfluidic technologies for cell-signaling and network studies as well as novel applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Generative AI models for learning flow maps of stochastic dynamical systems in bounded domains

Simulating stochastic differential equations (SDEs) in bounded domains, presents significant computational challenges due to particle exit phenomena, which requires accurate modeling of interior stochastic dynamics and boundary interactions. Despite the success of machine learning-based methods in learning SDEs, existing learning methods are not applicable to SDEs in bounded domains because they cannot accurately capture the particle exit dynamics. We present a unified hybrid data-driven approach that combines a conditional diffusion model with an exit prediction neural network to capture both interior stochastic dynamics and boundary exit phenomena. Our ML model consists of two major components: a neural network that learns exit probabilities using binary cross-entropy loss with rigorous convergence guarantees, and a training-free diffusion model that generates state transitions for non-exiting particles using closed-form score functions. The two components are integrated through a probabilistic sampling algorithm that determines particle exit at each time step and generates appropriate state transitions. Here, the performance of the proposed approach is demonstrated via three test cases: a one-dimensional simplified problem for theoretical verification, a two-dimensional advection-diffusion problem in a bounded domain, and a three-dimensional problem of interest to magnetically confined fusion plasmas.

Bounded domains↗
DOE OSTI2024

Simulation of gas mixture dynamics in a pipeline network using explicit staggered-grid discretization

Here we develop an explicit staggered finite difference discretization scheme for simulating the transport of highly heterogeneous gas mixtures through pipeline networks. This study is motivated by the proposed blending of hydrogen into natural gas pipelines to reduce end use carbon emissions while using existing pipeline systems throughout their planned lifetimes. Our computational method accommodates an arbitrary number of constituent gases with very different physical properties that may be injected into a network with significant spatiotemporal variation. In this setting, the gas flow physics are highly location- and time- dependent, so that local composition and nodal mixing must be accounted for. The resulting conservation laws are formulated in terms of pressure, partial densities and flows, and volumetric and mass fractions of the constituents. We include non-ideal equations of state that employ linear approximations of gas compressibility factors, so that the pressure dynamics propagate locally according to a variable wave speed that depends on mixture composition and density. We derive compatibility relationships for network edge boundary values that are more complex than for a homogeneous gas. The simulation method is evaluated on initial boundary value problems for a single pipe and a small network, is cross-validated with a lumped element simulation, and used to demonstrate a local monitoring and control policy for maintaining allowable concentration levels.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Hydrogen embrittlement of Zircaloy-4 fabricated by ultrasonic additive manufacturing

Ultrasonic additive manufacturing (UAM) was successfully applied to the zirconium material system to create a planar geometry. Following fabrication, SS-J3 type tensile specimens of the UAM and wrought Zircaloy-4 with a nominal gage section of 5×1.2×0.75 mm were machined for hydriding studies and mechanical testing. The SS-J3 type tensile specimens were gas-charged with hydrogen using a custom system that precisely controls hydrogen gas flow rate, hydrogen partial pressure, and temperature. To avoid altering the UAM material, the maximum process temperature was limited to 550°C. Using different initial hydrogen gas pressures and flow rates, various hydrogen contents (70–1755 wppm) were achieved for Zircaloy-4 specimens. Tensile testing shows that, regardless of hydrogen content, all UAM specimens measured yield strengths in the range of 557±16 MPa and ultimate tensile strengths of 660±4MPa. However, total elongation clearly decreased as a function of increasing hydrogen content. Further, at the lowest hydrogen content, the total plastic elongation measured 21.5%, and this value decreased to less than 2% when the hydrogen content was increased to 1000 wppm. Cross-sectional optical microscopy images revealed that the hydride distributions are randomly oriented. For the low hydrogen content specimen, these randomly distributed hydrides are isolated from each other. A sandwiched structure, consisting of high-density and low-density layers, was also observed for the specimens with hydrogen content between 200 and 400 wppm. As the hydrogen content increases, the hydrides diffuse into the low-density hydride layers to form a network across the whole specimen. Although the tensile properties of UAM and wrought Zircaloy-4 exhibit the same behavior with increasing hydrogen content, the distinctions in grain orientations led to differences in the hydride orientations at low and intermediate hydrogen concentrations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Giant anisotropic in-plane thermal conduction induced by Anomalous phonons in pentagonal PdSe 2

In two-dimensional materials, different atomic stacking induces anisotropic atomic interactions and phonon dispersions, leading to the anisotropy of in-plane thermal transport. In this work, we report an exceptional case in layered pentagonal PdSe 2 , where the bonding, force constants, and lattice constants are nearly-equal along the in-plane crystallographic axis directions, while the thermal conductivity is surprisingly much greater along b-axis than along a-axis with a ratio up to 1.8. Such strong anisotropy is not only unexpected in in-plane uniform structured materials, but also 2 comparable to the record high in-plane anisotropic thermal conductivity in the nonuniform structured material reported to date (the ratio is ~2.0 in TiS 3 ). By combining the inelastic X-ray scattering measurement and the first-principles calculations, we attribute such high anisotropy to the low-energy phonons along a-axis, particularly their lower group velocities and “avoided-crossing” behavior. The different buckling structures between a- (zigzag-type) and b-axis (flat-type) are mainly responsible for such unique phonon dynamics properties of PdSe 2 . This finding helps to discover materials with high anisotropic in-plane thermal conductivity in uniform structures and reveals new physics of anisotropy of in-plane thermal conduction. Due to the unique features in structure and thermal transport properties, PdSe 2 may serve as a new platform for designing novel devices to route heat flow precisely at the nanoscale.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Scale-Dependent Bedform Migration and Deformation in the Physical and Spectral Domains

Multi-scale bathymetries observed in laboratory channels and natural rivers have posed several challenges to the description of bedform geometry and kinematics. To better quantify scale-dependent migration velocity, high resolution spatio-temporal bed evolution data are analyzed using three independent methods: (i) a bedform tracking method that identifies individual bedforms in longitudinal bed elevation profiles, (ii) a cross-correlation based, image matching technique in consecutive bathymetry scans similar to Particle Image Velocimetry, and (iii) two-dimensional frequency-wavenumber spectra marking a dispersion relation between bedform period and length, allowing estimation of Fourier-based, scale-dependent migration velocities. The comparative studies show that the spectral and bedform tracking methods present small differences in the migration velocity of medium-large bedforms. However, significant deviations occur in the small bedforms, mostly due to bedform sheltering or amplified exposure to the flow, which also contribute to their enhanced scale-dependent deformation. The analysis also shows that some form of scale separation is needed to extract the large and slowly migrating bedforms as they are masked by secondary features when several bedform orders coexist. Bathymetric Image Velocimetry reasonably estimates the averaged migration velocity and exhibits potential to study relations between spanwise bed surface velocity and local bedform kinematics.

58 GEOSCIENCES↗
DOE OSTI2022

Hawaii Wave Surge Energy Converter (HAWSEC) OSU O.H. Hinsdale Basin

The following information and metadata applies to both the Phase I (Hydrodynamics) and Phase II (Full System Power Take-Off) zip folders which contain testing data from the OSU (Oregon State University) O.H. Hinsdale Wave Research Laboratory, from both OSU and the University of Hawaii at Manoa (UH). See zip folders provided further below in the downloads section. For experimental data of the full system, including PTO, see Phase II dataset. There are two main directories in each Phases's zip folder: "OSU_data" and "UH_data". The "OSU_data" directory contains data collected from their DAQ (data acquisition system), which includes all wave gauge observations, as well as body motions derived from their Qualisys motion tracking system. The organization of the directory follows OSU's convention. Detailed information on the instrument setup can be found under "OSU_data/docs/setup/instm_locations". The experiments conducted are documented in the "OSU_data/docs/daq_logs", which provides the trial number to the corresponding data located under "OSU_data/data" in several formats (e.g., ".mat" and ".txt"). Inside the trial directory, data is provided for each of the instruments defined in "OSU_data/docs/setup/instm_locations". The "UH_data" directory contains data collected from their DAQ. The data is stored in a ".tdms" file format. There are free plug-ins for Microsoft Excel and MathWorks MATLAB to read the ".tdms" format. Below are a few links providing methods to read in the data, but a Google search should identify alternatives sources if these no longer exist (valid as of January 2024): Excel: http://www.ni.com/example/27944/en/ MATLAB: https://www.mathworks.com/matlabcentral/fileexchange/30023-tdms-reader The Excel plugin is recommend to get a quick overview of the data. The UH data is organized by directory name, in which the sub-directories for each experiment contains a directory whose name defines the wave height and period for the experimental data within. For example, a directory name "H02_T0275" corresponds to an experiment with wave height 0.1m and a period of 2.75s. For random wave data, the gamma value is also included in the directory name. For example, a directory name "H02_T0225_G18" corresponds to an experiment with a significant wave height of 0.2m, a peak period of 2.25s, and a gamma value of 1.8, with each spectra being a TMA spectrum. For the free decay experiments, the directory name is defined by the initial angular displacement. For example, a directory name "ang05_run01" corresponds to an experiment with an initial angular displacement of 5 degrees. There is a dataset in the UH data for each corresponding experiment defined in the OSU DAQ logs. The ".tdms" data is output from the DAQ at fixed intervals. Therefore, if multiple files are contained within the folder, the data will need to be stitched together. Within the UH dataset, there are two input channels from the OSU DAQ providing a random square wave signal for time synchronization ("ENV-WHT-0010") and a high/low signal ("ENV-WHT-0012") to identify when the wave maker is active (+5V). The UH data is logged as a collection of channel outputs. Channels not in use for the OSU testing (either Phase I or Phase II) are marked "nan" below. If the sensor is disconnected, it will record noise throughout the experiment. Below are the channel definitions in terms of what they measure: GPS Time = time CYL-POS-0001 = position between flap and fixed reference CYL-LCA-0001 = force between flap and hydraulic cylinder REC-LPT-0001 = nan REC-HPT-0001 = nan REC-HPT-0002 = nan REC-HPT-0003 = nan HHT-HPT-0001 = pressure at exhaust ("head" only) REC-FQC-0001 = nan REC-FQC-0002 = nan HHT-FQC-0001 = flow at exhaust ("head" only) ENV-WHT-0001 = nan ENV-WHT-0002 = nan ENV-WHT-0003 = nan ENV-WHT-0010 = random signal from OSU DAQ ENV-WHT-0012 = high/low signal from OSU DAQ Also included is a calibration curve to convert the string pot data to flap pi...

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Precise relative magnitude measurement improves fracture characterization during hydraulic fracturing

SUMMARY Microseismic monitoring is an important technique to obtain detailed knowledge of in-situ fracture size and orientation during stimulation to maximize fluid flow throughout the rock volume and optimize production. Furthermore, considering that the frequency of earthquake magnitudes empirically follows a power law (i.e. Gutenberg–Richter), the accuracy of microseismic event magnitude distributions is potentially crucial for seismic risk management. In this study, we analyse microseismicity observed during four hydraulic fracture treatments of the legacy Cotton Valley experiment in 1997 at the Carthage gas field of East Texas, where fractures were activated at the base of the sand-shale Upper Cotton Valley formation. We perform waveform cross-correlation to detect similar event clusters, measure relative amplitude from aligned waveform pairs with a principal component analysis, then measure precise relative magnitudes. The new magnitudes significantly reduce the deviations between magnitude differences and relative amplitudes of event pairs. This subsequently reduces the magnitude differences between clusters located at different depths. Reduction in magnitude differences between clusters suggests that some attenuation-related biases could be effectively mitigated with relative magnitude measurements. The maximum likelihood method is applied to understand the magnitude frequency distributions and quantify the seismogenic index of the clusters. Statistical analyses with new magnitudes suggest that fractures that are more favourably oriented for shear failure have lower b-value and higher seismogenic index, suggesting higher potential for relatively larger earthquakes, rather than fractures subparallel to maximum horizontal principal stress orientation.

58 GEOSCIENCES↗
DOE OSTI2021

Two-Phase Turbulence Statistics from High Fidelity Dispersed Droplet Flow Simulations in a Pressurized Water Reactor (PWR) Sub-Channel with Mixing Vanes

In the dispersed flow film boiling regime (DFFB), which exists under post-LOCA (loss-of-coolant accident) conditions in pressurized water reactors (PWRs), there is a complex interplay between droplet dynamics and turbulence in the surrounding steam. Experiments have accredited particular significance to droplet collision with the spacer-grids and mixing vane structures and their consequent positive feedback to the heat transfer recorded in the immediate downstream vicinity. Enabled by high-performance computing (HPC) systems and a massively parallel finite element-based flow solver—PHASTA (Parallel Hierarchic Adaptive Stabilized Transient Analysis)—this work presents high fidelity interface capturing, two-phase, adiabatic simulations in a PWR sub-channel with spacer grids and mixing vanes. Selected flow conditions for the simulations are informed by the experimental data found in the literature, including the steam Reynolds number and collision Weber number (Wec={40,80}), and are characteristic of the DFFB regime. Data were collected from the simulations at an unprecedented resolution, which provides detailed insights into the continuous phase turbulence statistics, highlighting the effects of the presence of droplets and the comparative effect of different Weber numbers on turbulence in the surrounding steam. Further, axial evolution of droplet dynamics was analyzed through cross-sectionally averaged quantities, including droplet volume, surface area and Sauter mean diameter (SMD). The downstream SMD values agree well with the existing empirical correlations for the selected range of Wec. The high-resolution data repository from the simulations herein is expected to be of significance to guide model development for system-level thermal hydraulic codes.

Saini, Nadish↗
DOE OSTI2024

Characterization of 235 U, 238 U and 239 Pu fission ionization chamber foils by α and γ -ray spectrometry

The masses of thin, electroplated deposits of 235 U, 238 U and 239 Pu on titanium backings have been determined to better than 1.3% precision by employing three complementary measurement techniques: $\alpha$ spectrometry, 2$\pi$ gas counting, and $\gamma$-ray spectrometry. A dedicated $\alpha$ spectrometer was designed and calibrated to allow for repeatable, high-precision measurements of the absolute $\alpha$ activity of a sample. The isotopic composition of the sample was determined by fitting the peaks in the $\alpha$ spectrum, a technique which was demonstrated to be generally consistent with mass spectrometry. Independently, the total activity of the deposits were measured in a 2$\pi$ gas flow proportional counter. A detailed study of the effects of $\alpha$ particle straggling in the actinide deposit and $\alpha$ particle recoil off the sample backing material was performed. The masses of these actinide deposits were also measured via $\gamma$-ray spectrometry. Here, we found that all three systematically independent techniques produced consistent results, and could be combined into a weighted average with further reduced uncertainty. The mass ratios of these actinide deposits were determined to better than 1% precision via $\alpha$ spectrometry, enabling high-precision fission cross section ratio measurements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Unstructured mesh-based multi-physics reactor transient simulation with iMC code

In this study, we present a time-dependent Monte Carlo simulation result about a complex fuel-loaded light water reactor system with thermal-hydraulic feedback using the KAIST iMC code. The intra-pin temperature distribution is calculated based on Monte Carlo-tallied detailed power distribution and the finite element heat transfer method. The temperature effect of fuel, absorber, and coolant are adequately reflected into the dynamic response by the on-the-fly Doppler broadened cross-section correction. For the time-dependent Monte Carlo transport simulation, the predictor-corrector quasi-static Monte Carlo method is used for this study. We tested the coupled analysis framework with a 17-by-17 fuel assembly loaded with the centrally-shielded burnable absorber (CSBA) fuel element. The temperature feedback effect on the material cross-section and the moderator density incurred negative reactivity in response to the external positive reactivity insertion, showing the inherent stability of the reactor system. The pin-wise power and temperature change during the system transient are presented and discussed. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Plasma diagnostics and modeling of lithium-containing plasmas

Thin-film deposition from chemically reactive multi-component plasmas is complex, and the lack of electron collision cross-sections for even the most common metalorganic precursors and their fragments complicates their modeling based on fundamental plasma physics. This study focuses on understanding the plasma physics and chemistry in argon (Ar) plasmas containing lithium bis (trimethylsilyl) amide used to deposit Li x Si y thin films. These films are emerging as potential solid electrolytes for lithium-ion batteries, and the Li-to-Si ratio is a crucial parameter to enhance their ionic conductivity. We deposited Li x Si y films in an axial flow-through plasma reactor and studied the factors that determine the variation of the Li-to-Si ratio in films deposited at various points on a substrate spanning the entire reactor axis. While the Li-to-Si ratio is 1:2 in the precursor, the Li-to-Si ratio is as high as 3:1 in films deposited near the plasma entrance and decreases to 1:1 for films deposited downstream. Optical emission from the plasma is dominated by Li emission near the entrance, but Li emission disappears downstream, which we attribute to the complete consumption of the precursor. We hypothesized that the axially decreasing precursor concentration affects the electron energy distribution function in a way that causes different dissociation efficiencies for the production of Li and Si. We used Li line intensities to estimate the local precursor concentration and Ar line ratios to estimate the local reduced electric field to test this hypothesis. This analysis suggests that the mean electron energy increases along the reactor axis with decreasing precursor concentration. As a result, the decreasing Li-to-Si ratio with axially decreasing precursor concentration may be explained by Li release from the precursor having lower threshold energy than Si release.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Communication Network Layer State Estimation Measurement Model for a Cyber-Secure Smart Grid

Network communication has been proven to be a very important tool and a key factor in the recent development and progress of the power grid operation. It is also considered as the foundation for the smart grid because information and communication are integrated into electricity distribution to achieve reliable and accurate knowledge of the power grid. In previous years, absorbing energy from substations and delivering it to customers was the only type of interaction we knew between utility companies and customers. Presently, the growing connections of small distributed generation units caused by the cost reduction of most of the technologies used in generation and storage of electrical energy, along with the potential benefits of renewable energy have pushed many researchers to look into the improvement of information and communication technologies (ICT) in order to ensure a bidirectional flow of power and data. Moreover, the evolution of information and communication technologies and its applications to smart grid have converted the smart grid into a cyber-physical system where vulnerabilities and additional security challenges such as cyber-threats and cyber-attacks have emerged. Previously, we have demonstrated that using machine learning-based processing on data gathered from communication networks and the power grid was a promising solution for detecting cyber threats by implementing a co-simulation of cyber-security for cross-layer strategy. Since the majority of the challenges observed can only be solved in the network communication layer, we present in this work a physics-based state estimation model of the communication network system towards enhanced cyber-physical security of the smart grid. Information integration with the previously developed machine learning model is developed, providing a enhanced cyber-physical security application for the smart grid. Easy-to-implement model, without hard-to-derive parameters, highlight potential aspects of the model for real-life applications.

Mathieu, Reynold↗
DOE OSTI2022

Impact of nuclear data covariance libraries on uncertainty quantification of sodium cooled fast reactor simulation

Uncertainty quantification in core modelling to assess reactor performance is a regular practice to identify design safety limits. Best-estimate plus uncertainty (BEPU) methods are commonly used to propagate uncertainties to core responses from input parameters, characterized by variance-covariance libraries. This study evaluates the impact of different covariance libraries on the uncertainty associated with core safety parameters and their propagation through reactor transient modelling for a metallic core (ABR-1000). The COMAC-1.0, COMMARA-2.0, and COMMARA-2.1 covariance libraries are used to propagate nuclear data uncertainties to neutron feedback coefficients using the Generalized Perturbation Theory, then through transient safety analyses to establish the confidence interval for safety performance of the reactor using stochastic sampling. The considered transients include Unprotected Transient Over Power and Unprotected Loss Of Flow. Most influential nuclide and reaction pairs in different variance-covariance libraries are identified and their impact on the uncertainty of safety parameters are evaluated. Significant differences were observed between uncertainties from COMAC compared to COMMARA libraries. Uncertainties from COMMARA libraries were consistently larger than COMAC for all cases. For structure, coolant density feedback coefficients and control rod worth, total nuclear data uncertainties were approximately 50% larger for COMMARA than COMAC. Differences were also observed in identification of top contributors of uncertainty in terms of reaction-cross sections. For instance, largest contribution to total uncertainty of radial expansion coefficient from COMMARA-2.0 and 2.1 originates from Na{sup 23} P1 Elastic cross section (0.81%) but from Pu{sup 239} Fission (0.53%) with COM C-1.0. Such differences are evaluated in this research for all quantities of interest. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗