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 Monte Carlo calculation of the quark-liberating phase transition in lattice quantum chromodynamics is presented. The transition temperature as a function of the lattice coupling g does not scale according to the perturbative beta function for 6/g-squared less than 6.1. Finite-size scaling is used in analyzing the properties of the lattice system near the transition point.
We performed theoretical studies of liquid helium by applying state of the art simulation and finite-size scaling techniques. We calculated universal scaling functions for the specific heat and superfluid density for various confining geometries relevant for experiments such as the confined helium experiment and other ground based studies. We also studied microscopically how the substrate imposes a boundary condition on the superfluid order parameter as the superfluid film grows layer by layer. Using path-integral Monte Carlo, a quantum Monte Carlo simulation method, we investigated the rich phase diagram of helium monolayer, bilayer and multilayer on a substrate such as graphite. We find excellent agreement with the experimental results using no free parameters. Finally, we carried out preliminary calculations of transport coefficients such as the thermal conductivity for bulk or confined helium systems and of their scaling properties. All our studies provide theoretical support for various experimental studies in microgravity.
We carry out state of the art simulations of properties of confined liquid helium near the superfluid transition to a degree of accuracy which allows to make predictions for the outcome of fundamental physics experiments in microgravity. First we report our results for the finite-size scaling behavior of heat capacity of superfluids for cubic and parallel-plate geometry. This allows us to study the crossover from zero and two dimensions to three dimensions. Our calculated scaling functions are in good agreement with recently measured specific heat scaling functions for the above mentioned geometries. We also present our results of a quantum simulation of submonolayer of molecular hydrogen deposited on an ideal graphite substrate using path-integral quantum Monte Carlo simulation. We find that the monolayer phase diagram is rich and very similar to that of helium monolayer. We are able to uncover the main features of the complex monolayer phase diagram, such as the commensurate solid phases and the commensurate to incommensurate transition, in agreement with the experiments and to find some features which are missing from the experimental analysis.
Himalayan avalanches are examined for agreement with power-law behavior between frequency of outfall and rock volume. The data support a scaling state beyond the size ranges in which corresponding sandpile experiments have suffered from finite-size effects. Large-scale data (12 orders of magnitude in rock volume) additionally suggest a possible role for the numerically predicted 'snowball effects' of large rocks triggering more net dislodging events.
The scaling behavior of sandpile models is investigated analytically. First, it is shown that sandpile models contain a set of domain walls, referred to as troughs, which bound regions that can experience avalanches. It is further shown that the dynamics of the troughs is governed by a simple set of rules involving birth, death, and coalescence events. A simple trough model is then introduced, and it is proved that the model has a phase transition with the density of the troughs as an order parameter and that, in the thermodynamic limit, the trough density goes to zero at the transition point. Finally, it is shown that the observed scaling behavior is a consequence of finite-size effects.
An elecrostatic finite-sized particle simulation model is used to study the early time-scale phenomena asociated with an impulsively injected plasma, expanding radially, normal to a strong ambient magnetic field. Results of the simulation show the early formation of a radial polarization electric field due to ion-electron charge separation which causes electrons to drift in the azimuthal direction. Velocity shear within this motion gives rise to a Kelvin-Helmholtz (dioxotron) instability, creating a radial fluted pattern as the plasma expands. In the nonlinear stage of instability, azimuthal variations in the electric field are observed which cause electrons to drift across the magnetic field, thereby reducing the space charge created by energetic ion expansion. The results of a linear stability analysis reveal that for a decreasing amount of charge separation in the plasma, the number of unstable azimuthal modes at maximum growth rate increases. Also, as the thickness of the electron ring increases, a fewer number of unstable modes develops for a given amount of charge separation. It is concluded that electrons become unstable at the very early time scale, within an ion gyroperiod, and that unstable electrons lead to modification of the ion dynamics.
We analyze the fluctuations in the correlation exponents obtained for noncoding DNA sequences. We find prominent sample-to-sample variations as well as variations within a single sample in the scaling exponent. To determine if these fluctuations may result from finite system size, we generate correlated random sequences of comparable length and study the fluctuations in this control system. We find that the DNA exponent fluctuations are consistent with those obtained from the control sequences having long-range power-law correlations. Finally, we compare our exponents for the DNA sequences with the exponents obtained from power-spectrum analysis and correlation-function techniques, and demonstrate that the original "DNA-walk" method is intrinsically more accurate due to reduced noise.
The hydrodynamics of gas flow past a finite-sized gravitating central object is studied in the adiabatic approximation in two spatial dimensions. The flow morphology is shown to depend on the ratio (epsilon) of the axisymmetric accretion radius to the length scale of the density variation. For small to intermediate epsilon values the flows are highly time dependent, and for large epsilon values a quasi-steady state is obtained.
The paper reports on the results of hydrodynamical simulations of nonaxisymmetric gas flow past a finite-sized gravitating object. The asymmetry in the flow is assumed to be due to a transverse density gradient. The solutions are found to be highly time-dependent and exhibit phases in which a disk forms. The direction of circulation of the disk alternates between the direct and counter sense on time scales of approximately 1 hr. While the disk is present, the specific angular momentum is high, and the mass accretion rate, low. The angular momentum capture oscillates about a mean which secularly increases. Such temporal behavior may be relevant to those X-ray pulsing systems (e.g., Vela X-1) which exhibit fluctuations in the X-ray luminosity and pulse period on time scales much shorter than the orbital period of the binary system.
An explanation of the dynamical mechanism for apse-alignment of the eccentric Uranian rings is necessary before observations can be used to determine properties such as ring masses, particle sizes, and elasticities. The leading model relies on the ring self-gravity to accomplish this task, yet it yields equilibrium masses which are not in accord with Voyager radio measurements. We explore possible solutions such that the self-gravity and the collisional terms are both involved in the process of apse-alignment. We consider limits that correspond to a hot and a cold ring, and show that pressure terms may play a significant role in the equilibrium conditions for the narrow Uranian rings. In the cold ring case, where the scale height of the ring near periapse is comparable to the ring particle size, we introduce a new pressure correction pertaining to a region of the ring where the particles are locked in their relative positions and jammed against their neighbors, and the velocity dispersion is so low that the collisions are nearly elastic. In this case, we find a solution such that the ring self-gravity maintains apse-alignment against both differential precession (m = 1 mode) and the fluid pressure. We apply this model to the Uranian alpha ring, and show that, compared to the previous self-gravity model, the mass estimate for this ring increases by an order of magnitude. In the case of a hot ring, where the scale height can reach a value as much as fifty times larger than a particle size, we find velocity dispersion profiles that result in pressure forces which act in such a way as to alter the ring equilibrium conditions, again leading to a ring mass increase of an order of magnitude; however, such a velocity dispersion profile would require a different mechanism than is currently envisioned for establishing heating/cooling balance in a finite-sized, inelastic particle ring. Finally, we introduce an important correction to the model of Chiang and Goldreich.