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.
A biased random walk, Monte Carlo procedure is used to study the melting transition in small clusters of N atoms. Results show an abrupt change in physical properties near melting, but there is no evidence that the transition is either first or second order; rather the transition is more gradual.
An approximate analytic model, based on continuous electron slowing, has been used for survey calculations. Where more accuracy is required, a Monte Carlo technique is used which combines an analytic representation of Coulombic collisions with a random walk treatment of inelastic collisions. The calculated electron distributions have been incorporated into another code that evaluates both the excited atomic state densities within the plasma and the radiative flux emitted from the plasma.
The recent development of low-loss single-mode optical fiber waveguides for light has made possible a new class of inertial reference devices built on the principle of a closed loop interferometer. Light circulating through the loop in both directions experiences a relative phase delay proportional to rotation rate about the loop axis. This paper derives the phase delay and discusses signal detection, signal processing techniques and error sources. It is concluded that synchronous modulation and demodulation and an active gain control at the signal calculation level are required to eliminate drift errors. Potential performance is extraordinarily good; rotation rate sensitivity of a few milli-arc seconds per second and angular position random walk errors of an arc second per square root hour appear feasible.
It is shown that the global nature of trapped particle motion in phase space makes it entirely plausible that particles with pitch angles near 90 deg are trapped for extremely long times, probably indefinitely. The long observed lifetimes of trapped particles can immediately be accounted for by this former population. For particles with low pitch angles an 'overstability' develops which serves to produce rapid variations in pitch angle and, presumably, random walk into the loss cone.
A generalized method for obtaining individual level population ratios is used to obtain relative intensities of extreme ultraviolet Fe XV emission lines in the range 284-500 A, which are density dependent for electron densities in the tokamak regime or higher. Four lines in particular are found to attain quite high intensities in the high-density limit. The same calculation provides inelastic contributions to linewidths. The method connects level populations and level widths through total probabilities t(ij), related to 'taboo' probabilities of Markov chain theory. The t(ij) are here evaluated for a real atomic system, being therefore of potential interest to random-walk theorists who have been limited to idealized systems characterized by simplified transition schemes.
A model is constructed which views regolith evolution on asteroids as a stochastic process. Average values are shown to be poor descriptors of regolith depth. The utility of the average depth is not significantly increased by avoiding large craters or thick ejecta deposits, a procedure adopted in previous regolith studies. The statistical uncertainty associated with regolith depth severely limits the power of regolith models in predicting parent-body size for brecciated meteorites. A Monte Carlo algorithm was used to simulate the random walks and corresponding charged-particle irradiation histories of grains in regoliths. On rocky asteroids, only about 20 percent of the grains was exposed to solar cosmic ray ions. Results based on present-day conditions in the asteroid belt agree well with irradiation features observed in gas-rich meteorites. An origin during epochs of early solar system evolution is not required.
Mean and fluctuating gas velocities, liquid mass fluxes and drop sizes were in nonevaporating sprays. These results, as well as existing measurements in solid particle-laden jets, were used to evaluate models of these processes. The following models were considered: (1) a locally homogeneous flow (LHF) model, where slip between the phases was neglected; (2) a deterministic separated flow (DSF) model, where slip was considered but effects of particle dispersion by turbulence were ignored; and (3) a stochastic separated flow (SSF) model, where effects of interphase slip and turbulent dispersion were considered using random-walk computations for particle motion. The LHF and DSF models did not provide very satisfactory predictions over the present data base. In contrast, the SSF model performed reasonably well - including conditions in nonevaporating sprays where enhanced dispersion of particles by turbulence caused the spray to spread more rapidly than single-phase jets for comparable conditions. While these results are encouraging, uncertainties in initial conditions limit the reliability of the evaluation. Current work is seeking to eliminate this deficiency.
The structure of particle-laden jets and nonevaporating and evaporating sprays was measured in order to evaluate models of these processes. Three models are being evaluated: (1) a locally homogeneous flow model, where slip between the phases is neglected and the flow is assumed to be in local thermodynamic equilibrium; (2) a deterministic separated flow model, where slip and finite interphase transport rates are considered but effects of particle/drop dispersion by turbulence and effects of turbulence on interphase transport rates are ignored; and (3) a stochastic separated flow model, where effects of interphase slip, turbulent dispersion and turbulent fluctuations are considered using random sampling for turbulence properties in conjunction with random-walk computations for particle motion. All three models use a k-e-g turbulence model. All testing and data reduction are completed for the particle laden jets. Mean and fluctuating velocities of the continuous phase and mean mixture fraction were measured in the evaporating sprays.
Continuous-time models of oscillator phase noise x(t) usually have stationary nth differences, for some n. The covariance structure of such a model can be characterized in the time domain by the structure function: D sub n (t;gamma sub 1, gamma sub 2) = E delta (n) sub gamma sub 1 x(s+t) delta(n) sub gamma sub 2 x (s). Although formulas for the special case D sub 2 (0;gamma,gamma) (the Allan variance times 2 gamma(2)) exist for power-law spectral models, certain estimation problems require a more complete knowledge of (0). Exhibited is a much simpler function of one time variable, D(t), from which (0) can easily be obtained from the spectral density by uncomplicated integrations. Believing that D(t) is the simplest function of time that holds the same information as (0), D(t) is called the fundamental structure function. D(t) is computed for several power-law spectral models. Two examples are D(t) = K/t/(3) for random walk FM, D(t) = Kt(2) 1n/t/ for flicker FM. Then, to demonstrate its use, a BASIC program is given that computes means and variances of two Allan variance estimators, one of which incorporates a method of frequency drift estimation and removal.
A byproduct of stratospheric motion distinct from the normal downgradient transport represented by vertical eddy diffusion is studied. It is suggested that an air parcel executes a random walk over a range that could be several kilometers, and similar amplitudes could be present in wavelike oscillations. The photochemical loss rate should therefore be averaged over this range, which give an effect that is proportional to the curvature of the profile of loss rate, and is particularly large for fluorocarbon-11, CFCl3. Test computations are presented which indicate that the known discrepancy between observations and predictions for this compound can be at least partly reconciled for plausible values of the amplitude.
A passive hydrogen maser employing a smaller microwave cavity than previous masers is described. The volume of the entire resonator package, including a conventional source, hexapole state selector, and four magnetic shields, is estimated at 20 liters. The techniques for constructing the ceramic cavity, coating it with FEP, and attaching the endcaps to achieve a stable container (to allow its use as a single vacuum system) are described in detail. Frequency stability measurements were made of four such masers, against an NBS-4 cesium primary standard, against each other, and against eight commercial cesium clocks; frequency stability level is reported as 1.7 x 10 to the -12th tau exp -1/2 from 1 to 100,000 sec, and no drift above 1 x 10 to the -15th/day was found. Pressure and temperature effects on frequency stability were found at the 1 x 10 to the -14th level. Noise behavior beyond 200,000 sec is characterized as random walk. The potential use of a compact stable maser as a clock is suggested.
A numerical, grid-free algorithm is presented for one-dimensional reaction-diffusion model of laminar flame propagation in premixed gases. It is based on the random element method we developed for the analysis of diffusional processes. The effect of combustion is taken into account by applying the principle of fractional steps to separate the process of diffusion, modeled by the random walk of computational elements, from the exothermic effects of chemical reaction, monitoring their strength. The validity of the algorithm is demonstrated by application to flame propagation problems for which exact solutions exist. The flame speed evaluated by its use oscillates around the exact value at a relatively small amplitude, while the temperature and species concentration profiles are self-correcting in their convergence to the exact solution. A satisfactory resolution is obtained by the use of quite a small number of computational elements which automatically adjust their distribution of fit sharp gradients.
Measurements and analysis limited to the dilute portions of turbulent nonevaporating sprays injected into a still air environment were completed. Mean and fluctuating velocities and Reynolds stress were measured in the continuous phase. Liquid phase measurements included liquid mass fluxes, drop sizes and drop size and velocity correlation. Initial conditions needed for model evaluation were measured at a location as close to the injector exit as possible. The test sprays showed significant effects of slip and turbulent dispersion of the discrete phase. The measurements were used to evaluate three typical models of these processes: (1) a locally homogenous flow (LHF) model, where slip between the phases were neglected; (2) a deterministic separated flow (DSF) model, where slip was considered but effects of drop dispersion by turbulence were ignored; and (3) a stochastic separated flow (SSF) model, where effects of interphase slip and turbulent dispersion were considered using random-walk computations for drop motion. The LHF and DSF models did not provide very satisfactory predictions for the present measurements. In contrast, the SSF model performed reasonably well with no modifications in the prescription of eddy properties from its original calibration. Some effects of drops on turbulence properties were observed near the dense regions of the sprays.
Measurements and analysis limited to the dilute portions of turbulent evaporating sprays, injected into a still air environment were completed. Mean and fluctuating velocities and Reynolds stress were measured in the continuous phase. Liquid phase measurements included liquid mass fluxes, drop sizes and drop size and velocity correlation. Initial conditions needed for model evaluation were measured at a location as close to the injector exit as possible. The test sprays showed significant effects of slip and turbulent dispersion of the discrete phase. The measurements were used to evaluate three typical models of these processes: (1) a locally homogeneous flow (LHF) model, where slip between the phases were neglected; (2) a deterministic separated flow (DSF) model, where slip was considered but effects of drop dispersion by turbulence were ignored; and (3) a stochastic separated flow (SSF) model, where effects of interphase slip and turbulent dispersion were considered using random-walk computations for drop motion. For all three models, a k-epsilon model as used to find the properties of the continuous phase. The LHF and DSF models did not provide very satisfactory predictions for the present measurements. In contrast, the SSF model performed reasonably well--with no modifications in the prescription of eddy properties from its original calibration.
Measurements of mean and fluctuating velocities of both phases as well as particle mass fluxes were completed in turbulent, particle-laden jets containing monodisperse particles with well-defined initial and boundary conditions. The new measurements were used to evaluate a stochastic separated flow model of the process which treated effects of interphase slip and turbulent dispersion using random-walk computations for particle motion. The continuous phase was treated using a modified k-epsilon model allowing for direct contributions of interphase transport to both mean and turbulence properties. The model performed reasonably well over the new data base, with all empirical parameters fixed from earlier work. In contrast, simplified models ignoring either interphase slip or turbulent dispersion yielded poor agreement with the measurements.
Structure measurements were completed within the dilute portion of axisymmetric nonevaporating sprays (SMD of 30 and 87 microns) injected into a still air environment, including: mean and fluctuating gas velocities and Reynolds stress using laser-Doppler anemometry; mean liquid fluxes using isokinetic sampling; drop sizes using slide impaction; and drop sizes and velocities using multiflash photography. The new measurements were used to evaluate three representative models of sprays: (1) a locally homogeneous flow (LHF) model, where slip between the phases was neglected; (2) a deterministic separated flow (DSF) model, where slip was considered but effects of drop interaction with turbulent fluctuations were ignored; and (3) a stochastic separated flow (SSF) model, where effects of both interphase slip and turbulent fluctuations were considered using random sampling for turbulence properties in conjunction with random-walk computations for drop motion. The LHF and DSF models were unsatisfactory for present test conditions-both underestimating flow widths and the rate of spread of drops. In contrast, the SSF model provided reasonably accurate predictions, including effects of enhanced spreading rates of sprays due to drop dispersion by turbulence, with all empirical parameters fixed from earlier work.
A recent argument dismissing the applicability of large-cell renormalization schemes to systems whose dimensionality is larger than the upper marginal dimension is critically discussed. In this connection, new large-cell renormalization results for the random walk for a dimensionality of 3 and 4 are presented which indicate convergence to the correct results.
The random element method is a generalization of the random vortex method that was developed for the numerical modeling of momentum transport processes as expressed in terms of the Navier-Stokes equations. The method is based on the concept that random walk, as exemplified by Brownian motion, is the stochastic manifestation of diffusional processes. The algorithm based on this method is grid-free and does not require the diffusion equation to be discritized over a mesh, it is thus devoid of numerical diffusion associated with finite difference methods. Moreover, the algorithm is self-adaptive in space and explicit in time, resulting in an improved numerical resolution of gradients as well as a simple and efficient computational procedure. The method is applied here to an assortment of problems of diffusion of momentum and energy in one-dimension as well as heat conduction in two-dimensions in order to assess its validity and accuracy. The numerical solutions obtained are found to be in good agreement with exact solution except for a statistical error introduced by using a finite number of elements, the error can be reduced by increasing the number of elements or by using ensemble averaging over a number of solutions.