Search NASA⌕ Search

SEARCH · Search NASA

Results for “spatial stochastic 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.

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

At least 73 records · Page 4

Wake-Resolving Acoustic Tomography: Advances through Numerical Covariance Methods

Acoustic tomography offers path-integrated measurements of atmospheric velocity and temperature fluctuations with high spatial resolution. Classical implementations of time-dependent stochastic inversion rely on homogeneous, isotropic covariance models that are poorly suited to the anisotropic structure of wind turbine wakes. By directly estimating heterogeneous covariances from large-eddy simulations (LESs) into the time-dependent stochastic inversion operator, we relax implicit assumptions in the analytical models used historically. Retrievals using these LES-informed models improve agreement with true fields in variance, turbulent kinetic energy, and spectral content compared to analytical and precursor-based covariance models. The results indicate that LES-informed covariance models can enhance the accuracy of acoustic tomography retrievals in complex, anisotropic flows such as wind turbine wakes in some cases and highlight instances where analytical models still offer competitive performance, despite their simplifying assumptions.

17 WIND ENERGY↗

Critical Slowing Down in Time-to-Extinction: An Example of Critical Phenomena in Ecology

We study a model for two competing species that explicitly accounts for effects due to discreteness, stochasticity and spatial extension of populations. The two species are equally preferred by the environment and do better when surrounded by others of the same species. We observe that the final outcome depends on the initial densities (uniformly distributed in space) of the two species. The observed phase transition is a continuous one and key macroscopic quantities like the correlation length of clusters and the time-to-extinction diverge at a critical point. Away from the critical point, the dynamics can be described by a mean-field approximation. Close to the critical point, however, there is a crossover to power-law behavior because of the gross mismatch between the largest and smallest scales in the system. We have developed a theory based on surface effects, which is in good agreement with the observed behavior. The course-grained reaction-diffusion system obtained from the mean-field dynamics agrees well with the particle system.

Gandhi, Amar↗

Stochastic Models of Tropical Rain-Rate Fields

Because of the extreme variability of rain rate in space and time and the difficulties with remote sensing methods of measuring rain rates, accurate determination of rainfall over large areas and time periods has long been a problem for hydrologists, meteorogists, and climatologists. A number of statistical models of rain have been developed in order to investigate the impact of rain variability on satellite remote sensing methods, validation of satellite rain products, and generation of rain maps with accompanying error estimates. These models may be useful in examining 'sub-grid scale' issues in representing precipitation in numerical mdoels. A stochastic model will first be described which can generate time-dependent high-resolution spatial rain fields with space and time correlations similar to those seen in rain data, as well as representing the presence of areas with zero rain rate and log-normally distributed rain rates where there is rain. A simpler model derived from this, formulated in the spectral domain, seems to imply fractal-like rain statistics at small scales when fit to rain data.

Bell, Thomas L.↗

Visualizing Distributions from Multi-Return Lidar Data to Understand Forest Structure

Spatially distributed probability density functions (pdfs) are becoming relevant to the Earth scientists and ecologists because of stochastic models and new sensors that provide numerous realizations or data points per unit area. One source of these data is from multi-return airborne lidar, a type of laser that records multiple returns for each pulse of light sent towards the ground. Data from multi-return lidar is a vital tool in helping us understand the structure of forest canopies over large extents. This paper presents several new visualization tools that allow scientists to rapidly explore, interpret and discover characteristic distributions within the entire spatial field. The major contribution from-this work is a paradigm shift which allows ecologists to think of and analyze their data in terms of the distribution. This provides a way to reveal information on the modality and shape of the distribution previously not possible. The tools allow the scientists to depart from traditional parametric statistical analyses and to associate multimodal distribution characteristics to forest structures. Examples are given using data from High Island, southeast Alaska.

Kao, David L.↗

Multivariate space - time analysis of PRE-STORM precipitation

This paper presents the methodologies and results of the multivariate modeling and two-dimensional spectral and correlation analysis of PRE-STORM rainfall gauge data. Estimated parameters of the models for the specific spatial averages clearly indicate the eastward and southeastward wave propagation of rainfall fluctuations. A relationship between the coefficients of the diffusion equation and the parameters of the stochastic model of rainfall fluctuations is derived that leads directly to the exclusive use of rainfall data to estimate advection speed (about 12 m/s) as well as other coefficients of the diffusion equation of the corresponding fields. The statistical methodology developed here can be used for confirmation of physical models by comparison of the corresponding second-moment statistics of the observed and simulated data, for generating multiple samples of any size, for solving the inverse problem of the hydrodynamic equations, and for application in some other areas of meteorological and climatological data analysis and modeling.

Polyak, Ilya↗

Mechanical model of an earthquake fault

The dynamic behavior of a simple mechanical model of an earthquake fault is studied. This model, introduced originally by Burridge and Knopoff (1967), consists of an elastically coupled chain of masses in contact with a moving rough surface. The present version of the model retains the full Newtonian dynamics with inertial effects and contains no externally imposed stochasticity or spatial inhomogeneity. The only nonlinear feature is a velocity-weakening stick-slip friction force between the masses and the moving surface. This system is being driven persistently toward a slipping instability and therefore exhibits noisy sequences of earthquakelike events. These events are observed in numerical simulations, and many of their features can be predicted analytically.

Carlson, J. M.↗

High-LET Patterns of DSBs in DNA Loops, the HPRT Gene and Phosphorylation Foci

We present new results obtained with our model based on the track structure and chromatin geometry that predicts the DSB spatial and genomic distributions in a cell nucleus with the full genome represented. The model generates stochastic patterns of DSBs in the physical space of the nucleus filled with the realistic configuration of human chromosomes. The model was re-used to find the distribution of DSBs in a physical volume corresponding to a visible phosphorylation focus believed to be associated with a DSB. The data shows whether there must more than one DSB per foci due to finite size of the visible focus, even if a single DSB is radiochemically responsible for the phosphorylation of DNA in its vicinity. The same model can predict patterns of closely located DSBs in a given gene, or in a DNA loop, one of the large-scale chromatin structures. We demonstrated for the example of the HPRT gene, how different sorts of radiation lead to proximity effect in DSB locations, which is important for modeling gene deletions. The spectrum of intron deletions and total gene deletions was simulated for the HPRT gene. The same proximity effect of DSBs in a loop can hinder DSB restitutions, as parts of the loop between DSBs is deleted with a higher likelihood. The distributions of DSBs and deletions of DNA in a loop are presented.

Ponomarev, Artem L.↗

The Elastic Response of the Earth to Interannual Variations in Antarctic Precipitation

Measurements of elastic displacements of the bedrock surrounding large ice sheets have been proposed as a means to detect mass changes in these ice sheets. However, accumulation of glacial mass on the ice sheets is a noisy process, subject to large spatial and temporal variations in precipitation. We simulated the response of the Antarctic continent to a stochastic model of interannual precipitation variations and found that interannual variations in the elastic response of the earth are large when compared to the long-term mean of displacements produced by an assumed average ice mass imbalance of 10%. If, as some scientists predict, Antarctic ice mass changes in the future become dramatic, the long-term signal should be large enough to be detected by a few years of geodetic measurements, despite climatic noise.

Hager, Bradford H.↗

Image Discrimination Models With Stochastic Channel Selection

Many models of human image processing feature a large fixed number of channels representing cortical units varying in spatial position (visual field direction and eccentricity) and spatial frequency (radial frequency and orientation). The values of these parameters are usually sampled at fixed values selected to ensure adequate overlap considering the bandwidth and/or spread parameters, which are usually fixed. Even high levels of overlap does not always ensure that the performance of the model will vary smoothly with image translation or scale changes. Physiological measurements of bandwidth and/or spread parameters result in a broad distribution of estimated parameter values and the prediction of some psychophysical results are facilitated by the assumption that these parameters also take on a range of values. Selecting a sample of channels from a continuum of channels rather than using a fixed set can make model performance vary smoothly with changes in image position, scale, and orientation. It also facilitates the addition of spatial inhomogeneity, nonlinear feature channels, and focus of attention to channel models.

Ahumada, Albert J., Jr.↗

Combining model-based and model-free methods for stochastic control of distributed energy resources

Modern distribution systems are experiencing a fast transformation with the growing penetration of distributed energy resources (DERs). Along with the economic and environmental benefits of DERs, challenges arise to address the uncertainties caused by their inherent volatility. If properly coordinated, however, DERs have the potential to provide the controllability that grid operators need. Here, we propose a hierarchical control framework that combines the model-based and model-free methods for stochastic DER control in distribution systems. The upper-level scheduler considers a chance-constrained optimal power flow problem (model-based) that schedules DER setpoints to minimize the operational cost and maintain the operating reserve. The lower-level distributed DER controllers absorb real-time disturbances and uncertainties using the extremum seeking control (model-free) to achieve grid objectives. The combination of model-based and model-free methods allows us to take the advantages of both methods to effectively manage the uncertainty in distribution systems. The proposed work is demonstrated on the IEEE 13-node feeder.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Alternative Representations of Convective Processes in the NASA GEOS-5 AGCM

The gap in explicit resolution of phenomena between global climate models and cloud resolving models is shrinking at a steady pace with global integrations of several km in resolution now practical for at least some time scales. In moving toward finer resolution the long standing problem of convective parameterization is being examined along with the assumption of convective quasi-equilibrium and how this can be reconciled with the stochastic and intermittent nature of convection. In this context we examine the nature of parameterized convection in the NASA Goddard Earth Observing System (GEOS-5) Atmospheric General Circulation Model. Our analysis uses both coarse (2.5 degree) and fine scale 0.25 degree spatial resolution integrations. Two basic formulations are compared: the default option is the Relaxed Arakawa-Schubert (RAS) scheme which invokes a sequence of linearly entraining plumes and quasi-equilibrium closure. An optional modification of this method (the "Stochastic Tokioka" constraint) places a random lower limit on plume entrainment. An alternative representation is the Kain-Fritsch parameterization which was originally developed for mesoscale numerical modeling strategies. Here entrainment is determined by a crude buoyancy sorting approach that allows the plume spectrum to be more responsive to ambient vertical stratification of moisture. Diagnostics of the model behavior are referenced to recent observational evidence of continuous phase transition behavior. In particular we examine the relationship between column water vapor and probablility of convective presence and intensity. Sensitivity of the statistics of convective behavior to parcel mixing/entrainment formulations and parcel initial thermodynamics are considered. Observational statistics from A-Train and TRMM sensors provide validation of the model integrations.

Robertson, Franklin↗

SSRS (Stochastic Soaring Raptor Simulator)

SSRS (Stochastic Soaring Raptor Simulator) is a generalizable, probabilistic, and predictive tool for wind energy developers, ecologists, wildlife managers and industry consultants to estimate the potential for soaring raptors to interact with operating wind turbines, without the need for site-specific data collection. Rapid expansion of wind energy development across the world has exposed the risk of turbine collisions for birds and bats. The risk to obligate soaring raptors such as golden eagles is of particular concern due to their small population and influence on ecological balance. Golden eagles rely heavily on updrafts to subsidize their flight, putting them in direct conflict with operational wind turbines that utilize the same wind resource. Understanding the behavior of soaring raptors with varying atmospheric conditions is crucial for predicting and mitigating the risk of turbine collision. This software contains a predictive movement model that simulates individual flight paths of golden eagles during updraft-subsidized long-distance flight, including migration. For a given set of atmospheric conditions, the model simulates thousands of eagles at turbine-scale spatial resolution (50m) to produce a relative presence density map. The simulated eagles rely on updrafts to pursue uninterrupted directional flight with minimal energy expenditure, following fluid-flow principles. The simulator includes a stochastic model of eagle behavior and a systematic method of accounting for spatiotemporal variations in atmospheric conditions. This framework only requires publicly available atmospheric data to estimate orographic and thermal updrafts, ensuring general usability.

Sandhu, Rimple↗

Influence of spatial variability of hydraulic characteristics of soils on surface parameters obtained from remote sensing data in infrared and microwaves

The correct interpretation of thermal and hydraulic soil parameters infrared from remotely sensed data (thermal infrared, microwaves) implies a good understanding of the causes of their temporal and spatial variability. Given this necessity, the sensitivity of the surface variables (temperature, moisture) to the spatial variability of hydraulic soil properties is tested with a numerical model of heat and mass transfer between bare soil and atmosphere. The spatial variability of hydraulic soil properties is taken into account in terms of the scaling factor. For a given soil, the knowledge of its frequency distribution allows a stochastic use of the model. The results are treated statistically, and the part of the variability of soil surface parameters due to that of soil hydraulic properties is evaluated quantitatively.

Brunet, Y.↗

Modeling information flow in a computer processor with a multi-stage queuing model

In this paper, we introduce a nonlinear stochastic model to describe the propagation of information inside a computer processor. In this model, a computational task is divided into stages, and information can flow from one stage to another. The model is formulated as a spatially-extended, continuous-time Markov chain where space represents different stages. This model is equivalent to a spatially-extended version of the M/M/s queue. The main modeling feature is the throttling function which describes the processor slowdown when the amount of information falls below a certain threshold. We derive the stationary distribution for this stochastic model and develop a closure for a deterministic ODE system that approximates the evolution of the mean and variance of the stochastic model. In conclusion, we demonstrate the validity of the closure with numerical simulations.

97 MATHEMATICS AND COMPUTING↗

Stochastic model theory of broadband shock associated noise from supersonic jets

A method based on the work of Tam and Tanna (1982) for calculating the near field noise spectrum and the spatial distribution of broadband shock associated noise from supersonic jets is proposed. Multiple-scales expansion is used to decompose the quasi-periodic shock cells into time-independent waveguide modes of the jet flow, and the interaction of the instability waves with each of the waveguide modes is shown to generate unsteady disturbances which become part of the broadband shock associated noise when radiated to the far field. The observed broadband shock associated noise is composed of a superposition of the various distinct spectra of the different waveguide modes, and the multispectra can be easily identified in many of the existing far and near field noise measurements.

Tam, C. K. W.↗

Scalable computations for nonstationary Gaussian processes

Nonstationary Gaussian process models can capture complex spatially varying dependence structures in spatial datasets. However, the large number of observations in modern datasets makes fitting such models computationally intractable with conventional dense linear algebra. In addition, derivative-free or even first-order optimization methods can be very slow to converge when estimating many spatially varying parameters. In this paper, we present a computational framework which couples an algebraic block diagonal plus low-rank covariance matrix approximation with stochastic trace estimation to facilitate the efficient use of second-order solvers for maximum likelihood estimation of Gaussian process models with many parameters. We demonstrate the effectiveness of these methods by simultaneously fitting 192 parameters in the popular nonstationary model of Paciorek and Schervish using 107,600 sea surface temperature anomaly measurements.

97 MATHEMATICS AND COMPUTING↗

Derivation of heterogeneous material distributions and their sensitivity to HM-coupled two-phase flow models exemplified with the LASGIT experiment

Abstract Advective gas transport in bentonite, a possible buffer material in repositories for radioactive materials, is difficult to simulate in numerical continuum models, partly due to the complicated microstructure of bentonite. To generate reliable models of repositories nevertheless, spatially distributed heterogeneous material properties can be used to allow localization of gas flow. In this study, a pore-size-dependent stochastic approach of the gas entry pressure is derived from Mercury Intrusion Porosimetry, which is used to replicate measurements from the LASGIT experiment. In addition, three benchmark tests are simulated to investigate the dependence of heterogeneous distributions of material properties on the mesh discretization, the temporal dependence, and the coupling between the processes influenced by the heterogeneous parameters. The numerical modeling results of the LASGIT experiment show that the onset of gas flow into the system and the subsequent increase in pressure and stress can be well reproduced using heterogeneous distributions. Compared to a model with homogeneous material properties, heterogeneous distributions may allow the generation of dilatancy-controlled microfractures—an important feature with regard to the advective gas flow in bentonites. However, it can be observed that the heterogeneous distributions in LASGIT are less significant, as technical gaps or differences in material types could have a greater impact.

Environmental Sciences & Ecology↗

Environmental factors driving spatial heterogeneity in desert halophile microbial communities

Spatial heterogeneity in microbial communities is observed in all-natural ecosystems and can stem from both adaptations to local environmental conditions as well as stochastic processes. Extremophile microbial communities inhabiting evaporitic halite nodules (salt rocks) in the Atacama Desert, Chile, are a good model ecosystem for investigating factors leading to microbiome heterogeneity, due to their diverse taxonomic composition and the spatial segregation of individual nodules. We investigated the abiotic factors governing microbiome composition across different spatial scales, allowing for insight into the factors that govern halite colonization from regional desert-wide scales to micro-scales within individual nodules. We found that water availability and community drift account for microbiome assembly differently at different distance scales, with higher rates of cell dispersion at the smaller scales resulting in a more homogenous composition. This trend likely applies to other endoliths, and to non desert communities, where dispersion between communities is limited. At the intra-nodule scales, a light availability gradient was most important in determining the distribution of microbial taxa despite intermixing by water displacement via capillary action.

microbiome extremophile, desert, heterogeneity, Me↗