SEARCH · Engineering Papers
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.
In this paper we construct the family of Geometric quadratic stochastic operators defined on the countable sample space of nonnegative integers and investigate their trajectory behavior. Such operators can be reinterpreted in terms of of evolutionary operator of free population. We show that Geometric quadratic stochastic operators are regular transformations.
The parameter-induced stochastic resonance based on spectral entropy (PSRSE) method is introduced for the detection of a very weak signal in the presence of strong noise. The effect of stochastic resonance on the detection is optimized using parameters obtained in spectral entropy analysis. Upon processing employing the PSRSE method, the amplitude of the weak signal is enhanced and the noise power is reduced, so that the frequency of the signal can be estimated with greater precision through spectral analysis. While the improvement in the signal-to-noise ratio is similar to that obtained using the Duffing oscillator algorithm, the computational cost reduces from O(N{sup 2}) to O(N). The PSRSE approach is applied to the frequency measurement of a weak signal made by a vortex flow meter. The results are compared with those obtained applying the Duffing oscillator algorithm.
A control problem for stochastically forced nonlinear continuous-time systems is considered. We propose a method for construction of the regulator that provides a preassigned probabilistic distribution of random states in stochastic equilibrium. Geometric criteria of the controllability are obtained. Constructive technique for the specification of attainability sets is suggested.
In past work, we introduced model-based inverse methods, and applied them to problems in which the anomaly could be reasonably modeled by simple canonical shapes, such as rectangular solids. In these cases the parameters to be inverted would be length, width and height, as well as the occasional probe lift-off or rotation. We are now developing a formulation that allows more flexibility in modeling complex flaws. The idea consists of expanding the flaw in a sequence of basis functions, and then solving for the expansion coefficients of this sequence, which are modeled as independent random variables, uniformly distributed over their range of values. There are a number of applications of such modeling: 1. Connected cracks and multiple half-moons, which we have noted in a POD set. Ideally we would like to distinguish connected cracks from one long shallow crack. 2. Cracks of irregular profile and shape which have appeared in cold work holes during bolt-hole eddy-current inspection. One side of such cracks is much deeper than other. 3. L or C shaped crack profiles at the surface, examples of which have been seen in bolt-hole cracks. By formulating problems in a stochastic sense, we are able to leverage the stochastic global optimization algorithms in NLSE, which is resident in VIC-3D®, to answer questions of global minimization and to compute confidence bounds using the sensitivity coefficient that we get from NLSE. We will also address the issue of surrogate functions which are used during the inversion process, and how they contribute to the quality of the estimation of the bounds.
A quadratic stochastic operator (in short QSO) is usually used to present the time evolution of differing species in biology. Some quadratic stochastic operators have been studied by Lotka and Volterra. In the present paper, we first give a simple characterization of Volterra QSO in terms of absolutely continuity of discrete measures. Further, we introduce a notion of orthogonal preserving QSO, and describe such kind of operators defined on two dimensional simplex. It turns out that orthogonal preserving QSOs are permutations of Volterra QSO. The associativity of genetic algebras generated by orthogonal preserving QSO is studied too.
We develop a fully stochastic theory for coagulation and fragmentation (CF) in a finite system with a maximum cluster size constraint. The process is modeled using a high-dimensional master equation for the probabilities of cluster configurations. For certain realizations of total mass and maximum cluster sizes, we find exact analytical results for the expected equilibrium cluster distributions. If coagulation is fast relative to fragmentation and if the total system mass is indivisible by the mass of the largest allowed cluster, we find a mean cluster-size distribution that is strikingly broader than that predicted by the corresponding mass-action equations. Combinations of total mass and maximum cluster size under which equilibration is accelerated, eluding late-stage coarsening, are also delineated. Finally, we compute the mean time it takes particles to first assemble into a maximum-sized cluster. Through careful state-space enumeration, the scaling of mean assembly times is derived for all combinations of total mass and maximum cluster size. We find that CF accelerates assembly relative to monomer kinetic only in special cases. All of our results hold in the infinite system limit and can be only derived from a high-dimensional discrete stochastic model, highlighting how classical mass-action models of self-assembly can fail.
We investigate a class of Hamiltonian particle systems and their stochastic behaviors. Using both rigorous proof and numerical simulations, we show that the geometric configuration can qualitatively change key statistical characteristics of the particle system, which are expected to be retained by stochastic modifications. In particular, whether a particle system has an exponential mixing rate or a polynomial mixing rate depends on whether the geometric setting allows a slow particle being reached by adjacent fast particles.
This paper is concerned with the stochastic coupled fractional Ginzburg-Landau equation with additive noise. We first transform the stochastic coupled fractional Ginzburg-Landau equation into random equations whose solutions generate a random dynamical system. Then we prove the existence of random attractor for random dynamical system.
We investigate, both numerically and analytically, the diffusion properties of a stochastic sled sliding on a substrate, subject to a constant towing force. The problem is motivated by the growing interest in controlling transport of artificial microswimmers in 2D geometries at low Reynolds numbers. We simulated both symmetric and asymmetric towed sleds. Remarkable properties of their mobilities and diffusion constants include sidewise drifts and excess diffusion peaks. We interpret our numerical findings by making use of stochastic approximation techniques.
Hall cross structures in magnetic nanowires are commonly used for electrical detection of magnetization reversal in which a domain wall (DW) is conventionally nucleated by a local Oersted field. In this letter, we demonstrate DW nucleation in Co/Ni perpendicular magnetic anisotropy nanowire at the magnetic Hall cross junction. The DWs are nucleated by applying an in-plane pulsed current through the nanowire without the need of a local Oersted field. The change in Hall resistance, detected using anomalous Hall effect, is governed by the magnetic volume switched at the Hall junction, which can be tuned by varying the magnitude of the applied current density and pulse width. The nucleated DWs are driven simultaneously under the spin transfer torque effect when the applied current density is above a threshold. The possibility of multiple DW generation and variation in magnetic volume switched makes nucleation process stochastic in nature. The in-plane current induced stochastic nature of DW generation may find applications in random number generation.
The electron residual energy originated from the stochastic heating in under-dense field-ionized plasma is investigated here. Initially, the optical response of plasma is modeled by using two counter-propagating electromagnetic waves. In this case, the solution of motion equation of a single electron indicates that by including the ionization, the electron with higher residual energy compared with that without ionization could be obtained. In agreement with chaotic nature of the motion, it is found that the electron residual energy will be significantly changed by applying a minor change in the initial conditions. Extensive kinetic 1D-3V particle-in-cell simulations have been performed in order to resolve full plasma reactions. In this way, two different regimes of plasma behavior are observed by varying the pulse length. The results indicate that the amplitude of scattered fields in a proper long pulse length is high enough to act as a second counter-propagating wave and trigger the stochastic electron motion. On the contrary, the analyses of intensity spectrum reveal the fact that the dominant scattering mechanism tends to Thomson rather than Raman scattering by increasing the pulse length. A covariant formalism is used to describe the plasma heating so that it enables us to measure electron temperature inside and outside of the pulse region.
A dynamic model of the thermochemical process in the flow reactor is considered. We study an influence of the random disturbances on the stationary regime of this model. A phenomenon of noise-induced excitability is demonstrated. For the analysis of this phenomenon, a constructive technique based on the stochastic sensitivity functions and confidence domains is applied. It is shown how elaborated technique can be used for the probabilistic analysis of the generation of mixed-mode stochastic oscillations in the flow chemical reactor.
In this paper, we present a novel hybrid method to simulate discrete stochastic reaction-diffusion models arising in biochemical signaling pathways. We study moderately stiff systems, for which we can partition each reaction or diffusion channel into either a slow or fast subset, based on its propensity. Numerical approaches missing this distinction are often limited with respect to computational run time or approximation quality. We design an approximate scheme that remedies these pitfalls by using a new blending strategy of the well-established inhomogeneous stochastic simulation algorithm and the tau-leaping simulation method. The advantages of our hybrid simulation algorithm are demonstrated on three benchmarking systems, with special focus on approximation accuracy and efficiency.
Experimental and theoretical studies are presented on the design of perturbations that enhance desynchronization in populations of oscillators that are synchronized by periodic entrainment. A phase reduction approach is used to determine optimal perturbation timing based upon experimentally measured phase response curves. The effectiveness of the perturbation waveforms is tested experimentally in populations of periodically and stochastically synchronized chemical oscillators. The relevance of the approach to therapeutic methods for disrupting phase coherence in groups of stochastically synchronized neuronal oscillators is discussed.
Predicting phonon relaxation time from molecular dynamics (MD) requires a long simulation time to compute the mode energy auto-correlation function. Here, we present an alternative approach to infer the phonon life-time from an approximate form of the energy auto-correlation function. The method requires as an input a set of sampled equilibrium configurations. A stochastic sampling method is used to generate the equilibrium configurations. We consider a truncated Taylor series expansion of the phonon energy auto-correlation function. The different terms in the truncated correlation function are obtained using the stochastic sampling approach. The expansion terms, thus, obtained are in good agreement with the corresponding values obtained using MD. We then use the approximate function to compute the phonon relaxation time. The relaxation time computed using this method is compared with that obtained from the exact correlation function. The two values are in agreement with each other.
Stochastic modelling is critical for studying many biochemical processes in a cell, in particular when some reacting species have low population numbers. For many such cellular processes the spatial distribution of the molecular species plays a key role. The evolution of spatially heterogeneous biochemical systems with some species in low amounts is accurately described by the mesoscopic model of the Reaction-Diffusion Master Equation. The Inhomogeneous Stochastic Simulation Algorithm provides an exact strategy to numerically solve this model, but it is computationally very expensive on realistic applications. We propose a novel adaptive time-stepping scheme for the tau-leaping method for approximating the solution of the Reaction-Diffusion Master Equation. This technique combines effective strategies for variable time-stepping with path preservation to reduce the computational cost, while maintaining the desired accuracy. The numerical tests on various examples arising in applications show the improved efficiency achieved by the new adaptive method.
Self-assembly of proteins into amyloid aggregates is an important biological phenomenon associated with human diseases such as Alzheimer’s disease. Amyloid fibrils also have potential applications in nano-engineering of biomaterials. The kinetics of amyloid assembly show an exponential growth phase preceded by a lag phase, variable in duration as seen in bulk experiments and experiments that mimic the small volumes of cells. Here, to investigate the origins and the properties of the observed variability in the lag phase of amyloid assembly currently not accounted for by deterministic nucleation dependent mechanisms, we formulate a new stochastic minimal model that is capable of describing the characteristics of amyloid growth curves despite its simplicity. We then solve the stochastic differential equations of our model and give mathematical proof of a central limit theorem for the sample growth trajectories of the nucleated aggregation process. These results give an asymptotic description for our simple model, from which closed form analytical results capable of describing and predicting the variability of nucleated amyloid assembly were derived. We also demonstrate the application of our results to inform experiments in a conceptually friendly and clear fashion. Our model offers a new perspective and paves the way for a new and efficient approach on extracting vital information regarding the key initial events of amyloid formation.
This paper investigates the population growth of a certain species in which every generation reproduces thrice over a period of predefined time, under certain constraints of resources needed for survival of population. We study the survival period of a species by randomizing the reproduction probabilities within a window at same predefined ages and the resources are being produced by the working force of the population at a variable rate. This randomness in the reproduction rate makes the population growth stochastic in nature and one cannot predict the exact form of evolution. Hence we study the growth by running simulations for such a population and taking an ensemble averaged over 500 to 5000 such simulations as per the need. While the population reproduces in a stochastic manner, we have implemented a constraint on the amount of resources available for the population. This is important to make the simulations more realistic. The rate of resource production then is tuned to find the rate which suits the survival of the species. We also compute the mean life time of the species corresponding to different resource production rate. Study for these outcomes in the parameter space defined by the reproduction probabilities and rate of resource production is carried out.