Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Average model”

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

Port of New York and New Jersey Drayage Electrification Analysis

The National Renewable Energy Laboratory (NREL) evaluated the potential for drayage electrification in the Port of New York and New Jersey (PoNYNJ), with a focus on operators: Harbor Freight Transport (HF), Safeway Trucking (SWT), and International Motor Freight Inc (IMF). This report summarizes the data collection and electrification evaluation of all three drayage operators, includes detailed operational data, and identifies the performance requirements for battery electric tractors (BETs) and corresponding infrastructure operated within the context of PoNYNJ drayage operation. This report also details a methodology to evaluate opportunities, strategies, and challenges associated with future expansions of BETs in meeting PANYNJ emissions goals. The Port Authority has established a goal of achieving Net Zero carbon emissions by 2050 across all facilities, including from tenant and stakeholder sources such as drayage trucks. NREL used real-world performance data collected on the three PoNYNJ drayage operations, along with modeling and analysis tools to compare BET to diesel trucks. From March to July 2021, NREL collected 1Hz vehicle and engine data from 46 drayage trucks at the three operators totaling nearly 121,000 miles of operation, providing enough information to assess vehicle operations for electrification potential. A Future Automotive Systems Technology Simulator (FASTSim) electric truck powertrain model was validated using PoNYNJ data and scenarios were run to evaluate drayage truck electrification requirements over the real-world cycles. The first scenario examined BET viability with minimal changes to existing operations. This assumes the trucks charge when stopped for two hours or longer, have a functional battery size of 375 kWh, and can charge at 270 kilowatts (kW) average which are the specification of the commercially available Freightliner eCascadia. The second scenario looked at what operational, charging infrastructure, and BET technology changes would be needed to fully electrify. Finally, detailed analysis was run on charging rate structure to understand operational costs to the fleets. The studied drayage trucks averaged 5.1 MPG, spent roughly 9% of their energy at idle, and drove an average of 140 miles per day with a maximum daily distance of 573 miles. The FASTSim model results indicate a comparable BET would use 417 kWh of energy per day on average accounting for cargo weight, which is close to the full usable capacity of the eCascadia currently available on the market. Based on the daily average operating data, partial fleet electrification is possible with current technology. However, some specific days of operation would require over 1,600 kWh of energy due to longer distances traveled by the trucks and more intense operation. Trucks used for long distance and intense operation cannot be readily electrified with current technology without operational changes. Full adoption of BETs could reduce CO 2 emissions from these fleets by roughly 75% today, eliminating 76 metric tons of CO 2 (MTCO 2 ) per vehicle each year, which equates to 24,100 MTCO 2 per year for all three operators. Commercially available direct current fast chargers (DCFC) have charge rates up to 350 kW. Based on the average daily modeled energy use for each operator, current industrial rate structures, and the assumption of 350 kW peak charging, full drayage electrification would increase electricity consumption. In addition, peak demand usage would increase with unmanaged charging along with cost of electricity having a direct impact on cost per mile for electric vehicles. The resulting cost per mile for BETs along with comparable cost per mile for conventional diesel trucks are also examined at $\$$4.00 per gallon of diesel. It will be important for PANYNJ and the drayage operators within the PoNYNJ to consider these load impacts to their existing electrical infrastructure and devise operational strategies that avoid coincident charging of vehicles to mitigate demand charges. Despite these electricity cost increases, savings from reductions in diesel consumption will help offset the costs of this increased electricity consumption. However, prices of both electricity and diesel are subject to change based on various factors meaning the realized savings will vary over time. This shows BETs could be cost-competitive on an energy cost per mile basis for all scenarios while diesel is above $\$$3.00/gal. Further, if diesel prices dropped to the 15-year low of $2.33/gal, it would still be cost competitive to operate the EVs with electricity costs of 16.3 ¢/kWh or less.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

PANDEMIC: Occupancy driven predictive ventilation control to minimize energy consumption and infection risk

During the SARS-CoV-2 (COVID-19) pandemic, governments around the world have formulated policies requiring ventilation systems to operate at a higher outdoor fresh air flow rate for a sufficient time, which has led to a sharp increase in building energy consumption. Therefore, it is necessary to identify an energy-efficient ventilation strategy to reduce the risk of infection. In this study, we developed an occupant-number-based model predictive control (OBMPC) algorithm for building ventilation systems. First, we collected the occupancy and Heating, ventilation, and air conditioning system (HVAC) data from March to July 2021. Then, four different models (Auto regression moving average-based multilayer perceptron (ARMA_MLP), Recurrent neural networks (RNN), Long short-term memory networks (LSTM), and Nonhomogeneous Markov with change points detection (NH_Markov)) were used to predict the number of room occupants from 15 min to 24 h ahead with an interval output. We found that each model could predict the number of occupants with 85% accuracy using a one-person offset. Furthermore, the accuracy of 15 min of the ahead prediction could reach 95% with a one-person offset, but none of them could track abrupt changes. The occupancy prediction results were used to calculate the ventilation demand using the Wells-Riley equation, and the upper bound can maintain an infection risk lower than 2% for 93% of the day. This OBMPC model could reduce the coil load by 52.44% and shift the peak load by 3 h up to 5 kW compared with 24 × 7 h full outdoor air (OA) system when people wear masks in the space. The occupancy prediction uncertainty could cause a 9% to 26% difference in demand ventilation, a 0.3°C to 2.4°C difference in zone temperature, a 28.5% to 44.5% difference in outdoor airflow rate, and a 10.7% to 28.2% difference in coil load.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Anomalous Transport in Three-Dimensional Discrete Fracture Networks: Interplay Between Aperture Heterogeneity and Injection Modes

In this work, we study how the interplay between fracture aperture heterogeneity and tracer injection mode controls fluid flow and tracer transport in three-dimensional (3D) discrete fracture networks (DFNs). The direct 3-D DFN simulations show that tracer injection mode has substantial effects on tracer spreading across all levels of aperture heterogeneity. The key controlling factor for effective transport is the initial Lagrangian velocity distribution, which is determined by the interplay between injection mode and aperture heterogeneity. The fundamental difference between initial Lagrangian velocity distribution and domain-scale Eulerian velocity distribution plays a vital role in determining anomalous transport. We effectively capture the observed anomalous transport using an upscaled transport model that incorporates initial velocity distribution, stationary velocity distribution, velocity correlation length, and average advective tortuosity. With the upscaled transport model, we accurately capture the evolution of Lagrangian velocity distribution and predict longitudinal spreading in 3-D DFN.

58 GEOSCIENCES↗
DOE OSTI2021

Accurate Simulation of Both Sensitivity and Variability for Amazonian Photosynthesis: Is It Too Much to Ask?

Abstract Estimates of Amazon rainforest gross primary productivity (GPP) differ by a factor of 2 across a suite of three statistical and 18 process models. This wide spread contributes uncertainty to predictions of future climate. We compare the mean and variance of GPP from these models to that of GPP at six eddy covariance (EC) towers. Only one model's mean GPP across all sites falls within a 99% confidence interval for EC GPP, and only one model matches EC variance. The strength of model response to climate drivers is related to model ability to match the seasonal pattern of the EC GPP. Models with stronger seasonal swings in GPP have stronger responses to rain, light, and temperature than does EC GPP. The model to data comparison illustrates a trade‐off inherent to deterministic models between accurate simulation of a mean (average) and accurate responsiveness to drivers. The trade‐off exists because all deterministic models simplify processes and lack at least some consequential driver or interaction. If a model's sensitivities to included drivers and their interactions are accurate, then deterministically predicted outcomes have less variability than is realistic. If a GPP model has stronger responses to climate drivers than found in data, model predictions may match the observed variance and seasonal pattern but are likely to overpredict GPP response to climate change. High or realistic variability of model estimates relative to reference data indicate that the model is hypersensitive to one or more drivers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Design and Modeling of an 80 bar Oxy-Combustor for Direct Fired Supercritical CO2 Applications

As supercritical CO2 power cycles for fossil energy power generation continue to generate interest, there are significant issues and unanswered questions regarding injector design, flame stabilization, wall heat transfer, CO emissions, combustion dynamics and other combustion phenomenon. For natural gas, direct-fired cycles with carbon capture it is believed that Computational Fluid Dynamics (CFD) modeling will play an essential role in the combustor design process. To accurately model turbulent reacting flows at these unique conditions, experimental data is needed to validate CFD codes and sub-models and is currently lacking at conditions relevant for these cycles. This paper presents the conceptual design and CFD simulations of an experimental 80 bar oxy-combustion facility and test article currently under construction at NETL. The facility is targeted towards the testing of a single injector, direct-fired sCO2 combustor at the 100 kW thermal output level. While these conditions do not reflect the actual Allam cycle operating conditions (300 bar) they are viewed as a stepping stone in the model validation process at supercritical conditions. Reynolds Averaged Navier-Stokes as well as Large Eddy Simulation are used to model the turbulent combustion process and aid in the design of the combustor and injector. Process parameters including oxidizer preheat temperature and combustor flowrate are investigated.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Geometry-complete perceptron networks for 3D molecular graphs

Abstract Motivation The field of geometric deep learning has recently had a profound impact on several scientific domains such as protein structure prediction and design, leading to methodological advancements within and outside of the realm of traditional machine learning. Within this spirit, in this work, we introduce GCPNet, a new chirality-aware SE(3)-equivariant graph neural network designed for representation learning of 3D biomolecular graphs. We show that GCPNet, unlike previous representation learning methods for 3D biomolecules, is widely applicable to a variety of invariant or equivariant node-level, edge-level, and graph-level tasks on biomolecular structures while being able to (1) learn important chiral properties of 3D molecules and (2) detect external force fields. Results Across four distinct molecular-geometric tasks, we demonstrate that GCPNet’s predictions (1) for protein–ligand binding affinity achieve a statistically significant correlation of 0.608, more than 5%, greater than current state-of-the-art methods; (2) for protein structure ranking achieve statistically significant target-local and dataset-global correlations of 0.616 and 0.871, respectively; (3) for Newtownian many-body systems modeling achieve a task-averaged mean squared error less than 0.01, more than 15% better than current methods; and (4) for molecular chirality recognition achieve a state-of-the-art prediction accuracy of 98.7%, better than any other machine learning method to date. Availability and implementation The source code, data, and instructions to train new models or reproduce our results are freely available at https://github.com/BioinfoMachineLearning/GCPNet.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Jet quenching in the glasma phase: Medium-induced radiation

Inspired by the recent considerations of parton momentum broadening in the glasma phase, we study the medium-induced soft gluon radiation of jet partons at early times in heavy-ion collisions. The glasma state is assumed to be comprised of independent color domains with homogenous longitudinal fields that vary event by event, and we further complete this model with an event-averaging procedure accounting for the finite correlation length. Using this description, we evaluate the rate of medium-induced radiation from an energetic parton at midrapidity in the glasma phase. We mainly focus on SU(2) color fields for simplicity, also referring to the U(1) case and comparing with the Baier, Dokshitzer, Mueller, Peigné, Schiff, and Zakharov rate to gain further insight. Our results show that there is an intricate interplay of the synchrotron-like radiation in a single color domain with the destructive interference between different color domains, after the medium averaging is performed. Thus, we find that the emission rate is sensitive to the matter structure, decreasing for a glasma state populated by smaller color domains—i.e., for a glasma with a larger characteristic saturation scale. Our approach can be applied to more realistic backgrounds, and it sets the stage for the modeling of jet evolution in the early stages of heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

The influence of tillage and fertilizer on the flux and source of nitrous oxide with reference to atmospheric variation using laser spectroscopy

Nitrous oxide (N 2 O) is the third most important long-lived greenhouse gas and agriculture is the largest source of N 2 O emissions. Curbing N 2 O emissions requires understanding influences on the flux and sources of N 2 O. We measured flux and evaluated microbial sources of N 2 O using site preference ( S P ; the intramolecular distribution of 15 N in N 2 O) in flux chambers from a grassland tilling and agricultural fertilization experiments and atmosphere. We identified values greater than that of the average atmosphere to reflect nitrification and/or fungal denitrification and those lower than atmosphere as increased denitrification. Our spectroscopic approach was based on an extensive calibration with 18 standards that yielded S P accuracy and reproducibility of 0.7 ‰ and 1.0 ‰, respectively, without preconcentration. Chamber samples from the tilling experiment taken ~ monthly over a year showed a wide range in N 2 O flux (0–1.9 g N 2 O-N ha -1 d -1 ) and S P (- 1.8 to 25.1 ‰). Flux and S P were not influenced by tilling but responded to sampling date. Large fluxes occurred in October and May in no-till when soils were warm and moist and during a spring thaw, an event likely representing release of N 2 O accumulated under snow cover. These high fluxes could not be ascribed to a single microbial process as S P differed among chambers. However, the year-long S P and flux data for no-till showed a slight direct relationship suggesting that nitrification increased with flux. The comparative data in till showed an inverse relationship indicating that high flux events are driven by denitrification. Corn ( Zea mays ) showed high fluxes and S P values indicative of nitrification ~ 4 wk after fertilization with subsequent declines in S P indicating denitrification. Although there was no effect of fertilizer treatment on flux or S P in switchgrass ( Panicum virgatum) , high fluxes occurred ~1 month after fertilization. In both treatments, S P was indicative of denitrification in many instances, but evidence of nitrification/fungal denitrification also prevailed. At 2 m atmospheric N 2 O S P had a range of 31.1 ‰ and 14.6 ‰ in the grassland tilling and agricultural fertilization experiments, respectively. These data suggest the influence of soil microbial processes on atmospheric N 2 O and argue against the use of the global average atmospheric S P in isotopic modeling approaches.

Environmental Sciences & Ecology↗
DOE OSTI2026

Avoidance of disruptions on KSTAR due to vertical displacement events via novel real-time stability assessment

Disruption avoidance via the DECAF approach has been achieved on KSTAR using a novel real-time vertical stability assessment and a multiactuator feedback control strategy. The development of disruption avoidance strategies with reactor-relevant reliability is an urgent activity, enabling future fusion power plants. The stability metric employed is based on a new formulation of a vertical force gradient balance metric evaluated across the poloidal cross section of the plasma, with parameters tuned using historical data. Evaluation of this metric on a validation set of 400 recent KSTAR shots indicates >82% of Vertical displacement events can be avoided via feedback control. Essential to its calculation is the two-dimensional toroidal current density distribution in the plasma. Measurement of this profile faster than fully-converged equilibrium reconstructions can deliver is found to improve forecaster performance and is achieved with a surrogate model that takes as input magnetic diagnostic measurements and outputs the current profile on a basis comprising the top principal components of historical current profiles (from past equilibrium reconstructions). This method solves the non-uniqueness problem typically faced when reconstructing current profiles directly from diagnostics, while improving computational time and accuracy. On average, profiles produced by this model reach coefficients of determination of >0.99 with respect to those from equilibrium reconstructions. The avoidance actuators employed include poloidal field coils and an electron cyclotron current drive system. The multiactuator approach, as shown in this first demonstration, allows disruption avoidance while minimizing impact to operational performance. This ability, along with its flexibility and speed, makes this new approach an attractive option for avoiding these types of disruptions in reactors.

Tobin, Matthew [Columbia Univ., New York, NY (Unit↗
DOE OSTI2022

Implications of different income distributions for future residential energy demand in the U.S.

Abstract Future income distribution will affect energy demand and its interactions with various societal priorities. Most future model simulations assume a single average consumer and thus miss this important demand determinant. We quantify long-term implications of alternative future income distributions for state-level residential energy demand, investment, greenhouse gas, and pollutant emission patterns in the United States (U.S.) by incorporating income quintiles into the residential energy sector of the Global Change Analysis Model with 50-state disaggregation. We find that if the income distribution within each U.S. state becomes more egalitarian than present, what means that the difference on income between the richest and poorest decreases over time, residential energy demand could be 10% (4%–14% across states) higher in 2100. This increase of residential energy demand will directly reduce energy poverty, with a very modest increment on economywide CO 2 emissions (1%–2%). On the other hand, if U.S. states transition to a less equitable income distribution than present, with the difference between richest and poorest increasing over time, residential energy demand could be 19% (12%–26% across states) lower. While this study focuses on a single sector, we conclude that to improve understanding of synergies and tradeoffs across multiple societal goals such as energy access, emissions, and investments, future model simulations should explicitly consider subregional income distribution impacts.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

The Area Localized Coupled Model for Analytical Mean Flow Prediction in Arbitrary Wind Farm Geometries

This work introduces the area localized coupled (ALC) model, which extends the applicability of approaches that couple classical wake superposition models and atmospheric boundary layer models to wind farms with arbitrary layouts. Coupling wake and top–down boundary layer models is particularly challenging since the latter requires averaging over planform areas associated with turbine-specific regions of the flow that need to be specified. The ALC model uses Voronoi tessellation to define this local area around each turbine. A top–down description of a developing internal boundary layer is then applied over Voronoi cells upstream of each turbine to estimate the local mean velocity profile. Coupling between the velocity at hub-height based on this localized top–down model and a wake model is achieved by enforcing a minimum least-square-error in mean velocity in each cell. The wake model in the present implementation takes into account variations in wind farm inflow velocity and represents the wake profile behind each turbine as a super-Gaussian function that smoothly transitions between a top-hat shape in the region immediately following the turbine to a Gaussian profile downstream. Detailed comparisons to large-eddy simulation (LES) data from two different wind farms demonstrate the efficacy of the model in accurately predicting both wind farm power output and local turbine hub-height velocity for different wind farm geometries. These validations using data generated from two different LES codes demonstrate the model's versatility with respect to capturing results from different simulation setups and wind farm configurations.

49 EE - Wind and Water Power Program - Wind (EE-4W↗
DOE OSTI2022

An intermediate model for fitting triplet–triplet annihilation in phosphorescent organic light emitting diode materials

Triplet–triplet annihilation (TTA) is one of the primary contributors to efficiency roll-off and permanent material degradation in phosphorescent organic light-emitting diodes. Here, the two limiting case models typically used to quantify this quenching mechanism are multi-step Dexter and single-step Förster, which, respectively, assume ideal Fickian diffusion or perfect trapping of triplet excitons. For device-relevant guest doping levels (typically 5–12 vol. %), both significant diffusion of excitons and trapping due to spatial and energetic disorder exist, so neither conventional model fits experimental data well. We develop and validate an intermediate TTA model, which is a weighted average of the limiting cases of pure radiative decay (no TTA) and multi-step Dexter based TTA that returns an effective TTA rate constant and a parameter quantifying the portion of well-isolated excitons. Kinetic Monte-Carlo simulations and time-resolved photoluminescence measurements of an archetype host–guest system demonstrate that our intermediate model provides significantly improved fits with more realistic physical values, is more robust to variations in experimental conditions, and provides an analysis framework for the effects of trapping on TTA.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2022

Forecasting distributed energy resources adoption for power systems

Failing to incorporate accurate distributed energy resource penetration forecasts into long-term resource and transmission planning can lead to cost inefficiencies at best and system failures at worst. We have developed an open-source tool that employs an advanced Bass specification to calibrate and forecast technology adoption. The advanced specification includes geographic clustering, exogenously estimated market size, and dynamic time steps. Training on historical adoption of rooftop photovoltaics at the U.S. county-level and using detailed techno-economic estimates, our model achieves a two-year average mean-absolute-percentage-error of 19% in predicting system counts at the county-level, weighted by population. Model error was negatively correlated with market maturity - the error was 12% for counties in states with at least 28 W-per-capita of installed capacity. The advanced specification significantly reduces unweighted forecasting percent error compared to a conventional Bass specification: from 196% to 25% for capacity and from 226% to 22% for system count.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Transport in Stochastic Media with Random Chord Length Distributions

Thermal radiation transport computations in binary Markovian random mixtures rely almost exclusively on the Levermore-Pomraning (LP) model which is obtained by applying a heuristic closure to the ensemble averaged random medium transport equation. The validity of this model has been extensively tested by comparing numerical results over a broad parameter range (material types and mixing parameters) against benchmark solutions in planar geometry. The conditions under which the LP-model provides useful results and when it breaks down are now well established, but work to date has been largely restricted to homogeneous mixing statistics, i.e., the mean chord lengths of both materials are taken to be spatially constant. In recent work, this limitation was relaxed by allowing the mean chord lengths and, in a consistent fashion, the volume fractions in the LP-model to vary spatially and in a follow-up investigation benchmark solutions were obtained by ensemble averaging results over material realizations sampled from a nonhomogeneous Poisson process (NHPP). Numerical experiments in rod geometry with specifically linear and quadratic spatial dependence of chord lengths showed that the material averaged radiation intensities vary nonmonotonically with depth into the medium, in stark contrast to solutions obtained assuming uniform chord lengths. Moreover, depending on the local optical thickness and strength of scattering, the LP-model results showed locally more nuanced deviations from the benchmark solutions than was the case with constant chord lengths. These limited numerical investigations highlight the nontrivial qualitative and quantitative consequences of nonhomogeneous mixing statistics, in particular that closure approximations may not be uniformly valid or invalid over the problem domain.

42 ENGINEERING↗
DOE OSTI2024

Simulating energetic ions and enhanced fusion rates from ion-cyclotron resonance heating with a full-wave/Fokker–Planck model

Reproducing fast-ion enhanced fusion rates from ion-cyclotron resonance heating (ICRH) in tokamaks requires the self-consistent coupling of a full-wave solver and a Fokker–Planck solver, which evolves multiple simultaneously resonant ion species. We introduce a new self-consistent model that iterates the TORIC full-wave solver with the CQL3D Fokker–Planck solver using the integrated plasma simulator (IPS). This model evolves the bounce-averaged ion distribution functions in both parallel and perpendicular velocity-space with a quasilinear radio frequency (RF) diffusion operator valid in the ion finite Larmor radius (FLR) limit and the RF electric fields with the resultant non-Maxwellian FLR dielectric tensor. This produces non-Maxwellian ICRH simulations that are fully self-consistent, fast, and interoperable with integrated modeling frameworks, such as TRANSP/GACODE/IPS-FASTRAN. We demonstrate our model's capabilities by validating it against experimental data in Alcator C-Mod. We then perform the first RF heating simulations of SPARC using self-consistent non-Maxwellian ion distributions to investigate the potential to enhance fusion rates using ion cyclotron resonance heating generated fast ions.

Physics↗
DOE OSTI2025

Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident

This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Automatic detection of cataclysmic variables from SDSS images

Abstract Investigating rare and new objects have always been an important direction in astronomy. Cataclysmic variables (CVs) are ideal and natural celestial bodies for studying the accretion process of semi-detached binaries with accretion processes. However, the sample size of CVs must increase because a lager gap exists between the observational and the theoretical expanding CVs. Astronomy has entered the big data era and can provide massive images containing CV candidates. CVs as a type of faint celestial objects, are highly challenging to be identified directly from images using automatic manners. Deep learning has rapidly developed in intelligent image processing and has been widely applied in some astronomical fields with excellent detection results. YOLOX, as the latest YOLO framework, is advantageous in detecting small and dark targets. This work proposes an improved YOLOX-based framework according to the characteristics of CVs and Sloan Digital Sky Survey (SDSS) photometric images to train and verify the model to realise CV detection. We use the Convolutional Block Attention Module to increase the number of output features with the feature extraction network and adjust the feature fusion network to obtain fused features. Accordingly, the loss function is modified. Experimental results demonstrate that the improved model produces satisfactory results, with average accuracy (mean average Precision at 0.5) of 92.0%, Precision of 92.9%, Recall of 94.3%, and $F1-score$ of 93.6% on the test set. The proposed method can efficiently achieve the identification of CVs in test samples and search for CV candidates in unlabeled images. The image data vastly outnumber the spectra in the SDSS-released data. With supplementary follow-up observations or spectra, the proposed model can help astronomers in seeking and detecting CVs in a new manner to ensure that a more extensive CV catalog can be built. The proposed model may also be applied to the detection of other kinds of celestial objects.

Astronomy & Astrophysics↗
DOE OSTI2022

Verification and Validation Activities for the Multi-Fidelity Toolkit

The Multi-Fidelity Toolkit (MFTK) is a simulation tool being developed at Sandia National Laboratories for aerodynamic predictions of compressible flows over a range of physics fidelities and computational speeds. These models include the Reynolds-Averaged-Navier-Stokes (RANS) equations, the Euler equations, and modified Newtonian aerodynamics (MNA) equations, and they can be invoked independently or coupled with hierarchical Kriging to interpolate between high-fidelity simulations using lower-fidelity data. However, as with any new simulation capability, verification and validation are necessary to gather credibility evidence. This work describes formal code- and solution-verification activities as well as model validation with uncertainty considerations. Code verification is performed on the MNA model by comparing with an analytical solution for flat-plate and inclined-plate geometries. Solution-verification activities include grid-refinement studies of HIFiRE-1 wind tunnel measurements, which are used for validation, for all model fidelities. A thorough treatment of the validation comparison with prediction error and validation uncertainty is also presented.

42 ENGINEERING↗