SEARCH · Search NASA

Results for “Work Domain Analysis”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,081 records · Page 60

All sources NASA NTRS DOE OSTI

DOE OSTI2016

Characterization of 3D Cirrus Cloud and Radiation Fields Using ARS/AIRS/MODIS data and its Application to Climate Model

During the report period, we have made the following research accomplishments. First, we performed analysis for a number of MODIS scenes comprising of heavy dust events and ice clouds, covering regions of frequent dust outbreaks in East Asia, Middle East, and West Africa, as well as areas associated with long-range dust transports over the Equatorial Tropical Atlantic Ocean. These scenes contain both dust/aerosols and clouds. We collected suitable aerosol/ice-cloud data, correlated ice cloud and aerosol parameters by means of statistical analysis, and interpreted resulting correlation trends based on the physical principles governing cloud microphysics. Aerosol and cloud optical depths and cloud effective particle size inferred from MODIS for selected domains were analyzed from which the parameters including dust aerosol number concentration, ice cloud water path, and ice particle number concentration were subsequently derived. We illustrated that the Twomey (solar albedo) effect can be statistically quantified based on the slope of best-fit straight lines in the correlation study. Analysis of aerosol and cloud retrieval products revealed that for all cases, the region with a larger dust aerosol optical depth is always characterized by a smaller cloud particle size, consistent with the Twomey hypothesis for aerosol-cloud interactions. Second, we developed mean correlation curves with uncertainties associated with small ice-crystal concentration observations for the mean effective ice crystal size (De) and ice water content (IWC) by dividing the atmosphere into three characteristic regions: Tropics cirrus, Midlatitude cirrus, including a temperature classification to improve correlation, and Arctic ice clouds. We illustrated that De has a high correlation with IWC based on theoretical consideration and analysis of thousands of observed ice crystal data obtained from a number of ARM-DOE field campaigns and other experiments. The correlation has the form: ln(De) = a + b ln(IWC) + c ((ln(IWC))2, where a, b, and c are fitting coefficients and are functions of three regions. We demonstrated that this correlation can be effectively incorporated in GCMs and climate models that predict IWC - a significant advance in ice microphysics parameterization for interactive cloud-radiation analysis and feedback. Substantial July mean differences are shown in the OLR (W/m2) and precipitation (mm/day) patterns between UCLA GCM simulations based on Des determined from the De-IWC correlations and the control run using a fixed ice crystal size. Third, in order to improve the computation of spectral radiative transfer processes in the WRF model, we developed a consistent and efficient radiation scheme that can better resolve the spectral bands, determine the cloud optical properties, and provide more reliable and accurate radiative heating fields. In the newly developed radiation module, we have implemented in WRF a modified and improved version referred to as the Fu-Liou-Gu scheme, which includes a combination of delta-four-stream and delta-two-stream approximations for solar and IR flux calculations, respectively. This combination has been proven to be computationally efficient and at the same time to produce a high degree of accuracy. The incorporation of nongray gaseous absorption in multiple scattering atmospheres was based on the correlated k-distribution method. The solar and IR spectra are divided into 6 and 12 bands, respectively, according to the location of absorption bands of H2O, CO2, O3, CH4, N2O, and CFCs. We further included absorption by the water vapor continuum and a number of minor absorbers in the solar spectrum leading to an additional absorption of solar flux in a clear atmosphere on the order of 1-3 W/m2. Additionally, we incorporated the ice microphysics parameterization that includes an interactive mean effective ice crystal size in association with radiation parameterizations. The Fu-Liou-Gu scheme is an ideal tool for the simulation of radiative transfer and ice microphysics within the domain of WRF. It is particularly useful for studying direct and indirect aerosol radiative effects associated with ice cloud formation. The newly implemented radiation module has been demonstrated to work well in WRF and can be effectively used for studies related to cirrus cloud formation and evolution as well as aerosol-cloud-radiation interactions. With the newly implemented radiation scheme, the simulations of cloud cover and ice water path have been improved for cirrus clouds, with a more consistent comparison with the corresponding MODIS observations, especially for optically thin cirrus with an improvement of about 20% in the simulated mean ice water path.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2018

Developing Intensity-Duration-Frequency (IDF) Curves From Satellite-Based Precipitation: Methodology and Evaluation

Given the continuous advancement in the retrieval of precipitation from satellites, it is important to develop methods that incorporate satellite-based precipitation data sets in the design and planning of infrastructure. This is because in many regions around the world, in situ rainfall observations are sparse and have insufficient record length. A handful of studies examined the use of satellite-based precipitation to develop intensity-duration-frequency (IDF) curves; however, they have mostly focused on small spatial domains and relied on combining satellite-based with ground-based precipitation data sets. In this study, we explore this issue by providing a methodological framework with the potential to be applied in ungauged regions. This framework is based on accounting for the characteristics of satellite-based precipitation products, namely, adjustment of bias and transformation of areal to point rainfall. The latter method is based on previous studies on the reverse transformation (point to areal) commonly used to obtain catchment-scale IDF curves. The paper proceeds by applying this framework to develop IDF curves over the contiguous United States (CONUS); the data set used is Precipitation Estimation from Remotely Sensed Information Using Artificial Neural Networks – Climate Data Record (PERSIANN-CDR). IDFs are then evaluated against National Oceanic and Atmospheric Administration (NOAA) Atlas 14 to provide a quantitative estimate of their accuracy. Results of this work show that median errors are in the range of (17–22%), (6–12%), and (3–8%) for one-day, two-day and three-day IDFs, respectively, and return periods in the range (2–100) years. Furthermore, a considerable percentage of satellite-based IDFs lie within the confidence interval of NOAA Atlas 14.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2013

PROTOPLANETARY DISK MASSES FROM STARS TO BROWN DWARFS

We present SCUBA-2 850 {mu}m observations of seven very low mass stars (VLMS) and brown dwarfs (BDs). Three are in Taurus and four in the TW Hydrae Association (TWA), and all are classical T Tauri (cTT) analogs. We detect two of the three Taurus disks (one only marginally), but none of the TWA ones. For standard grains in cTT disks, our 3{sigma} limits correspond to a dust mass of 1.2 M{sub Circled-Plus} in Taurus and a mere 0.2 M{sub Circled-Plus} in the TWA (3-10 Multiplication-Sign deeper than previous work). We combine our data with other submillimeter/millimeter (sub-mm/mm) surveys of Taurus, {rho} Oph, and the TWA to investigate the trends in disk mass and grain growth during the cTT phase. Assuming a gas-to-dust mass ratio of 100:1 and fiducial surface density and temperature profiles guided by current data, we find the following. (1) The minimum disk outer radius required to explain the upper envelope of sub-mm/mm fluxes is {approx}100 AU for intermediate-mass stars, solar types, and VLMS, and {approx}20 AU for BDs. (2) While the upper envelope of apparent disk masses increases with M{sub *} from BDs to VLMS to solar-type stars, no such increase is observed from solar-type to intermediate-mass stars. We propose this is due to enhanced photoevaporation around intermediate stellar masses. (3) Many of the disks around Taurus and {rho} Oph intermediate-mass and solar-type stars evince an opacity index of {beta} {approx} 0-1, indicating significant grain growth. Of the only four VLMS/BDs in these regions with multi-wavelength measurements, three are consistent with considerable grain growth, though optically thick disks are not ruled out. (4) For the TWA VLMS (TWA 30A and B), combining our 850 {mu}m fluxes with the known accretion rates and ages suggests substantial grain growth by 10 Myr, comparable to that in the previously studied TWA cTTs Hen 3-600A and TW Hya. The degree of grain growth in the TWA BDs (2M1207A and SSPM1102) remains largely unknown. (5) A Bayesian analysis shows that the apparent disk-to-stellar mass ratio has a roughly constant mean of log{sub 10}[M{sub disk}/M{sub *}] Almost-Equal-To -2.4 all the way from intermediate-mass stars to VLMS/BDs, supporting previous qualitative suggestions that the ratio is {approx}1% throughout the stellar/BD domain. (6) Similar analysis shows that the disk mass in close solar-type Taurus binaries (sep <100 AU) is significantly lower than in singles (by a factor of 10), while that in wide solar-type Taurus binaries ({>=}100 AU) is closer to that in singles (lower by a factor of three). (7) We discuss the implications of these results for planet formation around VLMS/BDs, and for the observed dependence of accretion rate on stellar mass.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

The MCDB System for Management and Analysis of Petabyte-Scale Uncertain Data

Analysts working with very large data sets often use statistical models to “guess” at unknown, inaccurate, or missing information associated with the data. For example, a distant object viewed through an optical lens will have its position slightly shifted by imperfections in the lens. Thus, rather than considering the object’s observed position to be absolutely correct, it makes sense to take into account the lens’s imperfections to obtain a probabilistic guess as to the object’s true position. For another example, it might be important to associate some sort of error distribution with each of the individual sensors in an array of magnetometers. This error distribution may be complex and include spatially-driven covariances, because errors in nearby sensors are likely to be correlated (caused, for example, by the presence of some nearby, fixed metal object). This project is concerned with the design and implementation of a prototype data management system called the Monte Carlo Database System, or MCDB for short. MCDB allows an expert-level analyst or statistician to attach arbitrary stochastic models to very large data sets in order to “guess” the values for unknown or inaccurate data, such as the actual position of the observed object in the lens example above. When the resulting data set is analyzed, the underlying stochastic models are used to generate hundreds or thousands of possible data set instances, and each of those possible instances is analyzed separately by MCDB. Thus, MCDB does not just give a single answer to the analysis, but it actually gives an empirical distribution of query results that embody the underlying uncertainty, and can in turn be analyzed using standard statistical techniques. The stochastic models in MCDB are implemented as user-defined, external C++ libraries called Variable Generation functions (VG functions for short). Because the VG function interface is exceedingly general, it allows MCDB to be used in a very wide variety of application domains, in conjunction with virtually any stochastic model. Ultimately, we ended up calling the system that we developed over the course of the project SimSQL, to denote the fact that it allowed for stochastic simulation of data. The software developed during the lifetime of the project is available for download from the SimSQL project website: http://cmj4.web.rice.edu/SimSQL/SimSQL.html.

42 ENGINEERING↗

DOE OSTI2020

Thermal transport across nanoscale damage profile in sapphire irradiated by swift heavy ions

In this work, we studied the degradation of thermal conductivity in single crystal sapphire (α-Al 2 O 3 ) irradiated by 167 MeV Xe swift heavy ions (SHIs) over the multiple fluences in the range of 10 12 –10 14 ions/cm 2 . Thermal conductivity was measured primarily in the cross-plane direction using a noncontact ultrafast optical pump-probe technique called picosecond time domain thermoreflectance (TDTR). Multiple samples with variable ion fluences allowed us to probe distinct regions resulting from different regimes of microstructure evolution caused by electronic energy loss. By tuning the penetration depth of the thermal waves using different modulation frequencies, two regions with distinct conductivities were identified and the values of which were found to be consistent with phonon-mediated thermal transport models while the microstructure was confirmed by electron microscopy characterization. These damaged regions were determined to be a several micrometer thick ion track region and several tens of nanometer-thick amorphous layer present only above 5.0 × 10 13 ions/cm 2 . These results demonstrate the applicability of TDTR to resolve thermal transport behavior in SHI irradiated oxides having nonhomogeneous damage profile on a nanometer scale. The presented approach facilitates future studies aiming at resolving the impact of distinct damage resulting from electronic and nuclear stopping regimes under irradiation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRSDate not supplied

Using XR for Improving Scientific Discovery With Numerical Weather Models

Earth science (ES) digital twins will help us understand the complex interactions and interrelationships that make up our Earth system and the impacts of earth science phenomena on it. Our work addresses two underdeveloped areas in current ES digital twin work: improving the understanding and interaction with ES model outputs by using Virtual and Mixed Reality (XR) tools and improving the non-intuitive mapping of continuous ES natural phenomena to gridded reference frames in current numerical models. Traditionally, scientists working on ES view and analyze the results of calculated or measured observables with static 1-dimensional (1D), 2D or 3D plots displayed on flat computer screens or paper. Using such limited mediums, it can be very difficult to identify, track and understand the evolution of key features due to poor viewing angles and the nature of flat computer screens. In addition, numerical models, such as the NASA Goddard Earth Observing System (GEOS) ES model, are almost exclusively formulated, visualized and analyzed in an Eulerian reference frame with fixed grid points in space and time. However, ES phenomena such as convective clouds, hurricanes and wildfire smoke plumes are visualized and analyzed in a Lagrangian reference frame: therefore it is often difficult and unnatural to understand these phenomena in relation to each other, visualized either in an Eulerian or Lagrangian context. In 3D visualizations, data generally takes one of three forms: gridded (e.g., voxelized) data, where space is divided into regions; point clouds, where data is represented as a set of points; and meshes, where objects are rendered as surfaces composed of small polygons (usually triangles). A gridded, Eulerian reference frame has been the default representation for the 2D visual analysis of atmospheric data in part because the numerical methods used to generate atmospheric model data in the first place use a gridded approach, with equations defining the relationships between the physical variables in each of a grid's cells across successive timesteps. In our work, we are particularly interested in data from GEOS. Another reason why gridded representations tend to be used for visualizing data from such models is because trajectories are difficult to interpret from representations on 2D surfaces, due to line-of-sight ambiguity. Instead of a fixed grid from GEOS, we embed a trajectory model to simulate particles' movement throughout a GEOS run. We then ingest these particle trajectories as animated point clouds with a NASA open source XR toolkit, the Mixed Reality Exploration Toolkit (MRET), and merge GEOS data with ES phenomena data onto one combined visualization that the user can intuitively interact with. Efficient rendering of arbitrarily large point clouds is an ongoing challenge being addressed by the computer science community, with the GPU-based optimizations and efficient GPU memory utilization a common theme of recent advances, especially for XR, where sustained high frame rate is mandatory to save the user from suffering due to simulation sickness. In this work, we describe and evaluate our progress in choosing and implementing appropriate methods for rendering arbitrarily large point clouds within MRET for XR. While tracking the XR headset enables the immersion of a user within a 3D scene of a data visualization, tracking of XR handheld controllers or user’s hands enables us to implement intuitive user interactions with the visualized datasets. Conventional tools require a user working with an ES visualization to conduct many interactions to commit their intended selections or manipulations with a visualized dataset; for example to specify a set of points in 3D space. Doing so in a 2D flat screen interface has traditionally required specifying a set of points in three distinct 2D coordinate systems (XY, XZ, and YZ), which is cumbersome. In other scientific domains, it has been shown that specifying or selecting a location or volume in XR using handheld controllers or tracked hands allows for greater speed and accuracy. We anticipate the same will hold true for atmospheric data, and we will share initial results of measuring the utility of such an interface. Notably, as the data being visualized is generated by GEOS as a prediction based on initial conditions, an intended application of our tool is to serve as part of an iterative feedback loop. Through XR, a scientist will review and manipulate a GEOS model run, modifying the conditions as needed to do subsequent runs of GEOS. Thereby, XR-based improvements to speed and accuracy of 3D tagging of points minimizes the effort required by both the scientist and the computer cluster conducting the necessary calculations.

Thomas Grubb↗

DOE OSTI2013

Downscaling Solar Power Output to 4-Seconds for Use in Integration Studies: Preprint

High penetration renewable integration studies require solar power data with high spatial and temporal accuracy to quantify the impact of high frequency solar power ramps on the operation of the system. Our previous work concentrated on downscaling solar power from one hour to one minute by simulation. This method used clearness classifications to categorize temporal and spatial variability, and iterative methods to simulate intra-hour clearness variability. We determined that solar power ramp correlations between sites decrease with distance and the duration of the ramp, starting at around 0.6 for 30-minute ramps between sites that are less than 20 km apart. The sub-hour irradiance algorithm we developed has a noise floor that causes the correlations to approach ~0.005. Below one minute, the majority of the correlations of solar power ramps between sites less than 20 km apart are zero, and thus a new method to simulate intra-minute variability is needed. These intra-minute solar power ramps can be simulated using several methods, three of which we evaluate: a cubic spline fit to the one-minute solar power data; projection of the power spectral density toward the higher frequency domain; and average high frequency power spectral density from measured data. Each of these methods either under- or over-estimates the variability of intra-minute solar power ramps. We show that an optimized weighted linear sum of methods, dependent on the classification of temporal variability of the segment of one-minute solar power data, yields time series and ramp distributions similar to measured high-resolution solar irradiance data.

14 SOLAR ENERGY↗

DOE OSTI2023

HFTS-1 (Final Report)

Despite the long history of hydraulic fracturing, the optimal number of fracturing stages during multi-stage fracture stimulation in horizontal wells can be a complex subject. In addition to the increased expense of multistage fracturing in horizontal wells, increasing the number of fracturing stages does not always correlate with a commensurate increase in production. The problem is the application of a uniform fracture stimulation design to all stages with no consideration for geological variations along the wellbore. The result is an inefficient use and costly waste of energy and water. Optimization of the fracturing process requires an understanding of the cause-and-effect relationship between fracturing parameters and local geological properties at a given location along the wellbore. Realizing that the generalized rock mechanics theories and hypotheses are not truly applicable to fractured and laminated shales, quantifiable impacts of a shale’s geomechanical and depositional features are a prerequisite for design and implementation of optimized hydraulic fractures. The overarching goal of this project is to understand and define the relationships of shale geology and fracture dynamics using detailed field data that includes coring of the fracture domain. Finally, resource recovery from shale formations is estimated to be less than fifteen percent. Research proposed in this project will establish the foundations for investigating enhanced recovery techniques for increased resource recovery in existing fracture treated wells. Natural gas as an EOR fluid is in the initial stages of being used broadly in the oil and gas industry and as such, many aspects of the process need to be researched and addressed prior to widespread acceptance. In addition, re-fracturing techniques, particularly using cemented liners also need to be further investigated. This project explores several of these EOR techniques in greater detail, not just at the theoretical level, but through actual field test sites. This report delves into the research program, identified tasks, several studies, observations, results, and conclusions under the Hydraulic Fracturing Test Site (HFTS) project. This test project is subdivided into three distinct phases and each phase is discussed in detail under separate sections. The original test site for Phase 1 as well as the Phase 2 site are in Midland Basin in Texas and operated by Laredo Petroleum. In the first two phases, the program targeted the Wolfcamp formation in Midland Basin. Phase 3 of the study targeted the EagleFord play in South Texas with Devon Energy as the operator. Significant conclusions from this study include: 1) lateral and vertical fracture growth mapping as well as proppant transport behavior, 2) Use of field gas and its applicability for huff n puff EOR injection in the Midland Basin, 3) Development and validation of proppant log, both at core scale as well as in-fill well implementation by analyzing drilling mud return samples, 4) Applicability of liner refracturing including best practices for enhanced unit level productivity. This report will detail all these studies as well as technical results and conclusions from said studies. These include general subsurface characterization work, completion designs as implemented, studies looking at fracture geometry, high resolution microseismic study, advanced crosswell strain mapping studies including drainage characterization, core characterization results, proppant analysis (core and observation wells), geochemistry, etc. to name a few.

02 PETROLEUM↗

NASA NTRS2001

Accurate Time-Dependent Traveling-Wave Tube Model Developed for Computational Bit-Error-Rate Testing

The phenomenal growth of the satellite communications industry has created a large demand for traveling-wave tubes (TWT's) operating with unprecedented specifications requiring the design and production of many novel devices in record time. To achieve this, the TWT industry heavily relies on computational modeling. However, the TWT industry's computational modeling capabilities need to be improved because there are often discrepancies between measured TWT data and that predicted by conventional two-dimensional helical TWT interaction codes. This limits the analysis and design of novel devices or TWT's with parameters differing from what is conventionally manufactured. In addition, the inaccuracy of current computational tools limits achievable TWT performance because optimized designs require highly accurate models. To address these concerns, a fully three-dimensional, time-dependent, helical TWT interaction model was developed using the electromagnetic particle-in-cell code MAFIA (Solution of MAxwell's equations by the Finite-Integration-Algorithm). The model includes a short section of helical slow-wave circuit with excitation fed by radiofrequency input/output couplers, and an electron beam contained by periodic permanent magnet focusing. A cutaway view of several turns of the three-dimensional helical slow-wave circuit with input/output couplers is shown. This has been shown to be more accurate than conventionally used two-dimensional models. The growth of the communications industry has also imposed a demand for increased data rates for the transmission of large volumes of data. To achieve increased data rates, complex modulation and multiple access techniques are employed requiring minimum distortion of the signal as it is passed through the TWT. Thus, intersymbol interference (ISI) becomes a major consideration, as well as suspected causes such as reflections within the TWT. To experimentally investigate effects of the physical TWT on ISI would be prohibitively expensive, as it would require manufacturing numerous amplifiers, in addition to acquiring the required digital hardware. As an alternative, the time-domain TWT interaction model developed here provides the capability to establish a computational test bench where ISI or bit error rate can be simulated as a function of TWT operating parameters and component geometries. Intermodulation products, harmonic generation, and backward waves can also be monitored with the model for similar correlations. The advancements in computational capabilities and corresponding potential improvements in TWT performance may prove to be the enabling technologies for realizing unprecedented data rates for near real time transmission of the increasingly larger volumes of data demanded by planned commercial and Government satellite communications applications. This work is in support of the Cross Enterprise Technology Development Program in Headquarters' Advanced Technology & Mission Studies Division and the Air Force Office of Scientific Research Small Business Technology Transfer programs.

Kory, Carol L.↗

DOE OSTI2017

ISAVS: Interactive Scalable Analysis and Visualization System

Modern science is inundated with ever increasing data sizes as computational capabilities and image acquisition techniques continue to improve. For example, simulations are tackling ever larger domains with higher fidelity, and high-throughput microscopy techniques generate larger data that are fundamental to gather biologically and medically relevant insights. As the image sizes exceed memory, and even sometimes local disk space, each step in a scientific workflow is impacted. Current software solutions enable data exploration with limited interactivity for visualization and analytic tasks. Furthermore analysis on HPC systems often require complex hand-written parallel implementations of algorithms that suffer from poor portability and maintainability. We present a software infrastructure that simplifies end-to-end visualization and analysis of massive data. First, a hierarchical streaming data access layer enables interactive exploration of remote data, with fast data fetching to test analytics on subsets of the data. Second, a library simplifies the process of developing new analytics algorithms, allowing users to rapidly prototype new approaches and deploy them in an HPC setting. Third, a scalable runtime system automates mapping analysis algorithms to whatever computational hardware is available, reducing the complexity of developing scaling algorithms. Here, we demonstrate the usability and performance of our system using a use case from neuroscience: filtering, registration, and visualization of tera-scale microscopy data. We evaluate the performance of our system using a leadership-class supercomputer, Shaheen II.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2018

Automatic Imagery Data Analysis for Proactive Computer-Based Workflow Management during Nuclear Power Plant Outages

This report is being submitted for the task “Final Report” of DOE NEUP Project 15-8121 “Automatic Imagery Data Analysis for Proactive Computer-Based Workflow Management during Nuclear Power Plant Outages.” Typical nuclear power plant (NPP) outages always involve thousands of maintenance and refueling activities and a large number of workers in limited workspaces, while having tight schedules and zero-tolerance for accidents. During an outage, thousands of workers will be working in various workspaces across the NPP. High outage costs and expensive delays (approximately 1.5 million dollars of loss per day of delay) in NPP maintenance demand tight outage schedules. In packed workspaces, an automatic system that monitors human behaviors in real-time and provides insights about current and pending schedule deviations from the plan is critical for ensuring: 1) effective collaboration among workers and worker teams from different trades; 2) less waste of time and resources due to the lack of situational awareness; and 3) proactive outage project control. The overall goal of this project is to test the hypothesis that real-time imagery-based object tracking and spatial analysis, as well as human behavior modeling of outage participants, will significantly improve the efficiency of outage control while lowering the rates of accidents and incidents. Three objectives of this project are: 1) Establish real-time object tracking and spatiotemporal analysis methods that automatically assess the productivity of field activities and detect anomalous spatiotemporal relationships among activities that cause inefficiencies and risks; 2) Establish realtime human tracking and human factor modeling methods for automatically diagnosing unexpected actions of and interactions between outage participants, those which cause inefficient collaborations between Advanced Outage Control Center (AOCC), satellite outage centers, NPP workers, and maintenance service providers; and 3) Test the proposed automated object tracking, human behavior modeling, and spatiotemporal analysis methods in outage control case studies in order to characterize the effectiveness of automated imagery-data-driven methods in proactively improving the efficiency and safety of workflows in outage coordination and risk management. Recent studies of detailed human behavior monitoring on construction sites have examined the potential of applying advanced computer vision algorithms in detecting and tracking anomalous workers (i.e., workers who do not wear hard hats or safety vests) for ensuring job site safety. Some studies of human factor studies revealed the importance of modeling detailed interactions between individuals within and across teams in better understanding the impact of human in proactive project control. Other studies in the construction domain have developed computational simulation frameworks that formalize detailed spatiotemporal interactions between tasks in order to simulate the impacts of individual tasks on the performance of workflows. While integrated analyses that combine human factor assessment, as well as image processing and simulation, are in demand for effective decision-making, limited studies have examined the potential of such integrated analyses in NPP outage control. This research project has examined an automatic outage monitoring and control system that integrates human factor analysis, computer vision techniques, and simulation methods in order to enable engineers to better understand the interactions between humans, resources, and workflows during outage processes. The project aims at providing NPP maintenance agencies insights into more efficient use of limited resources in extending the life of a nuclear plant, as well as reducing waste while ensuring sufficient generation of electricity. This study is significant for the safety of nuclear plants, sustainable electricity generation for livable communities, and cost savings for maintaining electricity infrastructures in the United States.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1992

Research in computational fluid dynamics and analysis of algorithms

Recently, higher-order compact schemes have seen increasing use in the DNS (Direct Numerical Simulations) of the Navier-Stokes equations. Although they do not have the spatial resolution of spectral methods, they offer significant increases in accuracy over conventional second order methods. They can be used on any smooth grid, and do not have an overly restrictive CFL dependence as compared with the O(N(exp -2)) CFL dependence observed in Chebyshev spectral methods on finite domains. In addition, they are generally more robust and less costly than spectral methods. The issue of the relative cost of higher-order schemes (accuracy weighted against physical and numerical cost) is a far more complex issue, depending ultimately on what features of the solution are sought and how accurately they must be resolved. In any event, the further development of the underlying stability theory of these schemes is important. The approach of devising suitable boundary clusters and then testing them with various stability techniques (such as finding the norm) is entirely the wrong approach when dealing with high-order methods. Very seldom are high-order boundary closures stable, making them difficult to isolate. An alternative approach is to begin with a norm which satisfies all the stability criteria for the hyperbolic system, and look for the boundary closure forms which will match the norm exactly. This method was used recently by Strand to isolate stable boundary closure schemes for the explicit central fourth- and sixth-order schemes. The norm used was an energy norm mimicking the norm for the differential equations. Further research should be devoted to BC for high order schemes in order to make sure that the results obtained are reliable. The compact fourth order and sixth order finite difference scheme had been incorporated into a code to simulate flow past circular cylinders. This code will serve as a verification of the full spectral codes. A detailed stability analysis by Carpenter (from the fluid Mechanics Division) and Gottlieb gave analytic conditions for stability as well as asymptotic stability. This had been incorporated in the code in form of stable boundary conditions. Effects of the cylinder rotations had been studied. The results differ from the known theoretical results. We are in the middle of analyzing the results. A detailed analysis of the effects of the heating of the cylinder on the shedding frequency had been studied using the above schemes. It has been found that the shedding frequency decreases when the wire was heated. Experimental work is being carried out to affirm this result.

Gottlieb, David↗

DOE OSTI2019

Experimental Study of the 3 He and Neutron Spin Structure at Low Q 2 and The Upgrade of the JLab Polarized 3 He Target

In the scattering process off a nuclear or nucleon target, the Gerasimov-Drell-Hearn (GDH) sum rule for real photons (Q 2 =0 where Q 2 &equiv; -q 2 with q the photon's 4-momentum) relates static properties of the target particle's ground state to dynamic properties of all its excited states. On the other side of the Q 2 spectrum, the Bjorken sum rule holds in the Bjorken limit Q 2 → ∞. Bjorken sum rule relates the final structure functions of the proton and neutron to the nucleon axial coupling constant in weak decay. These two sum rules belong to domains where calculations are achievable but use different degrees of freedom: hadronic degrees of freedom at low Q 2 versus partonic degrees of freedom at intermediate Q 2 . Meanwhile, different methods have been used to connect the two sum rules at finite Q 2 values: Chiral Perturbation Theory is used to expand the GDH sum rule while Operator Product Expansion is used to expand the Bjorken sum rule. In recent decades, improvements in polarized beam and polarized target techniques have made it possible to test theoretical predictions in the intermediate Q 2 region. During the Jefferson Lab (JLab) Hall A E97110 experiment, a precise measurement of polarized cross sections was performed at 0.02 < Q 2 < 0.3 GeV 2 using a polarized 3 He target as an effective polarized neutron target. The measured data allowed us to test predictions of Chiral Perturbation Theory at very low Q 2 . Furthermore, an extrapolation to the real photon point Q 2 =0 tests the GDH sum rule on the neutron. In order to reach the small angles necessary for the low Q 2 range, a new septum magnet was installed in Hall A for this experiment. Unfortunately, the magnet was mis-wired during initial running. There were therefore two periods for this experiment: the first period had the defective magnet due to mis-wiring; while in the second period, the magnet had been fixed and was working properly. The cross sections and the asymmetries for the first period must be extracted using a difficult and unusual method employing focal plane variables in the spectrometer. In the work described in this thesis, the target-plane and focal-plane method will be established first and confirmed by elastic scattering cross section measurements from carbon foil and 3 He targets. Cross sections and asymmetries of inelastic scattering were then extracted using the same method, which constitute the main physics results of this thesis. Preliminary results on ITT , the integral of the polarized cross section, were also extracted and are presented in this thesis. In addition to the data analysis of E97110, this thesis includes instrumentation work performed on the polarized 3 He target at JLab. For the 12 GeV program of JLab, the polarized 3 He target is being upgraded to satisfy new experimental requirements. The new target is a convection-based cell with two transfer tubes rather than a single transfer tube as in the former 6 GeV design. The first stage of the target upgrade aims to produce a 40 cm long, 10 amg target that can withstand 30 μA of electron beam current with an in-beam polarization of 55%, doubling the figure-of-merit from the 6 GeV time. Several studies were conducted at the Test Lab of JLab on the new convection cell. Work and improvement done on the Pulse Nuclear Magnetic Resonance (NMR) polarimetry and its test results will be presented. As this thesis is written, the target is being installed in Hall C of JLab in preparation for a 6-month long running from November 2019 to May 2020.

Ton, Nguyen↗

DOE OSTI2022

Protection and Restoration Solutions to Reliable and Resilient Integration of Grid-connected PV Installations and Distributed Energy Resources: Design, Testbed, Proof of Work and Impact Studies (Final Report)

To gain a better understanding of the complex transients in a utility grid with a large number of solar PV installations coupled by PWM inverters during the protection and restoration period of the power grid, in this project a kW-level experimental system with dominating inverter-based resources and a hardware-in-loop simulation system with transmission and distribution models, as well as several SEL protective relays have been built and used. The major research findings of the project are summarized below in two perspectives: Experimental research: A self-organized ultra-high frequency solitary waveform is discovered and demonstrated in the hardware experimental system. Despite the familiarities, such a solitary waveform is distinct from any harmonics, transient, or resonance waves in that, it is 1) non-dispersive over time and space, 2) not associated with any active source or the linear superposition of sources, 3) half-cycle asymmetric, 4) not responding to filters or change of system characteristic resonance frequency, 5) ubiquitous as it occurs simultaneously everywhere in the system from DC supply, power lines, and the utility grid, and 6) explosive through tripping the protection or damaging susceptible devices or circuits. Analytical study: The nature of such a new waveform and an analytical explanation of its formation are studied based on related physics and non-linear science principles, with the following highlights: 1) such a new waveform follows the solutions to the Non-linear Schrödinger equation, so the classic linear perturbation theory is unable to explain or predict such a unique waveform as confirmed with the research team of RTDS; 2) the critical condition of occurrence of such a waveform is derived, which indicates that the breakings of solitary wave is a system synchronous issue between the utility grid in the 60Hz phasor domain and duty-ratio modulation of DC sources in the inverter switching frequency domain; 3) such a waveform carries energy mass so it could be detrimental; 4) such a waveform is deceiving as it is not readily detectable in the energy propagation direction by the primary protection equipment, while it is detrimental in the perpendicular direction of energy propagation, i.e. voltage direction. Thus, it has more likely challenges to the second primary equipment and devices, particularly at the weakest link and point such as aging insulation and inappropriate setting of susceptible devices. Therefore, such a waveform can be easily ignored in the aftermath investigation. The intellectual merits: The experimental discovery reveals certain unfamiliar transient phenomena that could challenge the integration of large-scale grid-connected solar PV installation during the group ride-through period of solar PV installations, commissioning of large-solar farms, or dramatic change of solar radiation conditions. Particularly, the research findings suggest that the occurrence of the unfamiliar transient phenomena is uniquely associated with the inverter-based solar PV installation, which is less likely to happen for rotary energy systems. The physics and non-linear science-based research work laid down an analytical path to the challenging transient stability problems including those that have been observed currently and future calls. Both experimental and analytical research results explain the limitation of many current research effects including some DoE research undertakings as well as possible solutions. The broader impacts: The research outcome of the project advances the understanding of the possible transient problems for the integration of inverter-based solar PV installation. It also highlights the theoretical and technical barriers for applying the conventional theory and technology to design and implement countermeasures against their adverse impacts of the large-scale solar PV installation and operation on grid reliability and security. More broadly, the research outcomes have shed some light on the open, fundamental challenges in the integration of large-scale solar photovoltaic energy into current and next-generation power grids across the country, and worldwide. The advanced physics and non-linear science-based analysis open up a path to harmonization of the renewables for societal energy needs via building a resilient and sustainable electric power infrastructure.

14 SOLAR ENERGY↗

DOE OSTI2014

Aerosol indirect effect on the grid-scale clouds in the two-way coupled WRF–CMAQ: model description, development, evaluation and regional analysis

This study implemented first, second and glaciation aerosol indirect effects (AIE) on resolved clouds in the two-way coupled Weather Research and Forecasting Community Multiscale Air Quality (WRF–CMAQ) modeling system by including parameterizations for both cloud drop and ice number concentrations on the basis of CMAQ-predicted aerosol distributions and WRF meteorological conditions. The performance of the newly developed WRF–CMAQ model, with alternate Community Atmospheric Model (CAM) and Rapid Radiative Transfer Model for GCMs (RRTMG) radiation schemes, was evaluated with observations from the Clouds and the See http://ceres.larc.nasa.gov/ . Earth's Radiant Energy System (CERES) satellite and surface monitoring networks (AQS, IMPROVE, CASTNET, STN, and PRISM) over the continental US (CONUS) (12 km resolution) and eastern Texas (4 km resolution) during August and September of 2006. The results at the Air Quality System (AQS) surface sites show that in August, the normalized mean bias (NMB) values for PM 2.5 over the eastern US (EUS) and the western US (WUS) are 5.3% (-0.1%) and 0.4% (-5.2%) for WRF–CMAQ/CAM (WRF–CMAQ/RRTMG), respectively. The evaluation of PM 2.5 chemical composition reveals that in August, WRF–CMAQ/CAM (WRF–CMAQ/RRTMG) consistently underestimated the observed SO 4 2- by -23.0% (-27.7%), -12.5% (-18.9%) and -7.9% (-14.8%) over the EUS at the Clean Air Status Trends Network (CASTNET), Interagency Monitoring of Protected Visual Environments (IMPROVE) and Speciated Trends Network (STN) sites, respectively. Both configurations (WRF–CMAQ/CAM, WRF–CMAQ/RRTMG) overestimated the observed mean organic carbon (OC), elemental carbon (EC) and and total carbon (TC) concentrations over the EUS in August at the IMPROVE sites. Both configurations generally underestimated the cloud field (shortwave cloud forcing, SWCF) over the CONUS in August due to the fact that the AIE on the subgrid convective clouds was not considered when the model simulations were run at the 12 km resolution. This is in agreement with the fact that both configurations captured SWCF and longwave cloud forcing (LWCF) very well for the 4 km simulation over eastern Texas, when all clouds were resolved by the finer resolution domain. The simulations of WRF–CMAQ/CAM and WRF–CMAQ/RRTMG show dramatic improvements for SWCF, LWCF, cloud optical depth (COD), cloud fractions and precipitation over the ocean relative to those of WRF default cases in August. The model performance in September is similar to that in August, except for a greater overestimation of PM 2.5 due to the overestimations of SO 4 2- , NH 4 + , NO 3 - , and TC over the EUS, less underestimation of clouds (SWCF) over the land areas due to the lower SWCF values, and fewer convective clouds in September. Finally, this work shows that inclusion of indirect aerosol effect treatments in WRF–CMAQ represents a significant advancement and milestone in air quality modeling and the development of integrated emissions control strategies for air quality management and climate change mitigation.

indirect aerosol forcing↗

DOE OSTI2017

Dynamic mapping of conical intersection seams: A general method for incorporating the geometric phase in adiabatic dynamics in polyatomic systems

The incorporation of the geometric phase in single-state adiabatic dynamics near a conical intersection (CI) seam has so far been restricted to molecular systems with high symmetry or simple model Hamiltonians. This is due to the fact that the ab initio determined derivative coupling (DC) in a multi-dimensional space is not curl-free, thus making its line integral path dependent. In a recent work [C. L. Malbon et al. , J. Chem. Phys. 145, 234111 (2016)], we proposed a new and general approach based on an ab initio determined diabatic representation consisting of only two electronic states, in which the DC is completely removable, so that its line integral is path independent in the simply connected domains that exclude the CI seam. Then with the CIs included, the line integral of the single-valued DC can be used to construct the complex geometry-dependent phase needed to exactly eliminate the double-valued character of the real-valued adiabatic electronic wavefunction. This geometry-dependent phase gives rise to a vector potential which, when included in the adiabatic representation, rigorously accounts for the geometric phase in a system with an arbitrary locus of the CI seam and an arbitrary number of internal coordinates. Here, we demonstrate this approach in a three-dimensional treatment of the tunneling facilitated dissociation of the S 1 state of phenol, which is affected by a C s symmetry allowed but otherwise accidental seam of CI. Here, since the space is three-dimensional rather than two-dimensional, the seam is a curve rather than a point. The nodal structure of the ground state vibronic wavefunction is shown to map out the seam of CI.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

AnisONet: A deep neural operator-based anisotropic permeability upscaler from pore to Darcy scale

Directional permeability variations, which govern directional fluid flow in porous media with anisotropy, are important to accurately predict flow behavior, reactive transport, and fluid–solid interactions for various processes such as enhanced geothermal systems, energy storage devices, and biological systems. However, the intricate architecture of porous media makes it difficult to predict directional permeabilities. In this work, we present a novel machine learning (ML) framework, AnisONet, built upon an integration of a convolutional neural network, Swin transformer, and the deep operator network architecture, designed to predict anisotropic permeability and upscale predictions to larger spatial domains. First, AnisONet was evaluated with three classes of two-dimensional (2D) porous media, including synthetic circular and elliptical grains and natural sandstone grains from micro-computed tomography images. A lattice Boltzmann model (LBM) was used to calculate directional permeabilities at every 10° angle, producing 19 data points per image of porous media. AnisONet is then trained to predict permeability as a function of rotation angle. AnisONet showed strong predictive capability of directional permeability. Second, we tested our model for five upscaling cases with a large image size in the finite-element method (FEM) for 2D Darcy flow with various permeability tensor construction methods. Overall, upscaled permeability tensors in FEM simulations produce a reasonably good match with LBM results, highlighting the importance of selecting appropriate tensor formation strategies for accurate permeability upscaling. AnisONet, as a directional permeability estimator, could be further developed for more complex geometries, with the potential to develop a foundational ML model for various applications in porous media.

42 ENGINEERING↗

DOE OSTI2015

Predictive Maturity of Multi-Scale Simulation Models for Fuel Performance

The project proposed to provide a Predictive Maturity Framework with its companion metrics that (1) introduce a formalized, quantitative means to communicate information between interested parties, (2) provide scientifically dependable means to claim completion of Validation and Uncertainty Quantification (VU) activities, and (3) guide the decision makers in the allocation of Nuclear Energy’s resources for code development and physical experiments. The project team proposed to develop this framework based on two complimentary criteria: (1) the extent of experimental evidence available for the calibration of simulation models and (2) the sophistication of the physics incorporated in simulation models. The proposed framework is capable of quantifying the interaction between the required number of physical experiments and degree of physics sophistication. The project team has developed this framework and implemented it with a multi-scale model for simulating creep of a core reactor cladding. The multi-scale model is composed of the viscoplastic self-consistent (VPSC) code at the meso-scale, which represents the visco-plastic behavior and changing properties of a highly anisotropic material and a Finite Element (FE) code at the macro-scale to represent the elastic behavior and apply the loading. The framework developed takes advantage of the transparency provided by partitioned analysis, where independent constituent codes are coupled in an iterative manner. This transparency allows model developers to better understand and remedy the source of biases and uncertainties, whether they stem from the constituents or the coupling interface by exploiting separate-effect experiments conducted within the constituent domain and integral-effect experiments conducted within the full-system domain. The project team has implemented this procedure with the multi- scale VPSC-FE model and demonstrated its ability to improve the predictive capability of the model. Within this framework, the project team has focused on optimizing resource allocation for improving numerical models through further code development and experimentation. Related to further code development, we have developed a code prioritization index (CPI) for coupled numerical models. CPI is implemented to effectively improve the predictive capability of the coupled model by increasing the sophistication of constituent codes. In relation to designing new experiments, we investigated the information gained by the addition of each new experiment used for calibration and bias correction of a simulation model. Additionally, the variability of ‘information gain’ through the design domain has been investigated in order to identify the experiment settings where maximum information gain occurs and thus guide the experimenters in the selection of the experiment settings. This idea was extended to evaluate the information gain from each experiment can be improved by intelligently selecting the experiments, leading to the development of the Batch Sequential Design (BSD) technique. Additionally, we evaluated the importance of sufficiently exploring the domain of applicability in experiment-based validation of high-consequence modeling and simulation by developing a new metric to quantify coverage. This metric has also been incorporated into the design of new experiments. Finally, we have proposed a data-aware calibration approach for the calibration of numerical models. This new method considers the complexity of a numerical model (the number of parameters to be calibrated, parameter uncertainty, and form of the model) and seeks to identify the number of experiments necessary to calibrate the model based on the level of sophistication of the physics. The final component in the project team’s work to improve model calibration and validation methods is the incorporation of robustness to non-probabilistic uncertainty in the input parameters. This is an improvement to model validation and uncertainty quantification stemming beyond the originally proposed scope of the project. We have introduced a new metric for incorporating the concept of robustness into experiment-based validation of numerical models. This project has accounted for the graduation of two Ph.D. students (Kendra Van Buren and Josh Hegenderfer) and two M.S. students (Matthew Egeberg and Parker Shields). One of the doctoral students is now working in the nuclear engineering field and the other one is a post-doctoral fellow at the Los Alamos National Laboratory. Additionally, two more Ph.D. students (Garrison Stevens and Tunc Kulaksiz) who are working towards graduation have been supported by this project.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗