Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spatial Temporal Model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 883 records · Page 49

Rift Valley Fever Risk Map Model and Seroprevalence in Selected Wild Ungulates and Camels from Kenya

Since the first isolation of Rift Valley fever virus (RVFV) in the 1930s, there have been multiple epizootics and epidemics in animals and humans in sub-Saharan Africa. Prospective climate-based models have recently been developed that flag areas at risk of RVFV transmission in endemic regions based on key environmental indicators that precede Rift Valley fever (RVF) epizootics and epidemics. Although the timing and locations of human case data from the 2006-2007 RVF outbreak in Kenya have been compared to risk zones flagged by the model, seroprevalence of RVF antibodies in wildlife has not yet been analyzed in light of temporal and spatial predictions of RVF activity. Primarily wild ungulate serum samples from periods before, during, and after the 2006-2007 RVF epizootic were analyzed for the presence of RVFV IgM and/or IgG antibody. Results show an increase in RVF seropositivity from samples collected in 2007 (31.8%), compared to antibody prevalence observed from 2000-2006 (3.3%). After the epizootic, average RVF seropositivity diminished to 5% in samples collected from 2008-2009. Overlaying maps of modeled RVF risk assessments with sampling locations indicated positive RVF serology in several species of wild ungulate in or near areas flagged as being at risk for RVF. Our results establish the need to continue and expand sero-surveillance of wildlife species Kenya and elsewhere in the Horn of Africa to further calibrate and improve the RVF risk model, and better understand the dynamics of RVFV transmission.

Seroprevalence↗

Mesoscale Convective Systems During SCSMEX: Simulations with a Regional Climate Model and a Cloud-Resolving Model

The South China Sea Monsoon Experiment (SCSMEX) was conducted in May-June 1998. One of its major objectives is to better understand the key physical processes for the onset and evolution of the summer monsoon over Southeast Asia and southern China (Lau et al. 2000). Multiple observation platforms (e.g., soundings, Doppler radar, ships, wind seafarers, radiometers, etc.) during SCSMEX provided a first attempt at investigating the detailed characteristics of convection and circulation changes, associated with monsoons over the South China Sea region. SCSMEX also provided precipitation derived from atmospheric budgets (Johnson and Ciesielski 2002) and comparison to those obtained from the Tropical Rainfall Measuring Mission (TRMM). In this paper, a regional climate model and a cloud-resolving model are used to perform multi-day integrations to understand the precipitation processes associated with the summer monsoon over Southeast Asia and southern China. The regional climate model is used to understand the soil - precipitation interaction and feedback associated with a flood event that occurred in and around China's Atlantic River during SCSMEX. Sensitivity tests on various land surface models, cumulus parameterization schemes (CASE), sea surface temperature (SST) variations and midlatitude influences are also performed to understand the processes associated with the onset of the monsoon over the S. China Sea during SCSMEX. Cloud-resolving models (CRMs) use more sophisticated and physically realistic parameterizations of cloud microphysical processes with very fine spatial and temporal resolution. One of the major characteristics of CRMs is an explicit interaction between clouds, radiation and the land/ocean surface. It is for this reason that GEWEX (Global Energy and Water Cycle Experiment) has formed the GCSS (GEWEX Cloud System Study) expressly for the purpose of improving the representation of the moist processes in large-scale models using CRMs. The Goddard Cumulus Ensemble (GCE) model is a CRM and is used to simulate convective systems associated with the onset of the South China Sea monsoon in 1998. The BRUCE model includes the same land surface model, cloud physics, and radiation scheme used in the regional climate model. A comparison between the results from the GCE model and regional climate model is performed.

Tao, W. K.↗

Multiphysics for nuclear energy applications using a cohesive computational framework

With the recent development of advanced numerical algorithms, software design, and low-cost high-performance computer hardware, reliance on coupled multiphysics to predict the behavior of complex physical systems is beginning to become standard practice. This is especially true in nuclear energy applications where strong nonlinear interdependencies exist between reactor physics, radiation transport, multi-scale nuclear fuels performance, thermal fluids, etc. Resolving these nonlinear dependencies requires choices in multiphysics software approaches. Two main multiphysics modeling and simulation approaches have emerged. The first is based upon "code coupling" where disparate physics codes of different software design, code languages, and spatial and temporal integration schemes are coupled together with relatively complex data passing interfaces. The second multiphysics software approach is to employ a "cohesive" framework where all physics applications are developed with a common software design, i.e., data structures, syntax, input format, integrated spatial and temporal discretization schemes, etc. In this paper we present the Multiphysics Object-Oriented Simulation Environment (MOOSE) development and runtime framework and describe the framework's cohesive modeling and simulation multiphysics approach. Then, a "cohesive-like" extension of the MOOSE framework is presented where MOOSE-based physics software applications are efficiently coupled to non-MOOSE (external) physics codes to form multiphysics applications using MOOSE's unique interface capabilities. Finally, several examples of MOOSE's cohesive and cohesive-like multiphysics applications will be demonstrated. These multiphysics demonstrations will incorporate both MOOSE-based applications and external codes, including Nek5000, RELAP-7, TRACE, BISON, and Pronghorn.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

OCO-2 Satellite-imposed Constraints on Terrestrial Biospheric CO2 Flux over South Asia

The spatiotemporal variability of terrestrial biospheric carbon dioxide (CO2) flux over South Asia is poorly understood. The inverse model “top-down” CO2 flux estimates which rely on atmospheric CO2 observational data are impeded by sparse in situ measurements over this region. The Orbiting Carbon Observatory 2 (OCO-2) provides much needed global retrievals of column-average CO2 dry-air mole fraction (XCO2) with the finest spatial and temporal resolution and highest sensitivity to surface fluxes available to-date. This study conducted global inverse model simulations as part of the second OCO-2 Multi-model Intercomparison Project assimilating version 9 of the OCO-2 XCO2 retrievals in land nadir (LN) and land glint (LG) observing modes and global in situ (IS) measurements. The four-dimensional variational assimilation system with the GEOS-Chem global chemical transport model was used to estimate CO2 fluxes from 2015 to 2018. We then assessed the spatiotemporal variability of optimized CO2 Net ecosystem Exchange (NEE) fluxes aggregated over the South Asia region. The most robust result found by assimilating OCO-2 observations was the constraints imposed on the seasonal cycle of terrestrial biospheric fluxes over South Asia. The seasonality of South Asian NEE estimated by assimilating OCO-2 or global IS data showed a larger seasonal cycle compared to the current understanding of NEE in this region (represented by the prior NEE used in the model). The satellite-in situ joint inversion (IS + LN + LG) led to land flux seasonal amplitude (absolute magnitude of the difference between peak and trough of monthly mean values over a year) of 4.0 PgC yr-1, compared to the prior model (1.5 PgC yr-1). Moreover, OCO-2 data imposed a phase shift in the seasonal cycle, resulting in a large CO2 source in April and a large uptake in September. Most of the features in the seasonal cycle imposed by OCO-2 data are in agreement with previous “top-down” studies assimilating regional aircraft observations. Nevertheless, we recommend a closer examination of the novel features highlighted by the OCO-2 column data over South Asia in future studies by conducting regional inverse modelling with denser in situ and vertically-resolved regional observations along with satellite data.

OCO-2↗

Quantifying the Known Unknown: Including Marine Sources of Greenhouse Gases in Climate Modeling

Researchers have recently estimated that Arctic submarine permafrost currently traps 60 billion tons of methane and contains 560 billion tons of organic carbon in seafloor sediments and soil, a giant pool of carbon with potentially large feedbacks on the climate system. Unlike terrestrial permafrost, the submarine permafrost system has remained a “known unknown” because of the difficulty in acquiring samples and measurements. Consequently, this potentially large carbon stock never yet considered in global climate models or policy discussions, represents a real wildcard in our understanding of Earth’s climate. This report summarizes our group’s effort at developing a numerical modeling framework designed to produce a first-of-its-kind estimate of Arctic methane gas releases from the marine sediments to the water column, and potentially to the atmosphere, where positive climate feedback may occur. Newly developed modeling capability supported by the Laboratory Directed Research and Development (LDRD) program at Sandia National Laboratories now gives us the ability to probabilistically map gas distribution and quantity in the seabed by using a hybrid approach of geospatial machine learning, and predictive numerical thermodynamic ensemble modeling. The novelty in this approach is its ability to produce maps of useful data in regions that are only sparsely sampled, a common challenge in the Arctic, and a major obstacle to progress in the past. By applying this model to the circum-Arctic continental shelves and integrating the flux of free gas from in situ methanogenesis and dissociating gas hydrates from the sediment column under climate forcing, we can provide the most reliable estimate of a spatially and temporally varying source term for greenhouse gas flux that can be used by global oceanographic circulation and Earth system models (such as DOE’s E3SM). The result will allow us to finally tackle the wildcard of the submarine permafrost carbon system, and better inform us about the severity of future national security threats that sustained climate change poses.

54 ENVIRONMENTAL SCIENCES↗

Atomization data requirements for rocket combustor modeling

The complex computer codes, which model liquid rocket combustors, require information about the distribution and atomization of these liquid reactants. The available information is, in general, of questionable validity and applicability. Authors and users of combustion codes are often unaware of, or underestimate the importance of, these deficiencies in atomization data. These deficiencies and their importance are examined. Results of analyses performed with a state-of-the-art rocket combustion code are presented which demonstrate the important effects of such atomization information as initial droplet sizes and size distribution on vaporization rate and losses. Also, the questionable aspects and inapplicability of the available atomization data are discussed. One important and often neglected or misunderstood aspect of atomization data is the differences between spatial (concentration) and flux (often called temporal) droplet size distributions. These are described, and a computer model constructed to assess the difference between concentration and flux droplet size distributions is described and results presented. Experimental data are also given to demonstrate this difference. Finally, experimental results are presented that demonstrate the very great, and often neglected effect, of the local gas velocity field on atomization.

Ferrenberg, A. J.↗

NMPC for Mode-Switching Operation of Reversible Solid Oxide Cell Systems

Solid oxide cells (SOCs) are a promising dual-mode technology that generates hydrogen through high-temperature water electrolysis and generates power through a fuel cell reaction that consumes hydrogen. Reversible operation of SOCs requires a transition between these two modes for hydrogen production setpoints as the demand and price of electricity fluctuate. Moreover, a well-functioning control system is important to avoid cell degradation during mode-switching operation. In this work, we apply nonlinear model predictive control (NMPC) to an SOC module and supporting equipment and compare NMPC performance to classical proportional integral (PI) control strategies, while ramping between the modes of hydrogen and power production. While both control methods provide similar performance in many metrics, NMPC significantly reduces cell thermal gradients and curvatures (mixed spatial-temporal partial derivatives) during mode switching. A dynamic process flowsheet of the reversible SOC system was developed in the open-source, equation-based IDAES modeling framework. Our IDAES dynamic simulation results show that NMPC can ramp the SOC system between hydrogen and power production targets within short mode-switching times. Moreover, NMPC can comply with operating limits in the SOC system more effectively than PI, and only NMPC can directly enforce user-specified limits for mixed spatial-temporal partial derivatives of temperature. This allows for management of the trade-off be-tween operating efficiency and cell degradation, which is dependent on these temperature curvatures.

Li, Mingrui↗

Lightning Strike Distance Distribution Beyond a Preexisting Lightning Area

The 45th Weather Squadron (45 WS) asked the Applied Meteorology Unit (AMU) to review the 30-year-old, lightning stand-off distances of 5 nautical miles (nmi) for applicability to today's operations. This was based on the realization that previous lightning strike distance studies did not match how 45 WS issues lightning warnings (Roeder, 2008). The previous lightning distance studies were from the point of origin of the lightning or the average starting location that would tend to be in the core of the thunderstorm. However, the 45 WS issues lightning warnings based on the edge of a preexisting lightning area. Before beginning the AMU project, it took several years to develop a method to calculate a distance distribution beyond a preexisting area (Roeder, 2015). The AMU pulled Lightning Detection and Ranging (LDAR) sensor data from 1/1/2013 to 12/31/2013. This dataset consisted of 37 million individual source data points from the LDAR sensors. Only sources within 50 km north, south, east or west of the LDAR grid center were included in the dataset. This limited the use of LDAR data to that with the greatest accuracy of source detection and increased data processing. Points were grouped into flashes based on spatial and temporal criteria. Based on the sensitivity analysis the AMU performed on the flash clustering algorithm, a time value of 0.3 seconds was found to model flashes adequately. Distance parameters were tested from 1,500 to 7,500 meters (m) in 500 m increments. Distance parameters of both 3,000 m and 4,000 m produced results in the plotting tool that were most representative of the physical behavior of lightning. Thus statistics were gathered for the most representative of these spatial and temporal criteria on the flash size and the polygon expansion distance in order to find the correct data distributions. The best fit curves for the LDAR polygon expansion frequency vs. distances for both the 3 kilometer (km) and 4 km distance threshold values were exponential decay functions and had R2 values of > 0.998, indicating good model fits. The equations of the best fit curves were then used to calculate a desired safety radius of 4 nmi for either 3 km or 4 km distance threshold criteria. The AMU analysis concludes the safe reduction of the 5 nmi lightning warning circles to 4 nmi should improve the operational impact by 36% if based on distance from the center of the property area being protected. If based on the edge of the property being protected, then the reduction is 4.5 nmi to 4 nmi and the operational impact is decreased by 21%. For the 6 nmi lightning warning circles, the recommended 4 nmi stand-off distance will result in a safe reduction of operational impact of 31% if based on the center of the area being protected, or 16% if based on the edge of the property.

Flash clustering algorithm↗

Solar Activity Modeling: From Subgranular Dynamical Scales to the Solar Cycles

Dynamical effects of solar magnetoconvection span a wide range spatial and temporal scales that extends from the interior to the corona and from fast turbulent motions to the global-Sun magnetic activity. To study the solar activity on short temporal scales (from minutes to hours), we use 3D radiative MHD simulations that allow us to investigate complex turbulent interactions that drive various phenomena, such as plasma eruptions, spontaneous formation of magnetic structures, funnel-like structures and magnetic loops in the corona, and others. In particular, we focus on multi-scale processes of energy exchange across the different layers, which contribute to the corona heating and eruptive dynamics, as well as interlinks between different layers of the solar interior and atmosphere. For modeling the global-scale activity we use the data assimilation approach that has demonstrated great potential for building reliable long-term forecasts of solar activity. In particular, it has been shown that the Ensemble Kalman Filter (EnKF) method applied to the Parker-Kleeorin-Ruzmakin dynamo model is capable of predicting solar activity up to one sunspot cycle ahead in time, as well as estimating the properties of the next cycle a few years before it begins. In this presentation, using the available magnetogram data, we discuss development of the methodology and forecast quality criteria (including forecast uncertainties and sources of errors). We demonstrate the influence of observational limitation on the prediction accuracy. We present the EnKF predictions of the upcoming Solar Cycle 25 based on both the sunspot number series and observed magnetic fields, and discuss the uncertainties and potential of the data assimilation approach for modeling and forecasting the solar activity.

Kitiashvili, I. N.↗

Machine Learned Empirical Numerical Integrator from Simulated Data

Recently, a number of state-of-the-art surrogate machine learning (ML) models have been designed for global weather and climate prediction, which have been trained using reanalysis data products. Reanalysis data products are constructed using numerical model simulations that combine numerical integration of partial differential equations and parameterization schemes. These products are typically only archived and made available using coarsened spatial and temporal resolutions. This study explores the impact of the numerical generation methods used to produce the training datasets and the temporal resolution of those datasets on machine learning surrogate models. Using the nonlinear vector autoregression (NVAR) machine as an explainable ML technique, simple dynamical systems are emulated with ML models trained on data produced by three classical numerical integration schemes. NVAR is validated as a skillful ML method, capable of producing accurate predictions and, more importantly, reconstructing both the underlying dynamics and the numerical integration scheme used to generate the training data. However, the machine fails to generalize predictions on unseen test data generated by different numerical integration schemes, despite the underlying dynamical system being the same. This result provides a word of caution for the growing field of machine learning emulation of weather and climate dynamics. Furthermore, we illustrate using NVAR that training on temporally coarsened data may increase the required complexity of ML models and potentially introduce new numerical challenges. Finally, we discover that empirical integration schemes with arbitrary time-stepping sizes can be constructed directly from the data, which implies a potential for the development of empirical numerical integration schemes.

54 ENVIRONMENTAL SCIENCES↗

Modeling detonation with CartaBlanca simulations

The accurate modeling of high explosive (HE) detonation and the con- sequent large solid deformation, failure, plastic flow, porosity growth, and shock wave propagation is important because simulations can capture spatial and temporal features that experimental diagnostics cannot capture. However, the simulation of the explosive event poses challenges to a computational scientist. These include the accurate modeling of ductile damage, crack formation, plastic deformation, as well as physical and nu- merical instabilities. The material response can be history-dependent and subject to large material deformation. Our research simulates the impact of a high explosive (Detasheet) onto a tantulum metal plate. We performed the simulations using CartaBlanca at different mesh resolutions. We decided the study would be impactful if we perform the simulations with the Material Point Method. Differences were observed at the different mesh resolutions in velocity and nodal stress magnitude, so increased mesh resolutions may be required. In addition, the Discontinuous Galerkin method would be needed to account for the large gas expansion.

97 MATHEMATICS AND COMPUTING↗

Understanding the Impacts of Soil Moisture Initial Conditions on NWP in the Context of Land-Atmosphere Coupling

The role of soil moisture in NWP has gained more attention in recent years, as studies have demonstrated impacts of land surface states on ambient weather from diurnal to seasonal scales. However, soil moisture initialization approaches in coupled models remain quite diverse in terms of their complexity and observational roots, while assessment using bulk forecast statistics can be simplistic and misleading. In this study, a suite of soil moisture initialization approaches is used to generate short-term coupled forecasts over the U.S. Southern Great Plains using NASA’s Land Information System (LIS) and NASA Unified WRF (NU-WRF) modeling systems. This includes a wide range of currently used initialization approaches, including soil moisture derived from ‘off-the-shelf’ products such as atmospheric models and land data assimilation systems, high-resolution land surface model spinups, and satellite-based soil moisture products from SMAP. Results indicate that the spread across initialization approaches can be quite large in terms of soil moisture conditions and spatial resolution, and that SMAP performs well in terms of heterogeneity and temporal dynamics when compared against high resolution land surface model and in-situ soil moisture estimates. Case studies are analyzed using the local land-atmosphere coupling (LoCo) framework that relies on integrated assessment of soil moisture, surface flux, boundary layer, and ambient weather, with results highlighting the critical role of inherent model background biases. In addition, simultaneous assessment of land vs. atmospheric initial conditions in an integrated, process-level fashion can help address the question of whether improvements in traditional NWP verification statistics are achieved for the right reasons.

Southern Great Plains↗

Temporal Large-Eddy Simulation

In 1999, Stolz and Adams unveiled a subgrid-scale model for LES based upon approximately inverting (defiltering) the spatial grid-filter operator and termed .the approximate deconvolution model (ADM). Subsequently, the utility and accuracy of the ADM were demonstrated in a posteriori analyses of flows as diverse as incompressible plane-channel flow and supersonic compression-ramp flow. In a prelude to the current paper, a parameterized temporal ADM (TADM) was developed and demonstrated in both a priori and a posteriori analyses for forced, viscous Burger's flow. The development of a time-filtered variant of the ADM was motivated-primarily by the desire for a unifying theoretical and computational context to encompass direct numerical simulation (DNS), large-eddy simulation (LES), and Reynolds averaged Navier-Stokes simulation (RANS). The resultant methodology was termed temporal LES (TLES). To permit exploration of the parameter space, however, previous analyses of the TADM were restricted to Burger's flow, and it has remained to demonstrate the TADM and TLES methodology for three-dimensional flow. For several reasons, plane-channel flow presents an ideal test case for the TADM. Among these reasons, channel flow is anisotropic, yet it lends itself to highly efficient and accurate spectral numerical methods. Moreover, channel-flow has been investigated extensively by DNS, and a highly accurate data base of Moser et.al. exists. In the present paper, we develop a fully anisotropic TADM model and demonstrate its utility in simulating incompressible plane-channel flow at nominal values of Re(sub tau) = 180 and Re(sub tau) = 590 by the TLES method. The TADM model is shown to perform nearly as well as the ADM at equivalent resolution, thereby establishing TLES as a viable alternative to LES. Moreover, as the current model is suboptimal is some respects, there is considerable room to improve TLES.

Pruett, C. D.↗

Comparison of Local Scale Measured and Modeled Brightness Temperatures and Snow Parameters from the CLPX 2003 by Means of a Dense Medium Radiative Transfer Theory Model

Microwave remote sensing offers distinct advantages for observing the cryosphere. Solar illumination is not required, and spatial and temporal coverage are excellent from polar-orbiting satellites. Passive microwave measurements are sensitive to the two most useful physical quantities for many hydrological applications: physical temperature and water content/state. Sensitivity to the latter is a direct result of the microwave sensitivity to the dielectric properties of natural media, including snow, ice, soil (frozen or thawed), and vegetation. These considerations are factors motivating the development of future cryospheric satellite remote sensing missions, continuing and improving on a 26-year microwave measurement legacy. Perhaps the biggest issues regarding the use of such satellite measurements involve how to relate parameter values at spatial scales as small as a hectare to observations with sensor footprints that may be up to 25 x 25 km. The NASA Cold-land Processes Field Experiment (CLPX) generated a dataset designed to enhance understanding of such scaling issues. CLPX observations were made in February (dry snow) and March (wet snow), 2003 in Colorado, USA, at scales ranging from plot scale to 25 x 25 km satellite footprints. Of interest here are passive microwave observations from ground-based, airborne, and satellite sensors, as well as meteorological and snowpack measurements that will enable studies of the effects of spatial heterogeneity of surface conditions on the observations. Prior to performing such scaling studies, an evaluation of snowpack forward modelling at the plot scale (least heterogeneous scale) is in order. This is the focus of this paper. Many forward models of snow signatures (brightness temperatures) have been developed over the years. It is now recognized that a dense medium radiative transfer (DMRT) treatment represents a high degree of physical fidelity for snow modeling, yet dense medium models are particularly sensitive to snowpack structural parameters such as grain size, density, and depth---parameters that may vary substantially within a snowpack. Microwave radiometric data and snow pit measurements collected at the Local-Scale Observation Site (LSOS) during the third Intensive Observation Period (IOP3) of the CLPX have been used to test the capabilities of a DMRT model using the Quasi Crystalline Approximation with Coherent Potential (QCA-CP). The radiometric measurements were made by the University of Tokyo s Ground Based Microwave Radiometer-7 (GBMR-7) system. We evaluate the degree to which a DMRT-based model can accurately reproduce the GBMR-7 brightness temperatures at different frequencies and incidence angles.

Tedescol, Marco↗

Extreme metrics from large ensembles: investigating the effects of ensemble size on their estimates

Abstract. We consider the problem of estimating the ensemble sizes required to characterize the forced component and the internal variability of a number of extreme metrics. While we exploit existing large ensembles, our perspective is that of a modeling center wanting to estimate a priori such sizes on the basis of an existing small ensemble (we assume the availability of only five members here). We therefore ask if such a small-size ensemble is sufficient to estimate accurately the population variance (i.e., the ensemble internal variability) and then apply a well-established formula that quantifies the expected error in the estimation of the population mean (i.e., the forced component) as a function of the sample size n, here taken to mean the ensemble size. We find that indeed we can anticipate errors in the estimation of the forced component for temperature and precipitation extremes as a function of n by plugging into the formula an estimate of the population variance derived on the basis of five members. For a range of spatial and temporal scales, forcing levels (we use simulations under Representative Concentration Pathway 8.5) and two models considered here as our proof of concept, it appears that an ensemble size of 20 or 25 members can provide estimates of the forced component for the extreme metrics considered that remain within small absolute and percentage errors. Additional members beyond 20 or 25 add only marginal precision to the estimate, and this remains true when statistical inference through extreme value analysis is used. We then ask about the ensemble size required to estimate the ensemble variance (a measure of internal variability) along the length of the simulation and – importantly – about the ensemble size required to detect significant changes in such variance along the simulation with increased external forcings. Using the F test, we find that estimates on the basis of only 5 or 10 ensemble members accurately represent the full ensemble variance even when the analysis is conducted at the grid-point scale. The detection of changes in the variance when comparing different times along the simulation, especially for the precipitation-based metrics, requires larger sizes but not larger than 15 or 20 members. While we recognize that there will always exist applications and metric definitions requiring larger statistical power and therefore ensemble sizes, our results suggest that for a wide range of analysis targets and scales an effective estimate of both forced component and internal variability can be achieved with sizes below 30 members. This invites consideration of the possibility of exploring additional sources of uncertainty, such as physics parameter settings, when designing ensemble simulations.

54 ENVIRONMENTAL SCIENCES↗

Development of Digital Twin Predictive Model for PWR Components: Updates on Multi Times Series Temperature Prediction Using Recurrent Neural Network, DMW Fatigue Tests, System Level Thermal-Mechanical-Stress Analysis

The long-term operation (LTO) of nuclear power plant (NPP) beyond their original design life of 40 years, can lead to more material damage associated with cyclic fatigue under thermal-mechanical loading cycles and associated long-term exposure of reactor material to the deleterious reactor-coolant environments. However, under this LTO condition the reactor components can still safely operate but may require more frequent Nondestructive Evaluation (NDE) of reactor components. Frequent NDE requirement may lead to frequent shutdown of the NPP. This in turn can lead to power outage and additional NDE-inspection-cost related economic loss. The economic loss can be minimized by reducing uncertainty in life estimation of safety-critical pressure boundary components and by implementing more digital approach such as by using upcoming digital-twin (DT) technology for predicting the structural states (e.g., time and location dependent inside/outside thickness temperature, stress, strain, plastic deformation, etc.) and associated fatigue life of a component in real time. Towards this goal Argonne National Laboratory (ANL) with the sponsorship of DOE Light Water Reactor Sustainability (LWRS) program is working on the development of a DT framework that can be used for real time environmental fatigue prediction of reactor components. The DT framework is based on limited experiment-data, Artificial-intelligence (AI) – Machine-Learning (ML) - Deep-Learning (DL) based techniques and Multiphysics-computational-mechanics such as finite element (FE) based modeling tools. Towards this overall goal, following are some of the major contributions made during the FY21: 1) Multiple 82/182 dissimilar metal weld (DMW) specimens (both solid-weld and joint-weld representing the actual reactor multi-metal nozzles) were fatigue tested. The resulting fatigue lives were compared to the NUREG-6909 based best-fit and design fatigue curves. Additionally, the results of 52/152 DMW fatigue specimens (which were recently tested at Republic of Korea under the sponsorship of International Nuclear Energy Research Initiative - INERI program) were compared to the NUREG-6909 based best-fit and design fatigue curves. From the comparison of 82/182 and 52/152 DMW test data with NUREG-6909 best-fit curve, most of the reported test data fall way away from the NUREG-6909 suggested best-fit or mean curve. The NUREG-6909 suggested best-fit curve is the best-fit curve of austenitic stainless steel and due to lack of enough data on Nickel-based welds, this is currently being used for predicting the life of Nickel-alloy-based welded components. However, the above observation may require higher scaling factor (e.g., ASME suggested factor of 20 on cycles rather than the current NUREG-6909 suggested factor of 12 on cycles) for scaling the austenitic-stainless-steel best-fit-curve for estimating the design or safe-life of a welded component. Accordingly, for example, if a DMW component experience a strain amplitude of 0.6% the PWR-water life of the component would be 52 cycles instead of 85 cycles. However, more DMW tests are required to further ascertain the above-mentioned observations. 2) A system level CAD and finite element model were developed which consists of reactor pressure vessel (RPV), part of steam generator (SG), part of pressurizer (PRZ), hot leg (HL), and surge line (SL). This is with detailed nozzle geometry and thermal-mechanical material properties of different metals to simulate realistic thermal-mechanical stress under connected system global thermal-mechanical boundary conditions. 3) Different system level heat transfer analyses were performed with estimation of relevant heat transfer coefficients. The resulting data were used in subsequent system level thermal-mechanical stress analysis and for generating spatial-temporal training and validation data for a system level digital-twin based temperature predictor. Transient heat transfer analyses were performed considering thermal boundary condition under design-basis (DB) loading and EDF (Électricité de France) data-based grid-load-following (EDF-GLF) loading cycles. 4) System level thermal-mechanical stress analysis was performed for identifying damage-prone hotspots and for future extension of the model for cyclic state prediction. From the system-level model simulation under DB loading cycle it is found that HL and the SL nozzle that connect to the HL can experience significant stress and strain and could be one of the weakest links in the overall reactor coolant system (RCS). 5) An AI/ML based DT model was developed for multi-time-series temperature prediction at any inside/outside thickness locations of PWR pressure boundary components. This is by using Recurrent-neural-network (RNN) and keras machine learning libraries. The RNN model was validated against two laboratory test-based data sets with one obtained through ANL’s in-air fatigue test system and other through PWR-water test loop. The experimentally validated DT model further validated against FE model results to predict thermal scarification related spatialtemporal temperatures at random locations of a component. The well validated DT model was then used for demonstrating spatial-temporal temperature prediction under 100+ years of reactor operation subjected to combined DB, EDF-GLF and randomized grid-load-following (RANDOMGLF) loading Cycles. The expert-elicitation DT model framework was developed assuming field/input/process measurements can be available from a few existing plant sensors and can readily be used by the NPP operators. The above temperature prediction model will feed to the next-step stress analysis model based on which the life of a component can be predicted in realtime, which is one of our future works.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The aluminium-26 distribution in a cosmological simulation of a Milky Way-type Galaxy

Context . The 1.8 MeV γ -rays corresponding to the decay of the radioactive isotope 26 Al (with a half-life of 0.72 Myr ) have been observed by the SPI detector on the INTEGRAL spacecraft and extensively used as a tracer of star formation and current nucleosynthetic activity in the Milky Way Galaxy. Further information is encoded in the observation related to the higher 26 Al content found in regions of the Galaxy with the highest line-of-sight (LoS) velocity relative to an observer located in the Solar System. However, this feature remains unexplained. Aims . We ran a cosmological “zoom-in” chemodynamical simulation of a Milky Way-type galaxy, including the production and decays of radioactive nuclei in a fully self-consistent way. We then analyzed the results to follow the evolution of 26 Al throughout the lifetime of the simulated galaxy to provide a new method for interpreting the 26 Al observations. Methods . We included the massive star sources of 26 Al in the Galaxy and its radioactive decay into a state-of-the-art galactic chemical evolution model, coupled with cosmological growth and hydrodynamics. This approach allowed us to follow the spatial and temporal evolution of the 26 Al content in the simulated galaxy. Results . Our results are in agreement with the observations with respect to the fact that gas particles in the simulation with relatively higher 26 Al content also have the highest LoS velocities. On the other hand, gas particles with relatively lower 26 Al content (i.e., not bright enough to be observed) generally display the lowest LoS velocities. However, this result is not conclusive because the overall rotational velocity of our simulated galaxy is higher than that observed for cold CO gas in the Milky Way Galaxy. Furthermore, we found no significant correlation between gas temperature, rotational velocity, and 26 Al content at any given radius. We also found the presence of transient 26 Al-rich spots at low LoS velocities and we show that one such spot had been captured by the INTEGRAL/SPI data. Based on our model, we present a prediction for the detection of 1.8 MeV γ -rays by the future COSI mission. We find that according to our model, the new instrument will be able to observe similar 26 Al-emission patterns to those seen by INTEGRAL/SPI.

Wehmeyer, B. (ORCID:0009000254149292)↗

The Effect of Spatially Varying Collision Frequency on the Development of the Rayleigh–Taylor Instability

Abstract The Rayleigh–Taylor (RT) instability is ubiquitously observed, yet has traditionally been studied using ideal fluid models. Collisionality can vary strongly across the fluid interface, and previous work demonstrates the necessity of kinetic models to completely capture dynamics in certain collisional regimes. Where previous kinetic simulations used spatially and temporally constant collision frequency, this work presents five-dimensional (two spatial, three velocity dimensions) continuum-kinetic simulations of the RT instability using a more realistic spatially varying collision frequency. Three cases of collisional variation are explored for two Atwood numbers: low to intermediate, intermediate to high, and low to high. The low-to-intermediate case exhibits no RT instability growth, while the intermediate-to-high case is similar to a fluid-limit kinetic case with interface widening biased toward the lower-collisionality region. A novel contribution of this work is the low-to-high collisionality case that shows significantly altered instability growth through an upward movement of the interface and damped spike growth due to increased free-streaming particle diffusion in the lower region. Contributions to the energy flux from the non-Maxwellian portions of the distribution function are not accessible to fluid models and are greatest in magnitude in the spike and regions of low collisionality. Increasing the Atwood number results in greater RT instability growth and reduced upward interface movement. Deviation of the distribution function from Maxwellian is inversely proportional to collision frequency and concentrated around the fluid interface. The linear phase of RT instability growth is well described by theoretical linear growth rates accounting for viscosity and diffusion.

79 ASTRONOMY AND ASTROPHYSICS↗