Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density modification”

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 379 records · Page 21

Development of an engineering model atmosphere for Mars

An engineering model atmosphere for Mars is being developed with many of the same features and capabilities for the highly successful Global Reference Atmospheric Model (GRAM) program for Earth's atmosphere. As an initial approach, the model is being built around the Martian atmosphere model computer subroutine (ATMOS) of Culp and Stewart (1984). In a longer-term program of research, additional refinements and modifications will be included. ATMOS includes parameterizations to stimulate the effects of solar activity, seasonal variation, diurnal variation magnitude, dust storm effects, and effects due to the orbital position of Mars. One of the current shortcomings of ATMOS is the neglect of surface variation effects. The longer-term period of research and model building is to address some of these problem areas and provide further improvements in the model (including improved representation of near-surface variations, improved latitude-longitude gradient representation, effects of the large annual variation in surface pressure because of differential condensation/sublimation of the CO2 atmosphere in the polar caps, and effects of Martian atmospheric wave perturbations on the magnitude of the expected density perturbation.

Justus, C. G.↗

A one-dimensional multispecies magnetohydrodynamic model of the dayside ionosphere of Venus

Using a modification of the one-dimensional multispecies 'one-major-ion' MHD model of Shinagawa et al. (1987), the behaviors of plasma and magnetic field in the dayside ionosphere of Venus was studied for both time-dependent and steady-state conditions. The present model is more complete than the one-major-ion model of Shinagawa et al., although a comparison of the results indicated that the one-major-ion treatment was a fairly good approximation. Two new cases are presented, including steady-state conditions for the magnetized ionosphere, and the inclusion of ion loss due to horizontal transport in the magnetized region. The resulting calculated profiles of the magnetic field and the electron density agree much better with the observations at high altitudes than those without the ion loss terms, indicating the importance of the horizontal transport processes in the ionosphere of Venus at high altitudes.

Shinagawa, H.↗

Airborne lidar and radiometric observations of PBL- and low clouds

Boundary layer- and low altitude clouds over open ocean and continent areas have been studied during several field campaigns since mid-1990 using the French airborne backscatter lidar LEANDRE in conjunction with on-board IR and visible radiometers. LEANDRE is an automatic system, and a modification of the instrumental parameters, when airborne, is computer controlled through an operator keyboard. The vertical range squared lidar signals and instrument status are displayed in real time on two dedicated monitors. The lidar is used either down- or up-looking while the aircraft is flying above or below clouds. A switching of the viewing configuration takes about a minute. The lidar measurements provide a high resolution description of cloud morphology and holes in cloud layers. The flights were conducted during various meteorological conditions on single or multilayer stratocumulus and cumulus decks. Analysis on a single shot basis of cloud top (or bottom) altitude and a plot of the corresponding histogram allows one to determine a probability density function (PDF). The preliminary results show the PDFs for cloud top are not Gaussian and symmetric about the mean value. The skewness varies with atmospheric conditions. An example of results recorded over the Atlantic ocean near Biarritz is displayed, showing: (1) the range squared lidar signals as a function of time (here 100 s corresponds to about 8 km, 60 shots are averaged on horizontal); the Planetary Boundary Layer (PBL) - up to 600 m - is observed at the beginning of the leg as well as on surface returns, giving an indication of the porosity; (2) the cloud top altitude variation between 2.4 to 2.8 km during the 150 to 320 s section; and (3) the corresponding PDF. Similar results are obtained on stratocumulus over land. Single shot measurements can be used also to determine an optical porosity at a small scale as well as a fractional cloudiness at a larger scale. A comparison of cloud top altitude retrieved from lidar and narrowbeam IR radiometer is conducted to study the scale integration problem. A good agreement within less than 100 m relies on spatial uniformity and an optically thick layer. In the presence of holes, a discrepancy is observed. This is illustrated in figure 2, displaying as a function of time (1) the lidar signals; (2) the target temperature (either clouds or sea surface) retreived from a narrowbeam IR radiometer, 17 C is the sea surface temperature on that day; and (3) the visible flux, linked to cloud albedo, measured by a pyranometer. In preparation of ASTEX, down- and up-looking measurements where conducted on stratocumulus clouds over the Atlantic Ocean near Quimper in Brittany. Depending on the flight pattern orientation with respect to the wind, the top and bottom cloud morphologies are different. Preliminary results are given on cloud morphology, cloud top PDFs, optical porosity, fractional cloudiness, and comparison of lidar and radiometric measurements.

Flamant, P. H.↗

Kasei Vallis of Mars: Dating the Interplay of Tectonics and Geomorphology

Crater density age dates on more than 250 small geomorphic surfaces in the Kasei Region of Mars show clusterings indicative of times of peak geomorphic and tectonic activity. Kasei Vallis is part of a 300 km wide channel system breaching a N-S trending ancient basement high (+50,000 crater age) separating the Chryse Basin from the Tharsis Volcanic Province of Mars. The basement high was covered by a least 3 groups of probable volcanic deposits. Major regional fracturing took place at age 4,000 to 5,000 and was immediately followed by deposition of regional volcanics of the Fesenkov Plains (age 3,000 to 4,200). Younger clusterings of dates in the 900 to 1,500 and 500 to 700 range represent only minor modification of the basic tectonic geomorphic landform. The data suggest that Kasei gap is a structurally controlled breach of a buried ridge by a rather brief episode of fluvial activity.

Wise, D. U.↗

Closed cycle MHD power generation experiments in the NASA Lewis Facility

Many modifications were made in the MHD facility. These include a redesign of the MHD duct interior, addition of mixing bars, increased electrical isolation of all the high temperature components from each other and from ground, and experimentation with various cesium seed vaporization and injection techniques. With the exception of the cesium system which needs further improvement the above modifications were quite successful and resulted in improvements in generator performance. The facility was run for a total of 400 hours in the past year, with 70 hours of this operation at temperatures of 2000 K or more with hot generator walls. With the exception of replacing one cracked brick in the MHD channel no repairs were required in the high temperature loop components for the duration of these tests. Uniform Faraday and Hall voltage profiles were obtained and the Faraday open circuit voltage varied from 90 to 100 percent of the ideal uBh. The magnitudes of the measured parameters are: Faraday open circuit voltage approximately 70 V, total Faraday current approximately 20 A, Hall voltage approximately 250 V, power output approximately 300 W, and power density .036 W/cu cm.

Sovie, R. J.↗

An Experimental System for Strategic Flight Path Management in Advanced Air Mobility

In the concept envisioned for Urban Air Mobility (UAM) operations, fleets of electric vertical takeoff and landing (eVTOL) vehicles would operate between vertiports distributed within a densely populated area. These operations would be largely independent from the existing air traffic control system and would place the responsibility for flight planning and aircraft separation on fleet operators. The fourth major level on the UAM Maturity Level scale, UML-4, relies on “collaborative and responsible” automation to enable operations in non-visual conditions with medium traffic density (hundreds of aircraft in one metropolitan region) and medium complexity. This level of service places many requirements on automation systems to assist the operators of these aircraft. NASA has developed the Autonomous Operations Planner (AOP), a reference prototype Flight Path Management automation system, and has modified AOP to support research of anticipated UML-4 operations. AOP creates a four-dimensional flight plan conforming to the constraints of these operations, evaluates and modifies the flight plan during flight as conditions and constraints evolve, and coordinates the flight plan with other airspace users and with service providers. This version of AOP has been integrated into the Sikorsky Autonomy Research Aircraft and used in a flight test activity. In this paper we discuss anticipated characteristics of UAM operations, modifications that were made to AOP to adapt to that environment or to support the flight test, and observations of software and aircraft performance during the flight test. The aircraft achieved four-dimensional conformance with the flight plan and AOP provided adequate planning in almost all cases. We discuss improvements that could be made to AOP to address deficiencies that were observed.

Autonomous Operations Planner↗

An Experimental System for Strategic Flight Path Management in Advanced Air Mobility

In the concept envisioned for Urban Air Mobility (UAM) operations, fleets of electric vertical takeoff and landing (eVTOL) vehicles would operate between vertiports distributed within a densely populated area. These operations would be largely independent from the existing air traffic control system and would place the responsibility for flight planning and aircraft separation on fleet operators. The fourth major level on the UAM Maturity Level scale, UML-4, relies on “collaborative and responsible” automation to enable operations in non-visual conditions with medium traffic density (hundreds of aircraft in one metropolitan region) and medium complexity. This level of service places many requirements on automation systems to assist the operators of these aircraft. NASA has developed the Autonomous Operations Planner (AOP), a reference prototype Flight Path Management automation system, and has modified AOP to support research of anticipated UML-4 operations. AOP creates a four-dimensional flight plan conforming to the constraints of these operations, evaluates and modifies the flight plan during flight as conditions and constraints evolve, and coordinates the flight plan with other airspace users and with service providers. This version of AOP has been integrated into the Sikorsky Autonomy Research Aircraft and used in a flight test activity. In this paper we discuss anticipated characteristics of UAM operations, modifications that were made to AOP to adapt to that environment or to support the flight test, and observations of software and aircraft performance during the flight test. The aircraft achieved four-dimensional conformance with the flight plan and AOP provided adequate planning in almost all cases. We discuss improvements that could be made to AOP to address deficiencies that were observed.

Autonomous Operations Planner↗

Improvements in the perturbation simulations of the global reference atmospheric model

The Global Reference Atmospheric Model (GRAM) program includes the capability for simulating pseudo-random perturbations in density, temperature, pressure, or wind components along a simulated reentry trajectory or other path through the atmosphere. Some concerns were expressed by GRAM users, however, that the mean-square perturbation gradients may be too large for small values of the vertical separation Delta z. The present GRAM perturbation simulations, based on a one-step autoregressive model, yield a power spectrum versus wavenumber k which is proportional to k sup -2 at high wavenumbers. This feature also produces mean-square perturbation differences which are directly proportional to Delta z, and mean-square perturbation gradients which are inversely proportional to Delta z. Thus, root-mean-square gradients, (Delta f/Delta z) sub rms, increase with decreasing Delta a as Delta z sup -1/2. A simple modification to GRAM is suggested which overcomes this problem, i.e., which produce root-mean-square gradient that remain bound as Delta z approaches zero. Possible applications of more sophisticated simulation approaches, e.g., second order autoregressive models, or fractal model techniques, were also explored briefly but found to yield improvements which appear too small to justify their considerable added complexity for use in the GRAM programs.

Justus, C. G.↗

Demonstration of Mars Crosslink Occultation Measurements for Future Small Spacecraft Constellations

Future planetary atmospheric profiling via radio occultations (RO) benefits significantly from increased received signal-to-noise ratio as well as geometrical coverage of links between two or more spacecraft in orbit around a target planet. These can be small spacecraft, possibly dedicated to quality radio-metrics with precision clock reference. Motivated by long-term research investigating the optimum SNR for the science performance, the optimum reference clock stability, the number and combination of wavelengths driving the number of transmitters and antennas, the optimum number and orbital spacecraft configuration, and design modifications to existing radio communication systems that would allow these new science objectives on future missions, three crosslink occultation experiments have been acquired between the Mars Odyssey and Mars Reconnaissance Orbiter spacecraft to probe the Martian atmosphere. While crosslink occultations between Earth orbiting satellites have long been used to profile the Earth’s atmosphere, this represents the first demonstration of crosslink occultation measurements at another planet. These measurements leverage the proximity link telecommunication payloads on each orbiter, which were designed to provide relay communication and navigation services to Mars landers and rovers. Analysis of the observed Doppler shift on each crosslink measurement reveals a clear signature of the Martian atmosphere, primarily the ionosphere. Inversion of the observed Doppler data yields vertical profiles of the Martian refractivity and electron density. The electron density profiles show the presence of two layers with peak densities and peak heights that are consistent with empirical models.

Asmar, Sami↗

Computer Aided Modeling to Determine the Effectiveness of Resistive Exercises as Countermeasures for Bone Mineral Density Loss

Due to the loss of gravitational loading, astronauts have a tendency to lose bone mineral density in their lumbar spine and lower extremities on orbit. NASA requires astronauts to perform exercises during space flight to help reduce the amount of demineralization. To test these exercises on earth, 17 week bed rest studies are conducted that consist of specific diet and exercise regimes. Developing a finite element model of these exercises will help to quantify the stress distribution imposed by of each of these exercises. To help develop this model, MRI images are acquired from individuals participating in the bed rest studies. The MRIs can be used to create a subject specific model of each individual for testing. The MRIs are processed in the Magnetic Resonance Imaging Data Transfer System program to develop a three-dimensional finite element model of the femur for evaluation. Modifications were made to the MRIDTS that simplified the model creation process. These modifications made it possible to construct two separate models of different portions of a bone simultaneously and then later connect them manually. This helped alleviate the warping problem associated with the drastic changes in geometry found in some body parts, such as the joints. The code was also modified to incorporate material properties of various bone components into the model. Interior meshing was also incorporated into the program to allow for both the cortical shell and the entire bone to be modeled. A prototype model of the right femur of an adult female is being constructed and tested to determine the feasibility of finite element analysis as a tool for evaluating exercise effectiveness. The model is being run through the ANSYS finite element program on the Alabama Super Computer Network. After the model is validated, models of bedrest subjects can be generated to investigate exercise countermeasures.

Murphy, Benjamin M.↗

Material Interactions with Solar Wind Ion Environments

Solar wind composition is dominated by hydrogen (approx.96%) and helium (approx.3 to 4%) with a minor fraction (less than or equal to 1%) of heavy ions. Hydrogen (helium) ions impact spacecraft surfaces with energies from 0.5 to 5 keV (1.8 to 21 keV) due to variations in solar wind velocity from 300 km/s to 1000 km/sec with extremes of a few 10 s keV during periods of extremely high solar wind velocity exceeding 1000 km/sec. Mean impact energies are typically on the order of approximately 1 keV and 4 keV for hydrogen ions and helium ions, respectively. These energies are typically of the peak of the energy dependent light ion sputter yields for hydrogen and helium on many metals. In addition, light ions with kilovolt energies have been shown to produce blister (or exfoliation) damage to metal surfaces due to formation of high pressure gas bubbles within the materials when exposed to ion fluences on the order of 10(exp 16 to (10(exp 17 ions/sq cm. A number of spacecraft designs for current and future missions include gossamer polymer structures with thin metallic reflection coatings to shield instruments from the Sun or solar sail propulsion systems for use in a variety of locations in the inner solar system from 0.5 to 1 AU. In addition, there is interest in designing spacecraft for solar physics missions requiring operations as close to the Sun as 0.16 to 0.2 AU. Integrity of the metallic coatings is critical in many of these applications since degradation will result in modification of material thermal properties or exposure of polymers to solar UV photons which can compromise mission requirements. This paper will evaluate the relative contributions of sputtering and blister formation to material degradation in solar wind environments over a range of radial distances from the Sun to demonstrate where solar wind environments become important for materials selection. We will first review the physics and results from laboratory measurements of light ion sputtering, blistering, and exfoliation of metallic surfaces to establish the order of magnitude ion fluence required for significant surface damage. Solar wind ion fluence environments will then be evaluated due to variations in solar wind conditions as a function of solar cycle for varying distances from the Sun using models for radial variations in solar wind ion number density, temperature, and velocity to determine where sputtering and blistering is most likely to be an issue. Finally, ion fluence statistics for varying radial distances from the Sun will be shown to establish the mission duration and radial distances from the Sun where missions will encounter sufficient ion fluence to exhibit damage to metallic surfaces.

Minow, Joseph I.↗

Study of non-equilibrium transport phenomena

Nonequilibrium phenomena due to real gas effects are very important features of low density hypersonic flows. The shock shape and emitted nonequilibrium radiation are identified as the bulk flow behavior parameters which are very sensitive to the nonequilibrium phenomena. These parameters can be measured in shock tubes, shock tunnels, and ballistic ranges and used to test the accuracy of computational fluid dynamic (CFD) codes. Since the CDF codes, by necessity, are based on multi-temperature models, it is also desirable to measure various temperatures, most importantly, the vibrational temperature. The CFD codes would require high temperature rate constants, which are not available at present. Experiments conducted at the NASA Electric Arc-driven Shock Tube (EAST) facility reveal that radiation from steel contaminants overwhelm the radiation from the test gas. For the measurement of radiation and the chemical parameters, further investigation and then appropriate modifications of the EAST facility are required.

Sharma, Surendra P.↗

High vacuum facility for hydrazine thruster testing

An ongoing modification is described of a large vacuum chamber to accommodate the ignition of an arcjet hydrazine thruster while maintaining a vacuum level of 1 x 10(exp -5) torr or less. The vacuum facility consists of a 20 ft stainless steel vacuum tank with an internal LN2 shroud, four 35 in. cryopumps and an 8 in. turbopump. To maintain a vacuum level of 1 x 10(exp -5) torr or less, 900 sq ft of liquid helium (LHe) shroud surface was installed to maintain the vacuum level and pumping requirements. A vacuum level of 1 x 10(exp -5) torr or less will allow the hydrazine thrust to exit the thruster nozzle and radiate into a space type environment so that the plume flow field can be analyzed and compared to the analytical model density distribution profile. Some other arcjet thruster characteristics measured are the electromagnetic interference (EMI) and exhaust contamination. This data is used to evaluate if the arcjet thruster with its high specific impulse in comparison to current chemical propulsion thruster can be used for the next generation of communication satellites.

Neary, Patrick F.↗

Computational Study of Field Initiated Surface Reactions for Synthesis of Diamond and Silicon

This project involves using quantum chemistry to simulate surface chemical reactions in the presence of an electric field for nanofabrication of diamond and silicon. A field delivered by a scanning tunneling microscope (STM) to a nanometer scale region of a surface affects chemical reaction potential energy surfaces (PES) to direct atomic scale surface modification to fabricate sub-nanometer structures. Our original hypothesis is that the applied voltage polarizes the charge distribution of the valence electrons and that these distorted molecular orbitals can be manipulated with the STM so as to change the relative stabilities of the electronic configurations over the reaction coordinates and thus the topology of the PES and reaction kinetics. Our objective is to investigate the effect of applied bias on surface reactions and the extent to which STM delivered fields can be used to direct surface chemical reactions on an atomic scale on diamond and silicon. To analyze the fundamentals of field induced chemistry and to investigate the application of this technique for the fabrication of nanostructures, we have employed methods capable of accurately describing molecular electronic structure. The methods we employ are density functional theory (DFT) quantum chemical (QC) methods. To determine the effect of applied bias on surface reactions we have calculated the QC PESs in various applied external fields for various reaction steps for depositing or etching diamond and silicon. We have chosen reactions which are thought to play a role in etching and the chemical vapor deposition growth of Si and diamond. The PESs of the elementary reaction steps involved are then calculated under the applied fields, which we vary in magnitude and configuration. We pay special attention to the change in the reaction barriers, and transition state locations, and search for low energy reaction channels which were inaccessible without the applied bias.

Musgrave, Charles Bruce↗

High Resolution Triple Axis X-Ray Diffraction Analysis of II-VI Semiconductor Crystals

The objective of this research program is to develop methods of structural analysis based on high resolution triple axis X-ray diffractometry (HRTXD) and to carry out detailed studies of defect distributions in crystals grown in both microgravity and ground-based environments. HRTXD represents a modification of the widely used double axis X-ray rocking curve method for the characterization of grown-in defects in nearly perfect crystals. In a double axis rocking curve experiment, the sample is illuminated by a monochromatic X-ray beam and the diffracted intensity is recorded by a fixed, wide-open detector. The intensity diffracted by the sample is then monitored as the sample is rotated through the Bragg reflection condition. The breadth of the peak, which is often reported as the full angular width at half the maximum intensity (FWHM), is used as an indicator of the amount of defects in the sample. This work has shown that high resolution triple axis X-ray diffraction is an effective tool for characterizing the defect structure in semiconductor crystals, particularly at high defect densities. Additionally, the technique is complimentary to X-ray topography for defect characterization in crystals.

Volz, H. M.↗

Hydrodynamic Instability and Thermal Coupling in a Dynamic Model of Liquid-Propellant Combustion

For liquid-propellant combustion, the Landau/Levich hydrodynamic models have been combined and extended to account for a dynamic dependence of the burning rate on the local pressure and temperature fields. Analysis of these extended models is greatly facilitated by exploiting the realistic smallness of the gas-to-liquid density ratio rho. Neglecting thermal coupling effects, an asymptotic expression was then derived for the cellular stability boundary A(sub p)(k) where A(sub p) is the pressure sensitivity of the burning rate and k is the disturbance wavenumber. The results explicitly indicate the stabilizing effects of gravity on long-wave disturbances, and those of viscosity and surface tension on short-wave perturbations, and the instability associated with intermediate wavenumbers for critical negative values of A(sub p). In the limit of weak gravity, hydrodynamic instability in liquid-propellant combustion becomes a long-wave, instability phenomenon, whereas at normal gravity, this instability is first manifested through O(1) wavenumbers. In addition, surface tension and viscosity (both liquid and gas) each produce comparable effects in the large-wavenumber regime, thereby providing important modifications to the previous analyses in which one or more of these effects was neglected. For A(sub p)= O, the Landau/Levich results are recovered in appropriate limiting cases, although this typically corresponds to a hydrodynamically unstable parameter regime for p << 1. In addition to the classical cellular form of hydrodynamic stability, there exists a pulsating form corresponding to the loss of stability of steady, planar burning to time-dependent perturbations. This occurs for negative values of the parameter A(sub p), and is thus absent from the original Landau/Levich models. In the extended model, however, there exists a stable band of negative pressure sensitivities bounded above by the Landau type of instability, and below by this pulsating form of hydrodynamic instability. Indeed, nonsteady modes of combustion have been observed at low pressures in hydroxylammonium nitrate (HAN)-based liquid propellants, which often exhibit negative pressure sensitivities. While nonsteady combustion may correspond to secondary and higher-order bifurcations above the cellular boundary, it may also be a manifestation of this pulsating type of hydrodynamic instability. In the present work, a nonzero temperature sensitivity is incorporated into our previous asymptotic analyses. This entails a coupling of the energy equation to the previous purely hydrodynamic problem, and leads to a significant modification of the pulsating boundary such that, for sufficiently large values of the temperature-sensitivity parameter, liquid-propellant combustion can become intrinsically unstable to this alternative form of hydrodynamic instability. For simplicity, further attention is confined here to the inviscid version of the problem since, despite the fact that viscous and surface-tension effects are comparable, the qualitative nature of the cellular boundary remains preserved in the zero-viscosity limit, as does the existence of the pulsating boundary. The mathematical model adopts the classical assumption that there is no distributed reaction in either the liquid or gas phases, but now the reaction sheet, representing either a pyrolysis reaction or an exothermic decomposition at the liquid/gas interface, is assumed to depend on local conditions there.

Margolis, S. B.↗

Modification of Mach 6 Freestream Flow by Pitot Probe Bow Shock as Measured by Line FLDI and SAFS

Measurements of the flow in front of flat-face Pitot probes in NASA Langley ResearchCenter’s 20-Inch Mach 6 Air Tunnel were made using a 35-point line focused laser differentialinterferometer (FLDI), a two-point FLDI, and a self-aligned focusing schlieren system. Thesesystems were used to make measurements of flow unsteadiness between the bow shock andPitot probe face, while also probing the freestream region just upstream of the bow shock.Shock standoff distances from the non-intrusive measurements match closely with analyticalsolutions using the probe face diameter. Acoustic wave angles and velocities in the freestreamcompare favorably with DNS simulations of a Mach 6 freestream flow, and wave velocities in thepost-shock region slow as they approach the probe face. Peak frequencies in the post-shockregion as measured by FLDI compare favorably to the estimates computed from the shockstandoff distance for most probe diameters. With increasing unit Reynolds number, fluctuationintensities in the post-shock region were attenuated below the intensities of the freestream, witha spectral, spatial, and probe size dependence. Pitot probe power spectral density estimatesare presented from the test campaign, which can be compared with FLDI and SAFS data toevaluate the applicability of the transfer functions derived from computational simulations.

Joshua M Weisberger↗

Dynamics of Sheared Granular Materials

This work focuses on the properties of sheared granular materials near the jamming transition. The project currently involves two aspects. The first of these is an experiment that is a prototype for a planned ISS (International Space Station) flight. The second is discrete element simulations (DES) that can give insight into the behavior one might expect in a reduced-g environment. The experimental arrangement consists of an annular channel that contains the granular material. One surface, say the upper surface, rotates so as to shear the material contained in the annulus. The lower surface controls the mean density/mean stress on the sample through an actuator or other control system. A novel feature under development is the ability to 'thermalize' the layer, i.e. create a larger amount of random motion in the material, by using the actuating system to provide vibrations as well control the mean volume of the annulus. The stress states of the system are determined by transducers on the non-rotating wall. These measure both shear and normal components of the stress on different size scales. Here, the idea is to characterize the system as the density varies through values spanning dense almost solid to relatively mobile granular states. This transition regime encompasses the regime usually thought of as the glass transition, and/or the jamming transition. Motivation for this experiment springs from ideas of a granular glass transition, a related jamming transition, and from recent experiments. In particular, we note recent experiments carried out by our group to characterize this type of transition and also to demonstrate/ characterize fluctuations in slowly sheared systems. These experiments give key insights into what one might expect in near-zero g. In particular, they show that the compressibility of granular systems diverges at a transition or critical point. It is this divergence, coupled to gravity, that makes it extremely difficult if not impossible to characterize the transition region in an earth-bound experiment. In the DE modeling, we analyze dynamics of a sheared granular system in Couette geometry in two (2D) and three (3D) space dimensions. Here, the idea is to both better understand what we might encounter in a reduced-g environment, and at a deeper level to deduce the physics of sheared systems in a density regime that has not been addressed by past experiments or simulations. One aspect of the simulations addresses sheared 2D system in zero-g environment. For low volume fractions, the expected dynamics of this type of system is relatively well understood. However, as the volume fraction is increased, the system undergoes a phase transition, as explained above. The DES concentrate on the evolution of the system as the solid volume fraction is slowly increased, and in particular on the behavior of very dense systems. For these configurations, the simulations show that polydispersity of the sheared particles is a crucial factor that determines the system response. Figures 1 and 2 below, that present the total force on each grain, show that even relatively small (10 %) nonuniformity of the size of the grains (expected in typical experiments) may lead to significant modifications of the system properties, such as velocity profiles, temperature, force propagation, and formation shear bands. The simulations are extended in a few other directions, in order to provide additional insight to the experimental system analyzed above. In one direction, both gravity, and driving due to vibrations are included. These simulations allow for predictions on the driving regime that is required in the experiments in order to analyze the jamming transition. Furthermore, direct comparison of experiments and DES will allow for verification of the modeling assumptions. We have also extended our modeling efforts to 3D. The (preliminary) results of these simulations of an annular system in zero-g environment will conclude the presentation.

Kondic, Lou↗