Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Observing System Simulation”

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 1,405 records · Page 78

A regional nuclear conflict would compromise global food security

A limited nuclear war between India and Pakistan could ignite fires large enough to emit more than 5 Tg of soot into the stratosphere. Climate model simulations have shown severe resulting climate perturbations with declines in global mean temperature by 1.8 °C and precipitation by 8%, for at least 5 y. Here we evaluate impacts for the global food system. Six harmonized state-of-the-art crop models show that global caloric production from maize, wheat, rice, and soybean falls by 13 (±1)%, 11 (±8)%, 3 (±5)%, and 17 (±2)% over 5 y. Total single-year losses of 12 (±4)% quadruple the largest observed historical anomaly and exceed impacts caused by historic droughts and volcanic eruptions. Colder temperatures drive losses more than changes in precipitation and solar radiation, leading to strongest impacts in temperate regions poleward of 30°N, including the United States, Europe, and China for 10 to 15 y. Integrated food trade network analyses show that domestic reserves and global trade can largely buffer the production anomaly in the first year. Persistent multiyear losses, however, would constrain domestic food availability and propagate to the Global South, especially to food-insecure countries. By year 5, maize and wheat availability would decrease by 13% globally and by more than 20% in 71 countries with a cumulative population of 1.3 billion people. In view of increasing instability in South Asia, this study shows that a regional conflict using <1% of the worldwide nuclear arsenal could have adverse consequences for global food security unmatched in modern history.

54 ENVIRONMENTAL SCIENCES↗

Preliminary Evaluation of Joint Electricity-Hydrogen Concept of Operations

An initial thermal power dispatch (TPD) concept of operations was evaluated that couples a nuclear power plant to a nearby hydrogen production plant. GSE Systems’ generic pressurized water reactor full-scope simulator was modified with a TPD model comprised of a thermal power extraction and delivery system. A prototype human-system interface (HSI) was developed to interact with the TPD model and allow participants to execute the basic operating scenarios for normal operations. Four retired operators performed the evaluation, and due to COVID-19 travel restrictions, the original in-person experimental design was restructured to support a remote participator evaluation using a web meeting platform. Data from operator feedback, observations from the research team, and quantitative survey responses revealed that the initial TPD concept of operations is feasible. The operators were comfortable with the engineered system and HSI and could manage it without adverse impacts to reactor power, plant safety, or equipment. Findings are discussed in terms of both the TPD system design and HSI performance.

human factors↗

A Finite-Time Ensemble Method for Mixed Layer Model Comparison

Here, this work evaluates the fidelity of various upper-ocean turbulence parameterizations subject to realistic monsoon forcing and presents a finite-time ensemble vector (EV) method to better manage the design and numerical principles of these parameterizations. The EV method emphasizes the dynamics of a turbulence closure multimodel ensemble and is applied to evaluate 10 different ocean surface boundary layer (OSBL) parameterizations within a single-column (SC) model against two boundary layer large-eddy simulations (LES). Both LES include realistic surface forcing, but one includes wind-driven shear turbulence only, while the other includes additional Stokes forcing through the wave-average equations that generate Langmuir turbulence. The finite-time EV framework focuses on what constitutes the local behavior of the mixed layer dynamical system and isolates the forcing and ocean state conditions where turbulence parameterizations most disagree. Identifying disagreement provides the potential to evaluate SC models comparatively against the LES. Observations collected during the 2018 monsoon onset in the Bay of Bengal provide a case study to evaluate models under realistic and variable forcing conditions. The case study results highlight two regimes where models disagree 1) during wind-driven deepening of the mixed layer and 2) under strong diurnal forcing.

54 ENVIRONMENTAL SCIENCES↗

Evaluating and Improving Convective Parameterization for GCMs Using ARM Observations. Final report

This technical report summarizes the achievements during the project period funded by the Atmospheric System Research (ASR) program. The general goal of the project is to improve the representation of atmospheric deep convection in global climate models. By using DOE Atmospheric Radiation Measurement (ARM) program observations and a hierarchy of models from cloud-resolving to global climate models, we evaluated and improved many aspects of convection parameterization schemes, such as the trigger function for convection and closure condition to determine the amount of convection. We also incorporated the stochasticity of convection into a convection scheme. These improvements, when implemented into the National Center for Atmospheric Research (NCAR) Community Earth System Model (CESM) and the DOE Exascale Energy Earth System Model (E3SM), led to improved simulation of precipitation characteristics, including precipitation intensity, Intertropical Convergence Zone, and Madden-Julian oscillation. The project also improved the representation of microphysical processes in convective clouds. It enables the interaction of cloud microphysics with aerosols in weather and climate models.

54 ENVIRONMENTAL SCIENCES↗

Physical Characterization of Interstellar Comet 2I/2019 Q4 (Borisov)

We present a study of interstellar comet 2I/2019 Q4 (Borisov) using both preperihelion and postperihelion observations spanning late September 2019 through late January 2020. The intrinsic brightness of the comet was observed to continuously decline throughout the timespan, likely due to the decreasing effective scattering cross-section as a result of volatile sublimation with a slope of −0.43 ± 0.02 km{sup 2} d{sup −1}. We witnessed no significant change in the slightly reddish color of the comet, with mean values of 〈g−r〉=0.68±0.04, 〈r−i〉=0.23±0.03, and the normalized reflectivity gradient across the g and i bands (S{sup ′})-bar (g,i)=(10.6±1.4)% per 10{sup 3} Å, all unremarkable in the context of solar system comets. Using the available astrometric observations, we confidently detect the nongravitational acceleration of the comet following a shallow heliocentric distance dependence of r{sub H}{sup −1±1}. Accordingly, we estimate that the nucleus is most likely ≲0.4 km in radius, and that a fraction of ≳0.2% of the total mass of the nucleus has been eroded due to the sublimation activity since the earliest observation of the comet in 2018 December by the time of perihelion. Our morphology simulation suggests that the dust ejection speed increased from ∼4 m s{sup −1} in 2019 September to ∼7 m s{sup −1} around perihelion for the optically dominant dust grains with β∼0.01, and that the observable dust grains are no smaller than micron size.

79 ASTRONOMY AND ASTROPHYSICS↗

Computationally efficient subglacial drainage modelling using Gaussian process emulators: GlaDS-GP v1.0

Subglacial drainage models represent water flow at the ice–bed interface through coupled distributed and channelized systems to determine water pressure, discharge, and drainage system geometry. While they are used to understand processes such as the relationship between surface melt and ice flow, the number of uncertain model parameters and the computational cost of running models makes it difficult to adequately explore the high-dimensional parameter space and evaluate uncertainty in model predictions. Here, we develop Gaussian process (GP) emulators that make fast predictions with associated uncertainty of subglacial drainage model outputs. Using a truncated principal component (PC) basis representation, we construct a GP emulator for diurnally averaged subglacial water pressure. We also explore emulation of scalar variables describing drainage efficiency and configuration. We train the emulators using ensembles of up to 512 simulations varying eight parameters of the Glacier Drainage System (GlaDS) model on a synthetic domain intended to represent an ice-sheet margin. The emulators make predictions ∼ 1000 times faster than GlaDS simulations, with errors <3 % for the water pressure field and ∼ 5 %–9 % for drainage efficiency and configuration. We apply the emulators to explore the eight-dimensional parameter space by computing variance-based parameter sensitivity indices, finding that three parameters (ice flow coefficient, bed bump aspect ratio, and the subglacial cavity system conductivity) explain 90 % of the variance in modelled water pressure in response to parameter changes. The GP emulator approach described here is well suited to integrating observational data with models to make calibrated, credible predictions of subglacial drainage.

58 GEOSCIENCES↗

Modeling the effects of active wake mixing on wake behavior through large-scale coherent structures

The use of active wake mixing (AWM) to mitigate downstream turbine wakes has created new opportunities for reducing power losses in wind farms. However, many current analytical or semi-empirical wake models do not capture the flow instabilities that are excited through the blade pitch actuation. In this work, we develop a framework, which accounts for the impacts of the large-scale coherent structures and turbulence on the mean flow, for modeling AWM. The framework uses a triple-decomposition approach for the unsteady flow field and models the mean flow and fine-scale turbulence with a parabolized Reynolds-averaged Navier–Stokes (RANS) system. The wave components are modeled using a simplified spatial linear stability formulation that captures the growth and evolution of the coherent structures. Comparisons with high-fidelity large eddy simulations (LESs) of the turbine wakes showed that this framework was able to capture the additional wake mixing and faster wake recovery in the far-wake regions for both the pulse and helix AWM strategies with minimal computational expense. In the near-wake region, some differences are observed in both the RANS velocity profiles and initial growth of the large-scale structures, which may be due to some simplifying assumptions used in the model.

17 WIND ENERGY↗

Microwave-Assisted Reactive CO 2 Capture with the SrCO 3 -Graphite System

This study details the initial development of a microwave heat-driven reactive CO 2 capture (RCC) process using the SrO/SrCO 3 cycle and graphitic carbon to absorb ppm-levels of CO 2 from humidified room temperature air and selectively convert it into CO. By combining first-principles density functional theory (DFT) simulations with thermogravimetric analysis (TGA) and X-ray diffraction (XRD) verification, it is demonstrated that SrO spontaneously absorbs moisture to generate Sr(OH) 2 followed by Sr(OH) 2 *1H 2 O, which then can spontaneously react with atmospheric levels of CO 2 to form SrCO 3 . The resulting SrCO 3 can then react with solid carbon at elevated temperatures to selectively produce CO and regenerate SrO. Graphitic carbon is an excellent microwave absorber that can quickly generate temperatures approaching 1000 °C and then cool to room temperature within minutes, potentially allowing for integration with intermittent electricity. It is shown that using microwaves to selectively heat a mixture of graphitic carbon and SrCO 3 produced CO with 85 ± 3% selectivity and stable performance over ten cycles. The rapid heating for release and room temperature CO 2 uptake processes appeared to prevent performance loss from particle sintering typically observed with traditional thermal systems. In conclusion, these results demonstrate a new RCC approach that adds to the growing number of potential technologies available for converting CO 2 to useful chemicals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Small to Full Height Scale Comparisons of Cesium Ion Exchange Performance with Crystalline Silicotitanate

The U.S. Department of Energy’s (DOE) Hanford Site houses 56 million gallons of high-level radioactive waste generated from plutonium production from 1944 to 1988. The supernatant waste, currently stored in underground tanks, is intended to be vitrified following filtration and 137 Cs removal at the Hanford Waste Treatment and Immobilization Plant (WTP) Low-Activity Waste (LAW) Vitrification Facility. The WTP Pretreatment Facility will not be operational for several years. The Tank Side Cesium Removal (TSCR) system is a technology demonstration that will remove cesium from tank waste supernate to support directly feeding LAW to the vitrification facility. The 137 Cs removal is important to meet the WTP LAW contract specification and ultimately for creating a contact-handled waste form. The waste acceptance criteria (WAC) limit for the WTP LAW Facility is <3.18E-5 Ci 137 Cs per mole Na. The TSCR system is skid mounted and employs two key technologies: dead end filtration for solids removal and ion exchange (IX) for cesium removal. Filtration is necessary to protect the functionality of the IX columns. The IX process utilizes sodium form crystalline silicotitanate (CST) IX media produced by UOP Honeywell, LLC (Des Plaines, IL) under the product name IONSIV R9140-B, 18 x 50 mesh, in a lead-lag-polish column configuration. Each column contains a CST bed height of 92 inches and a volume of approximately 157 gallons of CST IX media. The full-scale TSCR operation will run at a residence time of 1.9 bed volumes (BV)/h, which results in a superficial velocity of 7.3 cm/min and a flowrate of ~5 gal/min. Column testing at small (2.5% of the full bed height), medium (12% of the full bed height), and full-height scales has been previously conducted to evaluate process variables and scale up performance of Cs exchange onto the CST. Cesium load performances from various sieve cuts at the small scale indicated that a 30-mesh sieve cut be tested to determine if it better reflects the 12% and full-height column performances at the small scale. Two process flowrates were tested in the small-column configuration with <30-mesh CST and simulant solutions. These tests were compared to full-height column tests at the same residence times to assess CST particle size effect on column scaling. Table S.1 summarizes the observed column performance determined for the two flowrates juxtaposed to the previous work with 5.6 M Na simulant at higher scales (used as benchmarks). The WAC breakthroughs between the small and full-height tests at 1.3 BV/h differed by ~59 BVs. The WAC breakthroughs for small, medium, and full-height tests at 1.8 BV/h were consistent at 240 BVs. The 50% Cs breakthroughs were nominally equivalent for all column tests. The common transition zones and onsets of Cs breakthrough at the 1.8 BV/h tests indicated that the Cs mass transfers were equivalent and thus the <30-mesh CST Cs load performance at the small scale successfully modeled that of the full-height system. It is recommended that the <30-mesh CST be used in subsequent 10-mL CST bed tests.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

PATOKA: Simulating Electromagnetic Observables of Black Hole Accretion

Abstract The Event Horizon Telescope (EHT) has released analyses of reconstructed images of horizon-scale millimeter emission near the supermassive black hole at the center of the M87 galaxy. Parts of the analyses made use of a large library of synthetic black hole images and spectra, which were produced using numerical general relativistic magnetohydrodynamics fluid simulations and polarized ray tracing. In this article, we describe the PATOKA pipeline, which was used to generate the Illinois contribution to the EHT simulation library. We begin by describing the relevant accretion systems and radiative processes. We then describe the details of the three numerical codes we use, iharm , ipole , and igrmonty , paying particular attention to differences between the current generation of the codes and the originally published versions. Finally, we provide a brief overview of simulated data as produced by PATOKA and conclude with a discussion of limitations and future directions.

79 ASTRONOMY AND ASTROPHYSICS↗

Reliable extrapolation of deep neural operators informed by physics or sparse observations

Deep neural operators can learn nonlinear mappings between infinite-dimensional function spaces via deep neural networks. As promising surrogate solvers of partial differential equations (PDEs) for real-time prediction, deep neural operators such as deep operator networks (DeepONets) provide a new simulation paradigm in science and engineering. Pure data-driven neural operators and deep learning models, in general, are usually limited to interpolation scenarios, where new predictions utilize inputs within the support of the training set. However, in the inference stage of real-world applications, the input may lie outside the support, i.e., extrapolation is required, which may result to large errors and unavoidable failure of deep learning models. Here, we address this challenge of extrapolation for deep neural operators. First, we systematically investigate the extrapolation behavior of DeepONets by quantifying the extrapolation complexity, via the 2-Wasserstein distance between two function spaces and propose a new strategy of bias–variance trade-off for extrapolation with respect to model capacity. Subsequently, we develop a complete workflow, including extrapolation determination, and we propose five reliable learning methods that guarantee a safe prediction under extrapolation by requiring additional information—the governing PDEs of the system or sparse new observations. The proposed methods are based on either fine-tuning a pre-trained DeepONet or multifidelity learning. We demonstrate the effectiveness of the proposed framework for various types of parametric PDEs. Furthermore, our systematic comparisons provide practical guidelines for selecting a proper extrapolation method depending on the available information, desired accuracy, and required inference speed.

42 ENGINEERING↗

Simulation of compound flooding in Japan using a nationwide model

For this study, a high-resolution, summit-to-sea unstructured-grid model was used to simulate two compound flooding events in different geomorphic settings in Japan: the July 2012 flood in Kumamoto city and July 2018 flood in Okayama City, both caused by a torrential rainfall during seasonal Meiyu front but otherwise exhibiting quite different compound flood characteristics. The model shows good performance in simulating flooding extent; e.g., the Hit Rate for the inundation event in 2018 exceeds 0.9. Sensitivity tests were conducted to determine the roles/significance of levee structures along major rivers in coastal regions. Our results indicate that riverine and oceanic factors as well as levees significantly contribute to the observed flooding extents, and we are able to quantify the contributions from each factor at different locations. The model has been implemented operationally as the first nationwide 3-day compound flooding forecast system for Japan.

54 ENVIRONMENTAL SCIENCES↗

A Bayesian Framework for Milling Stability Prediction and Reverse Parameter Identification

This paper describes a physics-guided Bayesian framework for identifying the milling stability boundary and system parameters through iterative testing. Prior uncertainties for the parameters are identified through physical simulation and literature reviews, without physical testing of the actual milling system. Those uncertainties are then propagated to the stability map using a physics-based stability model, which is used to suggest a test point. The uncertainties are updated based on the new information acquired from the cutting test to form a new probability distribution, called the posterior. Finally, the posterior are compared to measured values for the stability boundary and system parameters to evaluate the approach. Based on experimental observations, the advantages and disadvantages of using a physics-guided model are discussed.

Cornelius, Aaron↗

Super‐Droplet Method to Simulate Lagrangian Microphysics of Nuclear Fallout in a Homogeneous Cloud

Abstract Nuclear detonations produce hazardous local and global particles or fallout. Predicting fallout size, chemical components, and location is necessary to inform officials and determine immediate guidance for the public. However, existing nuclear detonation fallout models prescribe the particle size distributions based on limited observations. In this work, we apply the super‐droplet method, which is a numerical modeling technique developed for cloud microphysics, to simulate size distributions of particles in a mushroom cloud formed post‐detonation of a nuclear device. We model fallout formation and evolution with homogeneous nucleation and condensation of a single species and a Monte Carlo coagulation algorithm. We verify the numerical methods representing coagulation and condensation processes against analytical test problems. Additionally, we explore several scenarios for the integral system mass and yield in equivalent kilotons (kt) of TNT (trinitrotoluene). The fallout size distribution median diameter d pg follows a scaling law based on the integral system mass m v 0 kg and yield Y kt: nm. We test the effect of cloud turbulence, enhanced nucleation and growth, and vapor volatility with a sensitivity study. The range in median diameter predictions for simulations of historical tests performed over the Pacific encompass the measurements of particles sampled from the cloud caps. Predicted median particle size ranges up to 217, 123, 86, and 35 nm for historical tests with yields of 0.2, 0.7, 2, and 10 Mt, respectively. This work can be expanded in many different directions to build a more predictive model for fallout formation.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Near field excited state imaging via stimulated electron energy gain spectroscopy of localized surface plasmon resonances in plasmonic nanorod antennas

Continuous wave (cw) photon stimulated electron energy loss and gain spectroscopy (sEELS and sEEGS) is used to image the near field of optically stimulated localized surface plasmon resonance (LSPR) modes in nanorod antennas. An optical delivery system equipped with a nanomanipulator and a fiber-coupled laser diode is used to simultaneously irradiate plasmonic nanostructures in a (scanning) transmission electron microscope. The nanorod length is varied such that the m = 1, 2, and 3 LSPR modes are resonant with the laser energy and the optically stimulated near field spectra and images of these modes are measured. Various nanorod orientations are also investigated to explore retardation effects. Optical and electron beam simulations are used to rationalize the observed patterns. As expected, the odd modes are optically bright and result in observed sEEG responses. The m = 2 dark mode does not produce a sEEG response, however, when tilted such that retardation effects are operative, the sEEG signal emerges. Thus, we demonstrate that cw sEEGS is an effective tool in imaging the near field of the full set of nanorod plasmon modes of either parity.

36 MATERIALS SCIENCE↗

Impact of spatiotemporal smoothing on the two-plasmon–decay instability

Higher levels of hot electrons from the two-plasmon–decay instability are observed when smoothing by spectral dispersion (SSD) is turned off in directly driven inertial confinement fusion experiments at the Omega Laser Facility. This finding is explained using a hot-spot model based on speckle statistics and simulation results from the laser–plasma simulation environment. The model accurately reproduces the relative increase in hot-electron activity at two different drive intensities, although it slightly overestimates the absolute number of hot electrons in all cases. Extrapolating from the current ≈ 360-GHz system while adhering to the logic of the hot-spot model suggests that larger SSD bandwidth should significantly mitigate hot-electron generation, and legacy 1-THz OMEGA experiments appear to support this conclusion. These results demonstrate that it is essential to account for laser speckles and spatiotemporal smoothing to obtain quantitative agreement with experiments. A compilation of hot-electron data from the past two decades reveals several other important points: (1) many prior experiments are more easily understood using recent results from multibeam absolute instability theory; and (2) experiments with ignition-scale conditions produce less hot electrons compared to OMEGA spherical experiments for a given vacuum overlapped intensity, which is a promising result for validating performance predictions based on hydrodynamic scaling relations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Pressure, temperature, and orientation dependent thermal conductivity of pentaerythritol tetranitrate (PETN)

We use reverse non-equilibrium molecular dynamics to determine the thermal conductivity tensor, κ, of tetragonal pentaerythritol tetranitrate (PETN). The most stable form under standard ambient conditions (PETN I, with space group P4¯2 1 c) is considered in the temperature and pressure intervals 200–500 K and 0–8 GPa, which covers approximately the stability range for this polymorph. We compute κ along the [100] and [001] directions, which is enough to construct the full thermal conductivity tensor for a system with tetragonal symmetry. In addition, we also determine κ along the [010] direction and confirm that κ 100 ~ κ 010 , with less than 4% average absolute error between the two quantities. We observe an anisotropic response for κ, with κ 100 > κ 001 across the whole (T,P) interval, and 37% difference at 300 K, 0 GPa. Furthermore, we provide analytical functions to interpolate κ(P,T) within the fitting interval and serve as input for continuum-scale simulations.

36 MATERIALS SCIENCE↗

Implications of phonon anisotropy on thermal conductivity of fluorite oxides

Abstract Fluorite oxides are attractive ionic compounds for a range of applications with critical thermal management requirements. In view of recent reports alluding to anisotropic thermal conductivity in this face-centered cubic crystalline systems, we perform a detailed analysis of the impact of direction-dependent phonon group velocities and lifetimes on the thermal transport of fluorite oxides. We demonstrate that the bulk thermal conductivity of this class of materials remains isotropic despite notable anisotropy in phonon lifetime and group velocity. However, breaking the symmetry of the phonon lifetime under external stimuli including boundary scattering present in nonequilibrium molecular dynamics simulations of finite size simulation cell gives rise to apparent thermal conductivity anisotropy. We observe that for accurate determination of thermal conductivity, it is important to consider phonon properties not only along high symmetry directions commonly measured in inelastic neutron or x-ray scattering experiments but also of those along lower symmetry. Our results suggests that certain low symmetry directions have a larger contribution to thermal conductivity compared to high symmetry ones.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗