Variance Reduction For A Discrete Velocity Gas

We extend a variance reduction technique developed by Baker and Hadjiconstantinou [1] to a discrete velocity gas. In our previous work, the collision integral was evaluated by importance sampling of collision partners [2]. Significant computational effort may be wasted by evaluating the collision integral in regions where the flow is in equilibrium. In the current approach, substantial computational savings are obtained by only solving for the deviations from equilibrium. In the near continuum regime, the deviations from equilibrium are small and low noise evaluation of the collision integral can be achieved with very coarse statistical sampling. Spatially homogenous relaxation of the Bobylev-Krook-Wu distribution [3,4], was used as a test case to verify that the method predicts the correct evolution of a highly non-equilibrium distribution to equilibrium. When variance reduction is not used, the noise causes the entropy to undershoot, but the method with variance reduction matches the analytic curve for the same number of collisions. We then extend the work to travelling shock waves and compare the accuracy and computational savings of the variance reduction method to DSMC over Mach numbers ranging from 1.2 to 10.Aerospace Engineering and Engineering Mechanic

Far Field Deposition Of Scoured Regolith Resulting From Lunar Landings

As a lunar lander approaches a dusty surface, the plume from the descent engine impinges on the ground, entraining loose regolith into a high velocity dust spray. Without the inhibition of a background atmosphere, the entrained regolith can travel many kilometers from the landing site. In this work, we simulate the flow field from the throat of the descent engine nozzle to where the dust grains impact the surface many kilometers away. The near field is either continuum or marginally rarefied and is simulated via a loosely coupled hybrid DSMC - Navier Stokes (DPLR) solver. Regions of two-phase and polydisperse granular flows are solved via DSMC. The far field deposition is obtained by using a staged calculation, where the first stages are in the near field where the flow is quasi-steady and the outer stages are unsteady. A realistic landing trajectory is approximated by a set of discrete hovering altitudes which range from 20m to 3m. The dust and gas motions are fully coupled using an interaction model that conserves mass, momentum, and energy statistically and inelastic collisions between dust particles are also accounted for. Simulations of a 4 engine configuration are also examined, and the erosion rates as well as near field particle fluxes are discussed.Astronom

Cosmic ray modulation and merged interaction regions

Beyond several AU, interactions among shocks and streams give rise to merged interaction regions in which the magnetic field is turbulent. The integral intensity of . 75 MeV/Nuc cosmic rays at Voyager is generally observed to decrease when a merged interaction region moves past the spacecraft and to increase during the passage of a rarefaction region. When the separation between interaction regions is relatively large, the cosmic ray intensity tends to increase on a scale of a few months. This was the case at Voyager 1 from July 1, 1983 to May 1, 1984, when the spacecraft moved from 16.7 to 19.6 AU. Changes in cosmic ray intensity were related to the magnetic field strength in a simple way. It is estimated that the diffusion coefficient in merged interaction regions at this distance is similar to 0.6 x 10 to the 22nd power sq cm/s

Chaotic string-capture by black hole

We consider a macroscopic charge-current carrying (cosmic) string in the background of a Schwarzschild black hole. The string is taken to be circular and is allowed to oscillate and to propagate in the direction perpendicular to its plane (that is parallel to the equatorial plane of the black hole). Nurmerical investigations indicate that the system is non-integrable, but the interaction with the gravitational field of the black hole anyway gives rise to various qualitatively simple processes like "adiabatic capture" and "string transmutation".Comment: 13 pages Latex + 3 figures (not included), Nordita 93/55

Shift Equivalence of Measures and the Intrinsic Structure of Shocks in the Asymmetric Simple Exclusion Process

We investigate properties of non-translation-invariant measures, describing particle systems on \bbz, which are asymptotic to different translation invariant measures on the left and on the right. Often the structure of the transition region can only be observed from a point of view which is random---in particular, configuration dependent. Two such measures will be called shift equivalent if they differ only by the choice of such a viewpoint. We introduce certain quantities, called translation sums, which, under some auxiliary conditions, characterize the equivalence classes. Our prime example is the asymmetric simple exclusion process, for which the measures in question describe the microscopic structure of shocks. In this case we compute explicitly the translation sums and find that shocks generated in different ways---in particular, via initial conditions in an infinite system or by boundary conditions in a finite system---are described by shift equivalent measures. We show also that when the shock in the infinite system is observed from the location of a second class particle, treating this particle either as a first class particle or as an empty site leads to shift equivalent shock measures.Comment: Plain TeX, 2 figures; [email protected], [email protected], [email protected], [email protected]

An east-west asymmetry in the solar wind velocity

East-west asymmetry in solar wind velocit

Sound radiation from a high speed axial flow fan due to the inlet turbulence quadrupole interaction

A formula is obtained for the total acoustic power spectra radiated out the front of the fan as a function of frequency. The formula involves the design parameters of the fan as well as the statistical properties of the incident turbulence. Numerical results are calculated for values of the parameters in the range of interest for quiet fans tested at the Lewis Research Center. As in the dipole analysis, when the turbulence correlation lengths become equal to the interblade spacing, the predicted spectra exhibit peaks around the blade passing frequency and its harmonics. There has recently been considerable conjecture about whether the stretching of turbulent eddies as they enter a stationary fan could result in the inlet turbulence being the dominant source of pure tones from nontranslating fans. The results of the current analysis show that, unless the turbulent eddies become quite elongated, this noise source contributes predominantly to the broadband spectrum

Cosmic ray modulation and turbulent interaction regions near 11 AU

When Voyager 2 was near 11 AU, the counting rate of nuclei approx 75 MeV/nucleon decreased during the interval from July, 1982 to November, 1982, and it increased thereafter until August, 1983. A decrease in cosmic ray flux was generally associated with the passage of an interaction region in which the magnetic field strength B was higher than that predicted by the spiral field model, B sub p. Several large enhancements in B/B sup p were associated with merged interaction regions which probably resulted from the interaction of two or more distinct flows. During the passage of interaction regions the cosmic ray intensity decreased at a rate proportional to (B/B sup p -1), and during the passage of rarefaction regions (where B/B sup p 1) the cosmic ray intensity increased at a constant rate. The general form of the cosmic ray intensity profile during this approx 13 month minicycle can be described by integrating these relations using the observed B(t). Latitudinal variations of the interaction regions and of the short-term cosmic ray variations were identified

A BPS Interpretation of Shape Invariance

We show that shape invariance appears when a quantum mechanical model is invariant under a centrally extended superalgebra endowed with an additional symmetry generator, which we dub the shift operator. The familiar mathematical and physical results of shape invariance then arise from the BPS structure associated with this shift operator. The shift operator also ensures that there is a one-to-one correspondence between the energy levels of such a model and the energies of the BPS-saturating states. These findings thus provide a more comprehensive algebraic setting for understanding shape invariance.Comment: 15 pages, 2 figures, LaTe

DSMC Modeling Of Gasdynamics, Radiation And Fine Particulates In Ionian Volcanic Jets

Aerospace Engineerin