Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spatial dependence”

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,693 records · Page 94

GERMcode: A Stochastic Model for Space Radiation Risk Assessment

A new computer model, the GCR Event-based Risk Model code (GERMcode), was developed to describe biophysical events from high-energy protons and high charge and energy (HZE) particles that have been studied at the NASA Space Radiation Laboratory (NSRL) for the purpose of simulating space radiation biological effects. In the GERMcode, the biophysical description of the passage of HZE particles in tissue and shielding materials is made with a stochastic approach that includes both particle track structure and nuclear interactions. The GERMcode accounts for the major nuclear interaction processes of importance for describing heavy ion beams, including nuclear fragmentation, elastic scattering, and knockout-cascade processes by using the quantum multiple scattering fragmentation (QMSFRG) model. The QMSFRG model has been shown to be in excellent agreement with available experimental data for nuclear fragmentation cross sections. For NSRL applications, the GERMcode evaluates a set of biophysical properties, such as the Poisson distribution of particles or delta-ray hits for a given cellular area and particle dose, the radial dose on tissue, and the frequency distribution of energy deposition in a DNA volume. By utilizing the ProE/Fishbowl ray-tracing analysis, the GERMcode will be used as a bi-directional radiation transport model for future spacecraft shielding analysis in support of Mars mission risk assessments. Recent radiobiological experiments suggest the need for new approaches to risk assessment that include time-dependent biological events due to the signaling times for activation and relaxation of biological processes in cells and tissue. Thus, the tracking of the temporal and spatial distribution of events in tissue is a major goal of the GERMcode in support of the simulation of biological processes important in GCR risk assessments. In order to validate our approach, basic radiobiological responses such as cell survival curves, mutation, chromosomal aberrations, and representative mouse tumor induction curves are implemented into the GERMcode. Extension of these descriptions to other endpoints related to non-targeted effects and biochemical pathway responses will be discussed.

Kim, Myung-Hee Y.↗

Using Principal Component Analysis (PCA) to Speed up Radiative Transfer (RT) Computations

Multiple scattering RT calculations time-consuming. Need a speed improvement of about 1000 (for OCO)! Solution: Make use of redundancies in spectra. Correlated-k (Lacis and Wang, Lacis and Oinas, Goody et al, Fu and Liou) Problem: Assume that spectral variation of atmospheric optical properties spatially correlated at all points along optical path. High accuracy (HI) and 2-stream (2S) calculations have high correlation. Single scattering (SS) computations highly scenario-dependent, but not time consuming. Perform SS and 2S calculations at every wavelength. Perform small number of HI computations. Need to compute correction factor B at every wavelength.

data sets↗

The Impact of ARM on Climate Modeling

Climate models are among humanity's most ambitious and elaborate creations. They are designed to simulate the interactions of the atmosphere, ocean, land surface, and cryosphere on time scales far beyond the limits of deterministic predictability, and including the effects of time-dependent external forcings. The processes involved include radiative transfer, fluid dynamics, microphysics, and some aspects of geochemistry, biology, and ecology. The models explicitly simulate processes on spatial scales ranging from the circumference of the Earth down to one hundred kilometers or smaller, and implicitly include the effects of processes on even smaller scales down to a micron or so. The atmospheric component of a climate model can be called an atmospheric global circulation model (AGCM). In an AGCM, calculations are done on a three-dimensional grid, which in some of today's climate models consists of several million grid cells. For each grid cell, about a dozen variables are time-stepped as the model integrates forward from its initial conditions. These so-called prognostic variables have special importance because they are the only things that a model remembers from one time step to the next; everything else is recreated on each time step by starting from the prognostic variables and the boundary conditions. The prognostic variables typically include information about the mass of dry air, the temperature, the wind components, water vapor, various condensed-water species, and at least a few chemical species such as ozone. A good way to understand how climate models work is to consider the lengthy and complex process used to develop one. Lets imagine that a new AGCM is to be created, starting from a blank piece of paper. The model may be intended for a particular class of applications, e.g., high-resolution simulations on time scales of a few decades. Before a single line of code is written, the conceptual foundation of the model must be designed through a creative envisioning that starts from the intended application and is based on current understanding of how the atmosphere works and the inventory of mathematical methods available.

atmospheric circulation↗

Multibaseline TanDEM-X Mangrove Height Estimation: The Selection of the Vertical Wavenumber

We generated a large-scale mangrove forest height map using multiple TanDEM-X (TDX) interferometric synthetic aperture radar (InSAR) acquisitions with various spatial baselines in order to improve the height estimation accuracy across a wide range of forest heights. The forest height inversion using InSAR data is strongly dependent upon the vertical wavenumber (i.e., perpendicular baseline). First, we investigated the role of the vertical wavenumber in forest height inversion from InSAR data using the sensitivity of the interferometric (volume) coherence to forest height. We used corrected but lower resolution and accuracy Shuttle Radar Topography Mission (SRTM) mangrove height maps as a priori information over Akanda and Pongara National Parks in Gabon to estimate lower and upper boundaries of the vertical wavenumber over test sites from the measured coherence-to-height sensitivity. Only TDX acquisitions within the boundaries of the vertical wavenumber were selected and combined for multibaseline mangrove height inversion. Mangrove forest height was obtained with multibaseline TDX acquisitions and was validated against the reference height derived from field measurement data providing improvements in multibaseline inversion over existing height estimates (i.e., SRTM height) and single-baseline inversions (multibaseline inversion: r(exp 2) = 0.98, root mean square error (RMSE) of 2.73 m; SRTM height: r(exp 2) = 0.86, RMSE = 7.21 m; single-baseline inversions: r(exp 2) = 0.08-0.97, RMSE = 3.86-11.10 m). As a result, to accurately estimate forest heights over a wide range (3-60 m), multibaseline InSAR acquisitions (at least three different baselines) are needed to exclude biases associated with the vertical wavenumber in forest height inversion.

vertical wavenumber.↗

Electrifying Road Trips to and from National Parks in the Western United States

We evaluate charging infrastructure requirement for electrified long-distance travels to/from national parks in the Western U.S. We conduct a high-resolution spatial and temporal analysis of road trips in the region, accounting for detailed travel volume and pattern, the type of electric vehicles, heterogeneous energy consumption rates depending on driving conditions, charging technology (e.g., power and speed), and charging behavior. The scope includes on-route (or way-point) charging as well as destination charging (in national parks). We estimate the locations and sizes of required charging stations in the entire Western U.S. to enable electrified road trips to/from national parks. We also examine environmental justice and energy equity aspects of projected charging infrastructure (for national parks) in the context of the overall vehicle electrification and charging infrastructure build-out.

ADVANCED PROPULSION SYSTEMS,ENERGY PLANNING, POLIC↗

Spatial profiles of photon chemical potential in near-field thermophotovoltaic cells

Emitted photons stemming from the radiative recombination of electron–hole pairs carry chemical potential in radiative energy converters. This luminescent effect can substantially alter the local net photogeneration in near-field thermophotovoltaic cells. Several assumptions involving the luminescent effect are commonly made in modeling photovoltaic devices; in particular, the photon chemical potential is assumed to be zero or a constant prescribed by the bias voltage. The significance of photon chemical potential depends upon the emitter temperature, the semiconductor properties, and the injection level. Hence, these assumptions are questionable in near-field thermophotovoltaic devices with nanoscale separation distances between the emitter and the cell. In the present work, an iterative solver that combines fluctuational electrodynamics with the drift-diffusion model is developed to tackle the coupled photon and charge transport problem, enabling the determination of the spatial profile of photon chemical potential beyond the detailed balance approach. Here, the difference between the results obtained by allowing the photon chemical potential to vary spatially and by assuming a constant value demonstrates the limitations of the conventional approaches. This study is critically important for performance evaluation of near-field thermophotovoltaic systems with nanoscale vacuum gaps.

30 DIRECT ENERGY CONVERSION↗

Nonoverlapping Localized Exponential Time Differencing Methods for Diffusion Problems

In this paper, we propose nonoverlapping localized exponential time differencing (ETD) methods for diffusion problems. The model time-dependent diffusion equation is first reformulated on subdomains based on the nonoverlapping domain decomposition, in which Neumann boundary conditions are imposed on the interfaces for the subdomain problems and Dirichlet type conditions are enforced to form a space-time interface problem. After spatial discretization by standard central finite differences and temporal integration with the first or second order ETD methods, the fully discrete interface problem is obtained. Such an interface problem is then solved iteratively either at each time step or over the whole time interval: the former involves the solution of stationary problems in each subdomain at each iteration while the latter involves the solution of time-dependent subdomain problems at each iteration. For both approaches, we prove that localized ETD solutions conserve mass exactly and converge in time to the exact space semidiscrete solution. Numerical experiments in two dimensions are also presented to illustrate the performance of the proposed methods.

Mass conservation↗

Effects of Planetary Thermal Structure on the Ascent and Cooling of Magma on Venus

Magellan radar images of the surface of Venus show a spatially broad distribution of volcanic features. Models of magmatic ascent processes to planetary surfaces indicate that the thermal structure of the interior significantly influences the rate of magmatic cooling and thus the amount of magma that can be transported to the surface before solidification. In order to understand which aspects of planetary thermal structure have the greatest influence on the cooling of buoyantly ascending magma, we have constructed magma cooling profiles for a plutonic ascent mechanism, and evaluated the profiles for variations in the surface and mantle temperature, surface temperature gradient, and thermal gradient curvature. Results show that, for a wide variety of thermal conditions, smaller and slower magma bodies are capable of reaching the surface on Venus compared to Earth, primarily due to the higher surface temperature of Venus. Little to no effect on the cooling and transport of magma are found to result from elevated mantle temperatures, elevation-dependent surface temperature variations, or details of the thermal gradient curvature. The enhanced tendency of magma to reach the surface on Venus may provide at least a partial explanation for the extensive spatial distribution of observed volcanism on the surface.

Sakimoto, Susan E. H.↗

A Binomial Stochastic Framework for Efficiently Modeling Discrete Statistics of Convective Populations

Abstract Understanding the coupling between convective clouds and the general circulation, as well as addressing the gray zone problem in convective parameterization, requires insight into the genesis and maintenance of spatial patterns in cumulus cloud populations. In this study, a simple toy model for recreating populations of interacting convective objects as distributed over a two‐dimensional Eulerian grid is formulated to this purpose. Key elements at the foundation of the model include i) a fully discrete formulation for capturing discrete behavior in convective properties at small population sample sizes, ii) object age‐dependence for representing life‐cycle effects, and iii) a prognostic number budget allowing for object interactions and co‐existence of multiple species. A primary goal is to optimize the computational efficiency of this system. To this purpose the object birth rate is represented stochastically through a spatially aware Bernoulli process. The same binomial stochastic operator is applied to horizontal advection of objects, conserving discreteness in object number. The applicability to atmospheric convection as well as behavior implied by the formulation is assessed. Various simple applications of the BiOMi model (Binomial Objects on Microgrids) are explored, suggesting that important convective behavior can be captured at low computational cost. This includes i) subsampling effects and associated powerlaw scaling in the convective gray zone, ii) stochastic predator‐prey behavior, iii) the downscale turbulent energy cascade, and iv) simple forms of spatial organization and convective memory. Consequences and opportunities for convective parameterization in next‐generation weather and climate models are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Estimation of Sea Ice Thickness Distributions through the Combination of Snow Depth and Satellite Laser Altimetry Data

Combinations of sea ice freeboard and snow depth measurements from satellite data have the potential to provide a means to derive global sea ice thickness values. However, large differences in spatial coverage and resolution between the measurements lead to uncertainties when combining the data. High resolution airborne laser altimeter retrievals of snow-ice freeboard and passive microwave retrievals of snow depth taken in March 2006 provide insight into the spatial variability of these quantities as well as optimal methods for combining high resolution satellite altimeter measurements with low resolution snow depth data. The aircraft measurements show a relationship between freeboard and snow depth for thin ice allowing the development of a method for estimating sea ice thickness from satellite laser altimetry data at their full spatial resolution. This method is used to estimate snow and ice thicknesses for the Arctic basin through the combination of freeboard data from ICESat, snow depth data over first-year ice from AMSR-E, and snow depth over multiyear ice from climatological data. Due to the non-linear dependence of heat flux on ice thickness, the impact on heat flux calculations when maintaining the full resolution of the ICESat data for ice thickness estimates is explored for typical winter conditions. Calculations of the basin-wide mean heat flux and ice growth rate using snow and ice thickness values at the 70 m spatial resolution of ICESat are found to be approximately one-third higher than those calculated from 25 km mean ice thickness values.

Kurtz, Nathan T.↗

Emergent disorder and mechanical memory in periodic metamaterials

Ordered mechanical systems typically have one or only a few stable rest configurations, and hence are not considered useful for encoding memory. Multistable and history-dependent responses usually emerge from quenched disorder, for example in amorphous solids or crumpled sheets. In contrast, due to geometric frustration, periodic magnetic systems can create their own disorder and espouse an extensive manifold of quasi-degenerate configurations. Inspired by the topological structure of frustrated artificial spin ices, we introduce an approach to design ordered, periodic mechanical metamaterials that exhibit an extensive set of spatially disordered states. While our design exploits the correspondence between frustration in magnetism and incompatibility in meta-mechanics, our mechanical systems encompass continuous degrees of freedom, and thus generalize their magnetic counterparts. We show how such systems exhibit non-Abelian and history-dependent responses, as their state can depend on the order in which external manipulations were applied. We demonstrate how this richness of the dynamics enables to recognize, from a static measurement of the final state, the sequence of operations that an extended system underwent. Thus, multistability and potential to perform computation emerge from geometric frustration in ordered mechanical lattices that create their own disorder.

36 MATERIALS SCIENCE↗

Time dependent wind fields

Two tasks were performed: (1) determination of the accuracy of Seasat scatterometer, altimeter, and scanning multichannel microwave radiometer measurements of wind speed; and (2) application of Seasat altimeter measurements of sea level to study the spatial and temporal variability of geostrophic flow in the Antarctic Circumpolar Current. The results of the first task have identified systematic errors in wind speeds estimated by all three satellite sensors. However, in all cases the errors are correctable and corrected wind speeds agree between the three sensors to better than 1 ms sup -1 in 96-day 2 deg. latitude by 6 deg. longitude averages. The second task has resulted in development of a new technique for using altimeter sea level measurements to study the temporal variability of large scale sea level variations. Application of the technique to the Antarctic Circumpolar Current yielded new information about the ocean circulation in this region of the ocean that is poorly sampled by conventional ship-based measurements.

Chelton, D. B.↗

An Adaptive Flow Solver for Air-Borne Vehicles Undergoing Time-Dependent Motions/Deformations

This report describes a concurrent Euler flow solver for flows around complex 3-D bodies. The solver is based on a cell-centered finite volume methodology on 3-D unstructured tetrahedral grids. In this algorithm, spatial discretization for the inviscid convective term is accomplished using an upwind scheme. A localized reconstruction is done for flow variables which is second order accurate. Evolution in time is accomplished using an explicit three-stage Runge-Kutta method which has second order temporal accuracy. This is adapted for concurrent execution using another proven methodology based on concurrent graph abstraction. This solver operates on heterogeneous network architectures. These architectures may include a broad variety of UNIX workstations and PCs running Windows NT, symmetric multiprocessors and distributed-memory multi-computers. The unstructured grid is generated using commercial grid generation tools. The grid is automatically partitioned using a concurrent algorithm based on heat diffusion. This results in memory requirements that are inversely proportional to the number of processors. The solver uses automatic granularity control and resource management techniques both to balance load and communication requirements, and deal with differing memory constraints. These ideas are again based on heat diffusion. Results are subsequently combined for visualization and analysis using commercial CFD tools. Flow simulation results are demonstrated for a constant section wing at subsonic, transonic, and a supersonic case. These results are compared with experimental data and numerical results of other researchers. Performance results are under way for a variety of network topologies.

Singh, Jatinder↗

Spatial non-locality of electronic correlations beyond GW approximation

The question of spatial locality of electronic correlations beyond GW approximation is one of the central issues of the famous combination of GW and dynamical mean field theory, GW+DMFT. In this study, the above question is addressed directly (for the first time) by performing calculations with and without assumption of locality of the corresponding diagrams. For this purpose we use sc(GW+G3W2) approach where the higher order part (G3W2) is evaluated with fully momentum dependent Green's function G and screened interaction W and with "local" variant, where the single site approximation is assumed for both G and W. For all three materials studied in this work (NiO, α-Ce, LiFeAs), we have found the spatial non-locality effects to be strong. For NiO and LiFeAs they, in fact, are decisive for the proper evaluation of vertex corrections. The results of this study have direct impact on our understanding of approximations made in practical implementations of GW+DMFT method, where all diagrams beyond GW (DMFT part) are assumed to be local. Taking into account the fact that the first diagrams beyond GW represent the most important contribution also in GW+DMFT calculations, we conclude that the basic assumption of GW+DMFT, namely the locality of diagrams evaluated in the DMFT part, is not as good as it is believed to be.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A spatially regularized detector for emergent/re-emergent disease outbreaks

Early detection of outbreaks caused by emergent pathogens, using epidemiological surveillance data i.e., daily case counts, is difficult. This is because the data tend to be noisy during the early epoch of the outbreak. In contrast, the spread-rate of the disease tends to be well-behaved, as it depends only on the mixing patterns of the population and the characteristics of the pathogen, neither of which behave erratically in space-time. In this report, we explore whether the spread-rate can be used for epidemiological surveillance, conditional on case count data. Estimating the spread-rate from case count data allows us to exploit exogenous information, e.g., incubation period distributions etc., which can considerably smooth out any erratic temporal behavior. Further, epidemiological dynamics are spatially correlated, and if case counts are available for multiple areal units e.g., counties, these correlations could potentially be used to suppress noise in the early epoch data. These exogenous information and structure are not exploited by conventional syndromic surveillance detectors to extract a well-behaved latent variable for monitoring purposes. The technical challenge lies in the estimation of the spread-rate field jointly over a collection of areal units; further, the spread-rate varies over time. We develop a method based on mean-field variational inference to approximately estimate the spread-rate field, using a Gaussian Random Field Model for spatial regularization. The method is tested on the estimation of spread-rate in the thirty-three counties of New Mexico and detect the arrival of the Fall 2020 COVID-19 wave in September 2020. We find that the method is scalable, but underestimates the uncertainty in the estimated spread-rate field. We detect the arrival of the Fall 2020 wave a week ahead of conventional syndromic surveillance algorithms, but our simplistic detection algorithm, based on simple anomaly detection, suffers from a high false positive rate, similar to conventional detectors.

59 BASIC BIOLOGICAL SCIENCES↗

Extraction of Profile Information from Cloud Contaminated Radiances

Clouds act to reduce the signal level and may produce noise dependence on the complexity of the cloud properties and the manner in which they are treated in the profile retrieval process. There are essentially three ways to extract profile information from cloud contaminated radiances: (1) cloud-clearing using spatially adjacent cloud contaminated radiance measurements, (2) retrieval based upon the assumption of opaque cloud conditions, and (3) retrieval or radiance assimilation using a physically correct cloud radiative transfer model which accounts for the absorption and scattering of the radiance observed. Cloud clearing extracts the radiance arising from the clear air portion of partly clouded fields of view permitting soundings to the surface or the assimilation of radiances as in the clear field of view case. However, the accuracy of the clear air radiance signal depends upon the cloud height and optical property uniformity across the two fields of view used in the cloud clearing process. The assumption of opaque clouds within the field of view permits relatively accurate profiles to be retrieved down to near cloud top levels, the accuracy near the cloud top level being dependent upon the actual microphysical properties of the cloud. The use of a physically correct cloud radiative transfer model enables accurate retrievals down to cloud top levels and below semi-transparent cloud layers (e.g., cirrus). It should also be possible to assimilate cloudy radiances directly into the model given a physically correct cloud radiative transfer model using geometric and microphysical cloud parameters retrieved from the radiance spectra as initial cloud variables in the radiance assimilation process. This presentation reviews the above three ways to extract profile information from cloud contaminated radiances. NPOESS Airborne Sounder Testbed-Interferometer radiance spectra and Aqua satellite AIRS radiance spectra are used to illustrate how cloudy radiances can be used in the profile retrieval process.

Smith, W. L.↗

Reciprocity relations in a biologically inspired nanomagnonic system with dipolar coupling

In this study, magnetosome chains in magnetotactic bacteria present ideal nanomagnonic model systems for studying collective resonance modes of dipolar-coupled single domain particles in relation to their spatial arrangement. Using microresonator-based ferromagnetic resonance (FMR) spectroscopy, electron microscopy, and micromagnetic modeling, we here provide insights into the complex magnonic activity within a single magnetosome chain. While the angular dependence of its FMR spectrum is dominated by twofold symmetry features due to the uniaxial anisotropy of linear chain segments, we also observed an unexpected behavior such as interrupted lines and flat bands due to the intricate geometrical details of this particular chain, such as a cross-like structural anomaly where a pair of particles is oriented perpendicular to the main axis of the chain and thus breaks the prevailing axial dipolar coupling symmetry. Such a cross junction formed by four particles exhibits interesting magnonic network properties. Notably, we observe reciprocity in the sense that the spectral response of one particle to an excitation of another one is identical to the response of the latter given an excitation of the former. Furthermore, we have identified that magnonic coupling between A and B can be facilitated via a dark state, as in magnonic stimulated Raman adiabatic passage, and that this dark-state coupling can be made non-reciprocal between A and B by breaking the symmetry of the spatial arrangement of the four particles.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Can Restoring Tidal Wetlands Reduce Estuarine Nuisance Flooding of Coasts Under Future Sea‐Level Rise?

Wetland restoration is an increasingly popular nature‐based method for flood risk mitigation in coastal communities. In this study, we present a novel method using hydrodynamic modeling and harmonic analysis to quantify wetlands' ability to reduce future nuisance flooding. The method leverages a hydrodynamic model calibrated to present day data and was run for a range of future sea‐level rise (SLR) and wetland restoration scenarios to quantify changes to tidal harmonic amplitudes and phases. The harmonic constituents are used to generate water surface elevations over a time period of interest (e.g., one year) and compared to critical exceedance thresholds such as levee elevations. Then, changes to nuisance flooding are calculated by counting the number of hours critical thresholds are exceeded under different SLR and wetland restoration scenarios. We applied the method to Coos Bay, Oregon, USA as a test case. We found restoration reduces the number of hours nuisance flooding occurs in downtown Coos Bay from 15 hr (present day conditions) to 0 hr (fully restored condition) under median SLR (82 cm by 2100). Restoration had spatially variable impacts on reducing peak flood elevations with minimal impacts near the estuary mouth and greatest impact 32 km inland. The effectiveness of restoration was heavily dependent on future SLR. Restoration was maximally effective in 2050 under all SLR scenarios, less effective in 2100 under median SLR, and not effective under high SLR. Modeling results suggest increased tidal prism and accommodation space are driving restoration‐associated reductions in tidal amplitudes.

Brand, Matthew W. [Louisiana State Univ., Baton Ro↗