Search NASASearch

SEARCH · Search NASA

Results for “Random walk model”

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

Physical interrelation between Fokker-Planck and random walk models with application to Coulomb interactions.

A model of the random walk is formulated to allow a simple computing procedure to replace the difficult problem of solution of the Fokker-Planck equation. The step sizes and probabilities of taking steps in the various directions are expressed in terms of Fokker-Planck coefficients. Application is made to many particle systems with Coulomb interactions. The relaxation of a highly peaked velocity distribution of particles to equilibrium conditions is illustrated.

Englert, G. W.

Neutral cometary atmospheres. I - An average random walk model for photodissociation in comets

A method for constructing physically realistic photochemical models, taking into account the isotropic ejection of dissociated molecular fragments as well as radiation pressure acceleration, has been developed using Monte Carlo techniques. The effect of the isotropic ejection, as opposed to the radial motion arbitrarily assumed in Haser's model, is adequately described by a simple average random walk model. It is shown that measured radial (Haser) scale lengths are in fact only lower limits to a range of possible true scale lengths for a given brightness profile, which explains the current discrepancies between observed scale lengths and those predicted by photochemistry.

Combi, M. R.

A polymer, random walk model for the size-distribution of large DNA fragments after high linear energy transfer radiation

DNA double-strand breaks (DSBs) produced by densely ionizing radiation are not located randomly in the genome: recent data indicate DSB clustering along chromosomes. Stochastic DSB clustering at large scales, from > 100 Mbp down to < 0.01 Mbp, is modeled using computer simulations and analytic equations. A random-walk, coarse-grained polymer model for chromatin is combined with a simple track structure model in Monte Carlo software called DNAbreak and is applied to data on alpha-particle irradiation of V-79 cells. The chromatin model neglects molecular details but systematically incorporates an increase in average spatial separation between two DNA loci as the number of base-pairs between the loci increases. Fragment-size distributions obtained using DNAbreak match data on large fragments about as well as distributions previously obtained with a less mechanistic approach. Dose-response relations, linear at small doses of high linear energy transfer (LET) radiation, are obtained. They are found to be non-linear when the dose becomes so large that there is a significant probability of overlapping or close juxtaposition, along one chromosome, for different DSB clusters from different tracks. The non-linearity is more evident for large fragments than for small. The DNAbreak results furnish an example of the RLC (randomly located clusters) analytic formalism, which generalizes the broken-stick fragment-size distribution of the random-breakage model that is often applied to low-LET data.

DNA/radiation effects

The rotation of photospheric magnetic fields: A random walk transport model

In an earlier study of solar differential rotation, we showed that the transport of magnetic flux across latitudes acts to establish quasi-stationary patterns, therby accounting for the observed rigid rotation of the large-scale photospheric field. In that paper, the effect of supergranular convection was represented by a continuum diffusion, limiting the applicability of the calculations to large spatial scales. Here we extend the model to scales comparable to that of the supergranulation itself by replacing the diffusive transport with a discrete random walk process. Rotation curves are derived by cross-correlating the simulated photospheric field maps for a variety of time lags and spatial resolutions. When the lag between maps is relatively short less than or approximately = 15 days), the midlatitude correlation functions show two distinct components: a broad feature associated with the large-scale unipolar patterns and a narrow feature originating from small magnetic structures encompossing from one to several supergranular cells. By fitting the broad component we obtain the rigid rotation profile of the patterns, whereas by fitting the narrow component, we recover the differential rate of the photospheric plasma itself. For time lags of 1 month or greater, only the broad feature associated with the long-lived patterns remains clearly identifiable in the simulations.

Wang, Y. -M.

Migration of lymphocytes on fibronectin-coated surfaces: temporal evolution of migratory parameters

Lymphocytes typically interact with implanted biomaterials through adsorbed exogenous proteins. To provide a more complete characterization of these interactions, analysis of lymphocyte migration on adsorbed extracellular matrix proteins must accompany the commonly performed adhesion studies. We report here a comparison of the migratory and adhesion behavior of Jurkat cells (a T lymphoblastoid cell line) on tissue culture treated and untreated polystyrene surfaces coated with various concentrations of fibronectin. The average speed of cell locomotion showed a biphasic response to substrate adhesiveness for cells migrating on untreated polystyrene and a monotonic decrease for cells migrating on tissue culture-treated polystyrene. A modified approach to the persistent random walk model was implemented to determine the time dependence of cell migration parameters. The random motility coefficient showed significant increases with time when cells migrated on tissue culture-treated polystyrene surfaces, while it remained relatively constant for experiments with untreated polystyrene plates. Finally, a cell migration computer model was developed to verify our modified persistent random walk analysis. Simulation results suggest that our experimental data were consistent with temporally increasing random motility coefficients.

NASA Program Fundamental Space Biology

Monte Carlo predictions of DNA fragment-size distributions for large sizes after HZE particle irradiation

DSBs (double-strand breaks) produced by densely ionizing space radiation are not located randomly in the genome: recent data indicate DSB clustering along chromosomes. DSB clustering at large scales, from >100 Mbp down to approximately 2 kbp, is modeled using a Monte-Carlo algorithm. A random-walk model of chromatin is combined with a track model, that predicts the radial distribution of energy from an ion, and the RLC (randomly-located-clusters) formalism, in software called DNAbreak. This model generalizes the random-breakage model, whose broken-stick fragment-size distribution is applicable to low-LET radiation. DSB induction due to track interaction with the DNA volume depends on the radiation quality parameter Q. This dose-independent parameter depends only weakly on LET. Multi-track, high-dose effects depend on the cluster intensity parameter lambda, proportional to fluence as defined by the RLC formalism. After lambda is determined by a numerical experiment, the model reduces to one adjustable parameter Q. The best numerical fits to the experimental data, determining Q, are obtained. The knowledge of lambda and Q allows us to give biophysically based extrapolations of high-dose DNA fragment-size data to low doses or to high LETs.

NASA Center JSC

Centric rings, acentric rings and excess acentric fragments based on a random-walk interphase chromosome model

Excess acentric fragments, consisting of acentric rings and acentric linear fragments, are among the most frequent kinds of chromosome-type aberrations produced by radiation. The frequency of acentric rings cannot be obtained directly by experiment but is estimated here from the ratio of acentric to centric rings, evaluated using a random-walk model for the organization of chromatin during interphase and an assumption that the probability of an exchange formation is proportional to the rate of collision between two DSB. This ratio is calculated to be 2.5 in low-LET irradiated human fibroblasts, significantly greater than the ratio if proximity effects are not considered. The calculated frequency of acentric rings is insufficient to account for all the observed excess acentric fragments. Assuming that the rest of the excess acentric fragments are due to incomplete exchanges, all possible recombinations between two DSB that result in acentric rings and acentric linear fragments have been identified. From the chromosome aberration data, the incompleteness parameter has been estimated. Intra-arm chromosome exchanges, either complete or incomplete, were estimated to account for more than 50% of the excess acentric fragments in human fibroblasts.

NASA Center JSC

Radiation-induced total-deletion mutations in the human hprt gene: a biophysical model based on random walk interphase chromatin geometry

PURPOSE: To develop a biophysical model that explains the sizes of radiation-induced hprt deletions. METHODS: Key assumptions: (1) Deletions are produced by two DSB that are closer than an interaction distance at the time of DSB induction; (2) Interphase chromatin is modelled by a biphasic random walk distribution; and (3) Misrejoining of DSB from two separate tracks dominates at low-LET and misrejoining of DSB from a single track dominates at high-LET. RESULTS: The size spectra for radiation-induced total deletions of the hprt gene are calculated. Comparing with the results of Yamada and coworkers for gamma-irradiated human fibroblasts the study finds that an interaction distance of 0.75 microm will fit both the absolute frequency and the size spectrum of the total deletions. It is also shown that high-LET radiations produce, relatively, more total deletions of sizes below 0.5 Mb. The model predicts an essential gene to be located between 2 and 3 Mb from the hprt locus towards the centromere. Using the same assumptions and parameters as for evaluating mutation frequencies, a frequency of intra-arm chromosome deletions is calculated that is in agreement with experimental data. CONCLUSIONS: Radiation-induced total-deletion mutations of the human hprt gene and intrachange chromosome aberrations share a common mechanism for their induction.

NASA Center JSC

Investigation for improving Global Positioning System (GPS) orbits using a discrete sequential estimator and stochastic models of selected physical processes

GEODYNII is a conventional batch least-squares differential corrector computer program with deterministic models of the physical environment. Conventional algorithms were used to process differenced phase and pseudorange data to determine eight-day Global Positioning system (GPS) orbits with several meter accuracy. However, random physical processes drive the errors whose magnitudes prevent improving the GPS orbit accuracy. To improve the orbit accuracy, these random processes should be modeled stochastically. The conventional batch least-squares algorithm cannot accommodate stochastic models, only a stochastic estimation algorithm is suitable, such as a sequential filter/smoother. Also, GEODYNII cannot currently model the correlation among data values. Differenced pseudorange, and especially differenced phase, are precise data types that can be used to improve the GPS orbit precision. To overcome these limitations and improve the accuracy of GPS orbits computed using GEODYNII, we proposed to develop a sequential stochastic filter/smoother processor by using GEODYNII as a type of trajectory preprocessor. Our proposed processor is now completed. It contains a correlated double difference range processing capability, first order Gauss Markov models for the solar radiation pressure scale coefficient and y-bias acceleration, and a random walk model for the tropospheric refraction correction. The development approach was to interface the standard GEODYNII output files (measurement partials and variationals) with software modules containing the stochastic estimator, the stochastic models, and a double differenced phase range processing routine. Thus, no modifications to the original GEODYNII software were required. A schematic of the development is shown. The observational data are edited in the preprocessor and the data are passed to GEODYNII as one of its standard data types. A reference orbit is determined using GEODYNII as a batch least-squares processor and the GEODYNII measurement partial (FTN90) and variational (FTN80, V-matrix) files are generated. These two files along with a control statement file and a satellite identification and mass file are passed to the filter/smoother to estimate time-varying parameter states at each epoch, improved satellite initial elements, and improved estimates of constant parameters.

Goad, Clyde C.

Fluctuation and rotation in diffusion-influenced monomolecular reactions

Models of partially absorbing boundaries for calculating diffusion-influenced chemical reactions are examined analytically. The limitations of diffusion-equation and random-walk models are indicated; the presence of an ambiguous length parameter in a standard radiation-boundary-condition (RBC) model is demonstrated; the RBC model for a boundary with oscillation (at a rate much faster than the rate of thermal collisions) between absorbing and nonabsorbing states is considered; an alternative stochastically switching absorber (SSA) model is introduced and shown to be physically consistent; the equivalence of one-dimensional RBC and SSA models under some conditions is shown; one time passage is calculated for an axially symmetric three-dimensional SSA model with anisotropic absorption; and the analogy between rotation and stochastic switching is explored.

Bashford, Donald

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

Analysis of Fibrous Felts for Flexible Ablators Using Synchrotron Hard X-Ray Micro-Tomography

We analyzed the material properties of low-density felts that are used as substrates for new-generation flexible and conformal carbon/phenolic ablators, and compared them with those of a rigid carbon fiber preform that is used to manufacture rigid carbon/phenolic ablators. Micro-tomography measurements were obtained using synchrotron X-rays, allowing the characterization of the materials microstructure at the scale of the fibers. Using the tomography voxels as computational grids, we computed tortuosity and room temperature conductivity. In addition we performed micro-scale simulations of the oxidation of carbon fibers using a random walk model for oxygen diffusion and a probability law to model surface reactions.

Ablation

Space infrared telescope pointing control system. Infrared telescope tracking in the presence of target motion

The use of charge-coupled-devices, or CCD's, has been documented by a number of sources as an effective means of providing a measurement of spacecraft attitude with respect to the stars. A method exists of defocussing and interpolation of the resulting shape of a star image over a small subsection of a large CCD array. This yields an increase in the accuracy of the device by better than an order of magnitude over the case when the star image is focussed upon a single CCD pixel. This research examines the effect that image motion has upon the overall precision of this star sensor when applied to an orbiting infrared observatory. While CCD's collect energy within the visible spectrum of light, the targets of scientific interest may well have no appreciable visible emissions. Image motion has the effect of smearing the image of the star in the direction of motion during a particular sampling interval. The presence of image motion is incorporated into a Kalman filter for the system, and it is shown that the addition of a gyro command term is adequate to compensate for the effect of image motion in the measurement. The updated gyro model is included in this analysis, but has natural frequencies faster than the projected star tracker sample rate for dim stars. The system state equations are reduced by modelling gyro drift as a white noise process. There exists a tradeoff in selected star tracker sample time between the CCD, which has improved noise characteristics as sample time increases, and the gyro, which will potentially drift further between long attitude updates. A sample time which minimizes pointing estimation error exists for the random drift gyro model as well as for a random walk gyro model.

Powell, J. D.

Radiation breakage of DNA: a model based on random-walk chromatin structure

Monte Carlo computer software, called DNAbreak, has recently been developed to analyze observed non-random clustering of DNA double strand breaks in chromatin after exposure to densely ionizing radiation. The software models coarse-grained configurations of chromatin and radiation tracks, small-scale details being suppressed in order to obtain statistical results for larger scales, up to the size of a whole chromosome. We here give an analytic counterpart of the numerical model, useful for benchmarks, for elucidating the numerical results, for analyzing the assumptions of a more general but less mechanistic "randomly-located-clusters" formalism, and, potentially, for speeding up the calculations. The equations characterize multi-track DNA fragment-size distributions in terms of one-track action; an important step in extrapolating high-dose laboratory results to the much lower doses of main interest in environmental or occupational risk estimation. The approach can utilize the experimental information on DNA fragment-size distributions to draw inferences about large-scale chromatin geometry during cell-cycle interphase.

Chromatin/radiation effects

Simulation of Inflated Pahoehoe Lava Flows

A new stochastic model simulates late-stage pahoehoe lobes where random processes dominate emplacement. The model prescribes probabilistic rules for determining where and when parcels of lava move within the lobe. Unlike a classical Brownian motion random walk, the model allows individual parcels to remain dormant, but fluid, for multiple time steps. The randomness of parcel volume transfers within the lobe interior as well as at the margins qualitatively reflects inflation processes observed in the field. The fraction of inflated volume to total volume increases with the total volume, with greater than 75% of the lobe volume contributed through inflation for typical lobes. The influence on planform shape and topographic cross-sectional profiles of total volume, source area and shape, topographic confinement, and sequential breakouts at the lobe margins, are all explored with the stochastic model. Each of these factors influences the overall lobe thickness and width. The model provides a means for assessing the relative importance of these processes through comparisons with field data. For the first time, Gaussian and parabolic functions are quantitatively fit to field measurements of pahoehoe lobes. Both functional forms provide adequate description of the cross-sectional flow shapes. When comparing simulated lobes to field data, sequential breakouts at the lobe margins are found to be an important process controlling the final topographic distribution of observed pahoehoe lobes.

modeling