Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiple time stepping”

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 199 records · Page 11

Restoring Distribution System Under Renewable Uncertainty Using Reinforcement Learning: Preprint

Distributed energy resources (DER) in distribution systems, including renewable generation, micro-turbine, and energy storage, can be used to restore critical loads following extreme events to increase grid resiliency. However, properly coordinating multiple DERs in the system for multi-step restoration process under renewable uncertainty and fuel availability is a complicated sequential optimal control problem. Due to its capability to handle system non-linearity and uncertainty, reinforcement learning (RL) stands out as a potentially powerful candidate in solving complex sequential control problems. Moreover, the offline training of RL provides excellent action readiness during online operation, making it suitable to problems such as load restoration, where in-time, correct and coordinated actions are needed. In this study, a distribution system prioritized load restoration based on a simplified single-bus system is studied: with imperfect renewable generation forecast, the performance of an RL controller is compared with that of a deterministic model predictive control (MPC). Our experiment results show that the RL controller is able to learn from experience, adapt to the imperfect forecast information and provide a more reliable restoration process when compared with the baseline MPC controller.

61 RADIATION PROTECTION AND DOSIMETRY↗

Restoring Distribution System Under Renewable Uncertainty Using Reinforcement Learning

Distributed energy resources (DER) in distribution systems, including renewable generation, micro-turbine, and energy storage, can be used to restore critical loads following extreme events to increase grid resiliency. However, properly coordinating multiple DERs in the system for multi-step restoration process under renewable uncertainty and fuel availability is a complicated sequential optimal control problem. Due to its capability to handle system non-linearity and uncertainty, reinforcement learning (RL) stands out as a potentially powerful candidate in solving complex sequential control problems. Moreover, the offline training of RL provides excellent action readiness during online operation, making it suitable to problems such as load restoration, where in-time, correct and coordinated actions are needed. In this study, a distribution system prioritized load restoration based on a simplified single-bus system is studied: with imperfect renewable generation forecast, the performance of an RL controller is compared with that of a deterministic model predictive control (MPC). Our experiment results show that the RL controller is able to learn from experience, adapt to the imperfect forecast information and provide a more reliable restoration process when compared with the baseline MPC controller.

61 RADIATION PROTECTION AND DOSIMETRY↗

Operator learning for predicting multiscale bubble growth dynamics

We report simulating and predicting multiscale problems that couple multiple physics and dynamics across many orders of spatiotemporal scales is a great challenge that has not been investigated systematically by deep neural networks (DNNs). Herein, we develop a framework based on operator regression, the so-called deep operator network (DeepONet), with the long-term objective to simplify multiscale modeling by avoiding the fragile and time-consuming “hand-shaking” interface algorithms for stitching together heterogeneous descriptions of multiscale phenomena. To this end, as a first step, we investigate if a DeepONet can learn the dynamics of different scale regimes, one at the deterministic macroscale and the other at the stochastic microscale regime with inherent thermal fluctuations. Specifically, we test the effectiveness and accuracy of the DeepONet in predicting multirate bubble growth dynamics, which is described by a Rayleigh–Plesset (R–P) equation at the macroscale and modeled as a stochastic nucleation and cavitation process at the microscale by dissipative particle dynamics (DPD). First, we generate data using the R–P equation for multirate bubble growth dynamics caused by randomly time-varying liquid pressures drawn from Gaussian random fields (GRFs). Our results show that properly trained DeepONets can accurately predict the macroscale bubble growth dynamics and can outperform long short-term memory networks. We also demonstrate that the DeepONet can extrapolate accurately outside the input distribution using only very few new measurements. Subsequently, we train the DeepONet with DPD data corresponding to stochastic bubble growth dynamics. Although the DPD data are noisy and we only collect sparse data points on the trajectories, the trained DeepONet model is able to predict accurately the mean bubble dynamics for time-varying GRF pressures. Taken together, our findings demonstrate that DeepONets can be employed to unify the macroscale and microscale models of the multirate bubble growth problem, hence providing new insight into the role of operator regression via DNNs in tackling realistic multiscale problems and in simplifying modeling with heterogeneous descriptions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reversible inhibition and reactivation of electron transfer in photosystem I

In photosystem I (PSI) complexes at room temperature light induced electron transfer can proceed down multiple protein-cofactor branches, and the process can be an order of magnitude faster on one branch compared to the other. One factor that might contribute to this branch asymmetry is a tryptophan amino acid at position 673 on the psaB protein. This amino acid is located between two of the pigments involved in electron transfer on the psaB protein branch. The corresponding residue on the psaA protein branch is a glycine amino acid. Microsecond time-resolved step-scan FTIR difference spectroscopy at 77 K has been used to study isolated PSI complexes from wild type and a mutant where the tryptophan residue was changed to phenylalanine. Photoaccumulated FTIR difference spectra indicate changes in the protein structure upon mutation. In the mutant we also find that the electron transfer processes is inhibited following long periods of repetitive flash illumination at room temperature. This is due to double protonation of the pigment involved in electron transfer. However, we show that we can restore electron transfer functionality by incubating the light-treated mutant PSI samples in the presence of newly added (not protonated) pigment.

59 BASIC BIOLOGICAL SCIENCES↗

Baltimore Social-Environmental Collaborative (BSEC) Doppler Lidar & Derived Products

This repository contains all processed Doppler‐lidar outputs from the PSU lidar deployed for the Baltimore Social‐Environmental Collaborative (BSEC) project. Vertical Stare Scans (fixed‐beam, vertical profiling): 1 Hz backscatter intensity (m⁻¹ sr⁻¹), signal‐to‐noise ratio (unitless), and Doppler vertical‐velocity (m s⁻¹) on ~30 m range gates, stored as CF-compliant NetCDF. Wind Profiles (horizontal‐wind retrieval): daily NetCDF outputs of retrieved horizontal wind speed (m s⁻¹) and direction (degrees), computed from the angled‐scan returns. Profile Statistics (summary statistics on the vertical velocity): 15 min windows (default) of mean, variance, skewness, kurtosis, high-frequency variance, etc., as a function of height; saved as CF-compliant NetCDF files. Boundary Layer Height (BLH) (fuzzy-logic output): 15 min BLH estimates (m), with lower/upper fuzzy bounds (m) and a quality flag (0–4) indicating data status (e.g., no data, good, below range, ran out of signal, cloud-topped). Cloud Base Height (Haar-gradient detection): 15 min estimates of cloud-base height (m) with a cloud-detection quality flag (0–3: none, low, moderate, high). All five product streams are organized by year and date under their own top-level folders (01_Vertical_Stare_Scans/ through 05_Cloud_Height/). Each folder contains a data_ /YYYY/ subdirectory with daily CF-compliant NetCDF outputs (96 windows per day at 15 min intervals). Global attributes in each file include creation history, version (2.0.0), institution, and source. Instrument & MeasurementsThe PSU Doppler Lidar samples aerosol backscatter (m⁻¹ sr⁻¹), signal-to-noise ratio, and radial velocity at ~1 Hz. Vertical stare scans point the beam straight up; after collecting angled scans through multiple elevation angles, the "Wind Profiles" product contains the fully retrieved horizontal wind speed and direction. Data were collected continuously at ~30 m range resolution, with a typical height ceiling of ~12 km. How to Use Open any NetCDF with Python's xarray, MATLAB, or similar CF-compliant tools. Stare scans and angled-scan retrievals (Wind Profiles) are CF-compliant daily NetCDF files. Profile-Statistics, BLH, and Cloud Height files are daily 15 min summaries (96 time steps per file). Inspect the included variables (e.g., vertical_velocity_variance, wind_speed, BLH, cloud_base_height) for your analyses. Use the quality flags (BLH_flag, cloud_flag) to filter out poor-quality retrievals. For more information or questions about processing methods, please contact:Nicholas E. Prince ⟨nec5299@psu.edu⟩Penn State Department of Meteorology & Atmospheric Science

Air Quality↗

Identifying the Gate-Opening Mechanism in the Flexible Metal–Organic Framework UTSA-300

Atomic-level understanding of the gate-opening phenomenon in flexible porous materials is an important step toward learning how to control, design, and engineer them for applications such as the separation of gases from complex mixtures. Here, we report such mechanistic insight through an in-depth study of the pressure-induced gate-opening phenomenon in our earlier reported metal–organic framework (MOF) Zn(dps) 2 (SiF 6 ) (dps = 4,4'-dipyridylsulfide), also called UTSA-300, using isotherm and calorimetry measurements, in situ infrared spectroscopy, and ab initio simulations. UTSA-300 is shown to selectively adsorb acetylene (C 2 H 2 ) over ethylene (C 2 H 4 ) and ethane (C 2 H 6 ) and undergoes an abrupt gate-opening phenomenon, making this framework a highly selective gas separator of this complex mixture. The selective adsorption is confirmed by pressure-dependent in situ infrared spectroscopy, which, for the first time, shows the presence of multiple C 2 H 2 species with varying strengths of bonding. A rare energetic feature at the gate-opening condition of the flexible MOF is observed in our differential heat energies, directly measured by calorimetry, showcasing the importance of this tool in adsorption property exploration of flexible frameworks and offering an energetic benchmark for further energy-based fundamental studies. Based on the agreement of this feature with ab initio-based adsorption energies of C 2 H 2 in the closed-pore structure UTSA-300a (“a” refers to the activated form), this feature is assigned to the weakening of the H-bond C–H···F formed between C 2 H 2 and fluorine of the MOF. Furthermore, our analysis identifies the weakening of this H-bond, the expansion of the closed-pore MOF upon successive C 2 H 2 coadsorption until its volume is close to that of the open-pore MOF, and the spontaneous gate opening to energetically favor C 2 H 2 adsorption in the open-pore structure as crucial steps in the gate-opening mechanism in this system.

36 MATERIALS SCIENCE↗

Implicit-explicit Runge-Kutta for radiation hydrodynamics I: Gray diffusion

Radiation hydrodynamics are a challenging multiscale and multiphysics set of equations. To capture the relevant physics of interest, one typically must time step on the hydrodynamics timescale, making explicit integration the obvious choice. On the other hand, the coupled radiation equations have a scaling such that implicit integration is effectively necessary in non-relativistic regimes. A first-order Lie-Trotter-like operator split is the most common time integration scheme used in practice, alternating between an explicit hydrodynamics step and an implicit radiation solve and energy deposition step. However, such a scheme is limited to first-order accuracy, and nonlinear coupling between the radiation and hydrodynamics equations makes a more general additive partitioning of the equations non-trivial. Here, we develop a new formulation and partitioning of radiation hydrodynamics with gray diffusion that allows us to apply (linearly) implicit-explicit Runge-Kutta time integration schemes. In conclusion, we prove conservation of total energy in the new framework, and demonstrate 2nd-order convergence in time on multiple radiative shock problems, achieving error 3–5 orders of magnitude smaller than the first-order Lie-Trotter operator split at the hydrodynamic CFL, even when Lie-Trotter applies a 3rd-order TVD Runge-Kutta scheme to the hydrodynamics equations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Calibrating the SPECTACULAR constitutive model using legacy Sandia data for two filled epoxy systems: 828/CTBN/DEA/GMB and 828/DEA/GMB

The SPECTACULAR model is a development extension of the Simplified Potential Energy Clock (SPEC) model. Both models are nonlinear viscoelastic constitutive models used to predict a wide range of time-dependent behaviors in epoxies and other glass-forming materials. This report documents the procedures used to generate SPECTACULAR calibrations for two particulate-filled epoxy systems, 828/CTBN/DEA/GMB and 828/DEA/GMB. No previous SPECTACULAR or SPEC calibration exists for 828/CTBN/DEA/GMB, while a legacy SPEC calibration exists for 828/DEA/GMB. To generate the SPECTACULAR calibrations, a step-by-step procedure was executed to determine parameters in groups with minimal coupling between parameter groups. This procedure has often been deployed to calibrate SPEC, therefore the resulting SPECTACULAR calibration is backwards compatible with SPEC (i.e. none of the extensions specific to SPECTACULAR are used). The calibration procedure used legacy Sandia experimental data stored on the Polymer Properties Database website. The experiments used for calibration included shear master curves, isofrequency temperature sweeps under oscillatory shear, the bulk modulus at room temperature, the thermal strain during a temperature sweep, and compression through yield at multiple temperatures below the glass transition temperature. Overall, the calibrated models fit the experimental data remarkably well. However, the glassy shear modulus varies depending on the experiment used to calibrate it. For instance, the shear master curve, isofrequency temperature sweep under oscillatory shear, and the Young's modulus in glassy compression yield values for the glassy shear modulus at the reference temperature that vary by as much as 15 %. Also, for 828/CTBN/DEA/GMB, the temperature dependence of the glassy shear modulus when fit to the Young's modulus at different temperatures is approximately four times larger than when it is determined from the isofrequency temperature sweep under oscillatory shear. For 828/DEA/GMB, the temperature dependence of the shear modulus determined from the isofrequency temperature sweep under oscillatory shear accurately predicts the Young's modulus at different temperatures. When choosing values for the shear modulus, fitting the glassy compression data was prioritized. The new and legacy calibrations for 828/DEA/GMB are similar and appear to have been calibrated from the same data. However, the new calibration improves the fit to the thermal strain data. In addition to the standard calibrations, development calibrations were produced that take advantage of development features of SPECTACULAR , including an updated equilibrium Helmholtz free energy that eliminates undesirable behavior found in previous work. In addition to the previously mentioned experimental data, the development calibrations require data for the heat capacity during a stress-free temperature sweep to calibrate thermal terms.

36 MATERIALS SCIENCE↗

An Indicator-based Approach to Sustainable Management of Natural Resources (Chapter 12)

Assessing the sustainability of natural resource management choices for agricultural and forest lands requires quantification of potential changes to a set of environmental and socioeconomic indicators selected to characterize reference scenarios relative to projected future scenarios. Correctly framing the questions with local stakeholders is a critical first step in the sustainability assessment, and the questions that can be addressed are often limited by data availability. Selecting and prioritizing indicators with stakeholders to address their needs and concerns improves the likelihood of investment in monitoring and evaluation of those indicators over time. Computational techniques for analyzing interactions between the selected indicators are inherently affected by the scales and formats of the assembled indicator datasets. Data analytics have the potential to improve understanding of the potential synergies and tradeoffs involved with meeting multiple environmental and socioeconomic goals simultaneously, but timely and appropriate indicator datasets are not always available—even in this new era of “big data.” Continued improvements in data science and data analytics are needed to broaden understanding and acceptance of problems and to provide valuable information for natural resource management. Advances in these areas will enable society to design future landscapes that meet multiple objectives, including the provisioning of agricultural and forest resources along with a variety of ecosystem services (e.g., clean water and healthy soils).

Parish, Esther↗

The phase-2 particle x-ray temporal diagnostic for simultaneous measurement of multiple x-ray and nuclear emission histories from OMEGA implosions (invited)

Electron-temperature (T e ) measurements in implosions provide valuable diagnostic information, as T e is negligibly affected by residual flows and other non-thermal effects unlike ion-temperature inferred from a fusion product spectrum. In OMEGA cryogenic implosions, measurement of T e (t) can be used to investigate effects related to time-resolved hot-spot energy balance. The newly implemented phase-2 Particle X-ray Temporal Diagnostic (PXTD) utilizes four fast-rise (~15 ps) scintillator-channels with distinct x-ray filtering. Titanium and stepped aluminum filtering were chosen to maximize detector sensitivity in the 10–20 keV range, as it has been shown that these x rays have similar density and temperature weighting to the emitted deuterium–tritium fusion neutrons (DTn) from OMEGA Cryo-DT implosions. High quality data have been collected from warm implosions at OMEGA. These data have been used to infer spatially integrated T e (t) with <10% uncertainty at peak emission. Nuclear and x-ray emission histories are measured with 10 ps relative timing uncertainty for x rays and DTn and 12 ps for x rays and deuterium- 3 He protons (D 3 Hep). A future upgrade to the system will enable spatially integrated T e (t) with 40 ps time-resolution from cryogenic DT implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High-resolution ion mobility based on traveling wave structures for lossless ion manipulation resolves hidden lipid features

Abstract High-resolution ion mobility (resolving power > 200) coupled with mass spectrometry (MS) is a powerful analytical tool for resolving isobars and isomers in complex samples. High-resolution ion mobility is capable of discerning additional structurally distinct features, which are not observed with conventional resolving power ion mobility (IM, resolving power ~ 50) techniques such as traveling wave IM and drift tube ion mobility (DTIM). DTIM in particular is considered to be the “gold standard” IM technique since collision cross section (CCS) values are directly obtained through a first-principles relationship, whereas traveling wave IM techniques require an additional calibration strategy to determine accurate CCS values. In this study, we aim to evaluate the separation capabilities of a traveling wave ion mobility structures for lossless ion manipulation platform integrated with mass spectrometry analysis (SLIM IM-MS) for both lipid isomer standards and complex lipid samples. A cross-platform investigation of seven subclass-specific lipid extracts examined by both DTIM-MS and SLIM IM-MS showed additional features were observed for all lipid extracts when examined under high resolving power IM conditions, with the number of CCS-aligned features that resolve into additional peaks from DTIM-MS to SLIM IM-MS analysis varying between 5 and 50%, depending on the specific lipid sub-class investigated. Lipid CCS values are obtained from SLIM IM ( TW(SLIM) CCS) through a two-step calibration procedure to align these measurements to within 2% average bias to reference values obtained via DTIM ( DT CCS). A total of 225 lipid features from seven lipid extracts are subsequently identified in the high resolving power IM analysis by a combination of accurate mass-to-charge, CCS, retention time, and linear mobility-mass correlations to curate a high-resolution IM lipid structural atlas. These results emphasize the high isomeric complexity present in lipidomic samples and underscore the need for multiple analytical stages of separation operated at high resolution. Graphical abstract

Reardon, Allison R. (ORCID:0000000165830134)↗

Producing Evaluation-quality 239 Pu Average Prompt Fission Neutron Multiplicities using a Correlated Fission Model

An evaluation of the average prompt fission neutron multiplicity, $\bar{ν}_p$, of 239 Pu(n,f) is shown. This evaluation includes (a) the correlated fission model CGMF, and (b) a detailed analysis of past and recently published experimental data. Using CGMF-calculated $\bar{ν}_p$ as prior enables to link, through the use of evaluated model input parameters, $\bar{ν}_p$ to other fission observables such as the prompt fission neutron spectrum (PFNS), preneutron emission fission yields as a function of mass, and the average total kinetic energy of the fragments. These evaluated parameters produce realistic predictions of many fission observables, while the evaluated $\bar{ν}_p$ agrees well (χ 2 ≈ 1) with data. Moreover, with the new evaluated $\bar{ν}_p$, the effective neutron multiplication factor of fast Pu ICSBEP critical assemblies are predicted with a mean bias of 58 pcm compared to 18 pcm with ENDF/B-VIII.0, when paired with a new 239 Pu PFNS and fission cross section. Due to these encouraging validation results, the evaluated $\bar{ν}_p$ is currently part of a release candidate for the 239 Pu ENDF/B-VIII.1 file. Hence, a correlated fission model was used for the first time for evaluating $\bar{ν}_p$ that is of evaluation quality. This is an important step towards consistent evaluations of prompt fission observables.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Foundations of an Electric Mobility Strategy for the City of Mexicali

The Foundations of an Electric Mobility Strategy for the city of Mexicali aligns with numerous energy, environmental, and transport plans and will help Mexicali meet multiple related goals. Mexicali’s energy mix, with 28% renewables, already enables plugin electric vehicles (PEVs) to reduce the mass of greenhouse gases (GHGs) per km driven 2/3 below that of their conventional counterparts. This GHG benefit will increase should Mexicali take steps to further increase their share of renewables in their electricity supply. Beyond increasing renewables, Mexicali could possibly deploy PEVs so that electric load is added in the right location (depending on further analysis of substations and feeders) and at the right time (between 21:00 and 11:00) in order to minimize grid upgrade costs. There are a handful of charge timing control mechanisms –at various stages of development– that Mexicali could implement. Transport electrification can facilitate mass transit by powering buses, trains, and small vehicles that get people from their homes or work to the transit stations and vice versa. Mexicali could utilize fleets as early PEV adopters in order to gain acceptance and add electric vehicle supply equipment (EVSE). Recommended prioritization of different types of fleets are suggested in this report: transit buses, school buses, airport ground support equipment (GSE), refuse trucks, taxis, shuttle buses, campus vehicles, delivery trucks, utility trucks, and finally semitrailers. There are a handful of policy options that Mexicali could use to incentivize fleets to purchase PEVs, including mandates, economic incentives, energy performance contracts, waivers to access restrictions, electricity discounts, and EVSE requirements in building codes. Mexicali’s taxi fleet was an early adopter of PEVs and had experienced some challenges—mostly related to the insufficient range of the taxis due to hot weather. In this report, we strategize ways to extend the range of the current electric taxis, including ways to make charging more convenient to the drivers, and more suggestions for appropriate vehicles to purchase in the future. This report also includes the groundwork of geotracking Mexicali’s taxi fleet so that more detailed recommendations can be made in the future. Once fleets have increased PEV acceptance and EVSE installations, the market will be ready to expand to private vehicle owners. In order to do this, more EVSE needs to be installed in the right locations. This mobility strategy lays out general local areas where EVSE could be well utilized, based on traffic patterns, land use, and demographic data. Mexicali could approach businesses within these areas that would likely make suitable hosts, based on how well they can profit from the additional business that EVSE would bring. Mexicali could then adopt a series of purchase incentives (including sales tax waivers or access to high-occupancy vehicles [HOV] lanes) that would encourage private vehicle owners to purchase PEVs. Purchase incentives run the risk of creating equity issues, which can be countered by promoting electrification in mass transit, creating more HOV lanes that have PEV exemptions, and installing EVSEs in underserved communities. Private PEV ownership will require a set of experts that Mexicali can help train, including PEV repair technicians, EVSE installation electricians, and first responders.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Development of a high-throughput method for processing sponge-stick samples to detect viable Bacillus anthracis spores

Since the national validation of the sponge-stick based method for detection of Bacillus anthracis spores in environmental samples, there have not been focused efforts to address the low throughput nature of the method, which processes only one sample at one time. Sample processing remains a serious bottleneck for rapidly analyzing large numbers of samples expected from a biological warfare attack. Therefore, we developed a high-throughput method to simultaneously process multiple sponge-stick samples to be better prepared for rapid response and recovery after wide area anthrax incidents. In this method, sponges are placed in 50 mL tubes containing 25 mL extraction buffer and shaken to release spores, after which the suspension is recovered for analysis. Here, we determined that an additional extraction step, conducted in the same tubes with 10 mL buffer, further increased spore recovery from sponge-stick by approximately 10 %. We determined that orbital shaking and multi-tube vortexing were both more effective than reciprocating shaking for recovering spores. We conducted simultaneous processing of up to 12 sponge-stick samples and demonstrated comparable spore recovery efficiencies to the traditional low-throughput stomacher-based method (approximately 60 % recovery at 10 2 -spore level and 75 % recovery at 10 4 -spore level for both methods in three replicate experiments, P > 0.05 for two-tailed t-tests for each experiment and spore level). We also demonstrated that our high-throughput method could be integrated with Rapid Viability-Polymerase Chain Reaction (RV-PCR) analysis and could detect levels as low as 40 spores per sponge even when challenged by a PCR particulate contaminant.

Anthrax↗

Progress on Demonstration of a MOOSE-Based Coupled Capability for Hot Channel Factors in Fast Reactors

Hot channel factors (HCFs) are computed values that account for the impact on predicted peak fuel, cladding, and coolant temperatures due to uncertainties in the as-built reactor’s material properties and geometry as well as uncertainties due to modeling approximations. Reduction in computed HCF values via reduction or elimination of modeling approximations may translate to significant economic savings if the reactor power can be raised due to the extra temperature margin gained. While limited historical datasets exist for sodium-cooled fast reactors (SFRs), there are no available HCF data for lead-cooled fast reactors (LFRs) outside of work generated previously within NEAMS. The computation of HCFs involves insights from reactor physics, thermal fluids and heat conduction calculations to determine how the peak temperatures respond to various uncertainties in the design. Due to the significant advantages for multi-physics coupling offered by the MOOSE framework, Griffin (MOOSE-based reactor physics code), MOOSE Heat Conduction Module, and Cardinal (MOOSE-wrapped multi-physics application which includes the NekRS thermal fluids code) are being coupled together using the MOOSE MultiApp System to develop a highfidelity multi-physics modeling capability for HCF simulations. This high-fidelity coupling workflow may also be beneficial for other fast reactor applications in the future. In previous work, Griffin and NekRS were individually assessed to ensure the necessary capabilities were in place. This work describes initial efforts to couple the codes (including folding in the MOOSE Heat Conduction Module) and determining the workflow for the perturbed calculations which will leverage the Stochastic Tools Module (STM). To our knowledge, this is the first coupling of Griffin and NekRS as well as the first exploratory use of Stochastic Tools Module for Cardinal. In this report, the neutronics code Griffin, the heat conduction solver in MOOSE, and the MOOSE-wrapped application containing NekRS (Cardinal) are linked together to demonstrate the coupled capability. Griffin and Cardinal are linked dynamically by specifying shared libraries. Different coupling hierarchies are tested for selecting the most appropriate coupling strategy. A coupling scheme is selected based on the efficiency of calculation and ease of data communication. Multiple tests are performed to choose suitable mesh structure, model configurations, scheme setup and boundary conditions to avoid loss of energy due to data interpolation between different modules or weak imposition of fluxes in finite element codes. Computational experiments are performed to study the tolerance control of each type of iteration to avoid false convergence. The coupled capability is demonstrated in both single pin and 7-pin models based on LFR materials and geometry. The study finds that the use of too large a time step size in the heat conduction module can lead to temperature oscillation even though the heat conduction equation does not have a time-derivative kernel, but only the time-dependent boundary condition. A 7-pin model without duct region achieved good convergence in the coupled calculation while a 7 pin model with duct region experienced data communication issues which need to be resolved.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Extended use of face masks during the COVID-19 pandemic - Thermal conditioning and spray-on surface disinfection

The current COVID-19 pandemic has resulted in globally constrained supplies for face masks and personal protective equipment (PPE). Production capacity is limited in many countries and the future course of the pandemic will likely continue with shortages for high quality masks and PPE in the foreseeable future. Therefore, expectations are that mask reuse, extended wear and similar approaches will enhance the availability of personal protective measures. Repeated thermal disinfection could be an important option and likely easier implemented in some situations, at least on the small scale, than UV illumination, irradiation or hydrogen peroxide vapor exposure. In this paper, an overview on thermal responses and ongoing filtration performance of multiple face mask types is provided. Most masks have adequate material properties to survive a few cycles (i.e. 30 min disinfection steps) of thermal exposure in the 75°C regime. Some are more easily affected, as seen by the fusing of plastic liner or warping, given that preferred conditioning temperatures are near the softening point for some of the plastics and fibers used in these masks. Hence adequate temperature control is equally important. As guidance, disinfectants sprayed via dilute solutions maintain a surface presence over extended time at 25 and 37°C. Some spray-on alcohol-based solutions containing disinfectants were gently applied to the top surface of masks. Neither moderate thermal aging (less than 24 h at 80 and 95°C) nor gentle application of surface disinfectant sprays resulted in measurable loss of mask filter performance. Subject to bio-medical concurrence (additional checks for virus kill efficiency) and the use of low risk non-toxic disinfectants, such strategies, either individually or combined, by offering additional anti-viral properties or short term refreshing, may complement reuse options of professional masks or the now ubiquitous custom-made face masks with their often unknown filtration effectiveness.

36 MATERIALS SCIENCE↗

Measuring and Extracting Activity from Time Series Data

This report summarizes the results of an LDRD focused on developing and demonstrating statistically rigorous methods for analyzing and comparing complex activities from remote sensing data. Identifying activity from remote sensing data, particularly those that play out over time and span multiple locations, often requires extensive manual effort because of the variety of features that describe the activity and the required domain expertise. Our results suggest that there are some hidden challenges in extracting and representing activities in sensor data. In particular, we found that the variability in the underlying behaviors can be difficult to overcome statistically, and the report identifies several examples of the issue. We discuss key lessons learned in the context of the project, and finally conclude with recommendations on next steps and future work.

97 MATHEMATICS AND COMPUTING↗

Detection, monitoring, and determination of location of changes in metallic structures using multimode acoustic signals

Acoustic transducers generate and receive acoustic signals at multiple locations along a surface of rigid structure, wherein longitudinal spacing between transducer locations define measurement zones. Acoustic signals with chosen amplitude-time-frequency characteristics excite multiple vibration modes in the structure within each zone. Small mechanical changes in inspection zones lead to scattering and attenuation of broadband acoustic signals, which are detectable as changes in received signal characteristics as part of a through-transmission technique. Additional use of short, narrowband pulse acoustic signals as part of a pulse-echo technique allows determination of the relative location of the mechanical change within each zone based on the differential delay profiles. For accurate acoustic modeling and simulation, the mesh size, time step, time delay, and time-window size are optimized. Frequency normalization of the Short-Time Fourier Transform of acoustic response output improves experiment-simulation cross-validation. Applications of the method to structures with arbitrarily complex geometries are also demonstrated.

Findikoglu, Alp Tugrul↗