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.

At least 19 records

First-passage problems: A probabilistic dynamic analysis for degraded structures

Structures subjected to random excitations with uncertain system parameters degraded by surrounding environments (a random time history) are studied. Methods are developed to determine the statistics of dynamic responses, such as the time-varying mean, the standard deviation, the autocorrelation functions, and the joint probability density function of any response and its derivative. Moreover, the first-passage problems with deterministic and stationary/evolutionary random barriers are evaluated. The time-varying (joint) mean crossing rate and the probability density function of the first-passage time for various random barriers are derived.

Shiao, Michael C.↗

A note on the first passage time problem

Kolmogoroff backward equation for analyzing first passage time problem for linear single degree of freedom vibratory system with linear viscous damping

Shinozuka, M.↗

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↗

Some heuristic procedures for analyzing random vibration of nonlinear oscillators.

The stationary response of a lightly damped nonlinear oscillator subjected to wideband random excitation can be examined as an example of thermal equilibrium. It may be assumed that the response consists of a series of free-vibration cycles with small random fluctuations in phase and amplitude. Certain statistical properties of the response can be estimated by averaging corresponding properties of the free vibration with respect to cycle amplitude distributions. Such heuristic procedures for determining the expected frequency and the autocorrelation function of the stationary response are outlined. Some additional results concerning first-passage problems for nonlinear oscillators are included.

Crandall, S. H.↗

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,↗

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↗

On the probability of cycle-slipping in first-order phase-locked loops.

The first-passage time boundary value problem for first-order phase-locked loops (PLL) is analyzed, and spectral representations are developed for the probability density function (pdf), the distribution function, and the moments of the first time to passage (or cycle-slip). For the sinusoidal PLL, an asymptotic formula, that is surprisingly accurate even at low loop SNR's and large frequency offsets, is obtained for the pdf of the time to cycle-slip, in terms of the mean time to slip.

La Frieda, J. R.↗

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↗

Airborne radar technology for windshear detection

The objectives and accomplishments of the two-and-a-half year effort to describe how returns from on-board Doppler radar are to be used to detect the presence of a wind shear are reported. The problem is modeled as one of first passage in terms of state variables, the state estimates are generated by a bank of extended Kalman filters working in parallel, and the decision strategy involves the use of a voting algorithm for a series of likelihood ratio tests. The performance issue for filtering is addressed in terms of error-covariance reduction and filter divergence, and the performance issue for detection is addressed in terms of using a probability measure transformation to derive theoretical expressions for the error probabilities of a false alarm and a miss.

Hibey, Joseph L.↗

Driven cavity simulation of turbomachine blade flows with vortex control

This paper presents a computational study of the three-dimensional flows in a rotating cavity with clearance between cavity walls and lid wall. The objectives of this study is to understand the interaction mechanism between tip leakage and blade passage flows and to assess the means to control the flow pattern and pressure losses. The classes of problems addressed include: passage geometry, passage loading including lid velocity and anti-vortex strength, and placement necessary to provide flow control. The computational model is first validated on generic flow problems and then applied to a specific blade passage configuration. Results of parametric studies for secondary flow pattern control are analyzed, and practical means of vortex control are discussed.

Athavale, M. M.↗

Limits on the isolation of stochastic vibration for microgravity space experiments

The limitations on the isolation of stochastic vibrations for microgravity space experiments are explored. These limitations result from the restricted interior space available for vibration isolation. A one-degree-of-freedom representation of the experiment-spacecraft system is used, and an ideal vibration actuator is assumed. A kinematic representation results, and the problem becomes one of finding the minimum acceleration trajectory within a pair of stochastic walls. The wall motion is characterized by an ergodic, stationary, zero-mean, Gaussian random process with known power spectral density. The geometry of the wall trajectories is defined in terms of their significant extrema and zero crossings. This geomemtry is used in defining a composite trajectory that has a mean square acceleratiaon lower than that on the optimal path satisfying the stochastic wall inequality constraints. The optimal control problem is solved on a return path yielding the mean square acceleration in terms of the distributions of significant maxima and first-passage time of the wall process. The methodology is applied to a microgravity isolation problem to find the lower bounds on root-mean-square acceleration given the disturbance power spectral density.

Knospe, C. R.↗

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↗

An implementation of a tree code on a SIMD, parallel computer

We describe a fast tree algorithm for gravitational N-body simulation on SIMD parallel computers. The tree construction uses fast, parallel sorts. The sorted lists are recursively divided along their x, y and z coordinates. This data structure is a completely balanced tree (i.e., each particle is paired with exactly one other particle) and maintains good spatial locality. An implementation of this tree-building algorithm on a 16k processor Maspar MP-1 performs well and constitutes only a small fraction (approximately 15%) of the entire cycle of finding the accelerations. Each node in the tree is treated as a monopole. The tree search and the summation of accelerations also perform well. During the tree search, node data that is needed from another processor is simply fetched. Roughly 55% of the tree search time is spent in communications between processors. We apply the code to two problems of astrophysical interest. The first is a simulation of the close passage of two gravitationally, interacting, disk galaxies using 65,636 particles. We also simulate the formation of structure in an expanding, model universe using 1,048,576 particles. Our code attains speeds comparable to one head of a Cray Y-MP, so single instruction, multiple data (SIMD) type computers can be used for these simulations. The cost/performance ratio for SIMD machines like the Maspar MP-1 make them an extremely attractive alternative to either vector processors or large multiple instruction, multiple data (MIMD) type parallel computers. With further optimizations (e.g., more careful load balancing), speeds in excess of today's vector processing computers should be possible.

Olson, Kevin M.↗

Unimolecular diffusion-mediated reactions with a nonrandom time-modulated absorbing barrier

A diffusion-reaction model with time-dependent reactivity is formulated and applied to unimolecular reactions. The model is solved exactly numerically and approximately analytically for the unreacted fraction as a function of time. It is shown that the approximate analytical solution is valid even when the system is far from equilibrium, and when the reactivity probability is more complicated than a square-wave function of time. A discussion is also given of an approach to problems of this type using a stochastically fluctuating reactivity, and the first-passage time for a particular example is derived.

Bashford, D.↗

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)↗

Dynamics of the Molten Contact Line

The purpose of this program is to develop a basic understanding of how a molten material front spreads over a solid that is below its melting point, arrests, and freezes. Our hope is that the work will contribute toward a scientific knowledge base for certain new applications involving molten droplet deposition, including the "printing" of arbitrary three-dimensional objects by precise deposition of individual molten microdrops that solidify after impact. Little information is available at this time on the capillarity-driven motion and arrest of molten contact line regions. Schiaffino and Sonin investigated the arrest of the contact line of a molten microcrystalline wax spreading over a subcooled solid "target" of the same material. They found that contact line arrest takes place at an apparent liquid contact angle that depends primarily on the Stefan number S=c(T(sub f) -T(sub t)/L based on the temperature difference between the fusion point and the target temperature, and proposed that contact line arrest occurs when the liquid's dynamic contact angle approaches the angle of attack of the solidification front just behind the contact line. They also showed, however, that the conventional continuum equations and boundary conditions have no meaningful solution for this angle. The solidification front angle is determined by the heat flux just behind the contact line, and the heat flux is singular at that point. By comparing experiments with numerical computations, Schiaffino and Sonin estimated that the conventional solidification model must break down within a distance of order 0.1 - 1 microns of the contact line. The physical mechanism for this breakdown is as yet undetermined, and no first-principles theory exists for the contact angle at arrest. Schiaffino and Sonin also presented a framework for understanding how to moderate Weber number molten droplet deposition in terms of similarity laws and experimentation. The study is based on experiments with three molten materials- molten wax on solid wax, water on ice, and mercury on frozen mercury- which between them span a considerable range of the deposition/solidification similarity parameters. Correlations are obtained for the spreading velocity, spreading time scales, the spreading factor (i.e. ratio of deposited drop's final footprint radius and the drop's initial radius), post-spreading liquid oscillation amplitudes and time scales, and bulk solidification time scales. Duthaler carried out an experimental and theoretical investigation of the relationship between the liquid's apparent contact angle and the Capillary number Ca=mu U/sigma based on contact line speed, for molten materials spreading over subcooled solids. This relationship is required for modeling of melt spreading. We have adapted Voinov's methodology to the molten contact line and formulated a theoretical model for the Ca vs. contact angle relationship, based Schiaffino and Sonin#s (1997a,b) wedge-like solidification front model. With the solidification front angle taken from Schiaffino and Sonin, the model is in good agreement with the experimental results for Ca vs. contact angle. Duthaler also extended the experimental investigation of droplet deposition and contact line freezing to more materials, including solder on glass, solder on solder, water on ice, and molten microcrystalline wax on wax. The latter also included tests on inclined targets. Deposition tests have also been done with molten octacosane (C28H58) on various targets. An important objective of our program has been the development of micron-scale sensors for measuring the transient temperature at a point on the substrate surface as a molten contact line moves over it. The expectation is that this temperature history will yield a better understanding of the thermal process in the contact line region. The sensors are of the thermistor type, either 2.5 microns or 1.5 microns square, microfabricated with silicon-based technology on either pure silicon or amorphous silicon dioxide chips. Each chip has 32 sensors distributed on its surface in arrays. The time response is better than 10 ms. At the time of writing, sensor calibration is in progress. Results on thermal transients during contact line passage will be discussed at the conference. While we expect that the data will provide information on the near-contact-line heat transfer process, we also foresee possible problems. First, the spatial resolution of the sensors may be insufficient to resolve the near-contact-line region. Second, the sensors protrude about 0.5 microns above the substrate surface, and may affect the contact line motion. Third, a sensor's temperature history depends on both the heat flux distribution into it from the fusion front and the thermal properties of the substrate below it and the solidified melt between it and the fusion front. The heat flux distribution in the contact line region must therefore be unfolded from computations of the overall system's transient thermal response.

Sonin, Ain A.↗

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↗