Search NASA⌕ Search

SEARCH · Search NASA

Results for “first passage problems”

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.

10 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Nucleation-mediated reshaping of facetted metallic nanocrystals: Breakdown of the classical free energy picture

Shape stability is key to avoiding degradation of performance for metallic nanocrystals synthesized with facetted non-equilibrium shapes to optimize properties for catalysis, plasmonics, and so on. Reshaping of facetted nanocrystals is controlled by the surface diffusion-mediated nucleation and growth of new outer layers of atoms. Here, Kinetic Monte Carlo (KMC) simulation of a realistic stochastic atomistic-level model is applied to precisely track the reshaping of Pd octahedra and nanocubes. Unexpectedly, separate constrained equilibrium Monte Carlo analysis of the free energy profile during reshaping reveals a fundamental failure of the classical nucleation theory (CNT) prediction for the reshaping barrier and rate. Why? Nucleation barriers can be relatively low for these processes, so the system is not locally equilibrated before crossing the barrier, as assumed in CNT. This claim is supported by an analysis of a first-passage problem for reshaping within a master equation framework for the model that reasonably captures the behavior in KMC simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Langevin and Fokker-Planck analyses for diffusion-mediated passing of circular and discorectangular species in two-dimensional channels

The propensity for pairs of diffusing species to pass each other within narrow channels or pores is of basic interest as a first-passage-type problem. It is also of relevance for solution-phase transport in nanoporous materials, and in particular for catalytic conversion reactions where high yield requires that product species can efficiently pass reactant species to exit the pores. Here, we analyze a two-dimensional model with nonoverlapping circular and discorectangular species confined to a rectangular channel, and where passing is mediated by Brownian dynamics in an implicit solvent. For narrower channels where passing is still possible, the discorectangle must align with the channel to pass the circular species. Behavior of the passing propensity, 𝑃, can be assessed by strongly damped Langevin simulations, or within an equivalent Fokker-Planck equation (FPE) formalism. The latter corresponds to a diffusion problem in a “higher-dimensional channel” with a constriction. We assess the variation of the passing propensity, 𝑃, for a broad range of channel width including its scaling just above the threshold where passing is sterically blocked. Analysis of 𝑃 versus the rotational diffusion coefficient 𝐷 𝑟 of the discorectangle reveals a significant decrease in 𝑃 for lower 𝐷 𝑟 for moderate channel width. This prompts a direct analysis of the regime where 𝐷 𝑟 → 0 , for which the FPE can be reduced to a three-dimensional diffusion problem, precise analysis of which is facilitated by adaptive-mesh finite element methods. The dependence of 𝑃 on the aspect ratio of the discorectangle is also assessed.

Rahman, Md Khaledur [Ames Laboratory (AMES), Ames,↗
DOE OSTI2022

Binding kinetics of harmonically confined random walkers

Diffusion-mediated binding of molecules under the influence of discrete spatially confining potentials is a commonly encountered scenario in systems subjected to explicit fields or implicit fields arising from tethering restraints. Here, in this work, we derive analytical expressions for the mean binding time of two random walkers geometrically confined by means of two harmonic potentials in one- and two-dimensional systems, which show excellent agreement with Brownian dynamics simulations. As a demonstration of its utility, we use this theory to maximize the communication speed in existing DNA walkers, obtaining quantitative agreement with previously reported experimental findings. The analytical expressions derived in this paper are broadly applicable to diverse systems, providing ways to characterize communication processes and optimize the rate of signal propagation for sensing and computing applications at the nanoscale.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Identification of hot water end-use process of electric water heaters from energy measurements

This paper presents an algorithm for the identification of parameters for a stochastic hot water end-use process that drives a homogeneous population of thermostatically controlled electric water heaters (EWH). Usually, only metered interval consumption data (kWh) is collected and the hot water end-use process is unobservable to utility and aggregators. However, the availability of EWHs for demand response (DR) is closely coupled with the hot water end-use process. In this context, the hot water end-use process is modeled as a two-state Markov chain (Use / No use), which causes the thermostatic ON-OFF switching process to behave as a Markov renewal process (MRP). A set of first passage-time problems is developed to obtain the moments of the transition probability densities of the MRP. These problems are addressed by establishing a system of coupled partial differential equations characterizing the temperature evolution of the EWH population. A key quantity in the methodology for estimating the parameters is the total time an EWH is ON within a period of interest. It is referred to as the total busy time. Total busy time in this approach is a random variable for which analytical expressions of the moments are developed as a function of the metered window length. The latter expressions become the basis of a hot water demand model identification algorithm which is validated using agent-based simulations of EWHs.

42 ENGINEERING↗
DOE OSTI2022

First-passage time statistics on surfaces of general shape: Surface PDE solvers using Generalized Moving Least Squares (GMLS)

Here, we develop numerical methods for computing statistics of stochastic processes on surfaces of general shape with drift-diffusion dynamics d X t = a (X t ) dt + b(X t ) d W t . We formulate descriptions of Brownian motion and general drift-diffusion processes on surfaces. We consider statistics of the form u (x) = E x [$∫^{τ}_{0}$ g (X t ) dt ] + E x [ f (X τ )] for a domain Ω and the exit stopping time τ = inf t { t >0 | X i Ω}, where f , g are general smooth functions. For computing these statistics, we develop high-order Generalized Moving Least Squares (GMLS) solvers for associated surface PDE boundary-value problems based on Backward- Kolmogorov equations. We focus particularly on the mean First Passage Times (FPTs) given by the case f = 0, g = 1 where u (x) = E x [τ]. We perform studies for a variety of shapes showing our methods converge with high-order accuracy both in capturing the geometry and the surface PDE solutions. We then perform studies showing how statistics are influenced by the surface geometry, drift dynamics, and spatially dependent diffusivities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Reservoir Sediment Management and Monitoring Database

Overview This dataset compiles dam sediment management and monitoring information from surveys, case studies, and journal articles. Additionally, features described by the National Inventory of Dams (i.e., presence of sluice gates) are included to indicate known infrastructure features that may address sediment releases. The location and description of records from downstream monitoring gages are catalogued in order to help with tracking conditions over time (e.g., before and after management actions, as operations change, etc.). The data help address national scale understanding of challenges and solutions related to the accumulation of sediment behind a dam as well as downstream passage. Sediment trapping causes problems as it reduces storage capacity, disrupts dam and reservoir function, impedes access for recreation, alters water quality/habitat conditions, and contributes to riverbank and coastal erosion within the reservoir. Data compilation from a variety of sources is a first step towards assessing system-wide efficacy of management solutions. This dataset was developed under the Water Power Technologies Office funded effort which began as a Seedling on Reservoir Sedimentation Data, and was supported by the Reservoir Sedimentation Modeling Framework and Data Analysis project. These projects have addressed challenges in describing sediment transport, trapping, and management at dams throughout the US. Methodology An outer join on dams/reservoirs with surveys and survey reports (documented in the RESSED database, USBR or USACE databases, project websites, etc.) with the National Inventory of Dams, based on the NIDID to determine dams with documented management and/or sluice gates. Additional dams with documented management activity were identified through review of technical articles from the past 25 years in Journal of Hydrology, Journal of Water Resources Planning and Management, Geomorphology, Journal of Hydraulic Engineering, Water, Journal of Cleaner Production, International Journal of Sediment Research, Nature Scientific Reports, Earth Surface Processes and Landforms, and Environmental Science and Pollution Research. Individual records were created for each survey or management activity documented. To evaluate downstream sediment monitoring records, the nhdPlusTools and dataRetrieval packages in R were used to find gages within 10km of each dam in the management database. Length of record and location of matched gages were retrieved for those parameters relevant to sediment concentration or total sediment discharge.

Hansen, Carly [ORNL] (ORCID:0000000193280838)↗
DOE OSTI2021

Conceptualizing Biogeochemical Reactions With an Ohm's Law Analogy

Abstract In studying problems like plant‐soil‐microbe interactions in environmental biogeochemistry and ecology, one usually has to quantify and model how substrates control the growth of, and interaction among, biological organisms (and abiotic factors, e.g., adsorptive mineral soil surfaces). To address these substrate‐consumer relationships, many substrate kinetics and growth rules have been developed, including the famous Monod kinetics for single‐substrate‐based growth and Liebig's law of the minimum for multiple‐nutrient‐colimited growth. However, the mechanistic basis that leads to these various concepts and mathematical formulations and the implications of their parameters are often quite uncertain. Here, we show that an analogy based on Ohm's law in electric circuit theory is able to unify many of these different concepts and mathematical formulations. In this Ohm's law analogy, a resistor is defined by a combination of consumers’ and substrates’ kinetic traits. In particular, the resistance is equal to the mean first passage time that has been used to derive the Michaelis‐Menten kinetics under substrate replete conditions for a single substrate as well as the predation rate of individual organisms. We further show that this analogy leads to important insights on various biogeochemical problems, such as (a) multiple‐nutrient‐colimited biological growth, (b) denitrification, (c) fermentation under aerobic conditions, (d) metabolic temperature sensitivity, and (e) the legitimacy of Monod kinetics for describing bacterial growth. We expect that our approach will help both modelers and nonmodelers to better understand and formulate hypotheses when studying certain aspects of environmental biogeochemistry and ecology.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Shock tube simulations for the three-layer Richtmyer–Meshkov instability with single-mode and multimode perturbations

While the canonical two-component, single-mode Richtmyer–Meshkov instability (RMI) has been extensively studied, relatively less work has focused on the effects of an additional intermediate-density middle layer. This work investigates such three-material RMI configurations at two Atwood number scenarios using the ARES hydrodynamics code. After validation against previous experimental and computational studies, setups corresponding to recent three-layer shock tube experiments are simulated. Cases with both single-mode and multimode perturbations are studied to quantify mixing across the interface between the materials with highest and intermediate density. In particular, this work is able to comprehensibly examine differences between two- and three-dimensional setups for the single-mode and multimode problems. Observations from previous two-layer investigations still apply in the three-layer setup, but over the time horizons considered, there appears to be insufficient nonlinear mode coupling to create significant differences between two- and three-dimensional simulations following the first passage of a shock. Finally, additional reshock simulations have additional nonlinear growth that does result in expected differences between two- and three-dimensional cases in this three-layer setup, but significant differences do not manifest during the time horizon studied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Compressible pairwise interaction extended point-particle model for force prediction of shock-particle bed interaction

We propose a pairwise influence framework for the complex unsteady compressible particle-laden flow problem by accounting for the scattered hydrodynamic waves emitting from neighboring particles in a Euler-Lagrange simulation. It has been observed from particle-resolved (PR) simulations of randomly dispersed particle beds under a loading shock that the compressible pseudoturbulence dominates the flow system even after the primary shock has passed, which causes fluctuations observed in the forces experienced by the particle. Moreover, the fact that each particle exists in the vicinity of a random arrangement of other particles modifies the time history of the drag force experienced by each particle during and after the passage of the shock. First, the scattering flow field due to an incoming shock interacting with a single sphere is constructed using an analysis of the flow in the acoustic limit. Then we examine the validity of the compressible Maxey-Riley-Gatignol force model by comparing the force prediction against a PR simulation of two interacting particles for various particle arrangements and incoming shock strength. Subsequently, the neighboring influences are stored as a library of maps that can be used readily in the calculation of the perturbation force. Lastly, the pairwise interaction assumption is evaluated by comparing the force predicted with the model with PR simulations of a randomly packed particle bed of 10% volume fraction for both water and air as the fluid medium for an incoming shock Mach number 1.22. With a considerably lower cost for the implementation of the model compared to PR simulations, it is verified that the model is reasonably accurate in pinpointing particles whose peak force is significantly larger or smaller than the mean drag but also to capture the prolonged fluctuations after the initial shock.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Early-time measure in eternal inflation

In a situation like eternal inflation, where our data is replicated at infinitely-many other space-time events, it is necessary to make a prior assumption about our location to extract predictions. The principle of mediocrity entails that we live at asymptotic late times, when the occupational probabilities of vacua has settled to a near-equilibrium distribution. In this paper we further develop the idea that we instead exist during the approach to equilibrium, much earlier than the exponentially-long mixing time. In this case we are most likely to reside in vacua that are easily accessed dynamically. Using first-passage statistics, we prove that vacua that maximize their space-time volume at early times have: 1. maximal ever-hitting probability; 2. minimal mean first-passage time; and 3. minimal decay rate. These requirements are succinctly captured by an early-time measure. The idea that we live at early times is a predictive guiding principle, with many phenomenological implications. First, our vacuum should lie deep in a funneled region, akin to folding energy landscapes of proteins. Second, optimal landscape regions are characterized by relatively short-lived vacua, with lifetime of order the de Sitter Page time. For our vacuum, this lifetime is ~ 10 130 years, which is consistent with the Standard Model estimate due to Higgs metastability. Third, the measure favors vacua with small, positive vacuum energy. This can address the cosmological constant problem, provided there are sufficiently many vacua in the entire ensemble of funnels. As a concrete example, we study the Bousso-Polchinski lattice of flux vacua, and find that the early-time measure favors lattices with the fewest number of flux dimensions. This favors compactifications with a large hierarchy between the lightest modulus and all other Kähler and complex structure moduli.

79 ASTRONOMY AND ASTROPHYSICS↗