Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “random walk”

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

Vortex Loops at the Superfluid Lambda Transition: An Exact Theory?

A vortex-loop theory of the superfluid lambda transition has been developed over the last decade, with many results in agreement with experiments. It is a very simple theory, consisting of just three basic equations. When it was first proposed the main uncertainty in the theory was the use Flory scaling to find the fractal dimension of the random-walking vortex loops. Recent developments in high-resolution Monte Carlo simulations have now made it possible to verify the accuracy of this Flory-scaling assumption. Although the loop theory is not yet rigorously proven to be exact, the Monte Carlo results show at the least that it is an extremely good approximation. Recent loop calculations of the critical Casimir effect in helium films in the superfluid phase T < Tc will be compared with similar perturbative RG calculations in the normal phase T > Tc; the two calculations are found to match very nicely right at Tc.

Williams, Gary A.↗

A Derivation of the Long-Term Degradation of a Pulsed Atomic Frequency Standard from a Control-Loop Model

The phase of a frequency standard that uses periodic interrogation and control of a local oscillator (LO) is degraded by a long-term random-walk component induced by downconversion of LO noise into the loop passband. The Dick formula for the noise level of this degradation is derived from an explicit solution of an LO control-loop model.

Dick formula local oscillator frequency phase freq↗

A Derivation of the Dick Effect from Control-Loop Models for Periodically Interrogated Passive Frequency Standards

The phase of a frequency standard that uses periodic interrogation and control of a local oscillator (LO) is degraded by a long-term random-walk component induced by downconversion of LO noise into the loop passband. The Dick formula for the noise level of this degradation can be derived from explicit solotions of two LO control-loop models. A summary of the derivations is given here.

Dick Effect Dick Formula local oscillator noise fr↗

Exploring the Connection Between Sampling Problems in Bayesian Inference and Statistical Mechanics

The Bayesian and statistical mechanical communities often share the same objective in their work - estimating and integrating probability distribution functions (pdfs) describing stochastic systems, models or processes. Frequently, these pdfs are complex functions of random variables exhibiting multiple, well separated local minima. Conventional strategies for sampling such pdfs are inefficient, sometimes leading to an apparent non-ergodic behavior. Several recently developed techniques for handling this problem have been successfully applied in statistical mechanics. In the multicanonical and Wang-Landau Monte Carlo (MC) methods, the correct pdfs are recovered from uniform sampling of the parameter space by iteratively establishing proper weighting factors connecting these distributions. Trivial generalizations allow for sampling from any chosen pdf. The closely related transition matrix method relies on estimating transition probabilities between different states. All these methods proved to generate estimates of pdfs with high statistical accuracy. In another MC technique, parallel tempering, several random walks, each corresponding to a different value of a parameter (e.g. "temperature"), are generated and occasionally exchanged using the Metropolis criterion. This method can be considered as a statistically correct version of simulated annealing. An alternative approach is to represent the set of independent variables as a Hamiltonian system. Considerab!e progress has been made in understanding how to ensure that the system obeys the equipartition theorem or, equivalently, that coupling between the variables is correctly described. Then a host of techniques developed for dynamical systems can be used. Among them, probably the most powerful is the Adaptive Biasing Force method, in which thermodynamic integration and biased sampling are combined to yield very efficient estimates of pdfs. The third class of methods deals with transitions between states described by rate constants. These problems are isomorphic with chemical kinetics problems. Recently, several efficient techniques for this purpose have been developed based on the approach originally proposed by Gillespie. Although the utility of the techniques mentioned above for Bayesian problems has not been determined, further research along these lines is warranted

Pohorille, Andrew↗

A Stable Clock Error Model Using Coupled First and Second Order Gauss-Markov Processes

Long data outages may occur in applications of global navigation satellite system technology to orbit determination for missions that spend significant fractions of their orbits above the navigation satellite constellation(s). Current clock error models based on the random walk idealization may not be suitable in these circumstances, since the covariance of the clock errors may become large enough to overflow flight computer arithmetic. A model that is stable, but which approximates the existing models over short time horizons is desirable. A coupled first- and second-order Gauss-Markov process is such a model.

Carpenter, Russell↗

Using LDPC Code Constraints to Aid Recovery of Symbol Timing

A method of utilizing information available in the constraints imposed by a low-density parity-check (LDPC) code has been proposed as a means of aiding the recovery of symbol timing in the reception of a binary-phase-shift-keying (BPSK) signal representing such a code in the presence of noise, timing error, and/or Doppler shift between the transmitter and the receiver. This method and the receiver architecture in which it would be implemented belong to a class of timing-recovery methods and corresponding receiver architectures characterized as pilotless in that they do not require transmission and reception of pilot signals. Acquisition and tracking of a signal of the type described above have traditionally been performed upstream of, and independently of, decoding and have typically involved utilization of a phase-locked loop (PLL). However, the LDPC decoding process, which is iterative, provides information that can be fed back to the timing-recovery receiver circuits to improve performance significantly over that attainable in the absence of such feedback. Prior methods of coupling LDPC decoding with timing recovery had focused on the use of output code words produced as the iterations progress. In contrast, in the present method, one exploits the information available from the metrics computed for the constraint nodes of an LDPC code during the decoding process. In addition, the method involves the use of a waveform model that captures, better than do the waveform models of the prior methods, distortions introduced by receiver timing errors and transmitter/ receiver motions. An LDPC code is commonly represented by use of a bipartite graph containing two sets of nodes. In the graph corresponding to an (n,k) code, the n variable nodes correspond to the code word symbols and the n-k constraint nodes represent the constraints that the code places on the variable nodes in order for them to form a valid code word. The decoding procedure involves iterative computation of values associated with these nodes. A constraint node represents a parity-check equation using a set of variable nodes as inputs. A valid decoded code word is obtained if all parity-check equations are satisfied. After each iteration, the metrics associated with each constraint node can be evaluated to determine the status of the associated parity check. Heretofore, normally, these metrics would be utilized only within the LDPC decoding process to assess whether or not variable nodes had converged to a codeword. In the present method, it is recognized that these metrics can be used to determine accuracy of the timing estimates used in acquiring the sampled data that constitute the input to the LDPC decoder. In fact, the number of constraints that are satisfied exhibits a peak near the optimal timing estimate. Coarse timing estimation (or first-stage estimation as described below) is found via a parametric search for this peak. The present method calls for a two-stage receiver architecture illustrated in the figure. The first stage would correct large time delays and frequency offsets; the second stage would track random walks and correct residual time and frequency offsets. In the first stage, constraint-node feedback from the LDPC decoder would be employed in a search algorithm in which the searches would be performed in successively narrower windows to find the correct time delay and/or frequency offset. The second stage would include a conventional first-order PLL with a decision-aided timing-error detector that would utilize, as its decision aid, decoded symbols from the LDPC decoder. The method has been tested by means of computational simulations in cases involving various timing and frequency errors. The results of the simulations ined in the ideal case of perfect timing in the receiver.

Jones, Christopher↗

A Monte-Carlo Model for the Formation of Radiation-induced Chromosomal Aberrations

Purpose: To simulate radiation-induced chromosome aberrations in mammalian cells (e.g., rings, translocations, and dicentrics) and to calculate their frequency distributions following exposure to DNA double strand breaks (DSBs) produced by high-LET ions. Methods: The interphase genome was assumed to be comprised of a collection of 2 kbp rigid-block monomers following the random-walk geometry. Additional details for the modeling of chromosomal structure, such as chromosomal domains and chromosomal loops, were included. A radial energy profile for heavy ion tracks was used to simulate the high-LET pattern of induced DSBs. The induced DSB pattern depended on the ion charge and kinetic energy, but always corresponded to the DSB yield of 25 DSBs/cell/Gy. The sum of all energy contributions from Poisson-distributed particle tracks was taken to account for all possible one-track and multi-track effects. The relevant output of the model was DNA fragments produced by DSBs. The DSBs, or breakpoints, were defined by (x, y, z, l) positions, where x, y, z were the Euclidian coordinates of a DSB, and where l was the relative position along the genome. Results: The code was used to carry out Monte Carlo simulations for DSB rejoinings at low doses. The resulting fragments were analyzed to estimate the frequencies of specific types of chromosomal aberrations. Histograms for relative frequencies of chromosomal aberrations and P.D.F.s (probability density functions) of a given aberration type were produced. The relative frequency of dicentrics to rings was compared to empirical data to calibrate rejoining probabilities. Of particular interest was the predicted distribution of ring sizes, irrespective of their frequencies relative to other aberrations. Simulated ring sizes were . 4 kbp, which are far too small to be observed experimentally (i.e., by microscopy) but which, nevertheless, are conjectured to exist. Other aberrations, for example, inversions, translocations, as well as multi-centrics were also recorded. Conclusion: High-LET DNA damage affects the frequencies of chromosomal aberrations. The ratio of rings to dicentrics is correct for the genomic size cut-offs corresponding to available experimental data. The present work predicts a relative abundance of small rings following irradiation by heavy ions.

Ponomarev, Artem L.↗

Modeling of the Ablation of Fibrous Materials in the Knudsen Regime

During atmospheric entry of planetary probes, the thermal protection system (TPS) of the probe is exposed to high temperatures under low pressures. In these conditions, carbonous TPS materials undergo gasification in the Knudsen regime leading to mass loss and wall recession called ablation. This work aims to improve the understanding of materiaVenvironment interactions through a study of the coupling between carbon dioxide transport in the Knudsen regime, heterogeneous oxidation of carbon, and sutface recession. A 3D Monte-Carlo simulation tool is used for this study. The fibrous architecture of the materiils, consisting of high porosity random array of carbon fibers, is numerically reproduced on a 3D Cartesian grid. Mass transport in the Knudsen regime from the boundary layer to the surface, and inside this porous material is simulated by random walk. A reaction probability is used to simulate the heterogeneous oxidation reaction. The surface recession is followed by front tracking using a simplified marching cube approach. The output data of the simulations are ablation velocity and dynamic evolution of the material porosity. A parametric study is carried out to analyze the material behavior as a function of Knudsen number for the porous media (length of the mean free path compared to the mean pore diameter) and the intrinsic reactivity of the carbon fibers. The results enable extrapolation of laboratory experimental data to actual entry conditions.

Lachaud, J.↗

TEACH (Train to Enable/Achieve Culturally Sensitive Healthcare)

Personnel from diverse ethnic and demographic backgrounds come together in both civilian and military healthcare systems, facing diagnoses that at one level are equalizers: coronary disease is coronary disease, breast cancer is breast cancer. Yet the expression of disease in individuals from different backgrounds, individual patient experience of disease as a particular illness, and interactions between patients and providers occurring in any given disease scenario, all vary enormously depending on the fortuity of the equation of "which patient happens to arrive in whose exam room." Previously, providers' absorption of lessons-learned depended on learning as an apprentice would when exposed over time to multiple populations. As a result, and because providers are often thrown into situations where communications falter through inadequate direct patient experience, diversity in medicine remains a training challenge. The questions then become: Can simulation and virtual training environments (VTEs) be deployed to short-track and standardize this sort of random-walk problem? Can we overcome the unevenness of training caused by some providers obtaining the valuable exposure to diverse populations, whereas others are left to "sink or swim"? This paper summarizes developing a computer-based VTE called TEACH (Training to Enable/Achieve Culturally Sensitive Healthcare). TEACH was developed to enhance healthcare providers' skills in delivering culturally sensitive care to African-American women with breast cancer. With an authoring system under development to ensure extensibility, TEACH allows users to role-play in clinical oncology settings with virtual characters who interact on the basis of different combinations of African American sub-cultural beliefs regarding breast cancer. The paper reports on the roll-out and evaluation of the degree to which these interactions allow providers to acquire, practice, and refine culturally appropriate communication skills and to achieve cultural and individual personalization of healthcare in their clinical practices.

Maulitz, Russell↗

Calculation of Dose Deposition in Nanovolumes and Simulation of gamma-H2AX Experiments

Monte-Carlo track structure simulations can accurately simulate experimental data: a) Frequency of target hits. b) Dose per event. c) Dose per ion. d) Radial dose. The dose is uniform in micrometers sized voxels; at the nanometer scale, the difference in energy deposition between high and low-LET radiations appears. The calculated 3D distribution of dose voxels, combined with chromosomes simulated by random walk is very similar to the distribution of DSB observed with gamma-H2AX experiments. This is further evidenced by applying a visualization threshold on dose.

Plante, Ianik↗

4U 1909+07: A Hidden Pearl

We present a detailed spectral and timing analysis of the High Mass X-ray Binary (HMXB) 4U 1909+07 with INTEGRAL and RXTE. 4U 1909+07 is a persistent accreting X-ray pulsar with a period of approximately 605 s. The period changes erratically consistent with a random walk expected for a wind accreting system. INTEGRAL detects the source with an average of 2.4 cps (corresponding to 15 mCrab), but sometimes exhibits flaring activity up to 50 cps (i.e. 300 mCrab). The strongly energy dependent pulse profile shows a double peaked structure at low energies and only a single narrow peak at energies above 20 keV. The phase averaged spectrum is well described by a powerlaw modified at higher energies by an exponential cutoff and photoelectric absorption at low energies. In addition at 6.4 keV a strong iron fluorescence line and at lower energies a black body component are present. We performed phase resolved spectroscopy to study the pulse phase dependence of the spectral parameters: while most spectral parameters are constant within uncertainties, the blackbody normalization and the cutoff folding energy vary strongly with phase.

Kreykenbohm, Ingo↗

Continuous Tuning and Calibration of Vibratory Gyroscopes

A method of control and operation of an inertial reference unit (IRU) based on vibratory gyroscopes provides for continuously repeated cycles of tuning and calibration. The method is intended especially for application to an IRU containing vibratory gyroscopes that are integral parts of microelectromechanical systems (MEMS) and that have cloverleaf designs, as described in several previous NASA Tech Briefs articles. The method provides for minimization of several measures of spurious gyroscope output, including zero-rate offset (ZRO), angle random walk (ARW), and rate drift. These benefits are afforded both at startup and thereafter during continuing operation, in the presence of unknown rotation rates and changes in temperature. A vibratory gyroscope contains a precision mechanically resonant structure containing two normal modes of vibration nominally degenerate in frequency and strongly coupled via a Coriolis term. In the case of the cloverleaf design MEMS gyro, these normal modes of vibration are plate rocking modes. The rocking motion of the plate is described by giving two angles, theta(sub 1) and theta(sub 2). A proof mass consisting of a post orthogonal to the plate ensures a high degree of Coriolis coupling of vibratory energy from one mode into the other under inertial rotation. The plate is driven and sensed capacitively across a few-microns-wide gap, and the normal mode frequencies can be tuned electrostatically by DC voltages applied across this gap. In order to sense rotation, the resonator plate is caused to rock in the theta(sub 1) direction, then any small motions in the theta(sub 2) direction are sensed, rebalanced, and interpreted as inertial rotation. In this scenario, the "drive" has been assigned to the theta(sub 1) direction, and the "sense" has been assigned to the theta(sub 2) direction.

Hayworth, Ken↗

Shape Shifting Satellites in Binary Near-Earth Asteroids: Do Meteoroid Impacts Play a Role in BYORP Orbital Evolution?

Less than catastrophic meteoroid impacts over 10(exp 5) years may change the shape of small rubble-pile satellites in binary NEAs, lengthening the average BYORP (binary Yarkovsky-Radzievskii-Paddack) rate of orbital evolution. An estimate of shape-shifting meteoroid fluxes give numbers close enough to causing random walks in the semimajor axis of binary systems to warrant further investigation

Rubincam, David Parry↗

Staring at 4U 1909+07 with Suzaku (Research Note)

We present an analysis of the neutron star High Mass X-ray Binary (HMXB) 4U 1909+07 mainly based on Suzaku data. We extend the pulse period evolution, which behaves in a random-walk like manner, indicative of direct wind accretion. Studying the spectral properties of 4U 1909+07 between 0.5 to 90keV we find that a power-law with an exponential cutoff can describe the data well, when additionally allowing for a blackbody or a partially covering absorber at low energies. We find no evidence for a cyclotron resonant scattering feature (CRSF), a feature seen in many other neutron star HMXBs sources. By performing pulse phase resolved spectroscopy we investigate the origin of the strong energy dependence of the pulse profile, which evolves from a broad two-peak profile at low energies to a profile with a single, narrow peak at energies above 20keV. Our data show that it is very likely that a higher folding energy in the high energy peak is responsible for this behavior. This in turn leads to the assumption that we observe the two magnetic poles and their respective accretion columns at different phases, and that these accretions column have slightly different physical conditions.

Fuerst, F.↗

Total Dose Effects on Bipolar Integrated Circuits at Low Temperature

Total dose damage in bipolar integrated circuits is investigated at low temperature, along with the temperature dependence of the electrical parameters of internal transistors. Bandgap narrowing causes the gain of npn transistors to decrease far more at low temperature compared to pnp transistors, due to the large difference in emitter doping concentration. When irradiations are done at temperatures of -140 deg C, no damage occurs until devices are warmed to temperatures above -50 deg C. After warm-up, subsequent cooling shows that damage is then present at low temperature. This can be explained by the very strong temperature dependence of dispersive transport in the continuous-time-random-walk model for hole transport. For linear integrated circuits, low temperature operation is affected by the strong temperature dependence of npn transistors along with the higher sensitivity of lateral and substrate pnp transistors to radiation damage.

linear integrated circuits↗

Using Markov Models of Fault Growth Physics and Environmental Stresses to Optimize Control Actions

A generalized Markov chain representation of fault dynamics is presented for the case that available modeling of fault growth physics and future environmental stresses can be represented by two independent stochastic process models. A contrived but representatively challenging example will be presented and analyzed, in which uncertainty in the modeling of fault growth physics is represented by a uniformly distributed dice throwing process, and a discrete random walk is used to represent uncertain modeling of future exogenous loading demands to be placed on the system. A finite horizon dynamic programming algorithm is used to solve for an optimal control policy over a finite time window for the case that stochastic models representing physics of failure and future environmental stresses are known, and the states of both stochastic processes are observable by implemented control routines. The fundamental limitations of optimization performed in the presence of uncertain modeling information are examined by comparing the outcomes obtained from simulations of an optimizing control policy with the outcomes that would be achievable if all modeling uncertainties were removed from the system.

Bole, Brian↗

Simulations of DSB Yields and Radiation-induced Chromosomal Aberrations in Human Cells Based on the Stochastic Track Structure Induced by HZE Particles

The formation of double‐strand breaks (DSBs) and chromosomal aberrations (CAs) is of great importance in radiation research and, specifically, in space applications. We are presenting a new particle track and DNA damage model, in which the particle stochastic track structure is combined with the random walk (RW) structure of chromosomes in a cell nucleus. The motivation for this effort stems from the fact that the model with the RW chromosomes, NASARTI (NASA radiation track image) previously relied on amorphous track structure, while the stochastic track structure model RITRACKS (Relativistic Ion Tracks) was focused on more microscopic targets than the entire genome. We have combined chromosomes simulated by RWs with stochastic track structure, which uses nanoscopic dose calculations performed with the Monte‐Carlo simulation by RITRACKS in a voxelized space. The new simulations produce the number of DSBs as function of dose and particle fluence for high‐energy particles, including iron, carbon and protons, using voxels of 20 nm dimension. The combined model also calculates yields of radiation‐induced CAs and unrejoined chromosome breaks in normal and repair deficient cells. The joined computational model is calibrated using the relative frequencies and distributions of chromosomal aberrations reported in the literature. The model considers fractionated deposition of energy to approximate dose rates of the space flight environment. The joined model also predicts of the yields and sizes of translocations, dicentrics, rings, and more complex‐type aberrations formed in the G0/G1 cell cycle phase during the first cell division after irradiation. We found that the main advantage of the joined model is our ability to simulate small doses: 0.05‐0.5 Gy. At such low doses, the stochastic track structure proved to be indispensable, as the action of individual delta‐rays becomes more important.

Ponomarev, Artem↗

Simulations of DSB Yields and Radiation-induced Chromosomal Aberrations in Human Cells Based on the Stochastic Track Structure iIduced by HZE Particles

The formation of double‐strand breaks (DSBs) and chromosomal aberrations (CAs) is of great importance in radiation research and, specifically, in space applications. We are presenting a new particle track and DNA damage model, in which the particle stochastic track structure is combined with the random walk (RW) structure of chromosomes in a cell nucleus. The motivation for this effort stems from the fact that the model with the RW chromosomes, NASARTI (NASA radiation track image) previously relied on amorphous track structure, while the stochastic track structure model RITRACKS (Relativistic Ion Tracks) was focused on more microscopic targets than the entire genome. We have combined chromosomes simulated by RWs with stochastic track structure, which uses nanoscopic dose calculations performed with the Monte‐Carlo simulation by RITRACKS in a voxelized space. The new simulations produce the number of DSBs as function of dose and particle fluence for high‐energy particles, including iron, carbon and protons, using voxels of 20 nm dimension. The combined model also calculates yields of radiation‐induced CAs and unrejoined chromosome breaks in normal and repair deficient cells. The joined computational model is calibrated using the relative frequencies and distributions of chromosomal aberrations reported in the literature. The model considers fractionated deposition of energy to approximate dose rates of the space flight environment. The joined model also predicts of the yields and sizes of translocations, dicentrics, rings, and more complex‐type aberrations formed in the G0/G1 cell cycle phase during the first cell division after irradiation. We found that the main advantage of the joined model is our ability to simulate small doses: 0.05‐0.5 Gy. At such low doses, the stochastic track structure proved to be indispensable, as the action of individual delta‐rays becomes more important.

Ponomarev, Artem↗