Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conditional random field”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS2001

Design and Calibration of Autonomous Coherent Doppler Lidar for Space Missions

Developed a new algorithm for the simulation of three dimensional homogeneous turbulent velocity fields. For typical atmospheric conditions it is impossible to produce a simulated velocity field that simultaneously satisfy a given spatial correlation and the corresponding spatial spectrum because of spectral aliasing. The new algorithms produce a turbulent velocity field which has accurate spatial correlations which is required for performance predictions from space-based systems. Developed a new algorithm for extracting the spatial statistics of the atmospheric velocity field using coherent Doppler lidar. The performance of the algorithm was compared with past methods and the new algorithm produces useful results for space-based data, which was not possible before. Developed new methods for verification of the errors in ground-based and space-based Doppler lidar wind measurements. These new methods do not require independent in situ data. This is an important issue for the verification of space-based Doppler lidar measurements of the global wind field. The performance of the new algorithm was compared with past results for both space-based and ground-based operation. The new algorithm has the best performance and is the only algorithm that performed satisfactory for spacebased operation. The performance of coherent Doppler lidar for a space missions with various scanning geometries was determined using computer simulation which contained the effects of random instrumental velocity errors, wind shear, wind variability along the range-gate and from shot-to-shot, and random variations in atmospheric aerosol backscatter over the measurement volume. The bias in the velocity estimates was small and the accuracy in the is typically less than 0.5 m/s for high signal conditions. For a large number of shot per velocity estimate, the threshold signal level for acceptable estimates is proportional to the number of shots to the minus one half power. This agrees with previous results determined for ground-based measurements. The use of multi-element optical detectors for autonomous operation of coherent Doppler lidar was shown to be a very promising technique. Optimal detector geometries were determined by computer simulation of performance: for ground-based testing with a fixed calibration target and for space-based operation using the random surface returns. The effects of refractive turbulence on ground-based calibration of coherent Doppler lidar was determined by computer simulations and compared with theoretical predictions. New techniques were required to correctly predict performance for the focused beam geometry commonly used for verification of space-based operation. An improved velocity estimator was evaluated for space-based applications were signal shot measurements are used to produce vector wind measurements. This permits more accurate measurements when the signal level is not known a priori or not available from multiple shot measurements. The average Doppler lidar signal spectrum including the effects of velocity turbulence was derived and calculated. This permits new estimation algorithms for turbulence based on spectral estimates. In situ atmospheric measurements were conducted and analyzed using an instrumented kite-platform. This work helps provide the required in situ data for verification of Doppler lidar velocity statistics.

Frehlich, Rod G.↗
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↗
NASA NTRS1988

A noise source identification technique using an inverse Helmholtz integral equation method

A technique is developed which utilizes numerical models and field pressure information to characterize acoustic fields and identify acoustic sources. The numerical models are based on boundary element numerical procedures. Either pressure, velocity, or passive boundary conditions, in the form of impedance boundary conditions, may be imposed on the numerical model. Alternatively, if no boundary information is known, a boundary condition can be left unspecified. Field pressure data may be specified to overdetermine the numerical problem. The problem is solved numerically for the complete sound field from which the acoustic sources may be determined. The model can then be used to idenfify acoustic intensity paths in the field. The solution can be modified and the model used to evaluate design alternatives. In this investigation the method is tested analytically and verified. In addition, the sensitivity of the method to random and bias error in the input data is demonstrated.

Gardner, B. K.↗
NASA NTRSDate not supplied

Quantifying Post-fire Vegetation Recovery on the Colorado Front Range

Forest composition and structure in the Colorado Front Range has been altered by changing wildfire regimes, particularly in response to increased fire severity. Subsequent reductions in post-fire tree regeneration can result in chronic impacts to ecological function and water quality. This project partnered with the U.S. Forest Service to estimate tree canopy recovery from 15 to 21 years after four Colorado Front Range fires that burned between 1996 and 2002— the Bobcat, Buffalo Creek, Hayman, and High Meadows fires—using Landsat 5 Thematic Mapper (TM), Landsat 7 Enhanced Thematic Mapper (ETM+), and Landsat 8 Operational Land Imager (OLI). First, we evaluated relationships between recovery metrics derived from spectral vegetation indices (i.e. post-fire vegetation, normalized burn ratio, percent unrecovered vegetation) and field-collected counts of post-fire tree seedling regeneration. Then, we used high-resolution imagery paired with Landsat timeseries spectral and geomorphometric data to map forest canopy cover 15-21 years post-fire. Finally, we assessed ecological drivers of recovery including topography, climate, soils, and fire severity using Random Forest. While we found that basic vegetation index metrics failed to distinguish field-measured tree regeneration, our data showed that both aspect and burn severity were important drivers of vegetation recovery. North-facing aspects tended to have higher pre-fire Normalized Difference Vegetation Index (NDVI) and thus higher burn severities. However, areas that burned at high severity experienced slower vegetation recovery, regardless of aspect suggesting that high severity fire can mute favorable site conditions on north-facing aspects (i.e. moister, more developed soils). While this work propels our understanding of variables that influence vegetative recovery, vegetation type, forest conversion, and watershed function, more work should focus on what vegetation indices identifying areas of recovery indicate across different forest types and burn severities.

NASA DEVELOP↗
NASA NTRS1998

Virtual and Experimental Visualization of Flows in Packed Beds of Spheres Simulating Porous Media Flows

A videotape presentation of flow in a packed bed of spheres is provided. The flow experiment consisted of three principal elements: (1) an oil tunnel 76.2 mm by 76.2 mm in cross section, (2) a packed bed of spheres in regular and irregular arrays, and (3) a flow characterization methodology, either (a) full flow field tracking (FFFT) or (b) computational fluid dynamic (CFD) simulation. The refraction indices of the oil and the test array of spheres were closely matched, and the flow was seeded with aluminum oxide particles. Planar laser light provided a two-dimensional projection of the flow field, and a traverse simulated a three-dimensional image of the entire flow field. Light focusing and reflection rendered the spheres black, permitting visualization of the planar circular interfaces in both the axial and transverse directions. Flows were observed near the wall-sphere interface and within the set of spheres. The CFD model required that a representative section of a packed bed be formed and gridded, enclosing and cutting six spheres so that symmetry conditions could be imposed at all cross-boundaries. Simulations had to be made with the flow direction at right angles to that used in the experiments, however, to take advantage of flow symmetry. Careful attention to detail was required for proper gridding. The flow field was three-dimensional and complex to describe, yet the most prominent finding was flow threads, as computed in the representative 'cube' of spheres with face symmetry and conclusively demonstrated experimentally herein. Random packing and bed voids tended to disrupt the laminar flow, creating vortices.

Hendricks, R. C.↗
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↗
NASA NTRS1993

Advanced Multi-Phase Flow CFD Model Development for Solid Rocket Motor Flowfield Analysis

It is known that the simulations of solid rocket motor internal flow field with AL-based propellants require complex multi-phase turbulent flow model. The objective of this study is to develop an advanced particulate multi-phase flow model which includes the effects of particle dynamics, chemical reaction and hot gas flow turbulence. The inclusion of particle agglomeration, particle/gas reaction and mass transfer, particle collision, coalescence and breakup mechanisms in modeling the particle dynamics will allow the proposed model to realistically simulate the flowfield inside a solid rocket motor. The Finite Difference Navier-Stokes numerical code FDNS is used to simulate the steady-state multi-phase particulate flow field for a 3-zone 2-D axisymmetric ASRM model and a 6-zone 3-D ASRM model at launch conditions. The 2-D model includes aft-end cavity and submerged nozzle. The 3-D model represents the whole ASRM geometry, including additional grain port area in the gas cavity and two inhibitors. FDNS is a pressure based finite difference Navier-Stokes flow solver with time-accurate adaptive second-order upwind schemes, standard and extended k-epsilon models with compressibility corrections, multi zone body-fitted formulations, and turbulence particle interaction model. Eulerian/Lagrangian multi-phase solution method is applied for multi-zone mesh. To simulate the chemical reaction, penalty function corrected efficient finite-rate chemistry integration method is used in FDNS. For the AL particle combustion rate, the Hermsen correlation is employed. To simulate the turbulent dispersion of particles, the Gaussian probability distribution with standard deviation equal to (2k/3)(exp 1/2) is used for the random turbulent velocity components. The computational results reveal that the flow field near the juncture of aft-end cavity and the submerged nozzle is very complex. The effects of the turbulent particles affect the flow field significantly and provide better prediction of the ASRM performance. The multi-phase flow analysis using the FDNS code in the present research can be used as a design tool for solid rocket motor applications.

Liaw, Paul↗
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↗