Search NASA⌕ Search

SEARCH · Search NASA

Results for “random processes”

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 559 records · Page 31

Asymptotic Equivalence of Probability Measures and Stochastic Processes

Let P{sub n} and Q{sub n} be two probability measures representing two different probabilistic models of some system (e.g., an n-particle equilibrium system, a set of random graphs with n vertices, or a stochastic process evolving over a time n) and let M{sub n} be a random variable representing a “macrostate” or “global observable” of that system. We provide sufficient conditions, based on the Radon–Nikodym derivative of P{sub n} and Q{sub n}, for the set of typical values of M{sub n} obtained relative to P{sub n} to be the same as the set of typical values obtained relative to Q{sub n} in the limit n→∞. This extends to general probability measures and stochastic processes the well-known thermodynamic-limit equivalence of the microcanonical and canonical ensembles, related mathematically to the asymptotic equivalence of conditional and exponentially-tilted measures. In this more general sense, two probability measures that are asymptotically equivalent predict the same typical or macroscopic properties of the system they are meant to model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Progress in Operational Analysis of Launch Vehicles in Nonstationary Flight

This paper presents recent results in an ongoing effort to understand and develop techniques to process launch vehicle data, which is extremely challenging for modal parameter identification. The primary source of difficulty is due to the nonstationary nature of the situation. The system is changing, the environment is not steady, and there is an active control system operating. Hence, the primary tool for producing clean operational results (significant data lengths and data averaging) is not available to the user. This work reported herein uses a correlation-based two step operational modal analysis approach to process the relevant data sets for understanding and development of processes. A significant drawback for such processing of short time histories is a series of beating phenomena due to the inability to average out random modal excitations. A recursive correlation process coupled to a new convergence metric (designed to mitigate the beating phenomena) is the object of this study. It has been found in limited studies that this process creates clean modal frequency estimates but numerically alters the damping.

James, George↗

On random fuzzy fractional partial integro-differential equations under Caputo generalized Hukuhara differentiability

This paper is devoted to the study of the solvability of random fuzzy fractional partial integro-differential equations under Caputo generalized Hukuhara differentiability. The notions of random fuzzy variables and fuzzy stochastic processes are developed for multivariable functions. The existence and uniqueness of two types of integral solutions generated from Darboux problem for nonlinear wave equations are proved using the successive approximations method and Gronwall’s inequality for stochastic processes. The continuous dependence on the data, the boundedness, and the stability with probability one of integral solutions are established to confirm the well posedness of our model. Some computational examples are presented to illustrate the theoretical results.

97 MATHEMATICS AND COMPUTING↗

Analysis of the Accuracy of Actuation Electronics in the Laser Interferometer Space Antenna Pathfinder

The Laser Interferometer Space Antenna Pathfinder (LPF) main observable, labelled∆g, is the differential force per unit mass acting on the two test masses in free fall condition, after the contribution of all non-gravitational forces have been compensated. At low frequencies the differential force is compensated by an applied electrostatic actuation force,which then must be subtracted from the measured acceleration to obtain∆g. Any inaccuracy in the actuation force contaminates the residual acceleration. This study investigates the accuracy of the electrostatic actuation system and its impact on the LPF main observable. It is shown that the inaccuracy is mainly caused by the rounding errors in the waveform processing and also by the random error caused by the analog to digital converter (ADC) random noise in the control loop. Both errors are one order of magnitude smaller than the resolution of the commanded voltages. We developed a simulator based on the LPF design to compute the close-to-reality actuation voltages and consequently, the resulting actuation forces. The simulator is applied in post-processing the LPF data.

M. Armano↗

Fast truncated SVD of sparse and dense matrices on graphics processors

We investigate the solution of low-rank matrix approximation problems using the truncated singular value decomposition (SVD). For this purpose, we develop and optimize graphics processing unit (GPU) implementations for the randomized SVD and a blocked variant of the Lanczos approach. Our work takes advantage of the fact that the two methods are composed of very similar linear algebra building blocks, which can be assembled using numerical kernels from existing high-performance linear algebra libraries. Furthermore, the experiments with several sparse matrices arising in representative real-world applications and synthetic dense test matrices reveal a performance advantage of the block Lanczos algorithm when targeting the same approximation accuracy.

Computer Science↗

Stochastic differential equation in a random environment

Solutions of the Itô stochastic differential equation in a random environment are considered. The random environment is formed by the generalized telegraph process. It is proved that the initial problem is equivalent to a system of two stochastic differential equations with nonrandom coefficients. The first equation is the Itô equation, and the initial process is its solution. The second equation is an equation with Poisson process, and its solution is a generalized telegraph process. The theorems of existence and uniqueness of strong and weak solutions are proved.

97 MATHEMATICS AND COMPUTING↗

Analysis of the Continuous Stellar Tracking Attitude Reference (CSTAR) attitude rate processor

The Continuous Stellar Tracking Attitude Reference (CSTAR) system is an in-house project for Space Station to provide high accuracy, drift free attitude and angular rate information for the GN&C system. The outputs of the solid state star trackers are processed to provide attitude information; rate data is then derived from the attitude. Rate derivation is based on discrete time polynomial approximation techniques. This gives simple algorithms which allow for interpolation by other users. Attitude rate is modeled as a constant with low amplitude, low frequency sinusoids superimposed. The rate processor is parameterized to account for the effects of random errors, sample rate, data processing rate and perturbation frequency. The baseline system may be characterized as follows: the three sigma attitude accuracy is 0.01 degrees, the three sigma rate accuracy is 0.0001 degrees per second, the sample rate is 100 Hertz, the sampled signal is bandlimited to 0.5 Hertz, and the data processing rate is 10 Hertz. The above system requires a differentiator of length 127. This will track rate perturbations of frequencies less than 0.01 Hertz with low systematic errors.

Uhde-Lacovara, J.↗

Matrix isolation technique for the study of some factors affecting the partitioning of trace elements

The factors that affect the preferred positions of cations in ionic solid solutions were investigated utilizing vibrational spectroscopy. Solid solutions of the sulfate and chromate ions codoped with La(+3) and Ca(+2) in a KBr host lattice were examined as a function of the polyvalent cation concentration. The cation-anion pairing process was found to be random for Ca(+2), whereas the formation of La(+3)-SO4(-2) ion pairs with a C2 sub v bonding geometry is highly preferential to any type of La(+3)-CrO4(-2) ion pair formation. The relative populations of ion pair site configurations are discussed in terms of an energy-entropy competition model which can be applied to the partition of trace elements during magmatic processes.

Grzybowski, J. M.↗

Meterwavelength Single-pulse Polarimetric Emission Survey. III. The Phenomenon of Nulling in Pulsars

A detailed analysis of nulling was conducted for the pulsars studied in the Meterwavelength Single-pulse Polarimetric Emission Survey. We characterized nulling in 36 pulsars including 17 pulsars where the phenomenon was reported for the first time. The most dominant nulls lasted for a short duration, less than five periods. Longer duration nulls extending to hundreds of periods were also seen in some cases. A careful analysis showed the presence of periodicities in the transition from the null to the burst states in 11 pulsars. In our earlier work, fluctuation spectrum analysis showed multiple periodicities in 6 of these 11 pulsars. We demonstrate that the longer periodicity in each case was associated with nulling. The shorter periodicities usually originate from subpulse drifting. The nulling periodicities were more aligned with the periodic amplitude modulation, indicating a possible common origin for both. The most prevalent nulls last for a single period and can be potentially explained using random variations affecting the plasma processes in the pulsar magnetosphere. On the other hand, longer-duration nulls require changes in the pair-production processes, which need an external triggering mechanism for the changes. The presence of periodic nulling puts an added constraint on the triggering mechanism, which also needs to be periodic.

79 ASTRONOMY AND ASTROPHYSICS↗

Variations in density and chemical composition at 120 km from chemical and dynamical processes.

Atmospheric parameters show both systematic and random patterns of behavior at 120 km altitude. Variations in density, temperature, pressure, and especially atomic oxygen concentration are important. The variations are particularly significant because constant boundary conditions at this altitude have commonly been assumed in atmospheric model making. It is surprising how well the assumption of constant boundary conditions has served over the past decade. However, their use has probably introduced erroneous concepts into the field of atmospheric structure, or at least held off the recognition and introduction of important new concepts. The variations at 120 km are caused by changes in energy input into the upper atmosphere, including internal gravity waves and tidal energy from below, and changes in the transport processes within the atmosphere.

Johnson, F. S.↗

Binary pseudo-random array standards for calibration of 3D optical surface profilers used for metrology with aspheric x-ray optics

High-accuracy surface metrology is vitally important in manufacturing ultra-high-quality free-form mirrors designed to manipulate x-ray light with nanometer-scale wavelengths. The current and potential capabilities of x˗ray mirror manufacturing are limited by inherent imperfections of the integrated metrology tools. Metrology tools are currently calibrated with super-polished flat test-standard/reference mirrors. This is acceptable for fabrication of slightly curved x-ray optics. However, for even moderately curved aspherical x-ray mirrors the flat-reference calibration is not sufficiently accurate. For micro-stitching interferometry developed for surface measurements with curved x-ray mirrors, the tool aberration errors are known to be transferred into the optical surface topography of x-ray mirrors. Our approach to improving metrology is to thoroughly calibrate the measuring tool and apply the results of the calibration to deconvolution of the measured data. Here we explore the application of a recently developed technique for calibrating the instrument transfer function (ITF) of 3D optical surface profilers to metrology with significantly curved x-ray optics. The technique, based on test standards patterned with two-dimensional (2D) binary pseudo-random arrays (BPRAs), employs the unique properties of the BPRA patterns in the spatial frequency domain. The inherent 2D power spectral density of the pattern has a deterministic white-noise-like character that allows direct determination of the ITF with uniform sensitivity over the entire spatial frequency range and field of view of an instrument. The high efficacy of the technique has been previously demonstrated in application to metrology with flat and slightly curved optics. Here, we concentrate on development of an efficient fabrication process for production of highly randomized (HR) BPRA test standards on flat and 500-mm spherical optical substrates. We also compare and discuss the results of the ITF calibration of an interferometric microscope when using the HR BPRA standards on flat and curved substrates.

Munechika, Keiko↗

Neural Algorithms for Low Power Implementation of Partial Differential Equations

The rise of low-power neuromorphic hardware has the potential to change high-performance computing; however much of the focus on brain-inspired hardware has been on machine learning applications. A low-power solution for solving partial differential equations could radically change how we approach large-scale computing in the future. The random walk is a fundamental stochastic process that underlies many numerical tasks in scientific computing applications. We consider here two neural algorithms that can be used to efficiently implement random walks on spiking neuromorphic hardware. The first method tracks the positions of individual walkers independently by using a modular code inspired by grid cells in the brain. The second method tracks the densities of random walkers at each spatial location directly. We present the scaling complexity of each of these methods and illustrate their ability to model random walkers under different probabilistic conditions. Finally, we present implementations of these algorithms on neuromorphic hardware.

97 MATHEMATICS AND COMPUTING↗

A high speed implementation of the random decrement algorithm

The algorithm is useful for measuring net system damping levels in stochastic processes and for the development of equivalent linearized system response models. The algorithm works by summing together all subrecords which occur after predefined threshold level is crossed. The random decrement signature is normally developed by scanning stored data and adding subrecords together. The high speed implementation of the random decrement algorithm exploits the digital character of sampled data and uses fixed record lengths of 2(n) samples to greatly speed up the process. The contributions to the random decrement signature of each data point was calculated only once and in the same sequence as the data were taken. A hardware implementation of the algorithm using random logic is diagrammed and the process is shown to be limited only by the record size and the threshold crossing frequency of the sampled data. With a hardware cycle time of 200 ns and 1024 point signature, a threshold crossing frequency of 5000 Hertz can be processed and a stably averaged signature presented in real time.

Kiraly, L. J.↗

A high-speed implementation of the random decrement algorithm

The algorithm is useful for measuring net system damping levels in stochastic processes and for the development of equivalent linearized system response models. The algorithm works by summing together all subrecords which occur after predefined threshold level is crossed. The random decrement signature is normally developed by scanning stored data and adding subrecords together. The high speed implementation of the random decrement algorithm exploits the digital character of sampled data and uses fixed record lengths of 2(n) samples to greatly speed up the process. The contributions to the random decrement signature of each data point was calculated only once and in the same sequence as the data were taken. A hardware implementation of the algorithm using random logic is diagrammed and the process is shown to be limited only by the record size and the threshold crossing frequency of the sampled data. With a hardware cycle time of 200 ns and 1024 point signature, a threshold crossing frequency of 5000 Hertz can be processed and a stably averaged signature presented in real time.

Kiraly, L. J.↗

The Nature of Grand Minima and Maxima from Fully Nonlinear Flux Transport Dynamos

We aim to investigate the nature and occurrence characteristics of grand solar minimum and maximum periods, which are observed in the solar proxy records such as {sup 10}Be and {sup 14}C, using a fully nonlinear Babcock–Leighton type flux transport dynamo including momentum and entropy equations. The differential rotation and meridional circulation are generated from the effect of turbulent Reynolds stress and are subjected to back-reaction from the magnetic field. To generate grand minimum- and maximum-like periods in our simulations, we used random fluctuations in the angular momentum transport process, namely the Λ-mechanism, and in the Babcock–Leighton mechanism. To characterize the nature and occurrences of the identified grand minima and maxima in our simulations, we used the waiting time distribution analyses, which reflect whether the underlying distribution arises from a random or a memory-bearing process. The results show that, in the majority of the cases, the distributions of grand minima and maxima reveal that the nature of these events originates from memoryless processes. We also found that in our simulations the meridional circulation speed tends to be smaller during grand maximum, while it is faster during grand minimum periods. The radial differential rotation tends to be larger during grand maxima, while it is smaller during grand minima. The latitudinal differential rotation, on the other hand, is found to be larger during grand minima.

79 ASTRONOMY AND ASTROPHYSICS↗

Snow Distribution Patterns Revisited: A Physics-Based and Machine Learning Hybrid Approach to Snow Distribution Mapping in the Sub-Arctic

Snowpack distribution in Arctic and alpine landscapes often occurs in repeating, year-to-year patterns due to local topographic, weather, and vegetation characteristics. Previous studies have suggested that with years of observational data, these snow distribution patterns can be statistically integrated into a snow process modeling workflow. Recent advances in snow hydrology and machine learning (ML) have increased our ability to predict snowpack distribution using in-situ observations, remote sensing data sets, and simple landscape characteristics that can be easily obtained for most environments. Here, we propose a hybrid approach to couple a ML snow distribution pattern (MLSDP) map with a physics-based, snow process model. We trained a random forest ML algorithm on tens of thousands of snow survey observations from a subarctic study area on the Seward Peninsula, Alaska, collected during peak snow water equivalent (SWE). We validated hybrid model outputs using in-situ snow depth and SWE observations, as well as a light detection and ranging data set and a distributed temperature profiling sensor data set. When the hybrid results were compared with the physics-based method, the hybrid method more accurately depicted the spatial patterns of the snowpack, areas of drifting snow, and years when no in-situ observations were used in the random forest ML training data set. The hybrid method also showed improvements in root mean squared error at 61% of locations where time-series estimations of snow depth were observed. These results can be applied to any physics-based model to improve the snow distribution patterning to reflect observed conditions in high latitude and high elevation cold region environments.

54 ENVIRONMENTAL SCIENCES↗

Dataset for scientific paper "Simulated plant‑mediated oxygen input has strong impacts on fine‑scale porewater biogeochemistry and weak impacts on integrated methane fluxes in coastal wetlands", a modeling study based on field observation at the tidal salt marshes of the Parker River Estuary, Massachusetts, United States

This dataset is the raw and processed data for the paper "Simulated plant ‑ mediated oxygen input has strong impacts on fine ‑ scale porewater biogeochemistry and weak impacts on integrated methane fluxes in coastal wetlands". This study investigated how plant-mediated oxygen input affects subsurface biogeochemical reactions of organic carbon degradation and the resulting methane emissions of coastal wetlands by model simulation. We used the subsurface geochemical simulator PFLOTRAN for the modeling, which produced the simulated changes in porewater chemical substances and methane emissions over 10 days under different scenarios of plant-mediated oxygen input.Specifically, this dataset contains: 1) the input files for PFLOTRAN of all simulation runs conducted in this study. Those files are with an extension of ".in", containing information of the biogeochemical reaction network (stoichiometry, reaction rate, Monod constants, etc), fluid flow rate and oxygen concentration in the fluid which together simulated the plant-mediated oxygen input, the configuration of artificial reactions that simulated the methane fluxes, etc. The PFLOTRAN input files are text files, which can be opened by NotePad, but running these input files will require proper installation of PFLOTRAN (instruction: https://documentation.pflotran.org/user_guide/how_to/installation/installation.html). 2) the raw and processed model output from PFLOTRAN of all simulation runs, and 3) the python scripts used to process the raw model output, including random allocation of root cells, converting raw data into organized formats, calculating the methane fluxes based on the model output, data visualization, etc. The raw and processed model output from PFLOTRAN are in .spydata format, which can be viewed with Python. and 3) the python scripts for data processing and analysis are programming scripts, which can be opened with Python.This modeling work, in particular the model parameterization of root density and initial conditions of porewater concentrations of biogeochemical substances, was based on field measurements at the salt marsh of the Upper Parker River Estuary, Massachusetts, United States.

54 ENVIRONMENTAL SCIENCES↗

Development of Collaborative Research Initiatives to Advance the Aerospace Sciences-via the Communications, Electronics, Information Systems Focus Group

The primary goal of the Adaptive Vision Laboratory Research project was to develop advanced computer vision systems for automatic target recognition. The approach used in this effort combined several machine learning paradigms including evolutionary learning algorithms, neural networks, and adaptive clustering techniques to develop the E-MOR.PH system. This system is capable of generating pattern recognition systems to solve a wide variety of complex recognition tasks. A series of simulation experiments were conducted using E-MORPH to solve problems in OCR, military target recognition, industrial inspection, and medical image analysis. The bulk of the funds provided through this grant were used to purchase computer hardware and software to support these computationally intensive simulations. The payoff from this effort is the reduced need for human involvement in the design and implementation of recognition systems. We have shown that the techniques used in E-MORPH are generic and readily transition to other problem domains. Specifically, E-MORPH is multi-phase evolutionary leaming system that evolves cooperative sets of features detectors and combines their response using an adaptive classifier to form a complete pattern recognition system. The system can operate on binary or grayscale images. In our most recent experiments, we used multi-resolution images that are formed by applying a Gabor wavelet transform to a set of grayscale input images. To begin the leaming process, candidate chips are extracted from the multi-resolution images to form a training set and a test set. A population of detector sets is randomly initialized to start the evolutionary process. Using a combination of evolutionary programming and genetic algorithms, the feature detectors are enhanced to solve a recognition problem. The design of E-MORPH and recognition results for a complex problem in medical image analysis are described at the end of this report. The specific task involves the identification of vertebrae in x-ray images of human spinal columns. This problem is extremely challenging because the individual vertebra exhibit variation in shape, scale, orientation, and contrast. E-MORPH generated several accurate recognition systems to solve this task. This dual use of this ATR technology clearly demonstrates the flexibility and power of our approach.

Knasel, T. Michael↗