Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Accounting”

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 703 records · Page 39

The Completion of a Geosynchronous Earth Orbit Survey with the Eugene Stansbery-Meter Class Autonomous Telescope

The Eugene Stansbery-Meter Class Autonomous Telescope (ES-MCAT) is the primary optical sensor used by the NASA Orbital Debris Program Office (ODPO) to statistically characterize the geosynchronous Earth orbit (GEO) debris environment and support future Orbital Debris Engineering Model (ORDEM) releases. The ES-MCAT completed its first optical survey of the GEO region from 2020 to 2022. The primary goal of this survey was to autonomously collect and process GEO data with calculated photometric and astrometric uncertainties. A pointing plan was developed to provide uniform sampling within the region of interest (ROI) while accounting for predicted downtime due to insufficient observing conditions. Detections are autonomously correlated to the Space Surveillance Network (SSN) catalog to determine if objects are correlated targets (CTs) or uncorrelated targets (UCTs), the latter of which are of interest for modeling the GEO orbital debris environment. To assess the size detection sensitivity over time and monitor the general performance of the telescope’s optics and software, the optical throughput of the system and limiting magnitudes are evaluated on a routine basis. While the telescope’s ability to operate autonomously and remotely allowed for the GEO survey to continue throughout the COVID-19 pandemic, travel restrictions hampered routine cleaning of the optics during this time, and the primary mirror degraded enough to require recoating. The mirror was removed in 2022, concluding the first GEO survey. The primary mirror received a new coating designed to be more robust against the harsh environment surrounding the ES-MCAT’s location on Ascension Island, accounting for experience gained during operations over the first GEO survey. In early 2023, the recoated primary mirror was reinstalled, and the second GEO survey was initiated. The primary goal of the second GEO survey is to characterize the evolving GEO debris environment with updated optics, software, and pointing strategies while allowing for the inclusion of non-GEO regimes or those that are outside of the ROI. While the pointing method implemented in the first survey allowed for adequate coverage of the ROI over two years, it has been improved to include pointings that avoid the Moon’s position and the galactic plane to reduce software processing time and maximize the detection capabilities of fainter objects. This method also accounts for the changing weather patterns throughout the year and reduces coverage gaps in the ROI. Provided the success of the first two-year GEO survey using autonomous operations, the ODPO is actively collaborating with the United States Space Force (USSF) to make the ES-MCAT a contributing sensor to the SSN. This paper presents results from the first GEO survey including magnitude distributions and orbital parameters for CTs and UCTs. Details are provided for the automated processing pipeline and the optical system throughput for the previous and current primary mirror coatings. In addition, an updated strategy for the second GEO survey to optimize coverage over the ROI is discussed, as are preliminary results from the ongoing second survey.

Corbin Cruz↗

The Transformation and Export of Organic Carbon Across an Arctic River-Delta-Ocean Continuum

The Arctic Ocean is surrounded by land that feeds highly seasonal rivers with water enriched in high concentrations of dissolved and particulate organic carbon (DOC and POC). Explicit estimates of the flux of organic carbon across the land-ocean interface are difficult to quantify and many interdependent processes makes source attribution difficult. A high-resolution 3-D biogeochemical model was built for the lower Yukon River and coastal ocean to estimate biogeochemical cycling across the land-ocean continuum. The model solves for complex reactions related to organic carbon transformation, including mechanistic photodegradation and multi-reactivity microbial processing, DOC to POC flocculation, and phytoplankton dynamics. The baseline DOC and POC flux out of the delta from April to September 2019, was 977 and 536 Gg C (~80% of the annual total), but only 50% of the DOC and 25% of the POC exited the plume across the 10 m isobath. Microbial breakdown of DOC accounted for a net loss of 168 Gg C (17% of delta export) within the plume and photodegradation accounted for a net loss of 46.6 Gg C DOC (5% of delta export) in 2019. Flocculation decreased the total organic carbon flux by only 6.4 Gg C (~1%), while POC sinking accounted for 63.3 Gg C (10%) settling in the plume. The loss of chromophoric dissolved organic matter (CDOM) due to photodegradation increased the light available for phytoplankton growth throughout the coastal ocean, demonstrating the secondary effects that organic carbon reactions can have on biological processes and the net coastal carbon flux.

Land-Ocean Continuum↗

On the Remote Impacts of Mid-Holocene Saharan Vegetation on South American Hydroclimate: A Modeling Intercomparison

Proxy reconstructions from the mid-Holocene (MH: 6,000 years ago) indicate an intensification of the West African Monsoon and a weakening of the South American Monsoon, primarily resulting from orbitally-driven insolation changes. However, model studies that account for MH orbital configurations and greenhouse gas concentrations can only partially reproduce these changes. Most model studies do not account for the remarkable vegetation changes that occurred during the MH, in particular over the Sahara, precluding realistic simulations of the period. Here, we study precipitation changes over northern Africa and South America using four fully coupled global climate models by accounting for the Saharan greening. Incorporating the Green Sahara amplifies orbitally-driven changes over both regions, and leads to an improvement in proxy-model agreement. Our work highlights the local and remote impacts of vegetation and the importance of considering vegetation changes in the Sahara when studying and modeling global climate.

Paleoclimate modeling↗

Earth-System-Model Evaluation of Cloud and Precipitation Occurrence for Supercooled and Warm Clouds Over the Southern Ocean's Macquarie Island

Over the remote Southern Ocean (SO), cloud feedbacks contribute substantially to Earth system model (ESM) radiative biases. The evolution of low Southern Ocean clouds (cloud-top heights < ∼ 3 km) is strongly modulated by precipitation and/or evaporation, which act as the primary sink of cloud condensate. Constraining precipitation processes in ESMs requires robust observations suitable for process-level evaluations. A year-long subset (April 2016–March 2017) of ground-based profiling instrumentation deployed during the Macquarie Island Cloud and Radiation Experiment (MICRE) field campaign (54.5∘ S, 158.9∘ E) combines a 95 GHz (W-band) Doppler cloud radar, two lidar ceilometers, and balloon-borne soundings to quantify the occurrence frequency of precipitation from the liquid-phase cloud base. Liquid-based clouds at Macquarie Island precipitate ∼ 70 % of the time, with deeper and colder clouds precipitating more frequently and at a higher intensity compared to thinner and warmer clouds. Supercooled cloud layers precipitate more readily than layers with cloud-top temperatures > 0 ∘C, regardless of the geometric thickness of the layer, and also evaporate more frequently. We further demonstrate an approach to employ these observational constraints for evaluation of a 9-year GISS-ModelE3 ESM simulation. Model output is processed through the Earth Model Column Collaboratory (EMC2) radar and lidar instrument simulator with the same instrument specifications as those deployed during MICRE, therefore accounting for instrument sensitivities and ensuring a coherent comparison. Relative to MICRE observations, the ESM produces a smaller cloud occurrence frequency, smaller precipitation occurrence frequency, and greater sub-cloud evaporation. The lower precipitation occurrence frequency by the ESM relative to MICRE contrasts with numerous studies that suggest a ubiquitous bias by ESMs to precipitate too frequently over the SO when compared with satellite-based observations, likely owing to sensitivity limitations of spaceborne instrumentation and different sampling methodologies for ground- versus space-based observations. Despite these deficiencies, the ESM reproduces the observed tendency for deeper and colder clouds to precipitate more frequently and at a higher intensity. The ESM also reproduces specific cloud regimes, including near-surface clouds that account for ∼ 25 % of liquid-based clouds during MICRE and optically thin, non-precipitating clouds that account for ∼ 27 % of clouds with bases higher than 250 m. We suggest that the demonstrated framework, which merges observations with appropriately constrained model output, is a valuable approach to evaluate processes responsible for cloud radiative feedbacks in ESMs.

cloud feedbacks↗

The Completion of a Geosynchronous Earth Orbit Survey with the Eugene Stansbery-Meter Class Autonomous Telescope

The Eugene Stansbery-Meter Class Autonomous Telescope (ES-MCAT) is the primary optical sensor used by the NASA Orbital Debris Program Office (ODPO) to statistically characterize the geosynchronous Earth orbit (GEO) debris environment and support future Orbital Debris Engineering Model (ORDEM) releases. The ES-MCAT completed its first optical survey of the GEO region from 2020 to 2022. The primary goal of this survey was to autonomously collect and process GEO data with calculated photometric and astrometric uncertainties. A pointing plan was developed to provide uniform sampling within the region of interest (ROI) while accounting for predicted downtime due to insufficient observing conditions. Detections are autonomously correlated to the Space Surveillance Network (SSN) catalog to determine if objects are correlated targets (CTs) or uncorrelated targets (UCTs), the latter of which are of interest for modeling the GEO orbital debris environment. To assess the size detection sensitivity over time and monitor the general performance of the telescope’s optics and software, the optical throughput of the system and limiting magnitudes are evaluated on a routine basis. While the telescope’s ability to operate autonomously and remotely allowed for the GEO survey to continue throughout the COVID-19 pandemic, travel restrictions hampered routine cleaning of the optics during this time, and the primary mirror degraded enough to require recoating. The mirror was removed in 2022, concluding the first GEO survey. The primary mirror received a new coating designed to be more robust against the harsh environment surrounding the ES-MCAT’s location on Ascension Island, accounting for experience gained during operations over the first GEO survey. In early 2023, the recoated primary mirror was reinstalled, and the second GEO survey was initiated. The primary goal of the second GEO survey is to characterize the evolving GEO debris environment with updated optics, software, and pointing strategies while allowing for the inclusion of non-GEO regimes or those that are outside of the ROI. While the pointing method implemented in the first survey allowed for adequate coverage of the ROI over two years, it has been improved to include pointings that avoid the Moon’s position and the galactic plane to reduce software processing time and maximize the detection capabilities of fainter objects. This method also accounts for the changing weather patterns throughout the year and reduces coverage gaps in the ROI. Provided the success of the first two-year GEO survey using autonomous operations, the ODPO is actively collaborating with the United States Space Force (USSF) to make the ES-MCAT a contributing sensor to the SSN. This paper presents results from the first GEO survey including magnitude distributions and orbital parameters for CTs and UCTs. Details are provided for the automated processing pipeline and the optical system throughput for the previous and current primary mirror coatings. In addition, an updated strategy for the second GEO survey to optimize coverage over the ROI is discussed, as are preliminary results from the ongoing second survey.

Corbin Cruz↗

Simulations of Cryogenic Line Chilldown with Advanced Sub-Grid Wall Boiling Models

A meso-scale model developed at MIT [2][11] for boiling processes in water was adapted for cryogens and demonstrated for chilldown of propellant transfer lines in both liquid nitrogen and hydrogen. The sub-grid boiling model accurately captures the contributions to heat transfer from the generation of bubble nuclei, growth, and interaction of the bubbles in the microlayer as well as quenching of the boiling surface following bubble departure. It was adapted for cryogenic fluids using thermodynamic scaling concepts taking into account non-dimensional pressures and temperatures that are scaled by the corresponding critical values for the fluid. The boiling model was demonstrated for line chilldown in liquid nitrogen (Darr et al. [6]). The predicted wall temperature at which quenching occurs was close to the test data while the slope of the temperature curve after quenching is initiated, showing a steeper variation than the test data. Simulations were also performed for liquid hydrogen by simulating experiments of Hartwig et al. [7]. Chilldown times in liquid hydrogen are much more rapid due to higher heat transfer in the film boiling regime and accounting for the higher turbulence levels were found to be important. Furthermore, at the much lower fluid temperatures in liquid hydrogen flows, accounting for the variable thermal properties of the solid material is critical and has a dramatic impact on the quench times. The need for additional studies to better understand the evolution of the quench front in liquid hydrogen was noted.

Computational Fluid Dynamics↗

An Urban Drainage Scheme for Large-Scale Flood Models

As flood modeling spatial resolutions get finer, physical processes normally neglected, such as urban drainage, must be accounted for. Here, we describe and evaluate an urban drainage scheme for large-scale flood models. The parameterization accounts for urban imperviousness, and water flow over streets and through a prescribed urban drainage network. A parameter sensitivity analysis is performed during three major extreme floods over Rio de Janeiro city, Brazil, at ∼200 m spatial resolution. Results show that, compared to a hypothetical case without urban drainage, representing a drainage network decreases urban flooding during selected extreme events across Rio de Janeiro by 31–53 %. Such a decrease is caused by an underground water storage of up to 2.5 billion m3 across the city during flood peaks. Underground water storage and transport smooth out and delay peak flows by a few hours over major rivers and channels draining the city. Simulations also indicate that the number of residents exposed to flooding drops by 60–80 %, from ∼5 million to 1–2 million, when an urban drainage system is considered during extreme events. Similar proportions are found for social infrastructure (i.e., schools and hospitals) exposed to flooding. Results reveal that racial minority and low-income populations could disproportionally be exposed to extreme floodings across the city. We conclude that representing urban drainage has a substantial impact on flood exposure and should be accounted for in fine resolution modeling. The proposed scheme is particularly useful in poorly monitored cities and where extreme floods are a frequent hazard yet to be tackled.

Urban drainage↗

Simulations of Cryogenic Line Chilldown with Advanced Sub-Grid Wall Boiling Models

A meso-scale model developed at MIT [2][11] for boiling processes in water was adapted for cryogens and demonstrated for chilldown of propellant transfer lines in both liquid nitrogen and hydrogen. The sub-grid boiling model accurately captures the contributions to heat transfer from the generation of bubble nuclei, growth, and interaction of the bubbles in the microlayer as well as quenching of the boiling surface following bubble departure. It was adapted for cryogenic fluids using thermodynamic scaling concepts taking into account non-dimensional pressures and temperatures that are scaled by the corresponding critical values for the fluid. The boiling model was demonstrated for line chilldown in liquid nitrogen (Darr et al. [6]). The predicted wall temperature at which quenching occurs was close to the test data while the slope of the temperature curve after quenching is initiated, showing a steeper variation than the test data. Simulations were also performed for liquid hydrogen by simulating experiments of Hartwig et al. [7]. Chilldown times in liquid hydrogen are much more rapid due to higher heat transfer in the film boiling regime and accounting for the higher turbulence levels were found to be important. Furthermore, at the much lower fluid temperatures in liquid hydrogen flows, accounting for the variable thermal properties of the solid material is critical and has a dramatic impact on the quench times. The need for additional studies to better understand the evolution of the quench front in liquid hydrogen was noted.

Computational Fluid Dynamics↗

NASA GPM Status and Future Activities

The joint U.S.-Japan Global Precipitation Measurement (GPM) mission is approaching a decade of operations, and continues to pursue research, dataset production, and outreach related to precipitation. Key activities over the last year were the release of an improved “Version 07” of all GPM precipitation and latent heating products, boosting the orbit of the GPM Core Observatory (GPM CO) to 435 km, and improving quality control on precipitation retrievals from the GPM constellation of passive microwave satellites. This presentation summarizes key improvements to the GPM products and provides some examples of the changes between Versions 06 and 07 in algorithm performance. One important operational change that affected Version 07 is that the scanning strategy for the Ka-band radar channel changed in May 2018; all products that depend on Ka were revised to accommodate this change. For example, in Version 07 the Goddard Profiling (GPROF) algorithm has implemented improvements in regions where orographic enhancement and suppression take place and where the surface is snowy/icy, and again covers radiometers reaching back to 1987. The Combined Radar Radiometer Algorithm (CORRA) now incorporates modified drop-size distribution constraints that substantially reduce bias. Revisions to the Convective-Stratiform Heating (CSH) algorithm employ new radiative transfer retrievals as well as accounting for terrain in the vertical coordinates. Each algorithm was adjusted to ensure continuity for each product across the boundary in 2014 between the predecessor Tropical Rainfall Measuring Mission (TRMM) and the GPM CO. The U.S. Science Team’s Integrated Multi-satellitE Retrievals for GPM (IMERG) was upgraded to account for distortions in the probability density function of regional precipitation rates due to weighted averaging in the Kalman filter used for “morphing” the passive microwave data. Maintaining the GPM CO orbital altitude in the the current very active solar cycle has been forcing the use of more fuel than planned and consequently shortening the forecasted life of the mission from the early 2030's to the late 2020's. It was considered vital to regain some of this lifetime to ensure overlap with the upcoming Atmosphere Observing System mission to provide crosscalibration of instruments. To accomplish this, the orbital altitude was raised from 400 to 435 km on 7-8 November 2023. Thereafter, the primary GPM CO algorithms had to be revised to account for the change in observing parameters. By meeting time this action should be complete. Recently, a screening algorithm based on auto-encoding was developed that uncovered 162 orbits (out of the many thousands of orbits across all years and all satellites) of passive microwave retrievals that had highly anomalous values. Removing these defective retrievals has improved the integrity of both the GPROF and IMERG records. However, the nature of the IMERG processing interacted sufficiently badly with the now-discovered anomalous orbits that it was necessary to completely reprocess the IMERG Final Run record, now labeled Version 07B. The presentation also considers major issues that require continued attention, including the use of machine learning algorithms and the operational challenge of swarms of “small”, perhaps short-lived satellites.

GPM↗

A Deep Neural Network for Achieving Spectrally Consistent and Seamless Infrared Radiance Measurements Across Geostationary Satellite Domains

The NASA Clouds and the Earth's Radiant Energy System (CERES) project provides the scientific community with observed top-of-atmosphere (TOA) shortwave and longwave fluxes for climate monitoring and climate model validation. To achieve this goal, CERES relies on TOA broadband fluxes derived from geostationary satellite (GEO) imagery to account for the diurnal flux variations between the CERES observation intervals. Consistent global flux derivation depends on accurate and consistent cloud retrievals. Scene-dependent spectral measurement inconsistency of the instruments that make up the contiguous ring of GEO observations (GEO-Ring), as well as limb darkening effects, can cause discontinuities in derived cloud properties and radiative fluxes at the boundaries of adjacent imager domains. Although the algorithms utilize radiative transfer models to account for instrument-band-dependent atmospheric correction and viewing zenith angle (VZA) dependency, small discontinuities may persist due to uncertainties inherent to the multiple imager-specific algorithms. Furthermore, while hyperspectral-instrument-based spectral band adjustment factors may effectively account for spectrally induced bias, they are less effective at reducing variance owed to the specific composition of the viewed scene, which is challenging to robustly characterize. As such, this article highlights the use of a deep neural network (DNN) to resolve spectral-and VZA-induced biases between GEO-Ring imagers. The DNN uses available infrared (IR) channels from the GEO instruments, along with viewing and solar illumination geometry, to estimate homogenized, VIIRS-like IR radiances for use in the GEO cloud algorithm. This approach is effective at mitigating scene-dependent spectral variance and VZA dependency, resulting in consistent radiance measurements across the GEO-Ring, thereby leading toward a more seamless global cloud assessment.

deep learning↗

Engineering Design Study of Laser Power Beaming for Applications on the Moon

Using a laser to send power to a photovoltaic receiver has been proposed to transmit electrical power on the moon, particularly for applications such as powering a rover in near-polar permanent-ly-shadowed regions (PSR) where solar power is not available, however detailed engineering design studies of the spacecraft for such applications has not previously been undertaken. In this work, we did an engineering design study of two applications of laser power beaming for near-term lunar ap-plications. In the first application studied, an application was studied to provide power during the lunar night for a global network of small landers spread across the lunar surface, at latitudes ranging from equatorial to near polar landing sites. We analyze a proposal to power such small landers from orbit, using a laser to beam power from an orbital power station to photovoltaic arrays on the landers that are tuned to the laser wavelength. A commercially-available high-power 1.07-µ diode-pumped fiber laser was chosen as the source. To minimize beam spread of the a spot at the required distance, a 1.5 meter optical element was provided, using a design based on the Kepler telescope. To provide power to landers at any location, three orbital stations (“beamcraft) are required, each carrying a 3-kW laser. (If the surface science platforms requiring power are only in polar locations, or converse-ly, only in near equatorial locations, only one orbital platform is needed.) The approach is seen to be feasible, and a systems analysis was completed, the concept of operations for the system out-lined and a design for the beamcraft put together. The second design study looked at surface-to-surface power beaming using the VSAT as the la-ser platform, for an application to provide power to the interior of a permanently shadowed lunar crater from a surface platform. To maximize the distance of beaming, taking into account possible surface irregularities and the short distance to the horizon of the moon, it is desirable to emplace the laser at an elevation above the surface. The Vertical Solar Array Technology (VSAT) is a NASA program developing a solar array mounted vertically on a 10-m tall mast, designed for emplacement on a Commercial Lunar Payload Services (CLPS) lander to provide 10-kW (BOL) power near the south polar region of the moon, with a target readiness date of 2028. We used this design as the starting platform and the power source for a laser power beaming station. By mounting the laser beam director at the top of the solar array mast, a viewing distance to power receivers up to 10 km is possible. Requirements for the system were to be able to provide 300 W of continuous usable power to users including CLPS landers, VIPER class rovers, or the proposed Lunar Terrain Vehicle. The specified requirement was to be able to transmit power to a distance of up to 10 km, over a design lifetime of 5 years, and fitting within a total system landed mass under 625 kg. Again, a 1.07-µ fiber laser is mounted on the deck of the lander, with laser output sent to the laser beam director by a fiber-optic cable. A 7 square meter deployable radiator keeps the laser within operating temperature limits. The beam director is based on the design of a prototype unit developed by the University of California Santa Barbara. The system beams power for 57% of the time, with 44% of the time idle (accounting for the time when the VSAT array is itself in shadow). 1595 Watts of optical power are output in the beam. Accounting for receiver efficiency and beam losses, this results in an output onto the 1.5-meter receiving photovoltaic array of 542 watts. Of this, 300 watts is directly available to the user, while 242 watts is directed to the batteries for use while the beam is not available.

laser power↗

Stochastic Analysis for Long Term Capital Structures, Systems, and Components Refurbishment and Replacement

As commercial Nuclear Power Plants (NPPs) pursue extended plant operation in the form of Second License Renewal (SLR), opportunities exist for these plants to provide capital investments to ensure long-term safe and economic performance. At the current time, several utilities have announced an intention to pursue extended operation for one or more of their NPPs via SLR . The goal of this research is to develop a risk-informed approach to evaluate and prioritize plant capital investments made in preparation for, and during the period of, extended plant operations to support decisions for NPP operations. Since the capital investments are influenced by various factors, such as markets, safety and regulatory, the decision-making process of NPP operations should take into account relevant factors for balancing risks, costs and profits. The traditional method of capital budgeting is based on the priority list of candidate projects using economic measures such as benefit-investment ratio, net present value (NPV) and internal rate of return. In the literatures, the problem of capital budgeting or the variant can be represented by an appropriate knapsack problem. The knapsack approach to capital budgeting takes as input as investment, along with the cost and profit of each project. The objective of capital budgeting is to find the combination of the binary decisions for every investment such that the overall profit is as large as possible. The output is a collection of projects to be carried out, and we refer this selected collection of projects as a project portfolio. One limitation of traditional optimization models for capital budgeting is that they do not account for risk/uncertainty in profit and cost streams associated with individual projects, they do not account for risk in resource availability in future years [1,2,3]. Projects can incur cost over-runs, especially when projects are large, performed infrequently, and when there is risk regarding technical viability, external contractors, and/or suppliers of requisite parts and materials. Occasionally, projects are performed ahead of schedule and with cost savings. Planned budgets for capital improvements can be cut and key personnel may be lost. Or, there may be surprise windfalls in budgets for maintenance activities due to decreased costs for “unplanned” maintenance. In these cases, how should we resolve capital budgeting when we have risk forecasts for costs, profits and budgets? One approach we proposed in this summary is to re-solve the optimization models based on assumed statistical distributions of given parameters. If these distributions were not available, a two-stage stochastic optimization approach can be used to provide priority lists to decision-makers to support better risk-informed decisions [4, 5]. In this summary, we will only focus on the first approach.

42 ENGINEERING↗

Is GW170817 a Multimessenger Neutron Star-Primordial Black Hole Merger?

We investigate the possibility of the gravitational-wave event GW170817 being a light, solar-mass black hole (BH) - neutron star (NS) merger. We explore two exotic scenarios involving primordial black holes (PBH) that could produce such an event, taking into account available observational information on NGC 4993. First, we entertain the possibility of dynamical NS-PBH binary formation where a solar-mass PBH and a NS form a binary through gravitational interaction. We find that while dynamical NS-PBH formation could account for the GW170817 event, the rate is highly dependent on unknown density contrast factors and could potentially be affected by galaxy mergers. We also find that PBH-PBH binaries would likely have a larger merger rate, assuming the density contrast boost factor of an order similar to the NS-PBH case. These exotic merger formations could provide new channels to account for the volumetric rate of compact-object mergers reported by LIGO/Virgo. Secondly, we consider the case where one of the NS's in a binary NS system is imploded by a microscopic PBH. We find that the predicted rate for NS implosion into a BH is very small, at least for the specific environment of NGC 4993. We point out that similar existing (e.g. GW190425 and GW190814) and future observations will shed additional light on these scenarios.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

INTEGRATION OF FLEX EQUIPMENT AND OPERATOR ACTIONS IN PLANT FORCE-ON-FORCE MODELS WITH DYNAMIC RISK ASSESSMENT

The overall operation and maintenance cost to protect nuclear power plants accounts for approximately 7% of the total cost of power generation, with labor accounting for half of this cost. In the current research, from interaction with utilities and other stakeholders, it was determined that physical security forces account for nearly 20% of the entire workforce at several nuclear power plants. Labor costs continue to rise in the U.S., so any measures to reduce the cost of operating a nuclear power plant will need to include a reduction in labor. The physical security pathway within the DOE’s Light Water Reactor Sustainability program aims to lower the cost of physical security through directed research into modeling and simulation, application of advanced sensors or deployment of advanced weapons. This report presents a modeling and simulation framework for integrating Diverse and Flexible Mitigation Capability (FLEX) portable equipment performance with Force on Force models of a plant’s physical security posture. The generic framework is described in detail, followed by a case study of modeling an adversarial attack aimed at causing a radiological release by sabotaging the plant’s power supply and its ultimate heat sink capabilities at a hypothetical nuclear power plant. Two different FLEX deployment strategies, series and parallel, are modeled with distinct timelines. The results of the adversarial attack modeled in a commercial Force on Force tool are integrated with the FLEX deployment model in INL’s dynamic modeling tool EMRALD. Monte Carlo simulation is used to model the distribution of the timeline in FLEX deployment strategies. The results demonstrate that, even in the extreme case of a successful adversarial attack, deployment of FLEX equipment can result in a significantly high likelihood of preventing radiological release. The modeling and simulation framework integrating FLEX equipment with Force on Force models enables the nuclear power plants to credit FLEX portable equipment in the plant security posture, resulting in an efficient and optimized physical security.

97 MATHEMATICS AND COMPUTING↗

Demonstrate Improved Ag Diffusion and Describe the Basis for Pd Penetration Modeling in SiC

In past work, an effective diffusivity coefficient was determined for Ag transport through the silicon carbide layer of a tristructural isotropic fuel particle. The effective diffusivity coefficient accounts for the microstructure of the silicon carbide and includes both bulk diffusion and grain boundary diffusion of Ag. In this report, the model has been improved by accounting for the enhanced concentration of vacancies in the bulk due to irradiation, which substantially influence bulk diffusivity at low temperatures. To improve the BISON model and make it fission rate dependent, effective diffusivity calculations have been performed that incorporate the radiation modified bulk diffusivity. The microstructure and irradiation-dependent effective diffusivity has also been implemented into BISON, and its predictions for Ag release from tristructural isotropic fuel have been successfully compared to AGR-1 post irradiation measurements. Moreover, a new feature has been developed in the Multiphysics Object-Oriented Simulation Environment (MOOSE) to account for different grain boundary types. The Ag diffusivity in 5 (210)/[001] grain boundaries has been computed and was found to be greater than in random high-angle grain boundaries. The presence of the fission product Pd can also have an important effect on the properties of the silicon carbide layer in tristructural isotropic particles. The penetration of Pd into the silicon carbide layer causes a corrosion reaction that can lead to the failure of the silicon carbide layer; however, this corrosion reaction is not well understood. To enable an improved understanding of the mechanism, ab-initio molecular dynamics simulations of Pd interaction with bulk silicon carbide have been performed. The improved understanding of the reaction will form a basis for future improvements to the BISON’s Pd penetration failure model

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Equitable Employment Access Assessed Through the Mobility Energy Productivity (MEP) Metric

This paper examines commuting options for an underserved neighborhood in Columbus, Ohio to a major employment center. The analysis is based on an emerging metric called the Mobility Energy Productivity (MEP) metric developed by the National Renewable Energy Laboratory (NREL) on behalf the Department of Energy (DOE). The purpose of the analysis is twofold. The first is to quantify relative attractiveness of commute modes between the two locations, using a perspective that includes travel time, energy and cost, while providing an equity lens to compare commute options between privately owned vehicles and pooled transportation options. The second objective is to apply MEP in a specific origin-destination (O-D) scenario, whereas previously it has been used primarily as an aggregate metropolitan-wide statistical measure. In so doing, parameters in MEP are further customized and the methodology is refined to account for unique aspects of this case study. Four commute options between the neighborhood and the industry employment based are analyzed: drive alone option, public transit express bus (historical), public transit normal route (current), and a proposed shuttle specific to the O-D pair. This analysis identified issues applying MEP that required further customization: (1) deprecation functions customized to modes other than driving, (2) accounting for first-mile last-mile travel times with transit, (3) accounting for transit frequency without resorting to full simulation. The results provide quantitative insights on the employment accessibility between these two locations, both across modes, and as equity of job accessibility for those who can and cannot operate a personal vehicle.

ADVANCED PROPULSION SYSTEMS↗

Efficient Modeling of Floating Wind Arrays Including Current Loads and Seabed Bathymetry: Preprint

Capabilities for modeling the effects of seabed bathymetry and current drag loads on a floating wind farm are now available in an open-source model for quasi-static analysis. In this model, mooring lines and dynamic cables are represented by a quasistatic solver that can quickly represent complex mooring/cabling arrangements and arrays of floating bodies. To account for seabed bathymetry, we expand the model to include a surface mesh that captures changes in water depth over a rectangular grid. We formulate modifications to the catenary equations that capture a mooring line's profile and tensions when contacting a slope seabed. To account for current drag loads on mooring lines and dynamic cables, we formulate a novel technique that rotates the reference frame so that the vector sum of the the weight and the current force are used in the catenary equations, while accounting for the seabed orientation. To complete the system, current drag loads on floating substructures are handled by inclusion of strip-theory drag calculations. These new capabilities are verified by comparing with results from the established offshore dynamics models MoorDyn and OrcaFlex in equivalent steady-state scenarios. The results show very good agreement for both sloped seabeds and current loads. With computation times of the quasi-static model typically under one second, the model additions are a useful capability toward rapidly evaluating a floating wind array's response to environmental loads under realistic site conditions.

bathymetry↗

A Survey of Multizone Airflow Simulation Software

H-Canyon simulates the building and process ventilation to ensure positive contamination control for any modifications as well as evaluating off-normal conditions. Simulations can also be used to identify air quality of personnel occupied areas. The models can be classified into two groups, single and multi-zone models. As the name states, a single zone model assumes that the building can be described by a single, well-mixed zone. Whereas large buildings have structures that must be described as multizone. Multizone airflow models deal with the difficulty of flows in a building by distinguishing the effects of internal flow restrictions. Several models have been developed that take internal partitions into account. CONTAM (currently used) is a multizone airflow program created by NIST to determine airflows including infiltration, exfiltration, wind pressures acting outside the building, buoyancy effects caused by the indoor and outdoor air temperature difference. COMIS is a multizone air infiltration, ventilation, and contaminant transport model that was developed at Lawrence Berkeley National Laboratory. COMIS deals with applied pressure distribution, airflow components, HVAC systems, schedules, and contaminant transport. EnergyPlus is designed by DOE. It supports flexible component level configuration of HVAC models and simulates sub-hourly time steps to control fast dynamic systems. The program features include providing accurate environmental performance data including energy consumption and carbon emissions, and reporting solar gains on surfaces. Heavent is used to model industrial ventilation systems to control airborne contaminants. It's designed based on the American Conference of Governmental Industrial Hygienists (ACGIH) Ventilation manual. DesignBuilder is the most recognized and advanced user interface to EnergyPlus. It gives access to all of the most frequently required simulation capabilities covering building fabric, thermal mass, glazing, shading renewables, and HVAC analysis. The objective of this study is to identify all relevant software available including commercially licensed, freely licensed, and open source, that can be used to simulate large, complex process facility ventilation systems. All discovered software tools are evaluated for salient features and compared to the software tool currently used to simulate the H-Canyon ventilation system. The design criteria of a ventilation system has two basic requirements: Supplying fresh air for the occupants and changing the air in the room sufficiently; therefore, the smells, fumes, and contaminants will be removed. In this study, first, whether the software simulates airflow (infiltration, natural ventilation, and interzonal) or considers them to be inputs is discussed. Second, if indoor contaminant levels are considered; and finally suggestions on how airflow analysis can be more appropriately accounted for when performing airflow simulation of a building. All of the airflow software reviewed for this study can account for constant infiltration and natural ventilation rates that are not affected by changes in indoor and outdoor conditions. Modelling infiltration in EnergyPlus or DesignBuilder are applicable for small buildings and not appropriate for larger buildings like the H-Canyon. The effect of wind on external pressures, infiltration, and natural and mechanical ventilation, can also be calculated using all of these software. These programs have control options to switch between natural to mechanical ventilation modes or both (hybrid). Interzonal airflow can be applied which means one zone can be specified to receive a fixed amount of air from another zone. The following ventilation modeling software is recommended for further testing: CONTAM, COMIS, Heavent. The next step in evaluating the software will consist of creating test facility simulations using evaluation copies of the software.

42 ENGINEERING↗