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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

Natural Hazard Risk Assessments at the Global Scale

Since 1990, natural hazards have led to over 1.6 million fatalities globally, and economic losses are estimated at an average of around USD 260–310 billion per year. The scientific and policy communities recognise the need to reduce these risks. As a result, the last decade has seen a rapid development of global models for assessing risk from natural hazards at the global scale. In this paper, we review the scientific literature on natural hazard risk assessments at the global scale, and we specifically examine whether and how they have examined future projections of hazard, exposure, and/or vulnerability. In doing so, we examine similarities and differences between the approaches taken across the different hazards, and we identify potential ways in which different hazard communities can learn from each other. For example, there are a number of global risk studies focusing on hydrological, climatological, and meteorological hazards that have included future projections and disaster risk reduction measures (in the case of floods), whereas fewer exist in the peer-reviewed literature for global studies related to geological hazards. On the other hand, studies of earthquake and tsunami risk are now using stochastic modelling approaches to allow for a fully probabilistic assessment of risk, which could benefit the modelling of risk from other hazards. Finally, we discuss opportunities for learning from methods and approaches being developed and applied to assess natural hazard risks at more continental or regional scales. Through this paper, we hope to encourage further dialogue on knowledge sharing between disciplines and communities working on different hazards and risk and at different spatial scales.

Philip J. Ward↗
NASA NTRS2019

Inter-annual variation of tropical cyclones simulated by GEOS-5 AGCM with modified convection scheme

The Goddard Earth Observing System version 5 (GEOS-5) global climate model with a 50-km horizontal resolution is forced by observed sea surface temperature (SST) to examine the fidelity of the seasonal-mean and inter-annual variation of tropical cyclones (TCs) in the western North Pacific (WNP) and the North Atlantic (NATL). The standard Relaxed Arakawa Schubert (RAS) deep convection scheme is modified to improve the representation of TCs, where the scheme implements a stochastic limit of the cumulus entrainment rate. The modification drives mid- and upper-tropospheric cooling and low- to mid-tropospheric drying in the background state, which tends to increase atmospheric instability. This enables the model to increase convective variability on an intra-seasonal timescale and improve the simulation of intense storms. Five-member ensemble runs with the modified RAS scheme for 12 years (1998–2009) exhibit realistic spatial distributions in the climatological-mean TC development area and their pathways over WNP and NATL. The GCM is able to reproduce the inter-annual variation of accumulated cyclone energy (ACE) by prescribing yearly varying observed SST even though the individual TC intensity is still underpredicted. A sensitivity of TC activity to El Niño–Southern Oscillation (ENSO) phase is also reproduced realistically over WNP in terms of the spatial pattern changes in the main development region and TC pathways. However, the model exhibits a notable deficiency in NATL in reproducing the observed inter-annual variation of TC activity and the sensitivity to the ENSO.

accumulated cyclone energy↗
DOE OSTI2020

Patchy snapshots of nuclear chain reactions

Stochastic fluctuations of the neutron population within a nuclear reactor are typically prevented by operating the core at a sufficiently high power. This regime, where the evolution of the neutron density is essentially deterministic, is key for automatic protection and safety systems to safely detect unwanted power excursions during an accident, and to rapidly initiate a reactor trip procedure in case it is needed. Recent works, supported by numerical simulations, have however reported that, for large reactors, the branching nature of the fission reactions might induce strongly non-Poissonian patterns in the neutron spatial distribution, and that stochastic fluctuations might still persist at reactor powers close to operating conditions (startup phase). An international program conducted by LANL, IRSN and CEA was therefore setup to experimentally detect and characterize such fluctuations and correlations. An experiment took place in 2017 at the Reactor Critical Facility (RCF) of the Rensselaer Polytechnic Institute (USA). In this paper we will report the main findings of this experimental program, supported by stochastic models and by the development of a dedicated high-fidelity Monte Carlo simulation code. We will in particular describe and explain the strong patchiness in neutron power distributions measured at the RCF, as well as a peculiar ‘blinking’ behavior of the core, and discuss the consequences of these findings on nuclear safety.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1991

Does ion tearing exist?

Theoretical models of reconnection in the geomagnetic tail and its relationship to the collisionless tearing mode are examined analytically. The focus is on attempts to remove the stabilizing effects of electron compressibility in the quasi-neutral sheet. It is shown that this goal can be achieved by introducing spatial diffusion of electrons across the magnetic-flux surfaces, but not by wave-turbulent pitch-angle scattering (Coroniti, 1980) or nonadiabatic stochastic first-invariant diffusion (Buechner and Zelenyi, 1988). Hence the ion tearing mode does not exist, but the electron tearing mode does.

Pellat, R.↗
DOE OSTI2026

Correlation-aware binning for small-angle neutron scattering via Gaussian-process inference

Binning in small-angle neutron scattering (SANS) is typically performed empirically, with fixed parameters chosen for convenience rather than statistical optimality. Such practices often fail to balance statistical precision and spatial resolution, leading to inconsistencies across instruments and datasets. Here we establish a correlation-aware framework that determines the optimal bin width from first principles by extending the classical Freedman–Diaconis (FD) rule to account for inter-bin correlations with a Gaussian process. In this formulation, the scattering intensity is treated as a smooth stochastic field whose statistical coherence is described by a covariance matrix. Analytical expressions of errors derived from this model yield closed-form criteria that separate the total deviation into contributions from counting noise, aliasing distortion and curvature-dependent correlation effects. Expressed in reduced variables, the resulting dimensionless error surface reveals a continuous transition from the uncorrelated FD regime to the correlation-dominated limit, providing a unified description of noise suppression and resolution control. Because the formulation depends only on the profile characteristics of scattering intensity I(Q), specifically its average intensity and first- and second-order derivatives, it applies generally to any SANS measurement regardless of sample, instrument or geometry. Experimental validation using small- and ultra-small-angle neutron scattering data confirms the predicted scaling behavior, demonstrating that correlation-aware inference systematically reduces mean-squared error and enables information-efficient reproducible data reduction across materials and instruments.

Tung, Chi-Huan [ORNL] (ORCID:0000000221972074)↗
DOE OSTI2022

A hierarchical stochastic modeling approach for representing point bar geometries and petrophysical property variations

The flow of fluids in point bars is affected by the existence of heterogeneities like shale drapes that are found on the surfaces of inclined heterolithic stratifications. In fact, these shale drapes can act as fluid flow baffles; therefore, developing a framework for modeling point bars and their associated heterogeneities is vital. In this study, a stochastic process-based modeling approach is presented for capturing the main point bar heterogeneities: accretion surfaces (i.e., the aerial heterogeneity) and inclined heterolithic stratifications (i.e., the vertical heterogeneity). The former was modeled using a sine-generation function and the latter, with a sigmoidal function after which they were combined into a 3D point bar model. To ensure proper modeling of petrophysical properties, we developed a more representative gridding scheme which generates curvilinear grids representative of the point bar geometry. This grid was then transformed into a rectilinear grid to allow for geostatistical simulation after which all petrophysical properties were mapped back into the original curvilinear grid. An essential element of this modeling approach is the stochastic representation of shale drapes at the interface between successive accretion surfaces. The workflow was tested using a real field dataset for the Cranfield field, Mississippi. The constructed model was then subjected to a flow simulation study mimicking a CO 2 storage scenario. Various sensitivities were simulated to evaluate the effect of heterogeneities on CO 2 flow within the point-bar. Results demonstrate the importance of representing point-bar related heterogeneity and the spatial distribution of shale drapes on CO 2 plume migration and storage.

58 GEOSCIENCES↗
DOE OSTI2023

Temperature uncertainty modelling with proxy structural data as geostatistical constraints for well siting: an example applied to Granite Springs Valley, NV, USA

Utilizing existing temperature and structural geology information around Granite Springs Valley, Nevada, we build 3D stochastic temperature models with the aims of evaluating the 3D uncertainty of temperature and choosing between candidate exploration well locations. The data used to support the modelling are measured temperatures and structural proxies from 3D geologic modelling (distance to fault, distance to fault intersections and terminations, Coulomb stress change and dilation tendency), the latter considered ‘secondary’ data. Two stochastic geostatistical techniques are explored for incorporating the structural proxies: cosimulation and local varying mean. With both the cosimulation and local varying mean methods, many equally-likely temperature models (i.e. realizations) are produced, from which temperature probability profiles are calculated at candidate well locations. To aid in choosing between the candidate locations, two quantities summarize the temperature probabilities: V prior and entropy. V prior quantifies the likelihood for economic temperatures at each candidate location, whereas entropy identifies where new information has the most potential to reduce uncertainty. In general, the cosimulation realizations have smoother spatial structure, and extrapolate high temperatures at candidate locations that are located along the direction of the longest spatial correlation, which are down dip from existing temperature logs. The smooth realizations result in tight temperature probability profiles that are easier to interpret, but they have unrealistic temperature reversals in some locations because of the dipping ellipsoid shape created and that the cosimulation technique does not enforce a conductive geothermal gradient as a baseline (i.e. linearly increasing temperature with depth). The local varying mean results produce realizations with more realistic geothermal gradients, with temperatures increasing downward since a depth-temperature relationship is included. However, because they have much noisier spatial nature compared to cosimulation, it is harder to interpret the temperature probability profiles. The different local varying mean results allow the geologist to determine which proxy (e.g. dilation v. distance to fault termination) should be used given the specific geothermal system. In general, V prior from local varying mean results identify locations that are close to high values for the structural proxies: areas with higher probabilities for higher temperatures. The entropy results identify where uncertainty is greatest and therefore new drilling information could be most useful. Though these techniques provide useful information, even when applied to areas of sparse data, our comparison of these two techniques demonstrates the need for new geothermal geostatistics techniques that combine the advantages of these two methods and that are tailored to the spatial uncertainty issues inherent in geothermal exploration.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Improving Solution Accuracy and Convergence for Stochastic Physics Parameterizations with Colored Noise

Stochastic parameterizations are used in numerical weather prediction and climate modeling to help capture the uncertainty in the simulations and improve their statistical properties. Convergence issues can arise when time integration methods originally developed for deterministic differential equations are applied naively to stochastic problems. In previous studies, it has been demonstrated that a correction term, known in stochastic analysis as the Itô correction, can help improve solution accuracy for various deterministic numerical schemes and ensure convergence to the physically relevant solution without substantial computational overhead. The usual formulation of the Itô correction is valid only when the stochasticity is represented by white noise. In this study, a generalized formulation of the Itô correction is derived for noises of any color. The formulation is applied to a test problem described by an advection–diffusion equation forced with a spectrum of fast processes. We present numerical results for cases with both constant and spatially varying advection velocities to show that, for the same time step sizes, the introduction of the generalized Itô correction helps to substantially reduce time integration error and significantly improve the convergence rate of the numerical solutions when the forcing term in the governing equation is rough (fast varying); alternatively, for the same target accuracy, the generalized Itô correction allows for the use of significantly longer time steps and, hence, helps to reduce the computational cost of the numerical simulation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Go with the Radflow: Thermal radiation experiments in the high energy-density regime [Slides]

Los Alamos National Laboratory is a nuclear weapons laboratory supporting our nation's defense. In support of the mission is a high energy-density physics program in which we design and execute experiments to study radiation-hydrodynamics phenomena and improve the predictive capability of our large-scale multi-physics software codes on our supercomputers. The Radflow Project is maturing a unique spectroscopic measurement that yields a spatially dependent supersonic radiation wavefront profile. The spectroscopic measurement is much more constraining than the typically used radiography, which measures the density changes due to the shock as the radiation wave goes subsonic. Recently, the spectroscopic measurement has been applied to advanced targets with a detailed heterogeneity representing a single realization of a stochastic medium. The Radflow Project also conducts opacity experiments to explain the current conundrum of the experimental iron opacity differing from theory, which is important in modeling the sun. More than that, though, we need to be able to answer these science questions and be able to validate our theory and codes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Evaluation of Two-Dimensional to One-Dimensional Site Response for Idaho National Laboratory

We perform two-dimensional (2D) site response analyses accounting for spatial variability of soil properties and subsurface geometry of the Eastern Snake River Plane (ESRP), and quantify their effects on ground surface motion relative to one-dimensional (1D) site response analyses at the Idaho National Laboratory (INL). We first present the development of random field idealizations of the repeated basalt lava flows, heterogeneously inter-layered with sediments, from seismic velocity data collected over four decades in the ESRP. Using realizations of the stochastic fields mapped on 2D deterministic finite element models, we perform 2D viscoelastic and equivalent-linear wave propagation simulations, and quantify the mean and variance of site response aggravation factors, defined as the response spectral ratio of 2D to 1D analyses on the ground surface. Results are shown to be insensitive to the constitutive material behavior considered here, for strains induced by rock outcrop peak ground acceleration (PGA) as high as 0.7g: viscoelastic and equivalent-linear analyses predict peak mean 2D/1D aggravation factor 1.05 at period T=0.075 sec (i.e. the 2D response spectrum is 5% higher than the corresponding 1D at that period, on average), which corresponds to the wavelength of the horizontal correlation length of the random field (50m). For periods longer than the fundamental period of the site (here, T 1 =0.3125 sec), the propagating wavelengths are too long to be affected by the 1D site response and the aggravation factor becomes equal to 1. The standard deviation of the natural logarithms of the 2D/1D aggravation factors is ~0.15 for periods shorter than the fundamental period of the site, and decays thereafter at a steady rate.

58 GEOSCIENCES↗
DOE OSTI2022

Fracton Hydrodynamics without Time-Reversal Symmetry

Here, we present an effective field theory for the nonlinear fluctuating hydrodynamics of a single conserved charge with or without time-reversal symmetry, based on the Martin-Siggia-Rose formalism. Applying this formalism to fluids with only charge and multipole conservation, and with broken time-reversal symmetry, we predict infinitely many new dynamical universality classes, including some with arbitrarily large upper critical dimensions. Using large scale simulations of classical Markov chains, we find numerical evidence for a breakdown of hydrodynamics in quadrupole-conserving models with broken time-reversal symmetry in one spatial dimension. Our framework can be applied to the hydrodynamics around stationary states of open systems, broadening the applicability of previously developed ideas and methods to a wide range of systems in driven and active matter.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1984

The spatial structure of terrain - A process signal in satellite digital images

Pattern recognition procedures applied to Landsat imagery carry an implicit assumption that the digital data are independently distributed. That assumption is incorrect over virtually any terrain. Deviations from independence occur because slopes follow a systematic pattern of variation arising from the slope-forming processes. That pattern can be identified using the stochastic process methodology of Box and Jenkins.Angles of adjacent slopes are autocorrelated and the bidirectional reflectance function transfers these systematic slope changes to the sensor. Imagery becomes autocorrelated through this transfer. Autocorrelation in the imagery can be removed through direct calculation from a digital elevation model or by use of stochastic process methodology. The latter has the advantage that the residuals are white noise; and it is applicable in any area, even where a D.E.M. is unavailable. The stochastic process signal can be used to study terrain processes.

Craig, R. G.↗
DOE OSTI2023

Compound Continental Risk of Multiple Extreme Floods in the United States

Abstract Understanding spatially correlated floods and modeling joint hazard associated with threshold exceedances across multiple locations is crucial for accurate estimation of continental‐scale portfolio risk. This work uses a non‐parametric copula‐based spatial simulator to analyze peak floods across the United States to derive the first‐of‐its‐kind continental portfolio risk estimates at the 10‐ and 100‐year return levels. We find significant interdependence in floods across the nation, revealing the recurring pattern of extreme events affecting the Northeast, Central, West, and Northwest United States in the same year. The stochastic simulator effectively manages high‐dimensional data and offers reliable uncertainty estimates for both spatially dependent floods and the aggregated flood losses at the continental level. El Niño‐Southern Oscillation and Atlantic Multidecadal Oscillation are identified as statistically significant tele‐connectors of aggregate loss. This research aims to advance the understanding of compound continental flood hazard and the potential large‐scale climate teleconnections that lead to such compound floods.

Geology↗
DOE OSTI2020

Grid-based diffusion Monte Carlo for fermions without the fixed-node approximation

A diffusion Monte Carlo algorithm is introduced that can determine the correct nodal structure of the wave function of a few-fermion system and its ground-state energy without an uncontrolled bias. This is achieved by confining signed random walkers to the points of a uniform infinite spatial grid, allowing them to meet and annihilate one another to establish the nodal structure without the fixed-node approximation. An imaginary-time propagator is derived rigorously from a discretized Hamiltonian, governing a non-Gaussian, sign-flipping, branching, and mutually annihilating random walk of particles. The accuracy of the resulting stochastic representations of a fermion wave function is limited only by the grid and imaginary-time resolutions and can be improved in a controlled manner. Here, the method is tested for a series of model problems including fermions in a harmonic trap as well as the He atom in its singlet or triplet ground state. For the latter case, the energies approach from above with increasing grid resolution and converge within 0.015 E h of the exact basis-set-limit value for the grid spacing of 0.08 a.u. with a statistical uncertainty of 10 –5 E h without an importance sampling or Jastrow factor.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1988

Use of LANDSAT images of vegetation cover to estimate effective hydraulic properties of soils

The estimation of the spatially variable surface moisture and heat fluxes of natural, semivegetated landscapes is difficult due to the highly random nature of the vegetation (e.g., plant species, density, and stress) and the soil (e.g., moisture content, and soil hydraulic conductivity). The solution to that problem lies, in part, in the use of satellite remotely sensed data, and in the preparation of those data in terms of the physical properties of the plant and soil. The work was focused on the development and testing of a stochastic geometric canopy-soil reflectance model, which can be applied to the physically-based interpretation of LANDSAT images. The model conceptualizes the landscape as a stochastic surface with bulk plant and soil reflective properties. The model is particularly suited for regional scale investigations where the quantification of the bulk landscape properties, such as fractional vegetation cover, is important on a pixel by pixel basis. A summary of the theoretical analysis and the preliminary testing of the model with actual aerial radiometric data is provided.

Eagleson, Peter S.↗
NASA NTRS2012

Upper-Ocean Heat Balance Processes and the Walker Circulation in CMIP5 Model Projections

Considerable uncertainty remains as to the importance of mechanisms governing decadal and longer variability of the Walker Circulation, its connection to the tropical climate system, and prospects for tropical climate change in the face of anthropogenic forcing. Most contemporary climate models suggest that in response to elevated CO2 and a warmer but more stratified atmosphere, the required upward mass flux in tropical convection will diminish along with the Walker component of the tropical mean circulation as well. Alternatively, there is also evidence to suggest that the shoaling and increased vertical stratification of the thermocline in the eastern Pacific will enable a muted SST increase there-- preserving or even enhancing some of the dynamical forcing for the Walker cell flow. Over the past decade there have been observational indications of an acceleration in near-surface easterlies, a strengthened Pacific zonal SST gradient, and globally-teleconnected dislocations in precipitation. But is this evidence in support of an ocean dynamical thermostat process posited to accompany anthropogenic forcing, or just residual decadal fluctuations associated with variations in warm and cold ENSO events and other stochastic forcing? From a modeling perspective we try to make headway on this question by examining zonal variations in surface energy fluxes and dynamics governing tropical upper ocean heat content evolution in the WCRP CMIP5 model projections. There is some diversity among model simulations; for example, the CCSM4 indicates net ocean warming over the IndoPacific region while the CSIRO model concentrates separate warming responses over the central Pacific and Indian Ocean regions. The models, as with observations, demonstrate strong local coupling between variations in column water vapor, downward surface longwave radiation and SST; but the spatial patterns of changes in the sign of this relationship differ among models and, for models as a whole, with observations. Our analysis focuses initially on probing the inter-model differences in energy fluxes / transports and Walker Circulation response to forcing. We then attempt to identify statistically the El Nino- / La Nina-related ocean heat content variability unique to each model and regress out the associated energy flux, ocean heat transport and Walker response on these shorter time scales for comparison to that of the anthropogenic signals.

Robertson, F. R.↗
DOE OSTI2023

Temperature Uncertainty Modeling with Proxy Structural Data as Geostatistical Constraints for Well Siting: An Example Applied to Granite Springs Valley, NV, USA

Utilizing existing temperature and structural information around Granite Springs Valley, Nevada, we build 3D stochastic temperature models with the aim of evaluating the 3D uncertainty of temperature and choosing between candidate exploration well locations . The data used to support the modeling are measured temperatures and structural proxies from 3D geologic modeling, the latter considered "secondary" data. Two stochastic geostatistical techniques are explored for incorporating the structural proxies: cosimulation and local varying mean. With both the cosimulation and local varying mean methods, many equally likely temperature models (i.e., realizations) are produced, from which temperature probability profiles are calculated at candidate well locations. To aid in choosing between the candidate locations, two quantities summarize the temperature probabilities: Vprior and entropy. Vprior quantifies the likelihood for economic temperatures at each candidate location, whereas entropy identifies where new information has the most potential to reduce uncertainty. In general, the cosimulation realizations have smoother spatial structure, and extrapolate high temperatures at candidate locations that are located along the direction of the longest spatial correlation, which are down dip from existing temperature logs. The smooth realizations result in tight temperature probability profiles that are easier to interpret, but they have unrealistic temperature reversals in some locations because the cosimulation technique does not enforce a conductive geothermal gradient as a baseline (i.e., linearly increasing temperature with depth). The local varying mean results produce realizations with more realistic geothermal gradients, with temperatures increasing downward since a depth-temperature relationship is included. However, because they have much noisier spatial nature compared to cosimulation, it is harder to interpret the temperature probability profiles. The different local varying mean results allow the geologist to determine which proxy (e.g., dilation versus distance to fault termination) should be used given the specific geothermal system. In general, Vprior from local varying mean results identify locations that are close to high values for the structural proxies: areas with highe r probabilities for higher temperatures. The entropy results identify where uncertainty is greatest and therefore new drilling information could be most useful. Though these techniques provide useful information, even when applied to areas of sparse data, our comp arison of these two techniques demonstrates the need for new geothermal geostatistics techniques that combine the advantages of these two methods and that are tailored to the spatial uncertainty issues inherent in geothermal exploration.

3D temperature modeling↗
DOE OSTI2021

Investigation of turbulent inflow specification in Euler–Lagrange simulations of mid-field spray

The process of atomization of a liquid jet by a parallel high-speed gas stream results in a spray, whose downstream development is of considerable interest to several applications. The round jet spray can be spatially divided into (i) a near-field (near-nozzle) region of liquid atomization and (ii) a downstream mid-field region of fully-dispersed droplets. In order to accurately model mid-field droplet dispersion, this work aims at developing a rigorous and robust injection model for Euler–Lagrange spray simulations. Results from experiments are used to obtain the relevant droplet number density, size distribution, and mean and standard deviation velocity distributions of the injection model, systematically in a step-by-step process. Two-phase large eddy simulations are performed by stochastically generating the Lagrangian droplets at the inlet of the mid-field region. Number flux, diameter distribution, mean velocity, and other time-averaged statistics at several downstream locations are shown to agree well with the corresponding experimental data.

42 ENGINEERING↗