Search NASA⌕ Search

SEARCH · Search NASA

Results for “conditional random fields”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Sensitivity Studies of the Inter-subchannel Mixing in Subchannel Analysis

In the PWR reactor core, the power profile is nonuniform and the inter-subchannel mixing is of great importance for flatting the temperature distribution of the coolant and for reactor core safety analysis. Subchannel codes are widely used in the reactor core safety analysis to assess the thermal hydraulic behavior of fuel assemblies. One of basic assumptions in subchannel analysis is that axial flow predominates in the channel flow. Under this assumption, the lateral momentum equation in subchannel codes can be simplified without bringing in significant errors. The inter-subchannel mixing leading to the exchange of mass, energy and momentum are modeled by adding source terms to the hydrodynamic equations or by separate models, and most of source terms and models are empirical correlations. With the advancement of mixing vane grid design, the lateral mixing has becoming a critical components in the evaluation of the inter-subchannel mixing in flow field. The accuracy of the subchannel analysis is strongly dependent on the modeling of inter-subchannel interactions between adjacent subchannels. So the inter-subchannel mixing models should be well considered and developed in subchannel analysis to better reflect their influence on flow field and temperature field in rod bundle. Generally speaking, the inter-subchannel mixing in subchannel analysis contains turbulent mixing and crossflow mixing in single phase flow. When it comes to two phase flow, the void drift will also contribute to inter-subchannel mixing [2]. Both the mechanisms and the generation conditions of each type of inter-subchannel mixing are different from each other. Turbulent mixing is a result of random turbulent flow and pressure fluctuations. Crossflow is caused by a net transverse pressure difference between adjacent subchannels. Lateral pressure differences can result from different subchannel hydraulic diameters, mass flux, heat flux distributions and changes in flow areas caused, for example, by wire wraps and spacer grids, respectively. Two-phase void drift is due to the tendency of the vapor phase to redistribute itself to a preferred equilibrium void distribution. This paper focuses on the sensitivity study of the inter-subchannel mixing with mixing vane grids for single phase flow in subchannel analysis. Both the turbulent mixing and crossflow mixing of the 5 x 5 rod bundle with mixing vane grids are predicted, and the results are compared and analyzed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Modeling and simulations of hydrodynamic shocks in a plasma flowing across randomized ICF scale laser beams

High-energy laser beams interacting with flowing plasmas can produce a plasma response that leads to deflection of the beam, beam bending. Such beams have usually a speckle structure generated by optical smoothing techniques that reduce the spatial and temporal coherence in the laser field pattern. The cumulative plasma response from laser speckles slows down the velocity of the incoming flow by momentum conservation. For slightly super-sonic flow the cumulative plasma response to the ponderomotive force exerted by the beam speckle ensemble is the strongest, such that slowing down the flow to subsonic velocities leads eventually to the generation of a shock around the cross section of the beam. This scenario has been predicted theoretically and is confirmed here by our hydrodynamic simulations in two dimensions with speckled beams and in one dimension with a reduced model. The conditions of shock generation are given in terms of the ponderomotive pressure, speckle size and the flow velocity. The nonlinear properties of the shocks are analyzed using Rankine–Hugoniot relations. According to linear theory, temporally smoothed laser beams exhibit a higher threshold for shock generation. Numerical simulations with beams that are smoothed by spectral dispersion compare well with the linear theory results, diverging from those produced by beams with only a random phase plates in the nonlinear regime. The conditions necessary for shock generation and their effects on the laser plasma coupling in inertial confinement fusion (ICF) experiments are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Efficient reconstruction and validation of heterogeneous microstructures for energy applications

The digital reconstruction of microstructures is necessary for simulations in fields ranging from geology to electrochemistry, but the state-of-the-art digital reconstruction techniques often compromise between resolution and field of view. It is challenging to retain detailed microstructure information in large-scale reconstructions. Here, this study investigates different aspects of the Yeong-Torquato algorithm based on correlation functions to make it more efficient. We achieve this goal by reducing the computational complexity of the chord-length distribution function and the two-point correlation function, applying the random sphere-packing method as the initial condition, and restricting potential voxel swaps to interfaces. In addition, a novel superposition parallel scheme is introduced to aid in searching for potential voxel swaps. The algorithm proposed is validated by comparing the pore-size distributions of reconstructed 3D custom battery electrodes from a sample dataset obtained from transmission X-ray microscopy. From a sample image with 200 x 200 pixels, the code can reconstruct a 300 x 300 x 300 structure in under 22 h and reconstruct a 400 x 400 x 400 structure in 43 h with eight cores.

42 ENGINEERING↗
DOE OSTI2015

Study of electron transport in a Hall thruster by axial–radial fully kinetic particle simulation

Electron transport across a magnetic field in a magnetic-layer-type Hall thruster was numerically investigated for the future predictive modeling of Hall thrusters. The discharge of a 1-kW-class magnetic-layer-type Hall thruster designed for high-specific-impulse operation was modeled using an r-z two-dimensional fully kinetic particle code with and without artificial electron-diffusion models. The thruster performance results showed that both electron transport models captured the experimental result within discrepancies less than 20% in thrust and discharge current for all the simulated operation conditions. The electron cross-field transport mechanism of the so-called anomalous diffusion was self-consistently observed in the simulation without artificial diffusion models; the effective electron mobility was two orders of magnitude higher than the value obtained using the classical diffusion theory. To account for the self-consistently observed anomalous transport, the oscillation of plasma properties was speculated. It was suggested that the enhanced random-walk diffusion due to the velocity oscillation of low-frequency electron flow could explain the observed anomalous diffusion within an order of magnitude. The dominant oscillation mode of the electron flow velocity was found to be 20 kHz, which was coupled to electrostatic oscillation excited by global ionization instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The direct and indirect drivers shaping RNA viral communities in grassland soil

Recent studies have revealed diverse RNA viral communities in soils. Yet, how environmental factors influence soil RNA viruses remains largely unknown. Here, we recovered RNA viral communities from 24 metatranscriptomes sequenced from grassland soils managed under a range of environmental conditions including 1) water content: 100% and 25% water holding capacity, 2) plant presence: planted with tall wheatgrass (Thinopyrum ponticum) and bare soil, 3) cultivar type: Alkar and Jose, and 4) soil depth: 0-5 cm and 15-25 cm. The recovered RNA viral communities were novel with nearly one-third of the RNA viral contigs uniquely detected in the studied grassland. The classified RNA viral contigs are mostly known as eukaryotic RNA viruses (74.7%) belonging to Phyla Duplornaviricota, Kitrinoviricota, Lenarviricota, and Pisuviricota. Eukaryotic RNA viruses of Family Mitoviridae as well as their natural hosts, Fungi, are one of the most dominant taxa. Consistent with the results of nonmetric multidimensional scaling analysis, the four environmental conditions (water content, plant presence, cultivar, and soil depth) significantly influence the assemblages of soil RNA viral communities as suggested by the correlation analysis and the random forest algorithm. The modularity analysis of the factor network and the structural equation modeling further support the hierarchical associations among the four environmental factors and the community factors representing the co-existing eukaryotic, prokaryotic, and RNA viral communities. The soil water content, plant presence, and type of cultivar demonstrate a significant positive impact on eukaryotic RNA viral richness directly as well as indirectly on eukaryotic RNA viral abundance via influencing the co-existing eukaryotic members in this soil. Our data also provide statistical support for the negative influence of soil depth on soil eukaryotic richness and abundances resulting in its indirect impact on soil eukaryotic RNA viral communities. This study provides field-relevant information on how environmental and community factors collectively shape soil RNA communities and contribute to ecological understanding of RNA viral survival under various environmental conditions and virus-host interactions in soil.

Wu, Ruonan↗
DOE OSTI2016

Investigating the effect of zirconium oxide microstructure on corrosion performance: a comparison between neutron, proton, and nonirradiated oxides - 2016-0069

Scanning precession electron diffraction in the transmission electron microscope has been used to simultaneously map the phase, orientation, and grain morphology of oxides formed on Zircaloy-2 after three and six cycles in a boiling water reactor in unprecedented detail. For comparison, a region of a preoxidized autoclave-formed oxide was also proton-irradiated at the Dalton Cumbrian Facility. The proton irradiation was observed to cause additional stabilization of the tetragonal phase that was attributed to the stabilizing effect of irradiation- induced defects in the oxide. In the reactor-formed oxides, no extra stabilization of the tetragonal grains was observed under neutron irradiation, as indicated by the similar tetragonal phase fraction and transformation twin-boundary distributions between the nonirradiated and reactor-formed oxides. It is suggested that the damage rate is too low in the newly formed oxide to cause significant stabilization of the tetragonal phase. This technique also reveals that the oxide formed under reactor conditions has a more heterogeneous microstructure, and the growth of well-oriented columnar monoclinic grains is significantly reduced compared with a nonirradiated oxide. High-angle annular dark-field scanning transmission electron microscopy also revealed the development of extensive networks of intergranular porosity and eventually grain de-cohesion in the reactor-formed oxides. These results suggest that the tetragonal-monoclinic transformation is not responsible for the accelerated corrosion exhibited under reactor conditions. It is proposed that the usual out-of-reactor oxide growth and nucleation processes are significantly modified under reactor conditions, resulting in a more heterogeneous and randomly oriented oxide microstructure with reduced columnar grain growth. It is suggested that this disordered oxide microstructure allows for the formation of extensive intergranular porosity that could lead to accelerated in-reactor corrosion. (authors)

36 MATERIALS SCIENCE↗
DOE OSTI2017

Field evaluation of remote wind sensing technologies: Shore-based and buoy mounted LIDAR systems

In developing a national energy strategy, the United States has a number of objectives, including increasing economic growth, improving environmental quality, and enhancing national energy security. Wind power contributes to these objectives through the deployment of clean, affordable and reliable domestic energy. To achieve U.S. wind generation objectives, the Wind and Water Power Program within the Department of Energy’s (DOE) Office of Energy Efficiency and Renewable Energy (EERE) instituted the U.S. Offshore Wind: Removing Market Barriers Program in FY 2011. Accurate and comprehensive information on offshore wind resource characteristics across a range of spatial and temporal scales is one market barrier that needs to be addressed through advanced research in remote sensing technologies. There is a pressing need for reliable offshore wind-speed measurements to assess the availability of the potential wind energy resource in terms of power production and to identify any frequently occurring spatial variability in the offshore wind resource that may impact the operational reliability and lifetime of wind turbines and their components and to provide a verification program to validate the “bankability” of the output of these alternative technologies for use by finance institutions for the financing of offshore wind farm construction. The application of emerging remote sensing technologies is viewed as a means to cost-effectively meet the data needs of the offshore wind industry. In particular, scanning and buoy mounted LIDAR have been proposed as a means to obtain accurate offshore wind data at multiple locations without the high cost and regulatory hurdles associated with the construction of offshore meteorological towers. However; before these remote sensing technologies can be accepted the validity of the measured data must be evaluated to ensure their accuracy. The proposed research will establish a unique coastal ocean test-bed in the Mid-Atlantic for the evaluation of LIDAR-based wind measurement systems to validate the accuracy of remotely measured wind data in marine applications. Specifically, the test-bed will be utilized to systematically evaluate the capability of emerging scanning LIDAR and buoy mounted vertically profiling LIDAR by: (1) Evaluating a fixed scanning LIDAR against land-based 50 and 60 meter high meteorological masts fitted with research quality cup-vane and/or sonic anemometers; (2) Evaluating a buoy mounted vertically profiling LIDAR fixed on land and floating in a sheltered bay against a co-located 60 meter high meteorological mast fitted with a research quality cup-vane and/or sonic anemometers and the fixed scanning LIDAR; and (3) Offshore field evaluation of both LIDAR platforms through a comparison of the fixed scanning LIDAR data and data obtained by the buoy mounted LIDAR located 10 miles offshore. The proposed research will systematically validate Light Detection and Ranging (LIDAR) based wind measurement systems and assess the temporal and spatial variability of the offshore wind resource in the Mid-Atlantic east of New Jersey. The goal of the proposed project is to address the technical and commercial challenges of the offshore wind energy industry by validating and assessing cost-effective, over ocean wind resource characterization technologies. The objective is to systematically evaluate the capability of both scanning and vertically profiling LIDARs to accurately measure 3D wind fields through comparison with fixed met masts and intercomparison among LIDAR platforms. Once validated, data collected by both buoy mounted vertically profiling LIDARs and shore-based, pulsed horizontally scanning LIDARs can be used to accurately assess offshore wind resources and to quantify the spatial and temporal variability in the offshore wind fields. One of the fundamental research questions to be addressed in phase 1 is the assessment of various measurement and data processing schemes to retrieve accurate wind vectors in the marine environment over large sampling ranges (10 to 12 km) and varying atmospheric aerosol levels. Atmospheric conditions and aerosol content within the coastal ocean region of the Mid-Atlantic seaboard of the US can vary significantly over short time periods in response to frontal passages and extratropical and tropical low pressure system passage offshore of the coast. Since aerosols provide the scattering medium for the determination of LIDAR Doppler shifts in the atmosphere the accuracy and range of LIDAR derived velocity measurements as a function of variation in aerosol content in the marine environment is a key research question to be addressed. In phase 1, it is desired to capture as much variation in atmospheric conditions and aerosol content as possible. To this end, collocated measurements of LIDAR and standard anemometer wind fields will be captured by the project PIs over all four seasons and during specific events (e.g., coastal low pressure system passage) in year 1. Additionally, since the meteorological masts are permanent structures, additional events can be captured over the three year duration of the field research project. All research instruments are owned by Fishermen’s Energy and made available to the PIs though a lease agreement as part of the DOE grant. Energy Fishermen’s Energy will be responsible for the operation and maintenance of the scanning LIDAR and met mast anemometers. On a daily basis, environmental data and systems performance indicators will be transmitted from each measurement station to the Fishermen’s project team consisting of both in-house personnel and equipment manufacturer engineers. Data sets include compiled LIDAR files as well as data sets from ancillary sensors. Diagnostic parameters to be monitored include standard deviations of measured values, battery levels and charging systems output, and the operational status. Once data have been confirmed as complete and reliable, files will be transferred to the Garrad Hassan (a subcontractor to Fishermen’s Energy) for incorporation in to the validation database, which is accessible to other scientific team members. Data collection times and durations will be determined by the PI and Co-PIs in consultation with instrument engineers to ensure the capture of data representative of the expected range of mid-Atlantic atmospheric conditions (e.g., temperature, moisture, coastal low pressure systems, tropical systems, rain, snow, fog). The collection and processing of the data is a function of site specific measurement requirements (Kelley et.al. 2007; Hannon et.al. 2008). To determine the optimal profiles of wind speed and direction from the LIDAR radial velocities as a function of azimuth angle, rigorous estimates of the bias and random error of each radial velocity estimate are required. Lockheed Martin Coherent Technologies, Inc., under contract with Fishermen’s Energy, will provide analyses of raw and processed data using various scan patterns to determine optimal performance settings for the pulsed scanning LIDAR. Once optimized, appropriate processing and analyses techniques will be evaluated by Garrad Hassan for use in validating the accuracy of the LIDAR wind field measurements against the standard anemometer measurements from the meteorological masts. The most attractive capability of the scanning LIDAR is the ability to provide high spatial resolution observations in a three-dimensional volume which provides superior statistical accuracy due to the large number of samples obtained. Each radial scan provides measurements in 100 range gates over a distance of 10 to 12 km at an update rate of 5 to 10 Hz and rotation of 2.5° per second. Each rotation at a fixed azimuth requires 2.4 minutes. Depending on the number of azimuths desired a complete scan can take up to 10 minutes or longer to complete. Once collected the radial velocities are processed to produce vector wind velocity estimates based on a set of data distributed in angle and range around points of interest, typically a standard grid within the radial wind map. To calculate wind vectors over a limited spatial area of interest for the comparison of data with other measurement platforms a localized least-squares approach has been applied by Hannon et.al. (2008) and a Variation Assimilation (VAR) processing technique has been applied by Chan and al Assimilation (VAR) processing technique has been applied by Chan and Shao (2006). Additionally, Kelley et.al. (2007) applied a “stare” technique that fixed a scanning LIDAR in both azimuth and elevation angles to measure over collocated sampling volumes of the LIDAR and a 3D sonic anemometer mounted to a fixed mast. Although the stare technique is limited to wind directions aligned with the sampling radial of the LIDAR, it does provide a direct comparison of sampling volumes. Each of the processing techniques described above (and possibly others) will be evaluated to determine the validity of the LIDAR derived wind fields in the marine environment. Numerical methods such as linear regression and comparison of probability density functions of wind fields measured by each instrument platform will be used to assess the processing techniques. Linear regression has the advantage of directly evaluating corresponding pairs of wind data measured by each instrument and can lend insight into deviations and bias between instruments as a function of wind speed. Assessment of the coherency between the probability density function of the wind measured by each instrument provides insight into processes that may not be accurately resolved by each instrument at specific frequencies. Once processed and assessed the most appropriate technique will be utilized to provide valid wind measurements from the pulsed scanning LIDAR. A detailed analysis of the measurement data from the LIDAR and the three meteorological towers and comparison of the coastal wind characteristics from the different systems will be performed by the CO_PI at NREL. This analysis will evaluate how the WT LIDAR performance and measurement of the wind characteristics vary with distance from the LIDAR and by atmospheric conditions; using the tower measurements at different distances (4.8, 9.6, and 19 km) from the LIDAR as a reference. The comparative analysis will include, to the extent possible, evaluation of parameters such as wind speed and direction distributions, wind shear, turbulence intensity and their variations by atmospheric conditions, month or season, and time of day.

17 WIND ENERGY↗
DOE OSTI2024

Mechanisms of Waterflood Inefficiency: Analysis of Geological, Petrophysical and Reservoir History, a Field Case Study of FWU (East Section)

The petroleum reservoir represents a complex heterogeneous system that requires thorough characterization prior to the implementation of any incremental recovery technique. One of the most commonly utilized and successful secondary recovery techniques is waterflooding. However, a lack of sufficient investigation into the inherent behavior and characteristics of the reservoir formation in situ can result in failure or suboptimal performance of waterflood operations. Therefore, a comprehensive understanding of the geological history, static and dynamic reservoir characteristics, and petrophysical data is essential for analyzing the mechanisms and causes of waterflood inefficiency and failure. In this study, waterflood inefficiency was observed in the Morrow B reservoir located in the Farnsworth Unit, situated in the northwestern shelf of the Anadarko Basin, Texas. To assess the potential mechanisms behind the inefficiency of waterflooding in the east half, geological, petrophysical, and reservoir engineering data, along with historical information, were integrated, reviewed, and analyzed. The integration and analysis of these datasets revealed that several factors contributed to the waterflood inefficiency. Firstly, the presence of abundant dispersed authigenic clays within the reservoir, worsened by low reservoir quality and high heterogeneity, led to unfavorable conditions for waterflood operations. The use of freshwater for flooding exacerbated the adverse effects of sensitive and migratory clays, further hampering the effectiveness of the waterflood. In addition to these factors, several reservoir engineering issues played a significant role in the inefficiency of waterflooding. These issues included inadequate perforation strategies due to the absence of detailed hydraulic flow units (HFUs) and rock typing, random placement of injectors, and uncontrolled injected fresh water. These external controlling parameters further contributed to the overall inefficiencies observed during waterflood operations in the east half of the reservoir. A detailed understanding of the mechanistic factors of inefficient waterflood operation will provide adequate insights into the development of the improved recovery technique for the field.

Morgan, Anthony (ORCID:0000000211519153)↗
DOE OSTI2015

Coagulation kinetics beyond mean field theory using an optimised Poisson representation

Binary particle coagulation can be modelled as the repeated random process of the combination of two particles to form a third. The kinetics may be represented by population rate equations based on a mean field assumption, according to which the rate of aggregation is taken to be proportional to the product of the mean populations of the two participants, but this can be a poor approximation when the mean populations are small. However, using the Poisson representation, it is possible to derive a set of rate equations that go beyond mean field theory, describing pseudo-populations that are continuous, noisy, and complex, but where averaging over the noise and initial conditions gives the mean of the physical population. Such an approach is explored for the simple case of a size-independent rate of coagulation between particles. Analytical results are compared with numerical computations and with results derived by other means. In the numerical work, we encounter instabilities that can be eliminated using a suitable “gauge” transformation of the problem [P. D. Drummond, Eur. Phys. J. B 38, 617 (2004)] which we show to be equivalent to the application of the Cameron-Martin-Girsanov formula describing a shift in a probability measure. The cost of such a procedure is to introduce additional statistical noise into the numerical results, but we identify an optimised gauge transformation where this difficulty is minimal for the main properties of interest. For more complicated systems, such an approach is likely to be computationally cheaper than Monte Carlo simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Influence of extreme fracture flow channels on the thermal performance of open-loop geothermal systems at commercial scale

Abstract Adequate stewardship of geothermal resources requires accurate forecasting of long-term thermal performance. In enhanced geothermal systems and other fracture-dominated reservoirs, predictive models commonly assume constant-aperture fractures, although spatial variations in aperture can greatly affect reservoir permeability, fluid flow distribution, and heat transport. Whereas previous authors have investigated the effects of theoretical random aperture distributions on thermal performance, here we further explore the influence of permeability heterogeneity considering field-constrained aperture distributions from a meso-scale field site in northern New York, USA. Using numerical models of coupled fluid flow and heat transport, we conduct thermal–hydraulic simulations for a hypothetical reservoir consisting of a relatively impervious porous matrix and a single, horizontal fracture. Our results indicate that in highly channelized fields, most well design configurations and operating conditions result in extreme rates of thermal drawdown (e.g., 50% drop in production well temperatures in under 2 years). However, some other scenarios that account for the risks of short-circuiting can potentially enhance heat extraction when mass flow rate is not excessively high, and the direction of geothermal extraction is not aligned with the most permeable features in the reservoir. Through a parametric approach, we illustrate that well separation distance and relative positioning play a major role in the long-term performance of highly channelized fields, and both can be used to help mitigate premature thermal breakthrough.

Rangel-Jurado, Nicolás↗
DOE OSTI2021

A novel random forest approach to revealing interactions and controls on chlorophyll concentration and bacterial communities during coastal phytoplankton blooms

Abstract Increasing occurrence of harmful algal blooms across the land–water interface poses significant risks to coastal ecosystem structure and human health. Defining significant drivers and their interactive impacts on blooms allows for more effective analysis and identification of specific conditions supporting phytoplankton growth. A novel iterative Random Forests (iRF) machine-learning model was developed and applied to two example cases along the California coast to identify key stable interactions: (1) phytoplankton abundance in response to various drivers due to coastal conditions and land-sea nutrient fluxes, (2) microbial community structure during algal blooms. In Example 1, watershed derived nutrients were identified as the least significant interacting variable associated with Monterey Bay phytoplankton abundance. In Example 2, through iRF analysis of field-based 16S OTU bacterial community and algae datasets, we independently found stable interactions of prokaryote abundance patterns associated with phytoplankton abundance that have been previously identified in laboratory-based studies. Our study represents the first iRF application to marine algal blooms that helps to identify ocean, microbial, and terrestrial conditions that are considered dominant causal factors on bloom dynamics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Phase field model of uranium carbide solidification through a combined KKS and orientation field approach

A Phase Field model is developed combining the Orientation Field approach to modeling solidification with the Kim, Kim, Suzuki (KKS) method of modeling binary alloys. These combined methods produce a model capable of simulating randomly oriented second phase dendrites with discrete control of the solid-liquid interface energy and thickness. The example system of carbon (C) in a liquid uranium (U) melt is used as a test for the model. The formation of uranium carbide (UC) within a liquid U melt is simulated for isothermal conditions and compares well with experiments.

Bair, Jacob L. (ORCID:0000000333273323)↗
DOE OSTI2022

High Fidelity CFD Simulations Supporting the KP-FHR

Kairos Power, LLC, is developing its version of the Fluoride-cooled High-temperature Reactor, the KP-FHR. The design uses a pebble bed core with fluoride salt as a coolant. The pebbles used in the KP-FHR have a diameter of 4 cm, with a shell fuel region where TRISO particles are embedded. A Pebble bed core design is adopted by several Gen IV reactors, They boast many benefits, such as fuel integrity, highly efficient heat transfer, and passive safety. However, it is challenging to accurately predict temperature and flow inside a pebble bed. Traditional approaches use the porous media model, which regards the pebble bed as a continuous medium, but with different temperature fields representing different levels, such as the fluid temperature, pebble surface temperature, and pebble center temperature. Empirical heat transfer correlations are adopted to calculate the heat transfer coefficient between different phases. However, empirical correlations are usually validated with experimental data, which usually lacks detail inside the pebble bed. The available experimental data is also generally at a high Reynolds number, which falls outside of the conditions of KP-FHR. Explicit computational fluid dynamics (CFD) simulations of randomly packed pebble beds have only become feasible recently. This is thanks to the rapid development of computational power and scalable algorithms. In this work, we used the Spectral Element Method (SEM) CFD code NekRS to simulate the randomly packed pebble bed in a cylindrical container. NekRS, which is the GPU variant of Nek5000, but refactored to utilize the computational power of GPUs using the OCCA library to run on hybrid architecture high performance computing systems. It was initially developed with the libParamunal library, but truncated and tuned for large-scale turbulence simulation. As a result, the SEM reaches higher precision with the same degrees of freedom by using a high-order Lagrange polynomial basis distributed on Gauss-Lobatto-Legendre quadrature inside each element, compared to lower-order methods, such the Finite Volume Method and Finite Element Method. The report is divided into five parts. We start with a general discussion of the pebble bed reactor, along with a specific investigation into the KP-FHR. The second part presents the numerical methodology. In the third part, we study a modular pebble bed with 1741 pebbles in a container of 7 pebble-diameter radius. Beyond LES simulations done by NekRS, we also leveraged the thermal radiation model in OpenFOAM to study heat transfer under no-forced-flow scenarios. Then, in the fourth part we simulated a pebble bed similar to the size of the Hermes Test Reactor. The total number of pebbles is in these simulations is 34,374. The container radius is 14 pebble-diameters. Finally, the report concludes in part five, with a discussion of future work.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2018

Predicting departure from nucleate boiling with advanced data- and physics-driven approaches

The subcooled and low-quality saturated flow boiling critical heat flux (CHF) corresponding to the departure from nucleate boiling (DNB) is one of the major limiting factors in the design and operation of pressurized water reactors (PWRs). Due to lack of agreement on its physical mechanisms, development of an accurate DNB-type CHF model has been elusive. The various steady-state (SS) predictive tools that have been proposed can be grouped into two categories: (1) data-driven methods—namely best-fit empirical correlations and look-up tables that result in relatively good agreement with specific experimental datasets but often fail to extend beyond their ranges of validity, and (2) the physics-based mechanistic models that rely on reasonable yet limited understanding of the underlying physics and apply constitutive relations to close the conservation equations. On one hand, with recent advances in computational capabilities and optimization techniques, machine learning (ML)-based methods constitute an alternative and a more advanced, data-driven approach. With enough data, they can be particularly useful in engineering fields where the physical phenomena are complex and challenging to model. On the other hand, a state-of-the-art physics-driven approach would not only provide more accurate prediction of SS CHF, but it would also help better understand and model transient applications. The data-driven regression models use deep feed-forward neural networks (NNs) and random forests to cross-validate with 1,865 CHF test cases, covering a wide range of flow conditions and channel geometries. The best-estimate ML-based predictors compare favorably with the widely used look-up table for annulus and plate, and sensitivity analysis has confirmed their effectiveness. The key advantage of ML-based methods is their online extensibility of applicability domain. The proposed physics-driven approach combines key assumptions and parameters in the relatively well-accepted theories of liquid sublayer dryout and near-wall bubble crowding. A more realistic understanding of local mechanisms has been modeled. The new model has been optimized and validated against 1,439 tube data, and it shows considerably improved performance when compared to a recent mechanistic model, demonstrating unbiased close agreement with measurements over a wide range of operating conditions. Future work will extend the improved physics-driven model to non-tube geometry applications and will couple the model with ML via a hybrid approach. The mechanistic SS work will also serve to improve understanding and to model transient CHF scenarios.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2015

Multiscale simulation of xenon diffusion and grain boundary segregation in UO₂

In light water reactor fuel, gaseous fission products segregate to grain boundaries, resulting in the nucleation and growth of large intergranular fission gas bubbles. The segregation rate is controlled by diffusion of fission gas atoms through the grains and interaction with the boundaries. Based on the mechanisms established from earlier density functional theory (DFT) and empirical potential calculations, diffusion models for xenon (Xe), uranium (U) vacancies and U interstitials in UO₂ have been derived for both intrinsic (no irradiation) and irradiation conditions. Segregation of Xe to grain boundaries is described by combining the bulk diffusion model with a model for the interaction between Xe atoms and three different grain boundaries in UO₂ (Σ5 tilt, Σ5 twist and a high angle random boundary), as derived from atomistic calculations. The present model does not attempt to capture nucleation or growth of fission gas bubbles at the grain boundaries. The point defect and Xe diffusion and segregation models are implemented in the MARMOT phase field code, which is used to calculate effective Xe and U diffusivities as well as to simulate Xe redistribution for a few simple microstructures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

On the Prediction of Aerosol-Cloud Interactions Within a Data-Driven Framework

Aerosol-cloud interactions (ACI) pose the largest uncertainty for climate projection. Among many challenges of understanding ACI, the question of whether ACI can be deterministically predicted has not been explicitly answered. Here we attempt to answer this question by predicting cloud droplet number concentration N c from aerosol number concentration N a and ambient conditions using a data-driven framework. We use aerosol properties, vertical velocity fluctuations, and meteorological states from the ACTIVATE field observations (2020–2022) as predictors to estimate N c . We show that the campaign-wide N c can be successfully predicted using machine learning models despite the strongly nonlinear and multi-scale nature of ACI. However, the observation-trained machine learning model fails to predict N c in individual cases while it successfully predicts N c of randomly selected data points that cover a broad spatiotemporal scale. This suggests that, within a data-driven framework, the N c prediction is uncertain at fine spatiotemporal scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The quality/cosmology tension for a post-inflation QCD axion

Abstract It is difficult to construct a post-inflation QCD axion model that solves the axion quality problem (and hence the Strong CP problem) without introducing a cosmological disaster. In a post-inflation axion model, the axion field value is randomized during the Peccei-Quinn phase transition, and axion domain walls form at the QCD phase transition. We emphasize that the gauge equivalence of all minima of the axion potential (i.e., domain wall number equals one) is insufficient to solve the cosmological domain wall problem. The axion string on which a domain wall ends must exist as an individual object (as opposed to a multi-string state), and it must be produced in the early universe. These conditions are often not satisfied in concrete models. Post-inflation axion models also face a potential problem from fractionally charged relics; solving this problem often leads to low-energy Landau poles for Standard Model gauge couplings, reintroducing the quality problem. We study several examples, finding that models that solve the quality problem face cosmological problems, and vice versa. This is not a no-go theorem; nonetheless, we argue that it is much more difficult than generally appreciated to find a viable post-inflation QCD axion model. Successful examples may have a nonstandard cosmological history (e.g., multiple types of cosmic axion strings of different tensions), undermining the widespread expectation that the post-inflation QCD axion scenario predicts a unique mass for axion dark matter.

Physics↗
DOE OSTI2017

Compositional phase diagram and microscopic mechanism of Ba 1-x Ca x Zr y Ti 1-y O 3 relaxor ferroelectrics

With extensive first-principles density-functional calculations, we construct a three-dimensional compositional phase diagram of Ba1-xCaxZryTi1-yO3 (BCZT) with the Ca and Zr content in the ranges of 0 ≤ xCa ≤ 0.2 and 0 ≤ yZr ≤ 1. Our calculations show that, when the Zr content increases, the difference in energy and difference in the structural parameters of the cubic, tetragonal, orthorhombic, and rhombohedral phases of BCZT are reduced. Eventually, all four phases merge into a multiphase with coexisting cubic structures (MPCCS) under Zr-rich conditions, indicating that BCZT undergoes phase transition from a normal ferroelectric (NFE) to a relaxor ferroelectric (RFE), consistent with experimental observations. The 3D diagram shows that the regions of merged and separated energy surfaces correspond to the regions of the RFE and NFE, respectively, which suggests that a MPCCS corresponds to a RFE. In addition, with the MPCCS model and Landau–Devonshire theory, we provide an interpretation of the high electromechanical properties of the BCZT relaxor ferroelectric and apply it to the classical local random field and micro–macro domain transition models.

Liu, Shi-Yu↗