Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “point contacts”

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 919 records · Page 51

The clumpy circumstellar medium around young supernova remnants

Each of the youngest supernova remnants known in the Milky Way, Cas A, Kepler's SNR, and Tycho's SNR, shows a different morphological structure caused by different conditions in the progenitor stars and their surroundings. In all three cases, however, the observed shells have a thickness of about 1/4 the radius, a sharp outer edge, and significant brightness irregularities. These features require that the circumstellar medium be highly clumped. To investigate the phenomenon, models of the expansion have been constructed using a one-dimensional spherical hydrodynamic code. As a supernova shock moves down the external density gradient of the star, material behind the shock begins to go into free expansion. Then as surrounding material is encountered a reverse shock moving back into the ejectum will be formed. Until the expansion has swept up about eight times the ejected mass when the situation can be considered as a point explosion in its surroundings, the dynamics are controlled by conditions between the shocks. The region is also where the synchrotron radio emission from relativistic electrons trapped in magnetic fields arises. Initial particles and fields are accelerated and amplified by eddy motion at the interface between the ejected and swept-up material and at the boundaries of clumps. Polarimetry shows that these SNR have a net radial orientation of their magnetic fields apparently from stretching by Rayleigh-Taylor instabilities at the contact surfaces. Without clumps the observed shell is much too narrow and steep on the inside.

Dickel, J. R.↗

A preliminary analysis of the Mariner 10 color ratio map of Mercury

A preliminary geological analysis of the Mariner 10 orange/UV color ratio map of Mercury is given, assuming a basaltic crust. Certain errors in the map are pointed out. The relationship between color and terrain are distinctly non-lunar. Rays and ejecta are bluer than average on Mercury, whereas they are redder on the Moon. This fact, along with the lack of the ferrous band in Mercury's spectral reflectance and smaller albedo contrasts, implies that the crust is low in Fe and Ti. There is no correlation between color boundaries and the smooth plains on Mercury, in contrast with the strong correlation between color and maria-highlands contacts on the Moon. The smooth plains are not Mercurian analogs of lunar maria, and a lunar-type of second wave melting did not occur. Ambiguous correlations between color and topography indicate that older, redder materials underlie younger, bluer rocks in many places on the planet, implying that the last stages of volcanism involved low-Fe lavas covering higher-Fe rocks. There is some evidence of late Fe-rich pyroclastic activity.

Rava, Barry↗

Nozomi Cis-Lunar Phase Orbit Determination

Japan's Institute of Space and Astronautical Science (ISAS) launched Nozomi, its first mission to the planet Mars using the newly developed M-V launch vehicle on July 3, 1998. Scientific objectives of the mission are to study the structure and dynamics of the Martian upper atmosphere and its interaction with the solar wind. Nozomi is a cooperative mission between ISAS and the National Aeronautics and Space Administration (NASA). The NASA contribution includes navigation and tracking services provided by the Jet Propulsion Laboratory (JPL). The spacecraft also serves as an engineering demonstration of basic technology for planetary exploration. One of the new technologies was a unique trajectory, developed by ISAS, which used solar gravitational perturbations at the weak stability boundary as an aid to achieve an Earth-Mars transfer orbit. This trajectory saves approximately 120 m/s of Delta V compared to direct hyperbolic insertion and is considered an enabling technology for the mission. Nozomi was the first spacecraft to employ this trajectory and provided on-orbit validation of the technique. The trajectory was achieved by initially placing the spacecraft in a highly elliptical cis-lunar phasing orbit. Six maneuvers were performed during this period to correct injection errors and target an outbound lunar swingby in September 1998. The gravity assist from the lunar swingby raised apogee to the vicinity of the weak stability boundary. After three more targeting maneuvers, Nozomi performed an inbound lunar swingby followed immediately by a powered Earth swingby in late December 1998. A 420 m/s Trans Mars Insertion (TMI) burn at the final Earth periapsis was intended to place the spacecraft on a heliocentric trajectory leading to Mars orbit insertion in October 1999. Orbit determination for Nozomi is performed in parallel by both ISAS and the Multi-Mission Navigation (MMNAV) group at JPL. This was an advantage for the mission because each group would generate solutions based on data collected from their respective tracking networks. Spacecraft events, such as sequence uplinks and maneuvers, were generally scheduled during passes at the Usuda tracking station in Japan. As a result, maneuver design and reconstruction was derived from MMNAV solutions based on JPL tracking data obtained immediately prior to or following maneuvers. Data was also exchanged between ISAS and MMNAV so orbit determination could be performed on joint data sets in support of critical targeting late in the cis-lunar phase. In this paper, information regarding the MMNAV orbit determination effort for the first six months of the mission is presented. The spacecraft trajectory is characterized first, followed by a discussion of the orbit determination estimation procedure and models. Results from selected orbit solutions are presented and compared against reconstructed trajectories. One area of emphasis in this paper is orbit determination in the vicinity of the weak stability boundary. Precise navigation was necessary to target the second lunar swingby and the powered Earth swingby. Delivery accuracy of 150 m was required for these critical encounters, but a number of factors contributed to the general degradation of orbit determination accuracy. This included the fact that the spacecraft was at apogee, at a range of 1.7 million km and moving at less than I km/sec perpendicular to the line of sight. Nozomi was also close to zero degrees declination where there are known limitations on orbit determination performance. Finally, S-band tracking data was acquired through the Nozomi backup low gain antenna. This antenna is offset from the axis of this spin stabilized spacecraft and superimposed large signatures in the Doppler and range data. These difficulties were overcome by combining long data arcs, spanning several maneuvers, with a high fidelity solar pressure model. The model included a physically accurate representation of the spacecraft structure and a high time resolution orientation model. Observation modeling included the removal of the spin induced Doppler bias, spin signature and per pass correction of range calibration errors applied for data leading up to critical events. As a result, all orbit determination goals were met. A second area of emphasis in this paper is the JPL tracking and orbit determination effort in support of the TMI maneuver. TMI occurred out of contact with ground stations and the JPL Goldstone tracking complex had the first pass following the bum. As a result, MMNAV had the responsibility to make a rapid assessment of the maneuver performance. MMNAV made the determination that a 100 m/s under bum had occurred and promptly informed ISAS via voice lines. ISAS immediately began preparations for a correction maneuver (TMIc), which had to be performed during the next Usuda pass. The near real time assessment by MMNAV provided accurate antenna frequency and pointing updates for the spacecraft acquisition at Usuda and the close coordination between the two agencies enabled the design and successful execution of the TMc maneuver. Propellant consumption during the correction burn dictated that the mission be redesigned. ISAS developed a new plan which adds 3 full solar orbits, two Earth swingbys and one lunar swingby with arrival at Mars in January 2004. The final Mars orbit will still enable the mission to achieve all of its science objectives.

Ryne, Mark↗

Diffusion in intact secondary cell wall models of plants at different equilibrium moisture content

Secondary plant cell walls are composed of carbohydrate and lignin polymers, and collectively represent a significant renewable resource. Leveraging these resources depends in part on a mechanistic understanding for diffusive processes within plant cell walls. Common wood protection treatments and biomass conversion processes to create biorefinery feedstocks feature ion or solvent diffusion within the cell wall. X-ray fluorescence microscopy experiments have determined that ionic diffusion rates are dependent on cell wall hydration as well as the ionic species through non-linear relationships. In this work, we use classical molecular dynamics simulations to map the diffusion behavior of different plant cell wall components (cellulose, hemicellulose, lignin), ions (Na + , K + , Cu2 + , Cl - ) and water within a model for an intact plant cell wall at various hydration states (3-30 wt% water). From these simulations, we analyze the contacts between different plant cell wall components with each other and their interaction with the ions. Generally, diffusion increases with increasing hydration, with lignin and hemicellulose components increasing diffusion by an order of magnitude over the tested hydration range. Ion diffusion depends on charge. Positively charged cations preferentially interact with hemicellulose components, which include negatively charged carboxylates. As a result, positive ions diffuse more slowly than negatively charged ions. Measured diffusion coefficients are largely observed to best fit piecewise linear trends, with an inflection point between 10 and 15% hydration. These observations shed light onto the molecular mechanisms for diffusive processes within secondary plant cell walls at atomic resolution.

09 BIOMASS FUELS↗

Development and evaluation of die materials for use in the growth of silicon ribbons by the inverted ribbon growth process, task 2, LSSA project

Silicon sessile drop experiments were performed on a variety of commercially available refractory carbides, nitrides, oxides, and borides to examine the potential of these materials for applications involving either direct contact with molten silicon or as substrates for CVD coatings in the fabrication of dies and crucibles for containing molten silicon. Simultaneous experiments were also conducted with CVD layers of SiC, Si3N4, and SiOxNy. Silicon nitride layers, deposited with NH3:SiH4 ratios ranging from 100:1 down to 5:1, were examined in sessile drop experiments to determine if the layers are degraded as a result of using lower reagent ratios. Preliminary experiments were undertaken on the stability of CVD Si3N4 near the melting point of silicon. Silicon ribbon segments were grown from vitreous carbon dies which had been coated with CVD Si3N4. Depending upon the purity of the die materials, ribbon resistivity values up to 40 Omega cm were obtained.

Duffy, M. T.↗

Critical point wetting drop tower experiment

The 100 m Drop Tower at NASA-Marshall was used to provide the step change in acceleration from 1.0 to 0.0005 g. An inter-fluid meniscus oscillates vertically within a cylindrical container when suddenly released from earth's gravity and taken into a microgravity environment. Oscillations damp out from energy dissipative mechanisms such as viscosity and interfacial friction. Damping of the oscillations by the later mechanism is affected by the nature of the interfacial junction between the fluid-fluid interface and the container wall. In earlier stages of the project, the meniscus shape which developed during microgravity conditions was applied to evaluations of wetting phenomena near the critical temperature. Variations in equilibrium contact angle against the container wall were expected to occur under critical wetting conditions. However, it became apparent that the meaningful phenomenon was the damping of interfacial oscillations. This latter concept makes up the bulk of this report. Perfluoromethyl cyclohexane and isopropanol in glass were the materials used for the experiment. The wetting condition of the fluids against the wall changes at the critical wetting transition temperature. This change in wetting causes a change in the damping characteristics of the interfacial excursions during oscillation and no measurable change in contact angle. The effect of contact line friction measured above and below the wetting transition temperature was to increase the period of vertical oscillation for the vapor-liquid interface when below the wetting transition temperature.

Kaukler, William F.↗

System for sealing flanged connections

A flange sealing system can include a loading member layer adapted to receive a sealing force thereon, at least one intermediate member receiving the sealing force from the loading member and adapted to distribute the load force to branch members in subsequent layers, and at least one contact member adapted to apply the sealing force evenly around a flange and press a surface of the flange and against at least one of a second flange surface or a cap to seal an opening associated with the flange via at least one of connection to the second flange surface or the cap. The system can be provided in the form of a first pipe receiving section and a second pipe receiving section for securing two flanges together with associated pipes with pivot points distributed about the first and second pipe receiving section and distributing even pressure around the flanges.

Baffes, Curtis M.↗

Multiple RGB ortho-mosaics and digital surface models in 2017 and 2018 across the Lower Montane site in the East River Watershed, Colorado

Aerial imagery was collected at the Lower Montane site (Pumphouse) in the East River Watershed, Colorado during the spring, summer, and fall seasons of 2017 and 2018 to improve the understanding of seasonal vegetation dynamics and their drivers. The datasets include Red-Green-Blue (RGB) ortho-mosaics and digital surface models (DSMs) inferred from the Unoccupied Aerial System (UAS) acquired aerial RGB imagery for June 3, June 19, July 7, and August 14, 2017, and for March 14, April 26, June 1, June 18, July 6, and August 7, 2018. Real-Time Kinematic Global Positioning System (RTK-GPS) surveyed Ground control points (GCPs) were used to increase the reconstruction accuracy. The reconstructed RGB mosaics and DSMs have been trimmed to cover a similar spatial domain. The accuracy of the RGB mosaics is considered high (~10 cm). DSM accuracy is highest (~10 cm) where sufficient GCPS are available, and more difficult to assess elsewhere (see reconstruction reports for uncertainty estimates). The dataset includes a total of 20 GeoTIFF (.tif) files, 10 PDF (.pdf) files, 3 data CSV (.csv) files, and 2 metadata CSV (.csv) files. Feel free to contact the authors with any questions or collaboration interests.This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗

Recent acceleration of Thwaites Glacier

The first velocity measurements for Thwaites Glacier were made by R. J. Allen in 1977. He compared features of Thwaites Glacier and Iceberg Tongue on aerial photography from 1947 and 1967 with 1972 Landsat images, and measured average annual displacements of 3.7 and 2.3 km/a. Using his photogrammetric experience and taking into consideration the lack of definable features and the poor control in the area, he estimated an average velocity of 2.0 to 2.9 km/a to be more accurate. In 1985, Lindstrom and Tyler also made velocity estimates for Thwaites Glacier. Using Landsat imagery from 1972 and 1983, their estimates of the velocities of 33 points ranged from 2.99 to 4.02 km/a, with an average of 3.6 km/a. The accuracy of their estimates is uncertain, however, because in the absence of fixed control points, they assumed that the velocities of icebergs in the fast ice were uniform. Using additional Landsat imagery in 1984 and 1990, accurate coregistration with the 1972 image was achieved based on fixed rock points. For the period 1972 to 1984, 25 points on the glacier surface ranged in average velocity from 2.47 to 2.76 km/a, with an overall average velocity of 2.62 +/- 0.02 km/a. For the period 1984 to 1990, 101 points ranged in velocity from 2.54 to 3.15 km/a, with an overall average of 2.84 km/a. During both time periods, the velocity pattern showed the same spatial relationship for three longitudinal paths. The 8-percent acceleration in a decade is significant. This recent acceleration may be associated with changes observed in this region since 1986. Fast ice melted and several icebergs calved from the base of the Iceberg Tongue and the terminus of Thwaites Glacier. However, as early as 1972, the Iceberg Tongue had very little contact with the glacier.

Ferrigno, J. G.↗

Theoretical and Experimental Investigation of the Stability of an Evaporating Constrained Vapor Bubble

The major accomplishments of the experimental portion of the research were documented in Ling Zheng's doctoral dissertation. Using Pentane, he obtained a considerable amount of data on the stability and heat transfer characteristics of an evaporating meniscus. The important points are that experimental equipment to obtain data on the stability and heat transfer characteristics of an evaporating meniscus were built and successfully operated. The data and subsequent analyses were accepted by the Journal of Heat Transfer for publication in 2004 [PU4]. The work was continued by a new graduate student using HFE-7000 [PU3] and then Pentane at lower heat fluxes. The Pentane results are being analyzed for publication. The experimental techniques are currently being used in our other NASA Grant. The oscillation of the contact line observed in the experiments involves evaporation (retraction part) and spreading. Since both processes occur with finite contact angles, it is important to derive a precise equation of the intermolecular forces (disjoining pressure) valid for non-zero contact angles. This theoretical derivation was accepted for publication by Journal of Fluid Mechanics [PU5]. The evaporation process near the contact line is complicated, and an idealized micro heat pipe has been proposed to help in elucidating the detailed evaporation process [manuscripts in preparation].

Wayner, P. C., Jr.↗

Mineralogy of SNC Meteorite EET79001 by Simultaneous Fitting of Moessbauer Backscatter Spectra

We have acquired M ssbauer spectra for SNC meteorite EET79001 with a MIMOS II backscatter M ssbauer spectrometer [1] similar to those now operating on Mars as part of the Mars Exploration Rover (MER) missions. We are working to compare the Fe mineralogical composition of martian meteorites with in-situ measurements on Mars. Our samples were hand picked from the >1 mm size fraction of saw fines on the basis of lithology, color, and grain size (Table 1). The chips were individually analyzed at approx.300K by placing them on a piece of plastic that was in turn supported by the contact ring of the instrument (oriented vertically). Tungsten foil was used to mask certain areas from analysis. As shown in Figure 1, a variety of spectra was obtained, each resulting from different relative contributions of the Fe-bearing minerals present in the sample. Because the nine samples are reasonably mixtures of the same Fe-bearing phases in variable proportions, the nine spectra were fit simultaneously (simfit) with a common model, adjusting parameters to a single minimum chi-squared convergence criterion [2]. The starting point for the fitting model and values of hyperfine parameters was the work of Solberg and Burns [3], who identified olivine, pyroxene, and ferrous glass as major, and ilmenite and a ferric phase as minor (<5%), Fe-bearing phases in EET79001.

Morris, Richard V.↗

Organic polymer compositions for removal of oxoanions from aqueous solutions

An organic polymer composition (OPC) comprising guanidinium-containing A units interconnected with B units of the formula *—R*)r, wherein each asterisk (*) in A units denotes a connection point with an asterisk in B units; R is a hydrocarbon linking group containing at least one carbon atom; r is an integer of 1, 2, or 3; and the composition necessarily includes an anionic species Xm− with a magnitude of charge m of at least 1, wherein the sum of negative charge provided by anionic species Xm− counterbalances the total positive charge provided by the A units. Also described herein is a method for removing one or more oxoanions from an aqueous source by (i) contacting the aqueous source with the above described OPC to result in absorption of the oxoanion into the OPC to produce an oxoanion-containing OPC; and (ii) removing the oxoanion-containing OPC from the aqueous source.

Jansone-Popova, Santa↗

Onset dynamics of air-water menisci on rock fracture surfaces

Information on menisci dynamics and equilibrium contact angles, θ e , is needed for modeling multiphase flow of fluids in geologic formations. The wettability of fracture faces is relevant to applications such as waterflooding in enhanced oil recovery and evaluation of caprock integrity for geologic carbon sequestration. We investigated the vertical capillary rise of air-water menisci on exposed fracture faces for a sandstone, a shale, and three granites. The experimental design approximated that employed in the Wilhelmy plate method. Replicate fracture faces were prepared by fracturing cylindrical rock cores using the Brazilian method and splitting the fractured cores apart along their fracture planes. Meniscus onset was imaged from a cross-sectional perspective using dynamic neutron radiography, and quantified with change point analysis. An existing model for meniscus onset fitted the experimental data very well (median R 2 = 0.96). Capillary rise was typified by t 2/3 scaling (where t is time), followed by a constant equilibrium height, z e . The average time taken for the menisci to reach z e was 0.39 s, with no differences between rock types. The menisci achieved a greater vertical extent on the fracture faces of the sedimentary rocks as compared to those of the granites. Apparent θ e values were calculated from the best fit estimates of z e using the Wilhelmy plate equation for a smooth and planar vertical surface. Despite the fractal roughness of some of the fracture faces, there was a significant linear relationship (R 2 = 0.80) between the apparent θ e values and intrinsic θ e values measured on polished surfaces of the same rocks using the sessile drop method. The regression equation indicated the apparent θ e values were consistently greater than the intrinsic values. This overestimation may be due to averaging of pixels at the air-water-solid interface, a required step in the image analysis. Further exploration of the range of applicability and measurement accuracy of this proposed new technique would be valuable.

58 GEOSCIENCES↗

Analysis Method for Non-Nominal First Acquisition

First this paper describes a method how the trajectory of the launcher can be modelled for the contingency analysis without having much information about the launch vehicle itself. From a dense sequence of state vectors a velocity profile is derived which is sufficiently accurate to enable the Flight Dynamics Team to integrate parts of the launcher trajectory on its own and to simulate contingency cases by modifying the velocity profile. Then the paper focuses on the thorough visibility analysis which has to follow the contingency case or burn performance simulations. In the ideal case it is possible to identify a ground station which is able to acquire the satellite independent from the burn performance. The correlations between the burn performance and the pointing at subsequent ground stations are derived with the aim of establishing simple guidelines which can be applied quickly and which significantly improve the chance of acquisition at subsequent ground stations. In the paper the method is applied to the Soyuz/Fregat launch with the MetOp satellite. Overall the paper shows that the launcher trajectory modelling with the simulation of contingency cases in connection with a ground station visibility analysis leads to a proper selection of ground stations and acquisition methods. In the MetOp case this ensured successful contact of all ground stations during the first hour after separation without having to rely on any early orbit determination result or state vector update.

Sieg, Detlef↗

Healable and conductive sulfur iodide for solid-state Li–S batteries

Solid-state Li-S batteries (SSLSBs) are made of low-cost and abundant materials free of supply chain concerns. Owing to their high theoretical energy densities, they are highly desirable for electric vehicles. However, the development of SSLSBs has been historically plagued by the insulating nature of sulfur and the poor interfacial contacts induced by its large volume change during cycling, impeding charge transfer among different solid components. We report an S 9.3 I molecular crystal with I 2 inserted in the crystalline sulfur structure, which shows a semiconductor-level electrical conductivity (approximately 5.9 × 10 -7 S cm -1 ) at 25 °C; an 11-order-of-magnitude increase over sulfur itself. Iodine introduces new states into the band gap of sulfur and promotes the formation of reactive polysulfides during electrochemical cycling. Further, the material features a low melting point of around 65 °C, which enables repairing of damaged interfaces due to cycling by periodical remelting of the cathode material. As a result, an Li-S 9.3 I battery demonstrates 400 stable cycles with a specific capacity retention of 87%. The design of this conductive, low-melting-point sulfur iodide material represents a substantial advancement in the chemistry of sulfur materials, and opens the door to the practical realization of SSLSBs.

25 ENERGY STORAGE↗

Prediction of hemiwicking dynamics in micropillar arrays

Dynamic hemiwicking behavior is observable in both nature and a wide range of industrial applications ranging from biomedical devices to thermal management. We present a semi-analytical modeling framework (without empirical fitting coefficients) to predict transient capillary-driven hemiwicking behavior of a liquid through a nano/microstructured surface, specifically a micropillar array. In our model framework, the liquid domain is discretized into micropillar unit cells to enable the time marching of the hemiwicking front. A simplified linear pressure drop is assumed along the hemiwicking length such that the local meniscus curvature, contact angle, and effective liquid height are determined at each time step in our transient model. This semi-analytical model is validated with experimental data from our own experiments and from published literature for different fluids. Our model predicts hemiwicking dynamics with <20% error over a broad range of micropillar geometries with height-to-pitch ratio ranging between ≈0.34 and 6.7 and diameter-to-pitch ratio in the range of ≈0.25–0.7 and without any fitting parameters. For lower diameter-to-pitch ratio data points related to sparse micropillar array arrangements, we suggest modifications to the semi-analytical model. This work sheds light on complex and dynamic solid–liquid–vapor interfacial interactions which could serve as a guide for the design of textured surfaces for wicking enhancement in multi-phase thermal and mass transport technologies and applications.

Mechanics↗

Hydrogen maser clocks in space for solid-Earth research and time-transfer applications: Experiment overview and evaluation of Russian miniature sapphire loaded cavity

The Observatoire Cantonal de Neuchatel (ON) is developing for ESTEC a compact H-maser for space use based upon a miniature sapphire loaded microwave cavity, a technique pioneered at VNIIFTRI. Various contacts between West-European parties, headed by ESA, and the Russian parties, headed by ESA, led to the proposal for flying two H-masers on Meteor 3M, a Russian meteorology satellite in low polar orbit. The experiment will include two masers, one provided by ON and the other by VNIIFTRI. T/F transfer and precise positioning will be performed by both a microwave link, using PRARE equipment, and an optical link, using LASSO-like equipment. The main objectives of the experiment are precise orbit determination and point positioning for geodetic/geophysical research, ultra-accurate time comparison and dissemination as well as in-orbit demonstration of operation and performance of H-masers. Within the scope of a preliminary space H-maser development phase performed for ESTEC at ON in preparation to the joint experiment, a Russian miniature sapphire loaded microwave cavity, on loan from VNIIFTRI, was evaluated in a full-size EFOS hydrogen maser built by ON. The experimental evaluation confirmed the theoretical expectation that with a hydrogen storage volume of only 0.65 liter an atomic quality factor of 1.5 x 10(exp 9) can be obtained for a -105 dBm output power. This represents a theoretical Allan deviation of 1.7 x 10(exp -15) averaged on a 1000 s time interval. From a full-size design to a compact one, therefore, the sacrifice in performance due to the reduction of the storage volume is very small.

Busca, G.↗

Aeolian Sand Transport in the Planetary Context: Respective Roles of Aerodynamic and Bed-Dilatancy Thresholds

The traditional view of aeolian sand transport generally estimates flux from the perspective of aerodynamic forces creating the airborne grain population, although it has been recognized that "reptation" causes a significant part of the total airborne flux; reptation involves both ballistic injection of grains into the air stream by the impact of saltating grains as well as the "nudging" of surface grains into a creeping motion. Whilst aerodynamic forces may initiate sand motion, it is proposed here that within a fully-matured grain cloud, flux is actually governed by two thresholds: an aerodynamic threshold, and a bed-dilatancy threshold. It is the latter which controls the reptation population, and its significance increases proportionally with transport energy. Because we only have experience with terrestrial sand transport, extrapolations of aeolian theory to Mars and Venus have adjusted only the aerodynamic factor, taking gravitational forces and atmospheric density as the prime variables in the aerodynamic equations, but neglecting reptation. The basis for our perspective on the importance of reptation and bed dilatancy is a set of experiments that were designed to simulate sand transport across the surface of a martian dune. Using a modified sporting crossbow in which a sand-impelling sabot replaced the bolt-firing mechanism, individual grains of sand were fired at loose sand targets with glancing angles typical of saltation impact; grains were projected at about 80 m/s to simulate velocities commensurate with those predicted for extreme martian aeolian conditions. The sabot impelling method permitted study of individual impacts without the masking effect of bed mobilization encountered in wind-tunnel studies. At these martian impact velocities, grains produced small craters formed by the ejection of several hundred grains from the bed. Unexpectedly, the craters were not elongated, despite glancing impact; the craters were very close to circular in planform. High-speed photography showed them to grow in both diameter and depth after the impactor had ricochetted from the crater site. The delayed response of the bed was "explosive" in nature, and created a miniature ejecta curtain spreading upward and outward for many centimeters for impact of 100-300 micron-diameter grains into similar material. Elastic energy deposited in the bed by the impacting grain creates a subsurface stress regime or "quasi-Boussinesq" compression field. Elastic recovery of the bed occurs by dilatancy; shear stresses suddenly convert the grains from closed to open packing, and grains are consequently able to eject themselves forcefully from the impact site. Random jostling of the grains causes radial homogenization of stress vectors and a resulting circular crater. There is a great temptation to draw parallels with cratering produced by meteorite impacts, but a rigorous search for common modelling ground between the two phenomena has not been conducted at this time. For every impact of an aerodynamically energized grain, there are several hundred grains ejected into the wind for the high-energy transport that might occur on Mars. Many of these grains will themselves become subject to the boundary layer's aerodynamic lift forces (their motion will not immediately die and add to the creep population), and these grains will become indistinguishable from those lifted entirely by aerodynamic forces. As each grain impacts the bed, it will eject even more grains into the flow. A cascading effect will take place, but because it must be finite in its growth, damping will occur as the number of grains set in motion causes mid-air collisions that prevent much of the impact energy from reaching the surface of the bed -thus creating a dynamic equilibrium in a high-density saltation cloud. It is apparent that for a given impact energy, the stress field permits a smaller volume of grains to convert to open packing as the size of the bed grains increases, or as the energy of the "percussive" grain decreases (by decrease in velocity or mass). Thus, the mass of the "repercussive" grain population that is ejected from the impact site becomes a function of the scale of the stress field in relation to the scale of the bed material (self-similarity being applicable if both bed size and energy are simultaneously adjusted). In other words, in a very high energy aeolian system where an aerodynamically raised grain can ballistically raise many more grains, the amount of material lifted into the wind becomes largely a function of a dilatancy threshold. If this threshold is exceeded, grains are repercussively injected into the saltation cloud. The "dilatancy threshold" may be defined in terms of the saltation percussive force required to convert the bed, through elastic response, from a closed to an open packing system. If open packing cannot be created, the grains cannot escape from the impact site, even though the elastic deformation and percussive force may be able to reorganize the grains with respect to one another. As the crossbow experiments showed, for an ever-increasing bed grain size, a point is reached when no material can be moved because the energy of the percussive grain is insufficient to dilate the relatively coarse bed. Although this seems to be stating the obvious -- that too little energy will not cause the bed to splash -- the consequences of exceeding the "splash threshold" by dilatancy are not so obvious for high-energy aeolian transport. It is noted that the force required to elastically dilate the bed has to overcome Coulombic grain attractions such as dipole-dipole coupling, dielectric, monopole, contact-induced dipole attractions, van der Waals forces, molecular monolayer capillary forces, as well as the mechanical interlocking frictional resistance of the grains. On Mars, it is predicted that the dilatancy threshold may be the prime control of grain flux. Additional information is contained in the original.

Marshall, J. R.↗