Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial representation”

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 91 records · Page 5

A high order Cartesian grid, finite volume method for elliptic interface problems

We present a higher-order finite volume method for solving elliptic PDEs with jump conditions on interfaces embedded in a 2D Cartesian grid. Second, fourth, and sixth order accuracy is demonstrated on a variety of tests including problems with high-contrast and spatially varying coefficients, large discontinuities in the source term, and complex interface geometries. We include a generalized truncation error analysis based on cell-centered Taylor series expansions, which then define stencils in terms of local discrete solution data and geometric information. In the process, we develop a simple method based on Green's theorem for computing exact geometric moments directly from an implicit function definition of the embedded interface. This approach produces stencils with a simple bilinear representation, where spatially-varying coefficients and jump conditions can be easily included and finite volume conservation can be enforced.

97 MATHEMATICS AND COMPUTING↗

Evaluation of the Ozone Fields in NASA's MERRA-2 Reanalysis

The assimilated ozone product from the Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2), produced at NASAs Global Modeling and Assimilation Office (GMAO) is summarized. The reanalysis begins in 1980 with the use of retrieved partial-column ozone concentrations from a series of Solar Backscatter Ultraviolet Radiometer (SBUV) instruments on NASA and NOAA spacecraft. Beginning in October 2004, retrieved ozone profiles from the Microwave Limb Sounder (MLS) and total column ozone from the Ozone Monitoring Instrument (OMI) on NASAs EOS Aura satellite are assimilated. While this change in data streams does lead to a discontinuity in the assimilated ozone fields in MERRA-2, making it not useful for studies in decadal (secular) trends in ozone, this choice was made to prioritize demonstrating the value NASAs high-quality research data in the reanalysis context. The MERRA-2 ozone is compared with independent satellite and ozonesonde data, focusing on the representation of the spatial and temporal variability of stratospheric and upper-tropospheric ozone. The comparisons show agreement within 10 (standard deviation of the difference) between MERRA-2 profiles and independent satellite data in most of the stratosphere. The agreement improves after 2004, when EOS Aura data are assimilated. The standard deviation of the differences between the lower-stratospheric and upper-tropospheric MERRA-2 ozone and ozonesondes is 11.2 and 24.5, respectively, with correlations of 0.8 and above. This is indicative of a realistic representation of the UTLS ozone variability in MERRA-2. After 2004, the upper tropospheric ozone in MERRA-2 shows a low bias compared to the sondes, but the covariance with independent observations is improved compared to earlier years. Case studies demonstrate the integrity of MERRA-2 analyses in representing important features such as tropopause folds.

Ozone↗

Evaluating the Impact of Model Resolutions and Cumulus Parameterization on Precipitation in NU-WRF: A Case Study in the Central Great Plains

Regional climate models are expected to exhibit improved skill at finer spatial resolutions due to improved representation of land surface heterogeneity. However, at spatial scales between 1 to 10 km (grey scales), these improvements are often illusive due to the competing benefits from spatial resolution and cumulus parameterization. This study provides insights into the impact of model resolution and cumulus parameterization on precipitation prediction in the Central Great Plains by using an object-based evaluation method. Our results show limited improvement solely from finer resolution but larger improvement without using the cumulus scheme at a 4-kmresolution. Compared to traditional evaluation methods, the object-based analysis shows that without the cumulus scheme the spatial properties of precipitation are better represented. In contrast, all model configurations show a dry bias in precipitation days and a tendency to produce widespread precipitation but with fewer hours with precipitationwhich indicates other shortcomings in the model.

Yuqi Zhang↗

MFA-DVR: direct volume rendering of MFA models

3D volume rendering is widely used to reveal insightful intrinsic patterns of volumetric datasets across many domains. However, the complex structures and varying scales of volumetric data can make efficiently generating high-quality volume rendering results a challenging task. Multivariate functional approximation (MFA) is a new data model that addresses some of the critical challenges: high-order evaluation of both value and derivative anywhere in the spatial domain, compact representation for large-scale volumetric data, and uniform representation of both structured and unstructured data. Here, in this paper, we present MFA-DVR, the first direct volume rendering pipeline utilizing the MFA model, for both structured and unstructured volumetric datasets. We demonstrate improved rendering quality using MFA-DVR on both synthetic and real datasets through a comparative study. We show that MFA-DVR not only generates more faithful volume rendering than using local filters but also performs faster on high-order interpolations on structured and unstructured datasets. MFA-DVR is implemented in the existing volume rendering pipeline of the Visualization Toolkit (VTK) to be accessible by the scientific visualization community.

97 MATHEMATICS AND COMPUTING↗

Acclimation of Phytoplankton Fe:C Ratios Dampens the Biogeochemical Response to Varying Atmospheric Deposition of Soluble Iron

Abstract Dissolved iron (dFe) plays an important role in regulating marine productivity. In high nutrient, low chlorophyll regions (>33% of the global ocean), iron is the primary growth limiting nutrient, and elsewhere iron can regulate nitrogen fixation by diazotrophs. The link between iron availability and carbon export is strongly dependent on the phytoplankton iron quotas or cellular Fe:C ratios. This ratio varies by more than an order of magnitude in the open ocean and is positively correlated with ambient dFe concentrations in field observations. Representing Fe:C ratios within models is necessary to investigate how ocean carbon cycling will interact with perturbations to iron cycling in a changing climate. The Community Earth System Model ocean component was modified to simulate dynamic, group‐specific, phytoplankton Fe:C that varies as a function of ambient iron concentration. The simulated Fe:C ratios improve the representation of the spatial trends in the observed Fe:C ratios. The acclimation of phytoplankton Fe:C ratios dampens the biogeochemical response to varying atmospheric deposition of soluble iron, compared to a fixed Fe:C ratio. However, varying atmospheric soluble iron supply has first order impacts on global carbon and nitrogen fluxes and on nutrient limitation spatial patterns. Our results suggest that pyrogenic Fe is a significant dFe source that rivals mineral dust inputs in some regions. Changes in dust flux and iron combustion sources (anthropogenic and wildfires) will modify atmospheric Fe inputs in the future. Accounting for dynamic phytoplankton iron quotas is critical for understanding ocean biogeochemistry and projecting its response to variations in atmospheric deposition.

Environmental Sciences & Ecology↗

City-Wide Distributed Roof-Top Photovoltaic System Adoption Forecast, Grid Impact Simulation, & Neighborhood Microgrid Contribution Assessment

The adoption of distributed photovoltaic (PV) systems grew significantly in recent years. Market projections anticipate future growth for both residential and commercial installations. To understand grid impacts associated with distributed PV, useful hosting capacity studies require accurate representations of the spatial distribution of PV adoptions. Prediction of PV locations and numbers depends on median income data, building use zoning maps, and permit records to understand existing trends and predict future adoption rates and locations throughout an entire city. Using the PV adoption data, advanced and realistic simulations were performed to capture the distributed PV impacts on the grid. Also, using graph theory community detection hundreds of neighborhood microgrids can be discovered for the entire city by identifying densely connected loads that are sparsely connected to other communities. Then, based on the PV adoption predictions, this work identified the contribution of PV within each of the newly discovered graph theory defined microgrid communities.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Nonlinear structural dynamics analysis using a modified modal method

The procedure for predicting the nonlinear dynamic response of structural components subjected to a step loading is presented. The procedure is a modified modal method that involves a change of dependent variables from the unknown nodal degrees of freedom of the finite element model of the structure to a smaller set of generalized coordinates. This change of dependent variables uses a combination of the nonlinear static solution and some selected vibration mode shapes. The vibration mode shapes correspond to the eigenvectors obtained by solving a standard free vibration eigenvalue problem wherein the stiffness matrix is expanded about the nonlinear static solution. A strategy is also presented for determining which and how many vibration mode shapes to include in the transformation. The effect of inaccurate representation of the spatial distribution of the applied load on the nonlinear dynamic response is discussed for two classes of structural behavior. Application of the procedure to structures which exhibit a stiffening behavior and to those with a softening behavior is presented.

Knight, N. F., Jr.↗

Benchmark solutions for the galactic heavy-ion transport equations with energy and spatial coupling

Nontrivial benchmark solutions are developed for the galactic heavy ion transport equations in the straightahead approximation with energy and spatial coupling. Analytical representations of the ion fluxes are obtained for a variety of sources with the assumption that the nuclear interaction parameters are energy independent. The method utilizes an analytical LaPlace transform inversion to yield a closed form representation that is computationally efficient. The flux profiles are then used to predict ion dose profiles, which are important for shield design studies.

Ganapol, Barry D.↗

Fast neural solution of a nonlinear wave equation

A neural algorithm for rapidly simulating a certain class of nonlinear wave phenomena using analog VLSI neural hardware is presented and applied to the Korteweg-de Vries partial differential equation. The corresponding neural architecture is obtained from a pseudospectral representation of the spatial dependence, along with a leap-frog scheme for the temporal evolution. Numerical simulations demonstrated the robustness of the proposed approach.

Toomarian, Nikzad↗

A Search for Giant Convection Cells on the Sun

Giant convection cells (slow, long-lived cellular flows in the Sun's convection zone with typical diameters greater than about 100,000 km) have been the subject of many searches over the last 30 years. If such flows exist, they should play a key role in structuring the Sun's large scale magnetic field and in driving the large scale axisymmetric flows: the differential rotation and meridional circulation. Detailed observations of the flows in these cells may also allow us to better predict future magnetic field configurations and the solar activity associated with them. Line-of-sight velocity data from the Michelson Doppler Interferometer on the ESA/NASA Solar and Heliospheric Observatory provides us with new opportunities to search for giant cells. This data is free of any atmospheric distortion and has been obtained continuously without any day/night gaps for more than two months at a time. These two-month datasets are important because giant cells are expected to have lifetimes somewhat longer that the Sun's 27 day rotation period. Any reappearance of a flow pattern after 27 days would be an important confirmation of the existence of these cells. The approach taken in this search is to separate the giant cell velocity signal from the other, stronger velocity signals by using a spherical harmonic representation of the spatial structures and a fourier decomposition of the temporal behavior. Any giant cell signal should be characterized by low spatial wavenumbers with temporal frequencies appropriate to the solar rotation of these patterns.

Hathaway, David H.↗

The Zodiacal Emission Spectrum as Determined by COBE and its Implications

We combine observations from the DIRBE and FIRAS instruments on the COBE satellite to derive an annually-averaged spectrum of the zodiacal cloud in the 10 to 1000 micron wavelength region. The spectrum exhibits a break at approx. 150 microns which indicates a sharp break in the dust size distribution at a radius of about 30 microns The spectrum can be fit with a single blackbody with a lambda(exp -2) emissivity law beyond 150 microns and a temperature of 240 K. We also used a more realistic characterization of the cloud to fit the spectrum, including a distribution of dust temperatures, representing different dust compositions and distances from the sun, as well as a realistic representation of the spatial distribution of the dust. We show that amorphous carbon and silicate dust with respective temperatures of 280 and 274 K at 1 AU, and size distributions with a break at grain radii of 14 and 32 microns, can provide a good fit to the average zodiacal dust spectrum. The total mass of the zodiacal cloud is 2 to 11 Eg (Eg=10(exp 18) g), depending on the grain composition. The lifetime of the cloud, against particle loss by Poynting- Robertson drag and the effects of solar wind, is about 10(exp 5) yr. The required replenishment rate is approx. 10(exp 14) g/yr. If this is provided by asteroid belt alone, the asteroids lifetime would be approx. 3 x 10(exp 10) yr. But comets and Kuiper belt objects may also contribute to the zodiacal cloud.

Fixsen, D. J.↗

Anisotropic Developments for Homogeneous Shear Flows

The general decomposition of the spectral correlation tensor R(sub ij)(k) by Cambon et al. (J. Fluid Mech., 202, 295; J. Fluid Mech., 337, 303) into directional and polarization components is applied to the representation of R(sub ij)(k) by spherically averaged quantities. The decomposition splits the deviatoric part H(sub ij)(k) of the spherical average of R(sub ij)(k) into directional and polarization components H(sub ij)(sup e)(k) and H(sub ij)(sup z)(k). A self-consistent representation of the spectral tensor in the limit of weak anisotropy is constructed in terms of these spherically averaged quantities. The directional polarization components must be treated independently: models that attempt the same representation of the spectral tensor using the spherical average H(sub ij)(k) alone prove to be inconsistent with Navier-Stokes dynamics. In particular, a spectral tensor consistent with a prescribed Reynolds stress is not unique. The degree of anisotropy permitted by this theory is restricted by realizability requirements. Since these requirements will be less severe in a more accurate theory, a preliminary account is given of how to generalize the formalism of spherical averages to higher expansion of the spectral tensor. Directionality is described by a conventional expansion in spherical harmonics, but polarization requires an expansion in tensorial spherical harmonics generated by irreducible representations of the spatial rotation group SO(exp 3). These expansions are considered in more detail in the special case of axial symmetry.

Cambon, Claude↗

Impact of Optimized land Surface Parameters on the Land-Atmosphere Coupling in WRF Simulations of Dry and Wet Extremes

Land-atmosphere (L-A) interactions play a critical role in determining the diurnal evolution of both planetary boundary layer (PBL) and land surface temperature and moisture budgets, as well as controlling feedbacks with clouds and precipitation that lead to the persistence of dry and wet regimes. Recent efforts to quantify the strength of L-A coupling in prediction models have produced diagnostics that integrate across both the land and PBL components of the system. In this study, we examine the impact of improved specification of land surface states, anomalies, and fluxes on coupled WRF forecasts during the summers of extreme dry (2006) and wet (2007) conditions in the U.S. Southern Great Plains. The improved land initialization and surface flux parameterizations are obtained through the use of a new optimization and uncertainty module in NASA's Land Information System (LIS-OPT), whereby parameter sets are calibrated in the Noah land surface model and classified according to the land cover and soil type mapping of the observations and the full domain. The impact of the calibrated parameters on the a) spin up of land surface states used as initial conditions, and b) heat and moisture fluxes of the coupled (LIS-WRF) simulations are then assessed in terms of ambient weather, PBL budgets, and precipitation along with L-A coupling diagnostics. In addition, the sensitivity of this approach to the period of calibration (dry, wet, normal) is investigated. Finally, tradeoffs of computational tractability and scientific validity (e.g.,. relating to the representation of the spatial dependence of parameters) and the feasibility of calibrating to multiple observational datasets are also discussed.

Kumar, Sujay↗

The Mexican Drought Atlas: Tree-Ring Reconstructions of the Soil Moisture Balance During the Late Pre-Hispanic, Colonial, and Modern Eras

Mexico has suffered a long history and prehistory of severe sustained drought. Drought over Mexico is modulated by ocean-atmospheric variability in the Atlantic and Pacific, raising the possibility for long-range seasonal climate forecasting, which could help mediate the economic and social impacts of future dry spells. The instrumental record of Mexican climate is very limited before 1920, but tree-ring chronologies developed from old-growth forests in Mexico can provide an excellent proxy representation of the spatial pattern and intensity of past moisture regimes useful for the analysis of climate dynamics and climate impacts. The Mexican Drought Atlas (MXDA) has been developed from an extensive network of 252 climate sensitive tree-ring chronologies in and near Mexico. The MXDA reconstructions extend from 1400 CE-2012 and were calibrated with the instrumental summer (JJA) self-calibrating Palmer Drought Severity Index (scPDSI) on a 0.5deg latitude/longitude grid extending over land areas from 14 to 34degN and 75-120degW using Ensemble Point-by-Point Regression (EPPR) for the 1944-1984 period. The grid point reconstructions were validated for the period 1920-1943 against instrumental gridded scPDSI values based on the fewer weather station observations available during that interval. The MXDA provides a new spatial perspective on the historical impacts of moisture extremes over Mexico during the past 600-years, including the Aztec Drought of One Rabbit in 1454, the drought of El Ano de Hambre in 1785-1786, and the drought that preceded the Mexican Revolution of 1909-1910. The El Nino/Southern Oscillation (ENSO) is the most important ocean-atmospheric forcing of moisture variability detected with the MXDA. In fact, the reconstructions suggest that the strongest central equatorial Pacific sea surface temperature (SST) teleconnection to the soil moisture balance over North America may reside in northern Mexico. This ENSO signal has stronger and more time-stable correlations than computed for either the Atlantic Multidecadal Oscillation or Pacific Decadal Oscillation. The extended Multivariate ENSO Index is most highly correlated with reconstructed scPDSI over northern Mexico, where warm events favor moist conditions during the winter, spring, and early summer. This ENSO teleconnection to northern Mexico has been strong over the past 150 years, but it has been comparatively weak and non-stationary in the MXDA over central and southern Mexico where eastern tropical Pacific and Caribbean/tropical Atlantic SSTs seem to be more important. The ENSO teleconnection to northern Mexico is weaker in the available instrumental PDSI, but analyses based on the millennium climate simulations with the Community Earth System Model suggest that the moisture balance during the winter, spring, and early summer over northern Mexico may indeed be particularly sensitive to ENSO forcing. Nationwide drought is predicted to become more common with anthropogenic climate change, but the MXDA reconstructions indicate that intense "All Mexico" droughts have been rare over the past 600 years and their frequency does not appear to have increased substantially in recent decades.

El Nino↗

Microstructural Topology as a Prescriptor for Quantum Coherence: Towards A Unified Framework for Decoherence in Superconducting Qubits

In superconducting quantum circuits, decoherence improvements are frequently obtained through process interventions that simultaneously modify surface chemistry, microstructural topology, and device geometry, leaving mechanistic attribution structurally underdetermined. Predictive materials engineering requires measurable structural statistics to be separated from geometry-dependent coupling coefficients into independently testable factors. We introduce the concept of classical and quantum microstructure. In that context, we formulate a channel-wise separable framework for decoherence in superconducting transmon qubits in which each loss channel is described by a reduced prescriptor. Here, a channel-specific microstructural state variable is determined independently of device geometry, and a geometry-dependent coupling functional is computable from field solutions without reference to surface chemistry. We derive this product form from a spatially resolved kernel representation and establish a perturbative separability criterion that defines the regime where independent variation of the variables is valid. The framework specifies five prescriptor classes for dominant loss pathways in transmon-class devices. Falsifiability is operationalized through a pre-committed 2x2 experimental protocol in which the variables must satisfy independent ratio checks within propagated uncertainty. A Minimum-Dataset Specification standardizes reporting for cross-laboratory inference. Part I establishes the conceptual and mathematical architecture; coordinated experimental validation is reserved for Part II.

Dravid, Vinayak P. [Northwestern U.]↗

Including a Cold Pool Representation in a Convection Parameterization and Simulating Its Impacts on the Spatial and Temporal Variability of the Precipitation in the NASA GEOS GCM

We developed and implemented a simple representation of a cold pool in the Grell-Freitas (GF) convection parameterization. The cold pool parameterization is based on the observation that convective-scale downdrafts produce a local ‘deficit’ of the moist static energy (MSE). This information is advected and becoming downwind available to trigger and intensify new convection. The cold pool is dissipated by a simple exponential decay using a lifetime of a few hours, or by interacting with the underneath surface by exchanging latent and sensible heat fluxes. Preliminary results show some improvement of the simulation of the diurnal cycle of the precipitation over the land, mainly during the nighttime.

Cold pools↗

Interactions Between Thresholds and Spatial Discretizations of Snow: Insights From Estimates of Wolverine Denning Habitat in the Colorado Rocky Mountains

Thresholds can be used to interpret environmental data in a way that is easily communicated and useful for decision making purposes. However, thresholds are often developed for specific data products and time periods, changing findings when the same threshold is applied to datasets or periods with different characteristics. Here, we test the impact of different spatial discretizations of snow on annual estimates of wolverine denning opportunities in the Colorado Rocky Mountains, defined using a snow water equivalent (SWE) threshold (0.20 m) and threshold date (15 May) from previous habitat assessments. Annual potential wolverine denning area (PWDA) was thresholded from a 36-year (1985 – 2020) snow reanalysis model with three different spatial discretizations: 1) 480 m grid cells (D480), 2) 90 m grid cells (D90), and 3) 480 m grid cells with implicit representations of subgrid snow spatial heterogeneity (S480). Relative to the D480 and S480 discretizations, D90 resolved shallower snow deposits on slopes between 3050 and 3350 m elevation, decreasing PWDA by 10%, on average. In years with warmer and/or drier winters, S480 discretizations with subgrid representations of snow heterogeneity increased PWDA, even within grid cells where mean 15 May SWE was less than the SWE threshold. These simulations increased PWDA by upwards of 30% in low snow years, as compared to the D480 and D90 simulations without subgrid snow heterogeneity. Despite PWDA sensitivity to different snow spatial discretizations, PWDA was controlled more by annual variations in winter precipitation and temperature. However, small changes to the SWE threshold (± 0.07 m) and threshold date (± 2 weeks) also affected PWDA by as much as 82%. Across these threshold ranges, PWDA was approximately 18% more sensitive to the SWE threshold than the threshold date. However, the sensitivity to the threshold date was larger in years with late spring snowfall, when PWDA depended on whether modeled SWE was thresholded before, during, or after spring snow accumulation. Our results demonstrate that snow thresholds are useful but may not always provide a complete picture of the annual variability in snow-adapted wildlife denning opportunities. Studies thresholding spatiotemporal datasets could be improved by including 1) information about the fidelity of thresholds across multiple spatial discretizations, and 2) uncertainties related to ranges of realistic thresholds.

Justin M. Pflug↗

Gaussian process analysis of electron energy loss spectroscopy data: multivariate reconstruction and kernel control

Abstract Advances in hyperspectral imaging including electron energy loss spectroscopy bring forth the challenges of exploratory and physics-based analysis of multidimensional data sets. The multivariate linear unmixing methods generally explore similarities in the energy dimension, but ignore correlations in the spatial domain. At the same time, Gaussian process (GP) explicitly incorporate spatial correlations in the form of kernel functions but is computationally intensive. Here, we implement a GP method operating on the full spatial domain and reduced representations in the energy domain. In this multivariate GP, the information between the components is shared via a common spatial kernel structure, while allowing for variability in the relative noise magnitude or image morphology. We explore the role of kernel constraints on the quality of the reconstruction, and suggest an approach for estimating them from the experimental data. We further show that spatial information contained in higher-order components can be reconstructed and spatially localized.

36 MATERIALS SCIENCE↗