Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “forcing fields”

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 1,009 records · Page 56

Enhanced Fair-Weather Electric Fields Soon After Sunrise

The typical fair weather electric field at the ground is between -100 and -300 V/m. At the NASA Kennedy Space Center and US Air Force Cape Canaveral Air Station (KSC) the electric field at the ground sometimes reaches -400 to -1200 V/m within an hour or two after sunrise on days that otherwise seem to be fair weather. We refer to the enhanced negative electric fields as the "sunrise enhancement." To investigate the sunrise enhancement at KSC we measured the electric field (E) in the first few hundred meters above the ground before and during several sunrise enhancements. From these E soundings we can infer the presence of charge layers and determine their thickness and charge density.

Marshall, T. C.↗

Assembly of Colloidal Aggregates by Electrohydrodynamic Flow: Kinetic Experiments and Scaling Analysis

Electric fields generate transverse flows near electrodes that sweep colloidal particles into densely packed assemblies. We interpret this behavior in terms of electrohydrodynamic motion stemming from distortions of the field by the particles that alter the body force distribution in the electrode charge polarization layer. A scaling analysis shows how the action of the applied electric field generates fluid motion that carries particles toward one another. The resulting fluid velocity is proportional to the square of the applied field and decreases inversely with frequency. Experimental measurements of the particle aggregation rate accord with the electrohydrodynamic theory over a wide range of voltages and frequencies.

Ristenpart, W. D.↗

Numerical Investigation of Two-Phase Flows With Charged Droplets in Electrostatic Field

A numerical method to solve two-phase turbulent flows with charged droplets in an electrostatic field is presented. The ensemble-averaged Navier-Stokes equations and the electrostatic potential equation are solved using a finite volume method. The transitional turbulence field is described using multiple-time-scale turbulence equations. The equations of motion of droplets are solved using a Lagrangian particle tracking scheme, and the inter-phase momentum exchange is described by the Particle-In-Cell scheme. The electrostatic force caused by an applied electrical potential is calculated using the electrostatic field obtained by solving a Laplacian equation and the force exerted by charged droplets is calculated using the Coulombic force equation. The method is applied to solve electro-hydrodynamic sprays. The calculated droplet velocity distributions for droplet dispersions occurring in a stagnant surrounding are in good agreement with the measured data. For droplet dispersions occurring in a two-phase flow, the droplet trajectories are influenced by aerodynamic forces, the Coulombic force, and the applied electrostatic potential field.

Kim, Sang-Wook↗

The Precise Orbit and the Challenge of Long Term Stability

The computation of a precise orbit reference is a fundamental component of the altimetric measurement. Since the dawn of the modern altimeter age, orbit accuracy has been determined by the quality of the GPS, SLR, and DORIS tracking systems, the fidelity of the measurement and force models, and the choice of parameterization for the orbit solutions, and whether a dynamic or a reduced-dynamic strategy is used to calculate the orbits. At the start of the TOPEX mission, the inaccuracies in the modeling of static gravity, dynamic ocean tides, and the nonconservative forces dominated the orbit error budget. Much of the error due to dynamic mismodeling can be compensated by reduced-dynamic tracking techniques depending on the measurement system strength. In the last decade, the launch of the GRACE mission has eliminated the static gravity field as a concern, and the background force models and the terrestrial reference frame have been systematically refined. GPS systems have realized many improvements, including better modeling of the forces on the GPS spacecraft, large increases in the ground tracking network, and improved modeling of the GPS measurements. DORIS systems have achieved improvements through the use of new antennae, more stable monumentation, and of satellite receivers that can track multiple beacons, and as well as through improved modeling of the nonconservative forces. Many of these improvements have been applied in the new reprocessed time series of orbits produced for the ERS satellites, Envisat, TOPEX/Poseidon and the Jason satellites, and as well as for the most recent Cryosat-2 and HY2A. We now face the challenge of maintaining a stable orbit reference for these altimetric satellites. Changes in the time-variable gravity field of the Earth and how these are modelled have been shown to affect the orbit evolution, and the calibration of the altimetric data with tide gauges. The accuracy of the reference frame realizations, and their projection into the future remains a source of error. Other sources of omission error include the geocenter for which no consensus model is as of yet applied. Although progress has been made in nonconservative force modeling through the use of detailed satellite-specific models, radiation pressure modeling, and atmospheric density modeling remain a potential source of orbit error. The longer term influence of variations in the solar and terrestrial radiation fields over annual and solar cycles remains principally untested. Also the long term variation in optical and thermal properties of the space vehicle surfaces would contribute to biases in the orbital frame if ignored. We review the status of altimetric precision orbit determination as exemplified by the recent computations undertaken by the different analysis centers for ERS, Envisat, TOPEX/Poseidon, Jason, Cryosat2 and HY2A, and we provide a perspective on the challenges for future missions such as the Jason-3, SENTINEL-3 and SWOT.

Lemoine, Frank G.↗

Global Patterns of Tectonism on Titan from Mountain Chains and Virgae

This research is based on the exploration of tectonic patterns on Titan from a global perspective. Several moons in the outer solar system display patterns of surface tectonic features that imply global stress fields driven or modified by global forces. Patterns such as these are seen in Europa's tidally induced fracture patterns, Enceladus's tiger stripes, and Ganymede's global expansion induced normal fault bands. Given its proximity to Saturn, as well as its eccentric orbit, tectonic features and global stresses may be present on Titan as well. Titan displays possible tectonic structures, such as mountain chains along its equator (Radebaugh et al. 2007), as well as the unexplored dark linear streaks termed virgae by the IAU. Imaged by Cassini with the RADAR instrument, mountain chains near the equator are observed with a predominante east-west orientation (Liu et al. 2012, Mitri et al. 2010). Orientations such as these can be explained by modifications in the global tidal stress field induced by global contraction followed by rotational spin-up. Also, due to Titan's eccentric orbit, its current rotation rate may be in an equilibrium between tidal spin-up near periapsis and spin-down near apoapsis (Barnes and Fortney 2003). Additional stress from rotational spin-up provides an asymmetry to the stress field. This, combined with an isotropic stress from radial contraction, favors the formation of equatorial mountain chains in an east-west direction. The virgae, which have been imaged by Cassini with both the Visual and Infrared Mapping Spectrometer (VIMS) and Imaging Science Subsystem (ISS) instruments, are located predominately near 30 degrees latitude in either hemisphere. Oriented with a pronounced elongation in the east-west direction, all observed virgae display similar characteristics: similar relative albedos as the surrounding terrain however darkened with an apparent neutral absorber, broken-linear or rounded sharp edges, and connected, angular elements with distinct, linear edges. Virgae imaged during northern latitude passes are oriented with their long dimensions toward Titan's antiSaturn point. If the virgae are of tectonic origin, for instance if the turn out to be i.e. grabens, they could serve as markers to Titan's global stress field. Using them in this way allows for a mapping of global tectonic patterns. These patterns will be tested for consistency against the various sources of global stress and orientations of mountain chains. By determining what drives Titan's tectonics globally, we will be able to place Titan within the context of the other outer planet icy satellites.

Cook, C.↗

Gamma-Ray Pulsar Light Curves in Vacuum and Force-Free Geometry

Recent studies have shown that gamma-ray pulsar light curves are very sensitive to the geometry of the pulsar magnetic field. Pulsar magnetic field geometries, such as the retarded vacuum dipole and force-free magnetospheres have distorted polar caps that are offset from the magnetic axis in the direction opposite to rotation. Since this effect is due to the sweepback of field lines near the light cylinder, offset polar caps are a generic property of pulsar magnetospheres and their effects should be included in gamma-ray pulsar light curve modeling. In slot gap models (having two-pole caustic geometry), the offset polar caps cause a strong azimuthal asymmetry of the particle acceleration around the magnetic axis. We have studied the effect of the offset polar caps in both retarded vacuum dipole and force-free geometry on the model high-energy pulse profiles. We find that, compared to the profiles derived from symmetric caps, the flux in the pulse peaks, which are caustics formed along the trailing magnetic field lines, increases significantly relative to the off-peak emission, formed along leading field lines. The enhanced contrast produces improved slot gap model fits to Fermi pulsar light curves like Vela, with vacuum dipole fits being more favorable.

Harding, Alice K.↗

Implementation of robotic force control with position accommodation

As the need for robotic manipulation in fields such as manufacturing and telerobotics increases, so does the need for effective methods of controlling the interaction forces between the manipulators and their environment. Position Accommodation (PA) is a form of robotic force control where the nominal path of the manipulator is modified in response to forces and torques sensed at the tool-tip of the manipulator. The response is tailored such that the manipulator emulates a mechanical impedance to its environment. PA falls under the category of position-based robotic force control, and may be viewed as a form of Impedance Control. The practical implementations are explored of PA into an 18 degree-of-freedom robotic testbed consisting of two PUMA 560 arms mounted on two 3 DOF positioning platforms. Single and dual-arm architectures for PA are presented along with some experimental results. Characteristics of position-based force control are discussed, along with some of the limitations of PA.

Ryan, Michael J.↗

Exhaust flow and propulsion characteristics of a pulsed MPD arc thruster

An experimental investigation of the near-field megawatt, single-shot exhaust is presented for the self-field and auxiliary field cases (0-2 T). Plasma impact pressure and number density are correlated to provide velocity profiles (30,000 to 70,000 m/sec), thrust (10 to 120 N), impulse (3 to 16 N-sec) and mass accounting. The data agree with Hugel's self-field theory for the case where thrust is produced entirely by electromagnetic force. The data show that the thrust monotonically increases with auxiliary field.

Michels, C. J.↗

Exhaust flow and propulsion characteristics of a pulsed MPD arc thruster.

Experimental investigation of the near-field, megawatt, single-shot exhaust for the self-field and auxiliary field cases (0 to 2 tesla) for a pulsed MPD arc thrustor. Plasma impact pressure and number density are correlated to provide velocity profiles (30,000 to 70,000 m/sec), thrust (10 to 120 N), impulse (3 to 16 N-sec) and mass accounting. The data agree with Hugel's self-field theory for the case where thrust is produced entirely by electromagnetic force. The data show that the thrust monotonically increases with auxiliary field.

Michels, C. J.↗

Qualifying the Sunpower M87N Cryocooler for Operation in the AMS-02 Magnetic Field

The Alpha Magnetic Spectrometer-02 (AMs-02) experiment uses a superfluid helium dewar to cool a large superconducting magnet. The outer vapor-cooled shields of the dewar are to be held at 80 K by four Sunpower M87N cryocoolers. These cryocoolers have magnetic components that might interact with the external applied field generated by the superconducting magnet, thereby degrading the cryocoolers' performance. Engineering models of the Sunpower M87 have been qualified for operation in a magnetic environment similar to the AMs-02 magnetic environment. Although there was no noticeable performance degradation at field levels that were comparable to AMs-02 field levels, there appears to be a small performance degradation at higher field levels. It was theorized that there were three possible issues related to these performance losses at high magnetic fields: i) induced piston rubbing on the cylinder wall due to forces and torques on the linear motor due to the applied magnetic fields; ii) Magnetic hysteretic and/or eddy current damping of the balancer due to its motion in the applied magnetic fields; iii) Inductance losses in motor due to the applied magnetic field. The experiments conducted at the Massachusetts Institute of Technology (MIT) cyclotron facility in June 2002 were designed to test these. Tests were performed over a range of field levels that were lower, comparable, and higher than the field levels that the cryocoolers will experience in the AMs-02 operating environment. This paper describes the experiments and the inferences derived from them.

Mustafi, Shuvo↗

Convection during Thermally Unstable Solidification of Pb-Sn in a Magnetic Field

Convection and macrosegregation in directionally solidified hypoeutectic Pb-38 wt pct Sn and hypereutectic Pb-64.5 wt pct Sn have been examined during upward and downward growth. Temperature fluctuations are observed along the length of the melt column during downward growth. With increasing Rayleigh number, these fluctuations change from none, to cyclic, to time periodic having multiple harmonics, and finally to random. At the higher convective driving force of 350 K temperature inversion, the transverse magnetic field decreased convective levels, strong random temperature fluctuations (flows) becoming smaller and periodic. The maximum field of 0.45 T was unable to completely eliminate convection. For the lower convective driving force of 150 K temperature inversion, the 0.05 T magnetic field decreased flows, and at 0.15 T, the field caused a dramatic decrease in the characteristic frequency of the temperature fluctuations, indicating a change in the nature of the flow, the waveform of the temperature fluctuations changing from sinusoidal to a pulsed wave. Temperature fluctuations and time delays between thermocouples were used to estimate flow velocities. Irrespective of the convection in the bulk melt (ahead of the mushy zone), longitudinal macrosegregation occurs only if the interdendritic melt mixes with the bulk melt.

Song, H.↗

Probabilistic Path Planning of Montgolfier Balloons in Strong, Uncertain Wind Fields

Lighter-than-air vehicles such as hot-air balloons have been proposed for exploring Saturn s moon Titan, as well as other bodies with significant atmospheres. For these vehicles to navigate effectively, it is critical to incorporate the effects of surrounding wind fields, especially as these winds will likely be strong relative to the control authority of the vehicle. Predictive models of these wind fields are available, and previous research has considered problems of planning paths subject to these predicted forces. However, such previous work has considered the wind fields as known a priori, whereas in practical applications, the actual wind vector field is not known exactly and may deviate significantly from the wind velocities estimated by the model. A probabilistic 3D path-planning algorithm was developed for balloons to use uncertain wind models to generate time-efficient paths. The nominal goal of the algorithm is to determine what altitude and what horizontal actuation, if any is available on the vehicle, to use to reach a particular goal location in the least expected time, utilizing advantageous winds. The solution also enables one to quickly evaluate the expected time-to-goal from any other location and to avoid regions of large uncertainty. This method is designed for balloons in wind fields but may be generalized for any buoyant vehicle operating in a vector field. To prepare the planning problem, the uncertainty in the wind field is modeled. Then, the problem of reaching a particular goal location is formulated as a Markov decision process (MDP) using a discretized space approach. Solving the MDP provides a policy of what actuation option (how much buoyancy change and, if applicable, horizontal actuation) should be selected at any given location to minimize the expected time-to-goal. The results provide expected time-to-goal values from any given location on the globe in addition to the action policy. This stochastic approach can also provide insights not accessible by deterministic methods; for example, one can evaluate variability and risk associated with different scenarios, rather than only viewing the expected outcome.

Wolf, Michael↗

Lorentz resonances and the structure of the Jovian ring

Charged dust orbiting through spatially periodic planetary magnetic fields will experience time-variable electromagnetic forces. When the forcing frequencies are nearly commensurate with the particle's orbital frequency, the particle undergoes large out-of-plane and radial excursions. Specific 'Lorentz' resonances, corresponding to particular spatial periodicities in the magnetic field, occur on either side of synchronous orbit. Lorentz resonance locations and strengths for the Jovian and Saturnian rings are described. The boundaries of the halo of the Jovian ring, and perhaps other ring structures, are near resonances.

Burns, J. A.↗

Thermal Field Imaging Using Ultrasound

It is often desirable to be able to determine the temperature field in the interiors of opaque fluids forced into convection by externally imposed temperature gradients. To measure the temperature at a point in an opaque fluid in the usual fashion requires insertion of a probe, and to determine the full field therefore requires either the ability to move this probe or the introduction of multiple probes. Neither of these solutions is particularly satisfactory, although they can lead to quite accurate measurements. As an alternative we have investigated the use of ultrasound as a relatively non-intrusive probe of the temperature field in convecting opaque fluids. The temperature dependence of the sound velocity can be sufficiently great to permit a determination of the temperature from timing the traversal of an ultrasound pulse across a chamber. In this paper we will present our results on convecting flows of transparent and opaque fluids. Our experimental cells consist of relatively narrow rectangular cavities made of thermally insulating materials on the sides, and metal top and bottom plates. The ultrasound transducer is powered by a pulser/receiver, the signal output of which goes to a very high speed signal averager. The average of several hundred to several thousand signals is then sent to a computer for storage and analysis. The experimental procedure is to establish a convective flow by imposing a vertical temperature gradient on the chamber, and then to measure, at several regularly spaced locations, the transit time for an ultrasound pulse to traverse the chamber horizontally (parallel to the convecting rolls) and return to the transducer. The transit time is related to the temperature of the fluid through which the sound pulse travels. Knowing the relationship between transit time and temperature (determined in a separate experiment), we can extract the average temperature across the chamber at that location. By changing the location of the transducer it is then possible to find the average temperature at different locations along the chamber, thereby determining the temperature profile along the system. (In the future we will construct an array of transducers. This will give us the capability to determine the temperature profile much more rapidly than at present, an important consideration if time-dependent phenomena are to be studied.) To validate our procedure we introduced encapsulated liquid crystal particles into glycerol. The liquid crystal particles' color varies depending on the temperature of the fluid. A photograph of the fluid through transparent sidewalls therefore gives a picture of the temperature field of the convecting fluid, independent of our ultrasound imaging. A representative result is shown in the Figure 1, which reveals a very satisfying correspondence between the two techniques. Therefore we have a great deal of confidence that the ultrasound imaging approach is indeed measuring the actual temperature profile of the fluid. The technique has also been applied to convecting liquid metal flows, and representative data will be presented from those experiments as well.

Andereck, D.↗

Key Observations and Discoveries from the Vector Electric Field Investigation on the C/NOFS Satellite

Launched in 2008 and operating for 7.5 years, the Air Force Communication /Navigation Outage Forecasting System (C/NOFS) satellite included the Vector Electric Field Investigation (VEFI) designed and built at NASA’s Goddard Space Flight Center. VEFI successfully met its objectives: 1) investigate the role of ambient electric fields in initiating nighttime density depletions and turbulence; 2) determine the quasi-DC electric fields associated with abrupt density depletions, and 3) quantify the spectrum of the irregularities associated with density depletions, providing many key observations and discoveries including: -- Global (low latitude) vector DC electric fields at 16 s/sec revealing variations with local time, longitude, and season between extremely low and moderate solar activity-- Reversed zonal E x B drifts below the F-region ledge at sunset and simultaneous observations of large scale Kelvin-Helmholtz instabilities as seeds of spread-F -- Large expanses of quasi-coherent kilometer-scale vector wave observations (electric field and density) below 450 km and their discovery as a source of scintillations using the C/NOFS GPS -- Reversed zonal DC electric fields and simultaneous observations of afternoon counter electrojet-- Enhanced zonal DC electric fields at sunrise and full vector plasma flow continuity at the terminator-- Vector electric field and density irregularities extending to meter-scales within equatorial plasma depletions -- Intense electric field structure within the equatorial ionosphere at night in the absence of density depletions -- First observations of the westward equatorial electrojet in post-midnight ionosphere and possible association with downward meridional drifts-- Spaceborne Dst observations and implications for ring current asymmetries-- Vector electric and magnetic field Poynting flux within depletions and TIDs-- Measurements of ionospheric reflectance and Alfvenic waves associated with TIDs-- Observations of Alfven resonators -- ULF magnetic field structure within the nightime equatorial ionosphere-- Discovery of Schumman resonances in space and implications of ELF radiowave propagation-- Vector observations of 50-60 Hz powerline radiation without harmonics-- Ion cyclotron resonance absorption lines associated with ELF hiss and the identification of ambient ions-- Parallel electric fields associated with lower hybrid waves driven by thunderstorm lightning-related sferics-- Discovery of Z-mode radiation in the equatorial ionosphere including its associated with density depletions-- Optical lightning detector waveforms and electric field sferics observed up to altitudes of 800 km

Robert F Pfaff↗

Magnetohydrodynamics (MHD) Aerocapture System for Enabling Faster-Larger Planetary Science & Human Exploration Missions

Since our completing the NIAC Phase I NIAC on this Advanced Aerocapture System, NASA Langley Research Center has funded or supported a number of studies and code enhancements through its Center Innovation Fund (CIF) and NASA’s NSTGRO and Internship Programs to mature the analysis capabilities and quantify the merits of the MHD Aerocapture System technology. These efforts have resulted in a plug and play analysis capability for assessing MHD aerocapture system performance for arrival at many planetary bodies of interest. Our efforts have especially focused on the potential mass savings for improving the capacity for science observations at Neptune and Triton. A re-cent Forbes article published “‘Orbital mechanics is probably going to decide for us whether we go to Uranus or Neptune because we need to flyby Jupiter,’ said Kunio Sayanagi at Hampton University, Virginia, who also worked on the Neptune Odyssey proposal…. Exactly when a mission can be sent to Uranus, or Neptune, depends on the relative position of Jupiter, which can help give a spacecraft a gravitational slingshot. That drastically shortens the cruise phase.” [1] Since shortening the cruise phase is important for these science missions, any mass savings enabled by the MHD Aerocapture System could be reallocated to increasing Thermal Protection System mass to allow faster arrival speeds and/or for onboarding additional payloads such as science instruments, batteries, or propellant for conducting more science for longer durations in the desired orbits. The analysis steps and codes for conducting trades and sizing vehicles for aerocapture are as follows: Step 1 is to conduct aeroheating analysis using LAURA of the selected entry vehicle shape to identify locations on the forebody where ionization and flow velocity are sufficient for producing Lo-rentz forces. LAURA is a multiblock structured grid finite-volume CFD solver developed at the NASA Langley Research Center. [2] LAURA has been used for aerothermal analysis support of the entry, de-scent and landing (EDL) phase of interplanetary missions over the last three decades [3-7]. Step 2 is to port the LAURA results into CFDWARP to calcu-late electrical and thermal conductivities of ionized flow for sizing MHD patch system and calculating Lorentz forces needed for controls analysis. CFDWARP is a CFD code that uses advanced nu-merical methods that enable the simulation of the full coupling between the aerodynamics, the magne-tohydrodynamics, and the non-neutral plasma sheaths. CFDWARP has the unique capability to simulate efficiently the non-neutral sheaths (near the electrodes) in coupled form with the quasi-neutral bulk MHD flow [8-11]. Step 3 is to link re-sults from LAURA and CFDWARP into POST2 for calculating entry trajectories and comparing MHD control results with other aerodynamic control strategies. The Program to Optimize Simulated Tra-jectories II (POST2) is a generalized point mass, discrete parameter targeting and optimization pro-gram. POST2 provides the capability to target and optimize point mass trajectories for multiple pow-ered or un-powered vehicles near an arbitrary rotat-ing, oblate planet [12]. Step 4: TPS sizing was per-formed using the Fully Implicit Ablation and Ther-mal-response code (FIAT) tool which computes the transient one-dimensional thermal response and surface thermochemistry of a multilayer stackup of thermal protection, bonding, and structural materi-als subject to aeroheating on one surface [13]. The sizing and margining methodology used was based on the approach documented by Mahzari and Milos [14] for the dual-layered heatshield for extreme entry environment technology (DL-HEEET) TPS concept. TPS analysis utilizes trajectory information from POST2. Using this step-wise plug and play MHD Aerocapture performance assessment process, our analysis targets a Neptune aerocapture trajectory that will place the spacecraft in an observation orbit for Triton. [15]. Magnetohydrodynamic (MHD) control of a 4.5-meter diameter MSL-style capsule resulted in TPS mass savings of nearly 2000 kg when using an MHD system mass of under 200 kg. The flight path for a vehicle using the MHD control strategy has a much lower heat rate and heat load compared to the conventional aerodynamic aerocapture strategies known as bank angle con-trolled (BAC) and direct force controlled (DFC). Both BAC and DFC have heat rates significantly greater than 1500 W/cm2 typically used as an upper limit for PICA. Thus, DL-HEEET TPS concept was required for the BAC and DFC control strategies. However, considering the more benign environ-ments for the MHD case, additional TPS concepts with improved mass efficiency were also assessed. PICA was considered for the MHD controlled strat-egy since the maximum heat rate was well within the limits (<1500 W/cm2) of PICA. TPS sizing re-sulted in a significant mass reduction. The PICA layer for this sizing case was about 7.8 cm. As a point of reference, the Mars 2020 mission, which used this same PICA concept, had a PICA thickness of 3.18 cm [16]. The trajectories used for the TPS sizing originat-ed from the POST2 simulations. The current, I, to an electromagnet configuration can be manipulated to allow for active control of the vehicle. Manipula-tion of the current, I, changes the magnetic field, B, which affects the Lorentz force and therefore the MHD drag force on the vehicle. Our analysis in-cluded both open-loop and close-loop control. Closed-loop control will enable improved overall performance when taking into account mission level uncertainties, such as interplanetary delivery errors and atmospheric modeling uncertainties. The open-loop and closed-loop MHD control cases do not dip as deep into the atmosphere as the aerodynamic cases. Three types of aerodynamic-only approach-es are investigated: bank angle modulation (BAM), director force control (DFC), and Drag Modulated. BAM and DFC make use of vehicle aerodynamic angles to steer the vehicle. Thus, changing the aer-odynamic forces acting on the vehicle for control, aerodynamic drag modulated case requires a vary-ing drag area to modulate the drag force. The MHD drag modulated case modulates MHD generated drag force that adds to the aerodynamic drag. This higher atmospheric activation of drag forces by the MHD patch results in significantly less heat flux on the vehicle. The MHD technology will enable shorter cruise times and deceleration of larger payloads for increasing the capacity for science at the Ice Giants or for returning astronauts to Earth from cislunar space or from Mars. The purpose of this presentation is to provide more details about this work and to highlight plans for further research and development including a flight demonstration.

R. W. Moses↗

A model for the formation of spokes in Saturn's rings

Evidence is mounting which implies that the generation, evolution, and motion of spokes require a dense plasma near the Saturn ring plane. It is presently suggested that spokes consist of charged, micron-sized dust particles which are elevated from the rings by radially moving, dense plasma columns due to meteor impacts on the ring. Electrostatic wall sheaths at the ring, and ring charging by electric fields sufficiently strong to overcome the gravitational force on small dust particles, arise from the dense plasma, which also increases the probability of dust particle excess electronic charge. The radial motion of the plasma column is due to an azimuthal polarization electric field due to the relative motion between the corotating plasma and the negatively charged dust particles moving at Keplerian speeds.

Goertz, C. K.↗

On the sensitivity of the N-body problem to small changes in initial conditions. III

The effects of increasing the softening parameter, epsilon, in the two body-forces are investigated in order to ascertain the relative importance for this instability of close encounters and the bulk mean field. Consideration is given to both the mean e-folding times for individual particle perturbations and the e-folding times for the total 3N-dimensional configuration space perturbation, which continues to grow exponentially until the typical perturbation of an individual particle becomes comparable to the size of the system. Close encounters and the bulk mean field are both important and apparently reinforce each other: the relatively smooth, systematic tidal force associated with the mean field tends to amplify divergences associated with 'kicks' induced by shorter range effects, which remain important even for large N.

Kandrup, Henry E.↗