Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Average model”

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

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

At least 1,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

The Use of CASES-97 Observations to Assess and Parameterize the Impact of Land-Surface Heterogeneity on Area-Average Surface Heat Fluxes for Large-Scale Coupled Atmosphere-Hydrology Models

To understand the effects of land-surface heterogeneity and the interactions between the land-surface and the planetary boundary layer at different scales, we develop a multiscale data set. This data set, based on the Cooperative Atmosphere-Surface Exchange Study (CASES97) observations, includes atmospheric, surface, and sub-surface observations obtained from a dense observation network covering a large region on the order of 100 km. We use this data set to drive three land-surface models (LSMs) to generate multi-scale (with three resolutions of 1, 5, and 10 kilometers) gridded surface heat flux maps for the CASES area. Upon validating these flux maps with measurements from surface station and aircraft, we utilize them to investigate several approaches for estimating the area-integrated surface heat flux for the CASES97 domain of 71x74 square kilometers, which is crucial for land surface model development/validation and area water and energy budget studies. This research is aimed at understanding the relative contribution of random turbulence versus organized mesoscale circulations to the area-integrated surface flux at the scale of 100 kilometers, and identifying the most important effective parameters for characterizing the subgrid-scale variability for large-scale atmosphere-hydrology models.

Chen, Fei↗
NASA NTRS2012

A HTAP Multi-Model Assessment of the Influence of Regional Anthropogenic Emission Reductions on Aerosol Direct Radiative Forcing and the Role of Intercontinental Transport

In this study, we assess changes of aerosol optical depth (AOD) and direct radiative forcing (DRF) in response to the reduction of anthropogenic emissions in four major pollution regions in the northern hemisphere by using results from 10 global chemical transport models in the framework of the Hemispheric Transport of Air Pollution (HTAP). The multi-model results show that on average, a 20% reduction of anthropogenic emissions in North America, Europe, East Asia and South Asia lowers the global mean AOD and DRF by about 9%, 4%, and 10% for sulfate, organic matter, and black carbon aerosol, respectively. The impacts of the regional emission reductions on AOD and DRF extend well beyond the source regions because of intercontinental transport. On an annual basis, intercontinental transport accounts for 10-30% of the overall AOD and DRF in a receptor region, with domestic emissions accounting for the remainder, depending on regions and species. While South Asia is most influenced by import of sulfate aerosol from Europe, North America is most influenced by import of black carbon from East Asia. Results show a large spread among models, highlighting the need to improve aerosol processes in models and evaluate and constrain models with observations.

Yu, Hongbin↗
NASA NTRS1977

Faulting parameters derived from computer simulation of earthquakes

Seismic source parameters, average displacement, rupture length, and strain energy release are investigated by computer simulation using a coupled massive block model of the sliding along an active fault. Average displacements and energy release vary considerably with the degree of heterogeneity in the friction and elastic parameters used in the model. Strain energy release is determined primarily by the product of dynamic friction, rupture length, and average displacement. Interrelationships among the faulting parameters are consistent with theoretical arguments and experimental data. The variation in the frequency of occurrence of simulation events with strain energy release is different from the variation in the frequency of naturally occurring events with seismic energy.

Cohen, S. C.↗
DOE OSTI2022

Sorghum Association Panel whole‐genome sequencing establishes cornerstone resource for dissecting genomic diversity

SUMMARY Association mapping panels represent foundational resources for understanding the genetic basis of phenotypic diversity and serve to advance plant breeding by exploring genetic variation across diverse accessions. We report the whole‐genome sequencing (WGS) of 400 sorghum ( Sorghum bicolor (L.) Moench) accessions from the Sorghum Association Panel (SAP) at an average coverage of 38× (25–72×), enabling the development of a high‐density genomic marker set of 43 983 694 variants including single‐nucleotide polymorphisms (approximately 38 million), insertions/deletions (indels) (approximately 5 million), and copy number variants (CNVs) (approximately 170 000). We observe slightly more deletions among indels and a much higher prevalence of deletions among CNVs compared to insertions. This new marker set enabled the identification of several novel putative genomic associations for plant height and tannin content, which were not identified when using previous lower‐density marker sets. WGS identified and scored variants in 5‐kb bins where available genotyping‐by‐sequencing (GBS) data captured no variants, with half of all bins in the genome falling into this category. The predictive ability of genomic best unbiased linear predictor (GBLUP) models was increased by an average of 30% by using WGS markers rather than GBS markers. We identified 18 selection peaks across subpopulations that formed due to evolutionary divergence during domestication, and we found six F st peaks resulting from comparisons between converted lines and breeding lines within the SAP that were distinct from the peaks associated with historic selection. This population has served and continues to serve as a significant public resource for sorghum research and demonstrates the value of improving upon existing genomic resources.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2009

Lunar Pole Illumination and Communications Maps Computed from GSSR Elevation Data

A Digital Elevation Model of the lunar south pole was produced using Goldstone Solar System RADAR (GSSR) data obtained in 2006.12 This model has 40-meter horizontal resolution and about 5-meter relative vertical accuracy. This Digital Elevation Model was used to compute average solar illumination and Earth visibility with 100 kilometers of the lunar south pole. The elevation data were converted into local terrain horizon masks, then converted into lunar-centric latitude and longitude coordinates. The horizon masks were compared to latitude, longitude regions bounding the maximum Sun and Earth motions relative to the moon. Estimates of Earth visibility were computed by integrating the area of the region bounding the Earth's motion that was below the horizon mask. Solar illumination and other metrics were computed similarly. Proposed lunar south pole base sites were examined in detail, with the best site showing yearly solar power availability of 92 percent and Direct-To-Earth (DTE) communication availability of about 50 percent. Similar analysis of the lunar south pole used an older GSSR Digital Elevation Model with 600-meter horizontal resolution. The paper also explores using a heliostat to reduce the photovoltaic power system mass and complexity.

lunar base↗
DOE OSTI2023

Adaptive Ensemble Refinement of Protein Structures in High Resolution Electron Microscopy Density Maps with Radical Augmented Molecular Dynamics Flexible Fitting

Recent advances in cryo-electron microscopy (cryo-EM) have enabled modeling macromolecular complexes that are essential components of the cellular machinery. The density maps derived from cryo-EM experiments are often integrated with manual, knowledge or artificial intelligence driven, and physics-guided computational methods to build, fit, and refine molecular structures. Going beyond a single stationary- structure determination scheme, it is becoming more common to interpret the experimental data with an ensemble of models, which contributes to an average observation. Hence, there is a need to decide on the quality of an ensemble of protein structures on-the-fly, while refining them against the density maps. Here, we introduce such an adaptive decision making scheme during the molecular dynamics flexible fitting (MDFF) of biomolecules. Using RADICAL-Cybertools, and the new RADICAL augmented MDFF implementation (R-MDFF) is examined in high-performance computing environments for refinement of two protein systems, Adenylate Kinase and Carbon Monoxide Dehydrogenase. For the test cases, use of multiple replicas in flexible fitting with adaptive decision making in R-MDFF improves the overall correlation to the density by 40% relative to the refinements of the brute-force MDFF. The improvements are particularly significant at high, 2 - 3 Å, map resolutions. More importantly, the ensemble model captures key features of biologically relevant molecular dynamics that is inaccessible to a single-model interpretation. Finally, the pipeline is applicable to systems of growing sizes, which is demonstrated using ensemble refinement of capsid proteins from Chimpanzee adenovirus. The overhead for decision making remaining low and robust to computing environments. The software is publicly available on GitHub and includes a short user guide to install the R-MDFF on different computing environments, from local Linux based workstations to High Performance Computing (HPC) environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1986

New models of Saturn's magnetic field using Pioneer 11 Vector Helium Magnetometer data

In a reanalysis of the Vector Helium Magnetometer data taken by Pioneer 11 during its Saturn encounter in 1979, using improvements in the data set and in the procedures, studies are made of a variety of models. The best is the P(11)84 model, an axisymmetric spherical harmonic model of Saturn's magnetic field within 8 Saturn radii of the planet. The appropriately weighted root mean square average of the difference between the observed and the modeled field is 1.13 percent. For the Voyager-based Z3 model of Connerney, Acuna, and Ness, this average difference from the Pioneer 11 data is 1.81 percent. The external source currents in the magnetopause, tail, bow shock, and perhaps ring currents vary with time and can only be crudely modeled. An algebraic formula is derived for calculating the L shells on which energetic charged particles drift in axisymmetric fields.

Davis, L., Jr.↗
DOE OSTI2026

Progressing Analysis of Variable Electric Rates (PAVER) Study

The Progressing Analysis of Variable Electric Rates (PAVER) study analyzed the impact of a range of time-varying electric rates on the performance of a regional electric grid and the resulting costs for participating and non-participating customers. This analysis leveraged and extended the work of PNNL’s Distribution System Operator with Transactive (DSO+T) study. Five different rate designs were included: a flat volumetric energy charge, a typical Time of Use (TOU) rate, a dynamic energy (DE) rate (based on wholesale locational marginal prices), a dynamic energy and capacity (DE+C) rate, and, finally, a Block and Swing (B&S) rate that billed customers based on their average load profile at constant pricing, but used the DE+C dynamic price for load deviations from their average profile. These rates were analyzed in a large-scale co-simulation of an entire regional grid with a customer population representative of the current state. A large fraction (80%) of residential and commercial customers were assumed to participate in these time-varying rates with automatically controlled HVAC, water heaters, electric vehicles, and batteries. This study assumed no industrial sector participation. The DE and DE+C rates saw system peak loads reduced by 6-7%, while the large participation in the TOU rate case saw a significant rebound effect and a resulting peak load increase of >5%. The impacts to the annual and peak system demand impacted system wholesale prices and the overall grid operating costs. This cost structure determined the revenue needed to be collected from customers by each rate design. Participating customers on the DE and DE+C rates (located in one of the modeled DSOs) saw reductions in average annual electricity bills of 11-17% with average increases in monthly bill variation of no more than 13%. At such high participation levels, TOU customers saw 10% higher average annual bills (due to system-wide rebound effects) and average increased monthly bill variation of 16%. Residential owners of large flexible loads (such as electric vehicles) saw larger bill savings (17-20%) when on a fully dynamic rate. The presence of on-site generation (such as rooftop solar) did not appear to appreciably change customer outcomes. Customers on the Block and Swing rate did see 6% lower monthly bill variation (as intended) than the flat rate case, but at the expense of appreciable bill savings, which were only 3%, comparable to the savings seen by non-participants. Given this finding we recommend that additional research be conducted into how best various bill protection mechanisms can balance minimizing customer bill variation with providing financial incentives commensurate with the flexibility customers provide. We also recommend that customer outcomes be explored across a range of regions using current actual customer and system cost data.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

A passive model for the evolution of subgrid-scale instabilities in turbulent flow regimes

In this work, a zero-dimensional model for the development of turbulent mixing layers driven by Rayleigh–Taylor (RT) or Richtmyer–Meshkov (RM) instabilities is presented, for application in Reynolds-Averaged Navier–Stokes (RANS) simulations. The model approximates the behavior of the turbulent mixing layer at a passively advected tracer particle that sits on the material interface of interest in the simulation, and is used to provide initial conditions for the RANS model once the mixing layer width has developed to a size that may be resolved on the computational mesh. The model relies on approximations to the Besnard–Harlow–Rauenzahn (BHR) RANS closure equations and is thus consistent the turbulent development of the mixing layer, unlike previous models of this type which have focused on the early-time laminar development of RT or RM instabilities. These instabilities will often grow through an initially laminar state before developing into turbulence, and criteria for switching from an existing laminar model to the proposed turbulent model are considered. The proposed model reproduces the mixing layer development of RT and RM instabilities seen in mixed width resolving BHR RANS simulations to a reasonable degree of accuracy, and relaxes the need for the RANS simulation to resolve the mixed width as soon as the mixing layer becomes turbulent, which may be impractical in cases where the mixing layer grows through a wide range of scales.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Bhatnagar-Gross-Krook subgrid model for neutrino quantum kinetics

Here, we present a new subgrid model for neutrino quantum kinetics, which is primarily designed to incorporate effects of collective neutrino oscillations into neutrino-radiation-hydrodynamic simulations for core-collapse supernovae and mergers of compact objects. We approximate the neutrino oscillation term in a quantum kinetic equation using the Bhatnagar-Gross-Krook (BGK) relaxation-time prescription, and the transport equation is directly applicable for classical neutrino transport schemes. The BGK model is motivated by recent theoretical indications that nonlinear phases of collective neutrino oscillations settle into quasisteady structures. We explicitly provide basic equations of the BGK subgrid model for both multiangle and moment-based neutrino transport to facilitate the implementation of the subgrid model in the existing neutrino transport schemes. We also show the capability of our BGK subgrid model by comparing it to fully quantum kinetic simulations for fast neutrino-flavor conversion. We find that the overall properties can be well reproduced in the subgrid model; the error of angular-averaged survival probability of neutrinos is within ~20%. By identifying the source of error, we also discuss perspectives to improve the accuracy of the subgrid model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Emergency Asset Positioning for Resilient Transmission Grid Operation

In this paper we study a problem of hurricane emergency preparedness via placement of emergency storage assets prior to its strike. We present a two-stage stochastic model for choosing locations and quantities of emergency storage to help support the power grid through a hurricane event. The expectation of losses in our two-stage model is estimated using the sample average approximation. We construct damage scenarios for sample average approximation using WIND Toolkit meteorological data and fragility curves of various electric grid components. We demonstrate the efficacy of our two-stage planning model by simulating operations during Hurricane Dolly on the 2000-bus transmission test system. Our model, coupled with our scenario selection strategy, is effective at mitigating loss of load when compared to a model without emergency storage assets placed prior to an extreme event.

emergency preparedness↗
NASA NTRS1995

Simulated interannual variability of the Greenland Sea deep water formation and its connection to surface forcing

A fully prognostic Arctic ice-ocean model is used to study the interannual variability of deepwater formation in the Greenland Sea Gyre based on the simulations for the Arctic ice-ocean system for the period 1955 and 1960 - 1985. The model uses monthly climatology for thermodynamic forcing components (such as air temperature and cloudiness), together with constant annual net precipitation and river runoff. The daily wind forcing is derived from analyzed sea level air pressures from the National Center for Atmospheric Research (NCAR). In summary, the model shows that the occurence of deep convection in the Greenland Sea Gyre is controlled by the extensive Fram Strait ice export and/or local wind conditions in the Greenland Sea. In the latter case the weakening of the local wind curl allows the Polar Front to move eastward. The movement of the Polar Front causes adverse ice conditions, often together with much larger than normal ice export from the Arctic, such as in 1968, which can block convection in the gyre. The density difference between upper and lower layers is investigated as an indication of water mass formation through convection, occurring as strong diffusion in the model. The model-simulated density difference between the average top 100 m and deep levels reveals that the period 1960 - 1985 had only a few distinct years with weak stratification, and, especially, the model predicts no deep convection since the nid-1970s. The common factor for the years of the weakest decrease of the model-predicted heat content of the upper 2000 m which can, to a high degree, be explained by local heat loss.

Haekkinen, Sirpa↗
DOE OSTI2023

Emergency Asset Positioning for Resilient Transmission Grid Operation

In this paper we study a problem of hurricane emergency preparedness via placement of emergency generation assets prior to its strike. We present a two-stage stochastic model for choosing locations and quantities of emergency storage to help support the power grid through a hurricane event. The expectation of losses in our two-stage model is estimated using the sample average approximation. We construct damage scenarios for sample average approximation using WIND Toolkit meteorological data and fragility curves of various electric grid components. We demonstrate the efficacy of our two-stage planning model by simulating operations during Hurricane Dolly on the 2000-bus transmission test system. Our model, coupled with our scenario selection strategy, is effective at mitigating loss of load when compared to a model without emergency generation assets placed prior to an extreme event.

emergency preparedness↗
DOE OSTI2023

Emergency Generation Asset Positioning for Resilient Transmission Grid Operation: Preprint

In this paper we study a problem of hurricane emergency preparedness via placement of emergency generation assets prior to its strike. We present a two-stage stochastic model for choosing locations and quantities of emergency storage to help support the power grid through a hurricane event. The expectation of losses in our two-stage model is estimated using the sample average approximation. We construct damage scenarios for sample average approximation using WIND Toolkit meteorological data and fragility curves of various electric grid components. We demonstrate the efficacy of our two-stage planning model by simulating operations during Hurricane Dolly on the 2000-bus transmission test system. Our model, coupled with our scenario selection strategy, is effective at mitigating loss of load when compared to a model without emergency generation assets placed prior to an extreme event.

emergency preparedness↗
NASA NTRS2022

VULCAN-CFD User Manual: Ver. 7.2.0

VULCAN-CFD offers a comprehensive set of capabilities to enable the simulation of continuum flowfields from subsonic to hypersonic conditions. The governing equations that are employed include allowances for both chemical and thermal nonequilibrium processes, coupled with a wide variety of turbulence models for both Reynolds-averaged and large eddy simulations. The software package can simulate two-dimensional, axisymmetric, or three-dimensional problems on structured multiblock meshes or fully unstructured meshes. A parabolic (i.e., space-marching) treatment can also be used for any subset of a structured mesh that can accommodate this solution strategy. The flow solver provides a significant level of geometric flexibility for structured grid simulations by allowing for arbitrary face-to-face C(0) continuous and non-C(0) continuous block interface connectivities. The unstructured grid paradigm allows for mixed element unstructured meshes that contain any combination of tetrahedral, prismatic, pyramidal, and hexahedral cell elements. The flow solver is also fully parallelized using MPI (Message Passing Interface) libraries in a data-parallel fashion, allowing for efficient simulations on modern High Performance Computing (HPC) systems. This document provides information related to the installation and execution of the VULCAN-CFD software package. A detailed description of the physical and numerical models available in the software are provided in the VULCAN-CFD Theory Manual.

VULCAN-CFD User Manual↗
NASA NTRS1994

Large-eddy simulation of flow through a plane, asymmetric diffuser

A challenge for traditional turbulence modeling, based on the Reynolds averaged Navier-Stokes equations, remains the accurate prediction of 'mild', adverse pressure-gradient driven separation from a smooth surface. With this study we want to explore the capability of large-eddy simulation to predict the separation which occurs on the deflected wall of an asymmetric, plane diffuser with opening angle of 10 deg. The flow through the plane diffuser exhibits some additional interesting physical phenomena which make it a challenging test case. In addition to 'mild' separation about halfway down the deflected ramp, the flow is characterized by a small backflow zone with stalled fluid in the rear part of the expanding section. The turbulent flow entering the diffuser is subject to combined adverse and radial pressure gradients stemming from the convex curvature. Finally the flow recovers into a developed, turbulent channel flow in the outlet section. Obi et al. provide measurements of mean flow, Reynolds stresses, and pressure recovery, which were obtained by means of LDV in a wind tunnel. The objective of this study is to investigate whether LES with the standard dynamic model is able to accurately predict the flow in the one-sided diffuser and to explore the resolution requirements and associated costs.

Kaltenbach, Hans-Jakob↗
NASA NTRS2005

Effects of High-frequency Wind Sampling on Simulated Mixed Layer Depth and Upper Ocean Temperature

Effects of high-frequency wind sampling on a near-global ocean model are studied by forcing the model with a 12 hourly averaged wind product and its 24 hourly subsamples in separate experiments. The differences in mixed layer depth and sea surface temperature resulting from these experiments are examined, and the underlying physical processes are investigated. The 24 hourly subsampling not only reduces the high-frequency variability of the wind but also affects the annual mean wind because of aliasing. While the former effect largely impacts mid- to high-latitude oceans, the latter primarily affects tropical and coastal oceans. At mid- to high-latitude regions the subsampled wind results in a shallower mixed layer and higher sea surface temperature because of reduced vertical mixing associated with weaker high-frequency wind. In tropical and coastal regions, however, the change in upper ocean structure due to the wind subsampling is primarily caused by the difference in advection resulting from aliased annual mean wind, which varies with the subsampling time. The results of the study indicate a need for more frequent sampling of satellite wind measurement and have implications for data assimilation in terms of identifying the nature of model errors.

mixed layered depth↗
NASA NTRSDate not supplied

Examination of L-Band Brightness Temperature Forecasts in the SMAP Level-4 Soil Moisture Analysis

The NASA Soil Moisture Active Passive (SMAP) mission [4] has been providing L-band (1.4 GHz) passive microwave brightness temperature (Tb) observations since April 2015. By assimilating the Tb observations into the NASA Catchment land surface model [5] using a spatially distributed ensemble Kalman filter, the NASA Global Modeling and Assimilation Office generates the SMAP Level-4 Soil Moisture (L4_SM) product, which provides global, 3-hourly, 9-km resolution estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture with ~2.5-day latency for use in research and applications [6]. The L4_SM product also includes estimates of soil temperature, land surface fluxes, and assimilation diagnostics such as the model forecast and observed Tb values [7]. The output from the L4_SM system is routinely monitored by the L4_SM team. Such monitoring provides valuable information; instances of unusually large Tb observation-minus-forecast (O-F) residuals can indicate events during which soil moisture conditions are poorly described in the land modeling system [7,8]. For example, repeated occurrences of very large Tb O-F values in central Australia were traced back to deficiencies in the gauge-based precipitation product used in the land modeling system through L4_SM Version 5; this discovery prompted the use, beginning in Version 6, of satellite- and gauge-based precipitation observations outside of North America [9]. Recently, a statistical analysis of the Tb O-F residuals revealed systematic errors in the “tau-omega” L-band radiative transfer model that converts the land-model simulated soil moisture and temperature into the forecast Tb prior to the L4_SM analysis. The L4_SM analysis is built on the ensemble Kalman filter and assumes unbiased forecast errors. The presence of systematic Tb forecast errors could thus adversely impact the quality of the analyzed soil moisture. In this paper, we examine the Tb O-F residuals of the latest Version 7 L4_SM data (Science Version ID Vv7030 and Vv7032) [10,11,12]. In the L4_SM Version 7 algorithm, key parameters of the L-band radiative transfer model, including soil roughness, scattering albedo, and a (seasonally varying) climatology of vegetation opacity, are obtained from the SMAP Level-2 Radiometer retrieval product (Version 5) [1]. Additionally, we also use ground measurements of surface soil moisture and soil temperature from the SMAP core validation sites as in situ reference of soil conditions [2,3]. As part of the L4_SM system calibration, the seasonally-varying bias between the model forecast Tb and the observed values is removed prior to the assimilation of the SMAP Tb observations [7]. That is, the L4_SM system is designed to only correct errors in synoptic-scale and interannual variations from the long-term mean seasonal cycle while maintaining the model’s (potentially erroneous) climatology. The Tb O-F residuals examined here are thus computed using SMAP Tb observations after they are rescaled to the mean seasonal cycle of Tb from the modeling system without SMAP data assimilation. Consequently, the long-term average of the Tb O-F residuals has a global mean of only 0.13 K and locally small values, ranging from -1 to 3 K. Despite the small time-average values of the Tb O-F residuals, the model forecast Tb was nevertheless found to exhibit undesirable systematic errors. At some locations, the time-average Tb O-F values strongly depend on surface soil moisture (SM). At the Yanco SMAP core validation site, for example, the Tb O-F residuals typically range from 5 to 15 K under dry soil moisture conditions (SM < 0.15 m3 m-3) yet are predominantly negative under wet soil moisture conditions (SM > 0.25 m3 m-3), with values ranging from 0 to -40 K. This results in soil moisture analysis increments that persistently make the soil drier under dry SM conditions and persistently make the soil wetter under wet SM conditions, suggesting an error in the dynamic range of the simulated Tb, soil moisture or soil temperature. In this paper, we describe the higher-order systematic Tb forecast errors in more detail, examine their impact on the L4_SM product quality, and explore potential avenues to further improve the L4_SM algorithm.

Rolf Reichle↗