SEARCH · Search NASA

Results for “spatially explicit modeling”

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.

SourceResearch connectionFrom yearThrough year

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

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Spatially Quasi-Periodic Water Waves of Infinite Depth

Abstract We formulate the two-dimensional gravity-capillary water wave equations in a spatially quasi-periodic setting and present a numerical study of solutions of the initial value problem. We propose a Fourier pseudo-spectral discretization of the equations of motion in which one-dimensional quasi-periodic functions are represented by two-dimensional periodic functions on a torus. We adopt a conformal mapping formulation and employ a quasi-periodic version of the Hilbert transform to determine the normal velocity of the free surface. Two methods of time-stepping the initial value problem are proposed, an explicit Runge–Kutta (ERK) method and an exponential time-differencing (ETD) scheme. The ETD approach makes use of the small-scale decomposition to eliminate stiffness due to surface tension. We perform a convergence study to compare the accuracy and efficiency of the methods on a traveling wave test problem. We also present an example of a periodic wave profile containing vertical tangent lines that is set in motion with a quasi-periodic velocity potential. As time evolves, each wave peak evolves differently, and only some of them overturn. Beyond water waves, we argue that spatial quasi-periodicity is a natural setting to study the dynamics of linear and nonlinear waves, offering a third option to the usual modeling assumption that solutions either evolve on a periodic domain or decay at infinity.

Wilkening, Jon (ORCID:0000000327827596)↗

DOE OSTI2023

Estimating National-Scale Wind Potential Using Spatially Explicit Turbine Layout Optimization

National renewable energy potential assessments play a broad and critical role in analysis of the clean energy transition by providing foundational estimates of developable clean resources. Common to all past wind potential assessments is an assumption that wraps the complexity of wind plant layout (arrangement of turbines) into a single metric known as capacity density, or rated power capacity per unit of land area. Quite often, a singular capacity density or rotor-diameter driven capacity density is used in wind potential assessments across broad geographies despite the complexities of local drivers. Here, we present a new wind technical potential assessment for the United States, leveraging a spatial optimization approach in lieu of the traditional uniform capacity density. The optimization approach is a spatially explicit method for determining the potential locations of individual wind turbines-taking into account the turbine configuration, plant economics and losses, wind resource, and siting considerations. Our approach accounts for the interactions between wind technology design, wind plant layout, and the vast array of regulatory, land use and infrastructure conflicts with wind development. Our results highlight the potential ability of larger turbines to enable increased wind capacity, up to a point, and increased generation when siting turbines in and around spatial constraints; moreover, they demonstrate and capture the LCOE benefit of relatively lower capacity densities and reduced wake losses when land is relatively abundant. These insights provide foundational knowledge for the wind sector as it develops and pursues future turbine models and as wind energy markets expand in zero-carbon futures. Further, when applied in capacity expansion models, supply curves developed by these methods can provide detailed local insights about where wind turbines might be deployed in and among known siting constraints for those regions where wind energy is determined to be economic, providing critical nuance to local decision-makers and stakeholders.

17 WIND ENERGY↗

DOE OSTI2026

Opportunities for bioenergy crops to support transitions from irrigated agriculture and conserve the U.S. High Plains Aquifer

This study investigates potential economic and groundwater driven transitions from irrigated maize production—the dominant irrigated cropping system in the High Plains Aquifer (HPA) region—to alternative crops such as sorghum and switchgrass, two common bioenergy feedstocks. Unsustainable groundwater extraction in the U.S. High Plains presents critical challenges including reduced irrigation capacities, diminished crop yields, lower land values, and escalating energy costs. Using a spatially explicit optimization framework combined with comprehensive economic and hydrogeological data, we evaluate optimal land-use strategies and irrigation system investments over a 30-year planning period. Results indicate significant regional variability in the future economic viability of irrigated agriculture, driven by differences in aquifer recharge rates, groundwater availability, and market conditions. Nebraska and parts of northern Texas can sustain irrigated maize profitability due to relatively favorable groundwater conditions and lower land rents, respectively. In contrast, many portions of Kansas and southern Texas are more likely to transition to dryland agriculture within two decades. Colorado and New Mexico show potential for significant adoption of switchgrass production as an alternative biomass-based energy crop. Overall, over the 30-year horizon, our model implies that approximately 23% of currently irrigated maize area across the HPA region may transition to dryland farming under business-as-usual conditions. This transition is complemented by a threefold increase in non-irrigated sorghum production, from 0.20 to 0.77 Mt yr−1, indicating that groundwater-driven shifts in agricultural production may support a larger regional base for bioenergy feedstocks. The study also reveals opportunities for producers to optimize economic returns and biomass production potential.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2020

PRISMS-PF: A general framework for phase-field modeling with a matrix-free finite element method

Abstract A new phase-field modeling framework with an emphasis on performance, flexibility, and ease of use is presented. Foremost among the strategies employed to fulfill these objectives are the use of a matrix-free finite element method and a modular, application-centric code structure. This approach is implemented in the new open-source PRISMS-PF framework. Its performance is enabled by the combination of a matrix-free variant of the finite element method with adaptive mesh refinement, explicit time integration, and multilevel parallelism. Benchmark testing with a particle growth problem shows PRISMS-PF with adaptive mesh refinement and higher-order elements to be up to 12 times faster than a finite difference code employing a second-order-accurate spatial discretization and first-order-accurate explicit time integration. Furthermore, for a two-dimensional solidification benchmark problem, the performance of PRISMS-PF meets or exceeds that of phase-field frameworks that focus on implicit/semi-implicit time stepping, even though the benchmark problem’s small computational size reduces the scalability advantage of explicit time-integration schemes. PRISMS-PF supports an arbitrary number of coupled governing equations. The code structure simplifies the modification of these governing equations by separating their definition from the implementation of the numerical methods used to solve them. As part of its modular design, the framework includes functionality for nucleation and polycrystalline systems available in any application to further broaden the phenomena that can be used to study. The versatility of this approach is demonstrated with examples from several common types of phase-field simulations, including coarsening subsequent to spinodal decomposition, solidification, precipitation, grain growth, and corrosion.

36 MATERIALS SCIENCE↗

DOE OSTI2023

A Vertically Resolved Analysis of Radiative Feedbacks on Moist Static Energy Variance in Tropical Cyclones

Abstract A vertically resolved moist static energy (MSE) variance budget framework is used to diagnose processes associated with the development of tropical cyclones (TCs) in a general circulation model (GCM) under realistic boundary conditions. Previous studies have shown that interactions between radiation and MSE promote TC development. Here, we examine the vertical contributions of radiation and its interactions with MSE by performing several mechanism-denial experiments in which synoptic-scale radiative interactions are suppressed either in the boundary layer or in the free troposphere. Partly suppressing radiative interactions results in a reduction in global TC frequency. However, the magnitude of reduction and structure of the feedback depend on the intensity and structure of the TCs in these mechanism-denial experiments, indicating that both the magnitude and the vertical location of radiative interactions can impact global TC frequency. Using instantaneous 6-hourly outputs, an explicit computation reveals distinct spatial patterns of the advection term: the vertical component is positive in the mid- to upper troposphere, which reflects an upward transport of MSE by deep convection, whereas the horizontal component is positive in the boundary layer. These results illustrate the impact of the vertical distribution of radiative interactions and vertically varied contribution of the advection term in the development of TCs.

Meteorology & Atmospheric Sciences↗

NASA NTRS1973

Martian thermal boundary layers: Subhourly variations induced by radiative-conductive heat transfer within the dust-laden atmosphere-ground system

Radiative-conductive heat transfer has been investigated for the ground-atmosphere system of the planet Mars. The basic goal was the quantitative determination of time dependent vertical distributions of temperature and static stability for Southern-Hemispheric summer season and middle and polar latitudes, for both dust-free and dust-laden atmospheric conditions. The numerical algorithm which models at high spatial and temporal resolution the thermal energy transports in the dual ground-atmosphere system, is based on solution of the applicable heating rate equation, including radiative and molecular-conductive heat transport terms. The two subsystems are coupled by an internal thermal boundary condition applied at the ground-atmosphere interface level. Initial data and input parameters are based on Mariner 4, 6, 7, and 9 measurements and the JPL Mars Scientific Model. Numerical experiments were run for dust-free and dust-laden conditions in the midlatitudes, as well as ice-free and ice-covered polar regions. Representative results and their interpretation are presented. Finally, the theoretical framework of the generalized problem with nonconservative Mie scattering and explicit thermal-convective heat transfer is formulated, and applicable solution algorithms are outlined.

Pallmann, A. J.↗

NASA NTRS2013

Infrared Transmission Spectroscopy of the Exoplanets HD 209458b and XO-1b Using the Wide Field Camera-3 on the Hubble Space Telescope

Exoplanetary transmission spectroscopy in the near-infrared using the Hubble Space Telescope (HST) NICMOS (Near Infrared Camera and Multi-Object Spectrometer) is currently ambiguous because different observational groups claim different results from the same data, depending on their analysis methodologies. Spatial scanning with HST/WFC3 (Wide Field Camera-3) provides an opportunity to resolve this ambiguity.We here report WFC3 spectroscopy of the giant planets HD 209458b and XO-1b in transit, using spatial scanning mode for maximum photon-collecting efficiency. We introduce an analysis technique that derives the exoplanetary transmission spectrum without the necessity of explicitly decorrelating instrumental effects, and achieves nearly photon-limited precision even at the high flux levels collected in spatial scan mode. Our errors are within 6 percent (XO-1) and 26 percent (HD 209458b) of the photon-limit at a resolving power of lambda divided by delta times lambda approximating 70, and are better than 0.01 percent per spectral channel. Both planets exhibit water absorption of approximately 200 ppm at the water peak near 1.38 m. Our result for XO-1b contradicts the much larger absorption derived from NICMOS spectroscopy. The weak water absorption we measure for HD209458b is reminiscent of the weakness of sodium absorption in the first transmission spectroscopy of an exoplanet atmosphere by Charbonneau et al. Model atmospheres having uniformly distributed extra opacity of 0.012 square centimeters per gram account approximately for both our water measurement and the sodium absorption. Our results for HD 209458b support the picture advocated by Pont et al. in which weak molecular absorptions are superposed on a transmission spectrum that is dominated by continuous opacity due to haze and or dust. However,the extra opacity needed for HD 209458b is grayer than for HD 189733b, with a weaker Rayleigh component.

planetary systems↗

DOE OSTI2024

Predicted range expansion of Prostephanus truncatus (Coleoptera: Bostrichidae) under projected climate change scenarios

Abstract The larger grain borer (Prostephanus truncatus [Horn] [Coleoptera: Bostrichidae]) is a wood-boring insect native to Central America and adapted to stored maize and cassava. It was accidentally introduced to Tanzania and became a pest across central Africa. Unlike many grain pests, P. truncatus populations can establish and move within forests. Consequently, novel infestations can occur without human influence. The objectives of our study were to (i) develop an updated current suitability projection for P. truncatus, (ii) assess its potential future distribution under different climate change scenarios, and (iii) identify climate variables that best inform the model. We used WALLACE and MaxEnt to predict potential global distribution by incorporating bioclimatic variables and occurrence records. Future models were projected for 2050 and 2070 with Representative Concentration Pathways (RCPs) 2.6 (low change) and 8.5 (high change). Distribution was most limited by high precipitation and cold temperatures. Globally, highly suitable areas (> 75%) primarily occurred along coastal and equatorial regions with novel areas in northern South America, India, southeastern Asia, Indonesia, and the Philippines, totaling 7% under current conditions. Highly suitable areas at RCPs 2.6 and 8.5 are estimated to increase to 12% and 15%, respectively, by 2050 and increase to 19% in 2070 under RCP 8.5. Centroids of highly suitable areas show distribution centers moving more inshore and away from the equator. Notably, the result is a range expansion, not a shift. Results can be used to decrease biosecurity risks through more spatially explicit and timely surveillance programs for targeting the exclusion of this pest.

Entomology↗

NASA NTRS1994

The Barberplaid Illusion

Mulligan showed that the perceived direction of a moving grating can be biased by the shape of the Gaussian window in which it is viewed. We sought to determine if a 2-D pattern with an unambiguous velocity would also show such biases. Observers viewed a drifting plaid (sum of two orthogonal 2.5 c/d sinusoidal gratings of 12% contrast, each with a TF of 4 Hz.) whose contrast was modulated spatially by a stationary, asymmetric 2-D Gaussian window (i.e. unequal standard deviations in the principal directions). The direction of plaid motion with respect to the orientation of the window's major axis (Delta Theta) was varied while all other motion parameters were held fixed. Observers reported the perceived plaid direction of motion by adjusting the orientation of a pointer. All five observers showed systematic biases in perceived plaid direction that depended on Delta Theta and the aspect ratio of the Gaussian window (lambda). For circular Gaussian windows Lambda = 1), plaid direction was veridically perceived. However, biases of up to 10 deg. were found for lambda = 2 and Delta Theta = 30 deg. These data present a challenge to models of motion perception which do not explicitly consider the integration of information across the visual field.

Beutter, B. R.↗

DOE OSTI2025

$2+1$ dimensional Floquet systems and lattice fermions: Exact bulk spectral equivalence

A connection has recently been proposed between periodically driven systems known as Floquet insulators in continuous time and static fermion theories in discrete time. This connection has been established in a (1+1) ( 1 + 1 ) -dimensional free theory, where an explicit mapping between the spectra of a Floquet insulator and a discrete-time Dirac fermion theory has been formulated. Here we investigate the potential of static discrete-time theories to capture Floquet physics in higher dimensions, where so-called anomalous Floquet topological insulators can emerge that feature chiral edge states despite having bulk bands with zero Chern number. Starting from a particular model of an anomalous Floquet system, we provide an example of a static discrete-time theory whose bulk spectrum is an exact analytic match for the Floquet spectrum. The spectra with open boundary conditions in a particular strip geometry also match up to finite-size corrections. However, the models differ in several important respects. The discrete-time theory is spatially anisotropic, so that the spectra do not agree for all lattice terminations, e.g. other strip geometries or on half spaces. This difference can be attributed to the fact that the static discrete-time model is quasi-one-dimensional in nature and therefore has a different bulk-boundary correspondence than the Floquet model.

Iadecola, Thomas (ORCID:0000000251456441)↗

NASA NTRS1985

Mixed models and reduction method for dynamic analysis of anisotropic shells

A time-domain computational procedure is presented for predicting the dynamic response of laminated anisotropic shells. The two key elements of the procedure are: (1) use of mixed finite element models having independent interpolation (shape) functions for stress resultants and generalized displacements for the spatial discretization of the shell, with the stress resultants allowed to be discontinuous at interelement boundaries; and (2) use of a dynamic reduction method, with the global approximation vectors consisting of the static solution and an orthogonal set of Lanczos vectors. The dynamic reduction is accomplished by means of successive application of the finite element method and the classical Rayleigh-Ritz technique. The finite element method is first used to generate the global approximation vectors. Then the Rayleigh-Ritz technique is used to generate a reduced system of ordinary differential equations in the amplitudes of these modes. The temporal integration of the reduced differential equations is performed by using an explicit half-station central difference scheme (Leap-frog method). The effectiveness of the proposed procedure is demonstrated by means of a numerical example and its advantages over reduction methods used with the displacement formulation are discussed.

Noor, A. K.↗

DOE OSTI2024

Large-scale spatially explicit analysis of carbon capture at cellulosic biorefineries

The large-scale production of cellulosic biofuels would involve spatially distributed systems including biomass fields, logistics networks and biorefineries. Better understanding of the interactions between landscape-related decisions and the design of biorefineries with carbon capture and storage (CCS) in a supply chain context is needed to enable efficient systems. Here we analyse the cost and greenhouse gas mitigation potential for cellulosic biofuel supply chains in the US Midwest using realistic spatially explicit land availability and crop productivity data and consider fuel conversion technologies with detailed CCS design for their associated CO 2 streams. Optimization methods identify trade-offs and design strategies leading to systems with attractive environmental and economic performance. Strategic and operational decisions depend on underlying spatial features and are sensitive to biofuel demand and CCS incentives. US CCS incentives neglect to motivate greenhouse gas mitigation from all supply chain emission sources, which leverage spatial interactions between CCS, electricity prices and the biomass landscape.

09 BIOMASS FUELS↗

NASA NTRS2005

The Use of Mascons to Resolve Time-Variable Gravity from GRACE

We have analyzed GRACE Level 1B data to resolve time-variable gravity using a local mascon approach. The spherical harmonic solutions released to date resolve the signal from surface hydrology over land areas at spatial scales of 750 to 1000 km over one month intervals [Wahr et al., 2004; Tapley et al., 2004]. In our local approach, we solve explicitly for the mass of water in surface blocks using only the KBRR data collected as GRACE overflies the region of interest. The local representation of gravity minimizes leakage of errors from other areas due to aliasing or mismodelling. In this paper, we report on the analysis of GRACE data from January 2003 through August 2004 over three regions: the Amazon, the Indian subcontinent, and the continental United States. We solve for mass change at 10-day intervals using 4 deg x 4 deg blocks. We give an overview of our latest results, and we present the results of error analyses, and comparisons to both hydrology models and in-situ data.

Lemoine, Frank G.↗

DOE OSTI2022

Modeling wave propagation in elastic solids via high-order accurate implicit-mesh discontinuous Galerkin methods

Here, a high-order accurate implicit-mesh discontinuous Galerkin framework for wave propagation in single-phase and bi-phase solids is presented. The framework belongs to the embedded-boundary techniques and its novelty regards the spatial discretization, which enables boundary and interface conditions to be enforced with high-order accuracy on curved embedded geometries. High-order accuracy is achieved via high-order quadrature rules for implicitly-defined domains and boundaries, whilst a cell-merging strategy addresses the presence of small cut cells. The framework is used to discretize the governing equations of elastodynamics, written using a first-order hyperbolic momentum-strain formulation, and an exact Riemann solver is employed to compute the numerical flux at the interface between dissimilar materials with general anisotropic properties. The space-discretized equations are then advanced in time using explicit high-order Runge–Kutta algorithms. Several two- and three-dimensional numerical tests including dynamic adaptive mesh refinement are presented to demonstrate the high-order accuracy and the capability of the method in the elastodynamic analysis of single- and bi-phases solids containing complex geometries.

42 ENGINEERING↗

DOE OSTI2026

A Multiscale Simulation Study of the Effects of Nantucket Island on Atmospheric Boundary Layer Flow

Multiscale atmospheric simulations were conducted to investigate how Nantucket Island modifies the marine atmospheric boundary layer (ABL) flow, how the mechanisms responsible are represented in mesoscale versus large-eddy simulation (LES) domains, and potential impacts on wind energy production downstream, in moderate and high wind speed cases. Here, to isolate the effects of the island, each case study was simulated (i) with and without island topography included and (ii) with observed versus spatially uniform sea surface temperature (SST) values, the latter removing the influence of spatial variability in SST, while also enabling the examination of how the island effects vary with near-surface thermodynamic stability over the surrounding ocean. Island effects were found to be strongest during daytime heating, which generates convective rolls that redistribute momentum vertically over the island and downstream. Colder SSTs and weaker winds strengthen these rolls and extend their influence. Flow statistics within a representative wind turbine rotor-swept area downstream show that while vertical momentum redistribution enhances near-surface wind speeds, those aloft decrease, resulting in minimal net change to rotor-equivalent wind speed and slightly reduced vertical shear. Comparisons of mesoscale and LES domains reveal that both capture the roughness-driven near-surface slowdown, but mesoscale models underestimate vertical momentum transport compared to LES, which explicitly resolves convective rolls. Increasing LES resolution beyond the grid spacing required to begin capturing the largest convective roll structures ∼ O(100) m yields minimal additional momentum transport but increases resolved turbulence kinetic energy, relevant for stress loading.

Atmosphere-land interaction↗

DOE OSTI2023

Cation–polymer interactions and local heterogeneity determine the relative order of alkali cation diffusion coefficients in PEGDA hydrogels

Current research efforts are focused on endowing polymer membranes with ion–ion selectivity by incorporating ion–polymer interactions into materials to bias the selective partitioning and or diffusivity of one species over another. However, little is known about the impact of such interactions on the mechanisms of ion transport. In this study, we probe the influence of cation–polymer interactions on cation, anion, and salt diffusivity in a model membrane material, poly(ethylene glycol) diacrylate (PEGDA) by modeling concentrated polyethylene oxide solutions via molecular dynamics simulations. These results are compared to published experimental data for LiCl, NaCl, and KCl diffusion in PEGDA. Experimentally, the order of salt and cation diffusion coefficients for LiCl, NaCl, and KCl deviate from the order in aqueous solutions. Here, simulations identify these deviations to arise from cation–polymer coordination in the membrane. Both the fraction of bound cations and the average binding lifetime increases with decreasing cation hydration free energy (moving down the alkali series), leading to different diffusivity trends in the membrane compared to solution. However, to recover the experimentally observed order of diffusivities cations and salt in our simulations, we needed to incorporate membrane heterogeneity explicitly via a polymer charge scaling procedure. Together, our results indicate that cation–polymer interactions, as well as spatial heterogeneity within the membrane, play a critical role in dictating the observed order of alkali cation and salt diffusion coefficients in membranes.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Critical Infrastructure Decision-Making under Long-Term Climate Hazard Uncertainty: The Need for an Integrated, Multidisciplinary Approach

U.S. critical infrastructure assets are often designed to operate for decades, and yet long-term planning practices have historically ignored climate change. With the current pace of changing operational conditions and severe weather hazards, research is needed to improve our ability to translate complex, uncertain risk assessment data into actionable inputs to improve decision-making for infrastructure planning. Decisions made today need to explicitly account for climate change – the chronic stressors, the evolution of severe weather events, and the wide-ranging uncertainties. If done well, decision making with climate in mind will result in increased resilience and decreased impacts to our lives, economies, and national security. We present a three-tier approach to create the research products needed in this space: bringing together climate projection data, severe weather event modeling, asset-level impacts, and contextspecific decision constraints and requirements. At each step, it is crucial to capture uncertainties and to communicate those uncertainties to decision-makers. While many components of the necessary research are mature (i.e., climate projection data), there has been little effort to develop proven tools for long-term planning in this space. The combination of chronic and acute stressors, spatial and temporal uncertainties, and interdependencies among infrastructure sectors coalesce into a complex decision space. By applying known methods from decision science and data analysis, we can work to demonstrate the value of an interdisciplinary approach to climate-hazard decision making for longterm infrastructure planning.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

NASA NTRS1988

Autocorrelation and regularization in digital images. I - Basic theory

Spatial structure occurs in remotely sensed images when the imaged scenes contain discrete objects that are identifiable in that their spectral properties are more homogeneous within than between them and other scene elements. The spatial structure introduced is manifest in statistical measures such as the autocovariance function and variogram associated with the scene, and it is possible to formulate these measures explicitly for scenes composed of simple objects of regular shapes. Digital images result from sensing scenes by an instrument with an associated point spread function (PSF). Since there is averaging over the PSF, the effect, termed regularization, induced in the image data by the instrument will influence the observable autocovariance and variogram functions of the image data. It is shown how the autocovariance or variogram of an image is a composition of the underlying scene covariance convolved with an overlap function, which is itself a convolution of the PSF. The functional form of this relationship provides an analytic basis for scene inference and eventual inversion of scene model parameters from image data.

Jupp, David L. B.↗