Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “three phase modeling”

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,351 records · Page 75

Preliminary design for a reverse Brayton cycle cryogenic cooler

A long life, single stage, reverse Brayton cycle cryogenic cooler is being developed for applications in space. The system is designed to provide 5 W of cooling at a temperature of 65 Kelvin with a total cycle input power of less than 200 watts. Key features of the approach include high speed, miniature turbomachines; an all metal, high performance, compact heat exchanger; and a simple, high frequency, three phase motor drive. In Phase 1, a preliminary design of the system was performed. Analyses and trade studies were used to establish the thermodynamic performance of the system and the performance specifications for individual components. Key mechanical features for components were defined and assembly layouts for the components and the system were prepared. Critical materials and processes were identified. Component and brassboard system level tests were conducted at cryogenic temperatures. The system met the cooling requirement of 5 W at 65 K. The system was also operated over a range of cooling loads from 0.5 W at 37 K to 10 W at 65 K. Input power to the system was higher than target values. The heat exchanger and inverter met or exceeded their respective performance targets. The compresssor/motor assembly was marginally below its performance target. The turboexpander met its aerodynamic efficiency target, but overall performance was below target because of excessive heat leak. The heat leak will be reduced to an acceptable value in the engineering model. The results of Phase 1 indicate that the 200 watt input power requirement can be met with state-of-the-art technology in a system which has very flexible integration requirements and negligible vibration levels.

Swift, Walter L.↗

Fundamentals of the Plasma Sail Concept: MHD and Kinetic Studies

The Mini-Magnetospheric Plasma Propulsion (M2P2), originally proposed by Winglee et al. [2000] predicts that a 15-km standoff distance (or 20-km cross-sectional dimension) of the magnetic bubble will provide for sufficient momentum transfer from the solar wind to accelerate a spacecraft to the unprecedented speeds of 50-80 km/s after an acceleration period of about three months. Such velocities will enable travel out of the solar system in period of about seven years-almost an order of magnitude improvement over present chemical based propulsion systems. However, for the parameters of the simulation of Winglee et al. [2000], a fluid model for the interaction of M2P2 with the solar wind is not valid. It is assumed in the MHD fluid model, normally applied to planetary magnetospheres, that the characteristic scale-size is much greater than the Larmor radius and ion skin depth of the solar wind. In the case of M2P2, the size of the magnetic bubble is actually less than or, comparable to, the scale of these characteristic parameters. Therefore, a kinetic approach, which addresses the small-scale physical mechanisms, must be used. We have adopted a two-component approach to determining a preliminary estimate of the momentum transfer to the plasma sail. The first component is a self-consistent MHD simulation of the small-scale expansion phase of the magnetic bubble. The fluid treatment is valid to roughly 5 km from the source and the steady-state MHD solution at the 5 km boundary was then used as initial conditions for the hybrid simulation. The hybrid simulations showed that the momentum transfer to the innermost regions of the plasma sail is negligible.

Khazanov, G.↗

Design and fabrication of a long-life Stirling cycle cooler for space application. Phase 3: Prototype model

A second-generation, Stirling-cycle cryocooler (cryogenic refrigerator) for space applications, with a cooling capacity of 5 watts at 65 K, was recently completed. The refrigerator, called the Prototype Model, was designed with a goal of 5 year life with no degradation in cooling performance. The free displacer and free piston of the refrigerator are driven directly by moving-magnet linear motors with the moving elements supported by active magnetic bearings. The use of clearance seals and the absence of outgassing material in the working volume of the refrigerator enable long-life operation with no deterioration in performance. Fiber-optic sensors detect the radial position of the shafts and provide a control signal for the magnetic bearings. The frequency, phase, stroke, and offset of the compressor and expander are controlled by signals from precision linear position sensors (LVDTs). The vibration generated by the compressor and expander is cancelled by an active counter balance which also uses a moving-magnet linear motor and magnetic bearings. The driving signal for the counter balance is derived from the compressor and expander position sensors which have wide bandwidth for suppression of harmonic vibrations. The efficiency of the three active members, which operate in a resonant mode, is enhanced by a magnetic spring in the expander and by gas springs in the compressor and counterbalance. The cooling was achieved with a total motor input power of 139 watts. The magnetic-bearing stiffness was significantly increased from the first-generation cooler to accommodate shuttle launch vibrations.

Keung, C.↗

SCATHA mission termination report

The SCATHA (Spacecraft Charging at High Altitudes) satellite was operated from the Consolidated Space Test Center in Sunnyvale, California from February 1979 to May 1991. It was a spin stabilized vehicle in a highly eccentric orbit that collected data on spacecraft charging. The purpose of such data gathering was to predict and/or model the effects of the Earth's magnetic field on synchronous and near synchronous satellites. During the majority of its lifetime, attitude precession maneuvers were done every 10-15 days to maintain solar panel orientation. Maneuver planning was difficult due to the structural characteristics of SCATHA. It is cylindrically shaped and has seven booms ranging in length from 2 to 50 meters. These precession maneuvers induced predictable nutation that damped out after a few days. Eventually fuel began running low due to these frequent maneuvers. Experiments that had required the spin axis be in the orbit plane had already been turned off or had collected all their data. To increase the vehicle lifetime, the spin axis was moved to ecliptic normal. While this stopped the need for frequent attitude maneuvering (only two per year required now), this movement of the spin axis caused nutation that would not damp out for the remainder of the mission. This phase of the mission, with the ecliptic normal orientation, lasted for approximately three years. Although nutation never damped, data gathering was uninterrupted. In late 1990, when SCATHA's transmitter became seriously degraded, the Air Force decided to turn SCATHA off. This would only be done after the satellite was made 'safe'. The most difficult part of making the vehicle safe was quickly purging the fuel. Several plans were considered. The selected plan was to perform a series of 20 degree attitude precession maneuvers (3 days apart to allow for the worst nutation to damp) until the fuel was depleted. Although this sounded simple, the actual execution proved difficult. This was due to a nearly complete lack of available telemetry data, large undamped motion of the long booms, inadequacies in attitude determination software, and an error in the fuel level calculation software. This paper discusses the various proposed termination plans and execution of the selected one. Attitude determination methodologies, nutation from maneuvers, and effects of the flexible booms on the termination mission are presented and analyzed from a satellite analyst point of view.

Stakkestad, Kjell↗

Crystal Growth of ZnSe by Physical Vapor Transport: A Modeling Study

Crystal growth from the vapor phase has various advantages over melt growth. The main advantage is from a lower processing temperature which makes the process more amenable in instances where the melting temperature of the crystal is high. Other benefits stem from the inherent purification mechanism in the process due to differences in the vapor pressures of the native elements and impurities, and the enhanced interfacial morphological stability during the growth process. Further, the implementation of Physical Vapor Transport (PVT) growth in closed ampoules affords experimental simplicity with minimal needs for complex process control which makes it an ideal candidate for space investigations in systems where gravity tends to have undesirable effects on the growth process. Bulk growth of wide band gap II-VI semiconductors by physical vapor transport has been developed and refined over the past several years at NASA MSFC. Results from a modeling study of PVT crystal growth of ZnSe arc reported in this paper. The PVI process is numerically investigated using both two-dimensional and fully three-dimensional formulation of the governing equations and associated boundary conditions. Both the incompressible Boussinesq approximation and the compressible model are tested to determine the influence of gravity on the process and to discern the differences between the two approaches. The influence of a residual gas is included in the models. The preliminary results show that both the incompressible and compressible approximations provide comparable results and the presence of a residual gas tends to measurably reduce the mass flux in the system. Detailed flow, thermal and concentration profiles will be provided in the final manuscript along with computed heat and mass transfer rates. Comparisons with the 1-D model will also be provided.

Ramachandran, Narayanan↗

Development of Optical Diagnostic Techniques for Microgravity Materials Processing

Materials processing including crystal growth, either under a gravity environment on ground or a microgravity environment in space, involves complicated phenomena of fluid motions in gas or liquid phases as well as interaction of various species. To obtain important physical insight, it is very necessary to provide gross-field optical diagnostics for monitoring various physical properties. Materials processing inhibits easy access by ordinary instruments and thus characterizing gross-field physical properties is very challenging. Typical properties of importance can be fluid velocity, temperature, and species concentration for fluids, and surface topology and defects for solids. Observing surface grow rate during crystal growth is also important. Material microstructures, i.e., integrity of crystal structures, is strongly influenced by the existence of thermally-induced flow as well as local nucleation of particles during solidification, which may act in many detrimental ways. In both ground-based and microgravity experiments, the nature of product property changes resulting from three-dimensional fluid or particle motions need be characterized. Gross-field diagnostics is thus required to identify their effects on product defects and process deficiencies. The quantitative visualization techniques can also be used for validation of numerical modeling. For optical nonintrusive gross-field diagnostic techniques, two approaches were developed as summer projects. One optical approach allows us to provide information of species concentration and temperature for monitoring in real time. The other approach, that is, the concept which is formulated for detection of surface topography measurement can provide unprecedented spatial resolution during crystal growth.

Cha, Soyoung Stephen↗

Fuzzy Logic Trajectory Design and Guidance for Terminal Area Energy Management

The second generation reusable launch vehicle will leverage many new technologies to make flight to low earth orbit safer and more cost effective. One important capability will be completely autonomous flight during reentry and landing, thus making it unnecessary to man the vehicle for cargo missions with stringent weight constraints. Implementation of sophisticated new guidance and control methods will enable the vehicle to return to earth under less than favorable conditions. The return to earth consists of three phases--Entry, Terminal Area Energy Management (TAEM), and Approach and Landing. The Space Shuttle is programmed to fly all three phases of flight automatically, and under normal circumstances the astronaut-pilot takes manual control only during the Approach and Landing phase. The automatic control algorithms used in the Shuttle for TAEM and Approach and Landing have been developed over the past 30 years. They are computationally efficient, and based on careful study of the spacecraft's flight dynamics, and heuristic reasoning. The gliding return trajectory is planned prior to the mission, and only minor adjustments are made during flight for perturbations in the vehicle energy state. With the advent of the X-33 and X-34 technology demonstration vehicles, several authors investigated implementing advanced control methods to provide autonomous real-time design of gliding return trajectories thus enhancing the ability of the vehicle to adjust to unusual energy states. The bulk of work published to date deals primarily with the approach and landing phase of flight where changes in heading angle are small, and range to the runway is monotonically decreasing. These benign flight conditions allow for model simplification and fairly straightforward optimization. This project focuses on the TAEM phase of flight where mathematically precise methods have produced limited results. Fuzzy Logic methods are used to make onboard autonomous gliding return trajectory design robust to a wider energy envelope, and the possibility of control surface failures, thus increasing the flexibility of unmanned gliding recovery and landing.

Burchett, Bradley↗

Runtime and Pressurization Analyses of Propellant Tanks

Multi-element unstructured CFD has been utilized at NASA SSC to carry out analyses of propellant tank systems in different modes of operation. The three regimes of interest at SSC include (a) tank chill down (b) tank pressurization and (c) runtime propellant draw-down and purge. While tank chill down is an important event that is best addressed with long time-scale heat transfer calculations, CFD can play a critical role in the tank pressurization and runtime modes of operation. In these situations, problems with contamination of the propellant by inclusion of the pressurant gas from the ullage causes a deterioration of the quality of the propellant delivered to the test article. CFD can be used to help quantify the mixing and propellant degradation. During tank pressurization under some circumstances, rapid mixing of relatively warm pressurant gas with cryogenic propellant can lead to rapid densification of the gas and loss of pressure in the tank. This phenomenon can cause serious problems during testing because of the resulting decrease in propellant flow rate. With proper physical models implemented, CFD can model the coupling between the propellant and pressurant including heat transfer and phase change effects and accurately capture the complex physics in the evolving flowfields. This holds the promise of allowing the specification of operational conditions and procedures that could minimize the undesirable mixing and heat transfer inherent in propellant tank operation. It should be noted that traditional CFD modeling is inadequate for such simulations because the fluids in the tank are in a range of different sub-critical and supercritical states and elaborate phase change and mixing rules have to be developed to accurately model the interaction between the ullage gas and the propellant. We show a typical run-time simulation of a spherical propellant tank, containing RP-1 in this case, being pressurized with room-temperature nitrogen at 540 R. Nitrogen, shown in blue on the right-hand side of the figures, enters the tank from the diffuser at the top of the figures and impinges on the RP-1, shown in red, while the propellant is being continuously drained at the rate of 1050 lbs/sec through a pipe at the bottom of the tank. The sequence of frames in Figure 1 shows the resultant velocity fields and mixing between nitrogen and RP-1 in a cross-section of the tank at different times. A vortex is seen to form in the incoming nitrogen stream that tends to entrain propellant, mixing it with the pressurant gas. The RP-1 mass fraction contours in Figure 1 are also indicative of the level of mixing and contamination of the propellant. The simulation is used to track the propagation of the pure propellant front as it is drawn toward the exit with the evolution of the mixing processes in the tank. The CFD simulation modeled a total of 10 seconds of run time. As is seen from Figure 1d, after 5.65 seconds the propellant front is nearing the drain pipe, especially near the center of the tank. Behind this pure propellant front is a mixed fluid of compromised quality that would require the test to end when it reaches the exit pipe. Such unsteady simulations provide an estimate of the time that a high-quality propellant supply to the test article can be guaranteed at the modeled mass flow rate. In the final paper, we will discuss simulations of the LOX and propellant tanks at NASA SSC being pressurized by an inert ullage. Detailed comparisons will be made between the CFD simulations and lower order models as well as with test data. Conditions leading to cryo collapse in the tank will also be identified.

Field, Robert E.↗

Fatigue Analysis of Notched Laminates: A Time-Efficient Macro-Mechanical Approach

A coupled transversely isotropic deformation and damage fatigue model is implemented within the finite element method and was utilized along with a static progressive damage model to predict the fatigue life, stiffness degradation as a function of number of cycles, and post-fatigue tension and compression response of notched, multidirectional laminates. Initially, the material parameters for the fatigue model were obtained utilizing micromechanics simulations and the provided [0], [90] and [plus or minus 45] experimental composite laminate S-N (stress-cycle) data. Within the fatigue damage model, the transverse and shear properties of the plies were degraded with an isotropic scalar damage variable. The damage in the longitudinal (fiber) ply direction was suppressed, and only the strength of the fiber was degraded as a function of fatigue cycles. A maximum strain criterion was used to capture the failure in each element, and once this criterion was satisfied, the longitudinal stiffness of the element was decreased by a factor of 10 (sup 4). The resulting, degraded properties were then used to calculate the new stress state. This procedure was repeated until final failure of the composite laminate was achieved or a specified number of cycles reached. For post-fatigue tension and compression behavior, four internal state variables were used to control the damage and failure. The predictive capability of the above-mentioned approach was assessed by performing blind predictions of the notched multidirectional IM7/977-3 composite laminates response under fatigue and post-fatigue tensile and compressive loading, followed by a recalibration phase. Although three different multidirectional laminates were analyzed in the course of this study, only detailed results (i.e., stiffness degradation and post-fatigue stress-train curves as well as damage evolution states for a single laminate ([30/60/90/minus 30/minus 60] (sub 2s)) are discussed in detail here.

Composite Fatigue analysis↗

A SUBLIME 3D Model for Cometary Coma Emission: The Hypervolatile-rich Comet C/2016 R2 (PanSTARRS)

The coma of comet C/2016 R2 (PanSTARRS) is one of the most chemically peculiar ever observed, in particular due to its extremely high CO/H2O and N(+2)/H2O ratios, and unusual trace volatile abundances. However, the complex shape of its CO emission lines, as well as uncertainties in the coma structure and excitation, has lead to ambiguities in the total CO production rate. We performed high-resolution, spatially, spectrally, and temporally resolved CO observations using the James Clerk Maxwell Telescope and Submillimeter Array to elucidate the outgassing behavior of C/2016 R2. Results are analyzed using a new, time-dependent, three-dimensional radiative transfer code (SUBlimating gases in LIME; SUBLIME, based on the open-source version of the LIne Modeling Engine), incorporating for the first time, accurate state-to-state collisional rate coefficients for the CO–CO system. The total CO production rate was found to be in the range of (3.8 − 7.6) × 10^(28) per s between 2018 January 13 and February 1 (at r(H) = 2.8–2.9 au), with a mean value of (5.3 ± 0.6) × 10^(28) per s. The emission is concentrated in a near-sunward jet, with a half-opening angle of ∼62° and an outflow velocity of 0.51 ± 0.01 km/s, compared to 0.25 ± 0.01 km/s in the ambient (and nightside) coma. Evidence was also found for an extended source of CO emission, possibly due to icy grain sublimation around 1.2 × 10^(5) km from the nucleus. Based on the coma molecular abundances, we propose that the nucleus ices of C/2016 R2 can be divided into a rapidly sublimating apolar phase, rich in CO, CO2, N2, and CH3OH, and a predominantly frozen (or less abundant), polar phase containing more H2O, CH4, H2CO, and HCN.

Martin Andrew Cordiner↗

Evolution of attached and detached slabs and their associated mantle dynamics

Over the two years of the NASA grant, this project has produced a significant amount of research results related to the plate subduction process and the surface crustal deformation at convergent boundaries (i.e., above subduction zones). While some research objectives are completely accomplished, other research tasks remain active and continue to be investigated at present. A steady state analytic thermal model for subducting slabs was used to examine the torques acting on a descending slab. It is found that gravitational torque vanishes when a slab is dipping either vertically or horizontally, unlike previous studies indicating that the magnitude of gravitational torque decreases as dip angle increases. Subsequently, a new time-dependent, analytic thermal model for a subducting slab was developed. The new model enables us to study transient phenomena associated with plate subduction analytically. On the basis of this model, the nature of slab dip angles was evaluated. Slab dip angles are found to be transient features. As they penetrate into the mantle and increase their lengths, the associated gravitational torque also increases resulting in a downward pulling of the slab to the steeper dip angle. This is especially true once a slab penetrates the olivine-spinel phase boundary at about 400 km depth. However, if the phase transformation does not follow the equilibrium condition, the gravitational torque may have a different behavior. This problem was investigated. Except for fast descending slabs, non-equilibrium phase transformation can only slow down the transient increase of slab dip angles discussed earlier. Its effect is not sufficiently strong to reverse the downward pulling for most of the slabs. However, when slabs subducting at 10 cm/yr or faster, a sufficient amount of metastable olivine can exist beneath 400 km. Because of its low density compared with the surrounding spinel, an upward buoyancy is produced resulting in an upward bending of the slab and possibly an upward rotation of the slab such that smaller dip angles are formed. Seismic studies of the Japanese Slab seem to support this interpretation. The development of oroclinal geometries at convergent boundaries was also examined to study plate obduction which is an important ingredient to the initiation of plate subduction. Although the study suggests that surface features are better modeled by block models, the large scale deformation can be adequately studied by viscous models. Such a model is now under development to complete our original objective to study the initiation of plate subduction. Finally, a three-dimensional, finite element, spherical convective model is developed to study dynamic plate subductions. The model development is now complete and it is being tested to ensure its proper operation. The model is able to generate convection results with a viscosity contrast of about 100. Our research continues to push the viscosity contrast to a level that is appropriate for a subducting slab.

Hsui, Albert T.↗

Orbiter Aerodynamic Acceleration Flight Measurements in the Rarefied-Flow Transition Regime

Acceleration data taken from the Orbital Acceleration Research Experiment (OARE) during reentry on STS-62 have been analyzed using calibration factors taken on orbit. This is the first Orbiter mission which collected OARE data during the Orbiter reentry phase. The data examined include the flight regime from orbital altitudes down to about 90 km which covers the free-molecule-flow regime and the upper altitude fringes of the rarefied-flow transition into the hypersonic continuum. Ancillary flight data on Orbiter position, orientation, velocity, and rotation rates have been used in models to transform the measured accelerations to the Orbiter center-of-gravity, from which aerodynamic accelerations along the Orbiter body axes have been calculated. Residual offsets introduced in the measurements by unmodeled Orbiter forces are identified and discussed. Direct comparisons are made between the OARE flight data and an independent micro-gravity accelerometer experiment, the High Resolution Accelerometer Package (HiRAP), which also obtained flight data on reentry during the mission down to about 95 km. The resulting OARE aerodynamic acceleration measurements along the Orbiter's body axis, aid the normal to axial acceleration ratio in the free-molecule-flow and transition-flow regimes are presented and compared with numerical simulations from three direct simulation Monte Carlo codes.

Blanchard, Robert C.↗

Solid Surface Combustion Experiment Completes a Series of Eight Successful Flights

The Solid Surface Combustion Experiment (SSCE) was the first combustion experiment to fly in the space shuttle and the first such experiment in the NASA spaceflight program since Skylab. SSCE was actually a series of experiments designed to begin to characterize flame spreading over solid fuels in microgravity and the differences of this flame spreading from normal gravity behavior. These experiments should lead to a better understanding of the physical processes involved--increasing our understanding of fire behavior, both in space and on Earth. SSCE results will help researchers evaluate spacecraft fire hazards. These experiments were conceived by the principal investigator, Professor Robert A. Altenkirch, Dean of Engineering at Washington State University. In the first five flights, the fuel sample--ashless filter paper instrumented with three thermocouples--was mounted in a sealed chamber filled with a 50-percent or 35-percent mixture of oxygen in nitrogen at pressures of 1.0, 1.5, and 2.0 atm. In the next three flights, a polymethyl methacrylate (plexiglass) fuel was instrumented with three thermocouples and tested in a 70-percent or 50-percent mixture of oxygen and nitrogen at pressures of 1.0 and 2.0 atm. SSCE is a self-contained, battery-operated experiment that can be flown either in the shuttle middeck or in the Spacelab module. More information about the hardware configuration have been published. This past year, the final two of eight flights were completed on STS-64 and STS-63. The NASA Lewis Research Center designed and built the SSCE payload and performed engineering, testing, scientific, and flight operations support. The SSCE project was supported in some way by nearly every major sector of Lewis' organization. Professor Altenkirch developed a numerical simulation of the flame-spreading process from first principles (of fluid mechanics, heat transfer, and reaction kinetics). The spread rates, flame shape, and thermodynamic data from the SSCE flights are being compared directly with the results of the computational model. Results from the eight flights will be used to formulate an improved solid-phase pyrolysis model. In addition, some results of the flights have been published and presented at international combustion symposiums. Additional solid fuel combustion experiments are being investigated for future tests with the existing hardware.

Source record↗

A Martian general circulation experiment with large topography

A three-layer general circulation model of the Martian atmosphere is described, and the assumptions governing the model are discussed. The simulated, zonally averaged circulation is found to have only limited sensitivity to differences between this model and an earlier general circulation model; this circulation compares reasonably well with observations. It is also found that the meridional mass flow produced by the seasonal condensation of CO2 in the winter polar region has a major influence on the circulation; owing to the weak influence of atmospheric heat transport, however, the mass flow is governed almost entirely by radiation. Quasi-barotropic stationary waves, which are forced kinematically by the topography and which resemble topographically forced terrestrial planetary waves, are generated by the model in the winter hemisphere region of strong eastward flow, while baroclinic stationary waves are thermally forced by topography in the tropics and summer subtropics. It is also concluded that transient baroclinically unstable waves, of somewhat lower dominant wavenumber than those found on the earth, are generated in winter midlatitudes; their amplitudes, wavenumbers, and phase speeds closely agree with what has been deduced from the Viking lander observations.

Pollack, J. B.↗

An application of miniscale experiments on Earth to refine microgravity analysis of adiabatic multiphase flow in space

Adiabatic two-phase flow is of interest to the design of multiphase fluid and thermal management systems for spacecraft. This paper presents original data and unifies existing data for capillary tubes as a step toward assessing existing multiphase flow analysis and engineering software. Comparisons of theory with these data once again confirm the broad accuracy of the theory. Due to the simplicity and low cost of the capillary tube experiments, which were performed on earth, we were able to closely examine for the first time a flow situation that had not previously been examined appreciably by aircraft tests. This is the situation of a slug flow at high quality, near transition to annular flow. Our comparison of software calculations with these data revealed overprediction of pipeline pressure drop by up to a factor of three. In turn, this finding motivated a reexamination of the existing theory, and then development of a new analytical and is in far better agreement with the data. This sequence of discovery illustrates the role of inexpensive miniscale modeling on earth to anticipate microgravity behavior in space and to complete and help define needs for aircraft tests.

Rothe, Paul H.↗

Using Global Simulations to Relate the Three-Part Structure of Coronal Mass Ejections to in Situ Signatures

White-light observations of coronal mass ejections (CMEs) often show the classic "three-part" structure consisting of (1) a bright front; (2) a dark cavity; and (3) a bright, compact core. It has proven difficult to unambiguously associate these features with in situ measurements of interplanetary CMEs (ICMEs), in all but a few cases. In this study we use a global MHD model to simulate the eruption and evolution of a CME out to 0.25 AU, allowing us to continuously track these features from the Sun and through the solar wind. Our results support the generally held view that the interplanetary flux rope corresponds to the dark cavity. We find that the bright front merges with solar wind material swept up by the ICME. Thus, the sheath material found ahead of fast ejecta is in fact composed from both ambient solar wind material, as well the bright front. We also note that, in this simulation, the bright front is formed from the overlying streamer configuration from within which the CME erupted and is not itself coronal material swept up during the early phase of the eruption. The conclusions reached in this study are undoubtedly sensitive to the initial configuration and mechanism used to initiate the CME, and thus care should be taken when using them to interpret specific observations. On the other hand, they provide a unique, unbroken connection between remote solar and interplanetary observations. Ultimately, detailed comparisons between observations and simulation results may be able to constrain or even rule out some mechanisms of CME initiation.

Riley, Pete↗

Using Global Simulations to Relate the Three-Part Structure of Coronal Mass Ejections to In Situ Signatures

White-light observations of coronal mass ejections (CMEs) often show the classic "three-part" structure consisting of (1) a bright front; (2) a dark cavity; and (3) a bright, compact core. It has proven difficult to unambiguously associate these features with in situ measurements of interplanetary CMEs (ICMEs), in all but a few cases. In this study we use a global MHD model to simulate the eruption and evolution of a CME out to 0.25 AU, allowing us to continuously track these features from the Sun and through the solar wind. Our results support the generally held view that the interplanetary flux rope corresponds to the dark cavity. We find that the bright front merges with solar wind material swept up by the ICME. Thus, the sheath material found ahead of fast ejecta is in fact composed from both ambient solar wind material, as well the bright front. We also note that, in this simulation, the bright front is formed from the overlying streamer configuration from within which the CME erupted and is not itself coronal material swept up during the early phase of the eruption. The conclusions reached in this study are undoubtedly sensitive to the initial configuration and mechanism used to initiate the CME, and thus care should be taken when using them to interpret specific observations. On the other hand, they provide a unique, unbroken connection between remote solar and interplanetary observations. Ultimately, detailed comparisons between observations and simulation results may be able to constrain or even rule out some mechanisms of CME initiation

Riley, Pete↗

Investigation of Mars Rotational Dynamics Using Earth-based Radio Tracking of Mars Landers

The development of space geodetic techniques over the past two decades has made it possible to measure the rotational dynamics of the Earth at the milliarcsecond level, improving our geophysical models of the Earth 's interior and the interactions between the solid Earth and its atmosphere. We have found that the rotational dynamics of Mars can be determined to nearly the same level of accuracy by acquiring Earth-based two-way radio tracking observations of three or more landers globally distributed on the surface of Mars. Our results indicate that the precession and long-term obliquity changes of the Mars pole direction can be determined to an angular accuracy corresponding to about 15 cm/yr at the planet's surface. In addition, periodic nutations of the pole and seasonal variations in the spin rate of the planet can be determined to 10 cm or less. Measuring the rotation of Mars at this accuracy would greatly improve the determination of the planet' s moment of inertia and would resolve the size of a planetary fluid core, providing a valuable constraint on Mars interior models. Detecting seasonal variations in the spin rate of Mars would provide global constraints on atmospheric angular momentum changes due to sublimation of the Mars CO2 polar ice caps. Finally, observation of quasisecular changes in Mars obliquity would have significant implications for understanding long-term climatic change. The key to achieving these accuracies is a globally distributed network of Mars landers with stable, phase-coherent radio transponders. By simultaneously acquiring coherent two-way carrier phase observations between a single Earth tracking station and multiple Mars landers, Earth media errors are essentially eliminated, providing an extremely sensitive measure of changes in the differential path lengths between the Earth tracking station and the Mars landers due to Mars rotation. Time variability of the instrumental phase delay through the radio transponder may represent the limiting error source for this technique. Calibration of the transponder stability to about 0.1 ns or less, over a single tracking arc of up to 12 hr, is sufficient to provide the decimeter-level determination of Mars orientation parameters quoted above. We will provide a detailed description of the multilander tracking technique and the requirements it imposes on both the lander radio system and the Earth-based ground-tracking system. This concept is currently part of the strawman science plan for the Mars Environmental Survey (MESUR) mission and complements many of the other MESUR science goals.

Edwards, C. D., Jr.↗