Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “embedded boundary”

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 325 records · Page 18

A first-order global model of Late Cenozoic climatic change: Orbital forcing as a pacemaker of the ice ages

The development of a theory of the evolution of the climate of the earth over millions of years can be subdivided into three fundamental, nested, problems: (1) to establish by equilibrium climate models (e.g., general circulation models) the diagnostic relations, valid at any time, between the fast-response climate variables (i.e., the 'weather statistics') and both the prescribed external radiative forcing and the prescribed distribution of the slow response variables (e.g., the ice sheets and shelves, the deep ocean state, and the atmospheric CO2 concentration); (2) to construct, by an essentially inductive process, a model of the time-dependent evolution of the slow-response climatic variables over time scales longer than the damping times of these variables but shorter than the time scale of tectonic changes in the boundary conditions (e.g., altered geography and elevation of the continents, slow outgassing, and weathering) and ultra-slow astronomical changes such as in the solar radiative output; and (3) to determine the nature of these ultra-slow processes and their effects on the evolution of the equilibrium state of the climatic system about which the above time-dependent variations occur. All three problems are discussed in the context of the theory of the Quaternary climate, which will be incomplete unless it is embedded in a more general theory for the fuller Cenozoic that can accommodate the onset of the ice-age fluctuations. We construct a simple mathematical model for the Late Cenozoic climatic changes based on the hypothesis that forced and free variations of the concentration of atmospheric greenhouse gases (notably CO2), coupled with changes in the deep ocean state and ice mass, under the additional 'pacemaking' influence of earth-orbital forcing, are primary determinants of the climate state over this period. Our goal is to illustrate how a single model governing both very long term variations and higher frequency oscillatory variations in the Pleistocene can be formulated with relatively few adjustable parameters.

Saltzman, Barry↗

Adjoint Algorithm for CAD-Based Shape Optimization Using a Cartesian Method

Adjoint solutions of the governing flow equations are becoming increasingly important for the development of efficient analysis and optimization algorithms. A well-known use of the adjoint method is gradient-based shape optimization. Given an objective function that defines some measure of performance, such as the lift and drag functionals, its gradient is computed at a cost that is essentially independent of the number of design variables (geometric parameters that control the shape). More recently, emerging adjoint applications focus on the analysis problem, where the adjoint solution is used to drive mesh adaptation, as well as to provide estimates of functional error bounds and corrections. The attractive feature of this approach is that the mesh-adaptation procedure targets a specific functional, thereby localizing the mesh refinement and reducing computational cost. Our focus is on the development of adjoint-based optimization techniques for a Cartesian method with embedded boundaries.12 In contrast t o implementations on structured and unstructured grids, Cartesian methods decouple the surface discretization from the volume mesh. This feature makes Cartesian methods well suited for the automated analysis of complex geometry problems, and consequently a promising approach to aerodynamic optimization. Melvin et developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the Euler equations. In both approaches, a boundary condition is introduced to approximate the effects of the evolving surface shape that results in accurate gradient computation. Central to automated shape optimization algorithms is the issue of geometry modeling and control. The need to optimize complex, "real-life" geometry provides a strong incentive for the use of parametric-CAD systems within the optimization procedure. In previous work, we presented an effective optimization framework that incorporates a direct-CAD interface. In this work, we enhance the capabilities of this framework with efficient gradient computations using the discrete adjoint method. We present details of the adjoint numerical implementation, which reuses the domain decomposition, multigrid, and time-marching schemes of the flow solver. Furthermore, we explain and demonstrate the use of CAD in conjunction with the Cartesian adjoint approach. The final paper will contain a number of complex geometry, industrially relevant examples with many design variables to demonstrate the effectiveness of the adjoint method on Cartesian meshes.

Nemec, Marian↗

Minnowbrook VI: 2009 Workshop on Flow Physics and Control for Internal and External Aerodynamics

Topics covered include: Flow Physics and control for Internal and External Aerodynamics (not in TOC...starts on pg13); Breaking CFD Bottlenecks in Gas-Turbine Flow-Path Design; Streamwise Vortices on the Convex Surfaces of Circular Cylinders and Turbomachinery Blading; DNS and Embedded DNS as Tools for Investigating Unsteady Heat Transfer Phenomena in Turbines; Cavitation, Flow Structure and Turbulence in the Tip Region of a Rotor Blade; Development and Application of Plasma Actuators for Active Control of High-Speed and High Reynolds Number Flows; Active Flow Control of Lifting Surface With Flap-Current Activities and Future Directions; Closed-Loop Control of Vortex Formation in Separated Flows; Global Instability on Laminar Separation Bubbles-Revisited; Very Large-Scale Motions in Smooth and Rough Wall Boundary Layers; Instability of a Supersonic Boundary-Layer With Localized Roughness; Active Control of Open Cavities; Amplitude Scaling of Active Separation Control; U.S. Air Force Research Laboratory's Need for Flow Physics and Control With Applications Involving Aero-Optics and Weapon Bay Cavities; Some Issues Related to Integrating Active Flow Control With Flight Control; Active Flow Control Strategies Using Surface Pressure Measurements; Reduction of Unsteady Forcing in a Vaned, Contra-Rotating Transonic Turbine Configuration; Active Flow Control Stator With Coanda Surface; Controlling Separation in Turbomachines; Flow Control on Low-Pressure Turbine Airfoils Using Vortex Generator Jets; Reduced Order Modeling Incompressible Flows; Study and Control of Flow Past Disk, and Circular and Rectangular Cylinders Aligned in the Flow; Periodic Forcing of a Turbulent Axisymmetric Wake; Control of Vortex Breakdown in Critical Swirl Regime Using Azimuthal Forcing; External and Turbomachinery Flow Control Working Group; Boundary Layers, Transitions and Separation; Efficiency Considerations in Low Pressure Turbines; Summary of Conference; and Final Plenary Session Transcript.

LaGraff, John E.↗

The onset and development of Kelvin-Helmholtz instability at the Venus ionopause

The paper investigates the Kelvin-Helmholtz instability at the Venus ionopause resulting from the flow of the (shocked) solar wind tangential to the ionopause for the case where the interplanetary field is oriented normal to the direction of flow. It is found that gravity stabilizes the long wavelength perturbations, and the finite thickness of the boundary layer stabilizes short wavelength modes. The magnetic 'gyroviscosity' due to finite Larmor radius effects either destabilizes the boundary or stabilizes it according to whether the solar wind electric field points away from or toward the ionosphere. For solar wind and ionosphere plasma parameters consistent with Pioneer Venus observations, it is found that the instabilities with the greatest growth rates (shortest growth times) have wavelengths of 50-150 km and growth times of 0.5 to several seconds. In addition, it is found how distortion of the ionopause by Kelvin-Helmholtz instability might lead to the formation of magnetic 'flux ropes' inside the ionosphere as well as ionospheric 'bubbles' embedded in the solar wind.

Wolff, R. S.↗

Filamentary current structures in the postnoon sector: Observations from UARS

During an intense geomagentic storm (Kp 7+) that began at approximately 1830 UT on October 1, 1991, the Upper Atmosphere Research Satellite (UARS) satellite encountered the dayside postnoon auroral oval. On two consecutive crossings of the northern hemisphere between 2040 and 2240 UT, the vector magnetometer detected region 1 and 2 Birkeland and ionospheric currents in the postnoon sector. Low-energy electron events were observed near 1400 MLT within a narrow portion of the region 1 current system. Simultaneous magnetic field measurements revealed the presence of intense (approximately 20 microA/sq m) bipolar filament current structures embedded in the auroral oval. The upward-directed currents were associated with the more concentrated region of precipitating electrons. Ions associated with the more intense flux of low-energy electrons exhibited a dispersion signature typical of an ion velocity filter. The dispersion, aligned along the orbit, exhibited higher-energy ions at lower latitudes and earlier local times. The colocation of filament currents and ion dispersion signatures at such late postnoon local times is not consistent with typical E x B 'cusp' dispersions. These features more likely result from dayside boundary wave phenomena.

Bythrow, P. F.↗

Development of BEM for ceramic composites

It is evident that for proper micromechanical analysis of ceramic composites, one needs to use a numerical method that is capable of idealizing the individual fibers or individual bundles of fibers embedded within a three-dimensional ceramic matrix. The analysis must be able to account for high stress or temperature gradients from diffusion of stress or temperature from the fiber to the ceramic matrix and allow for interaction between the fibers through the ceramic matrix. The analysis must be sophisticated enough to deal with the failure of fibers described by a series of increasingly sophisticated constitutive models. Finally, the analysis must deal with micromechanical modeling of the composite under nonlinear thermal and dynamic loading. This report details progress made towards the development of a boundary element code designed for the micromechanical studies of an advanced ceramic composite. Additional effort has been made in generalizing the implementation to allow the program to be applicable to real problems in the aerospace industry.

Henry, D. P.↗

Control of Thermal Deflection, Panel Flutter and Acoustic Fatigue at Elevated Temperatures Using Shape Memory Alloys

The High Speed Civil Transport (HSCT) will have to be designed to withstand high aerodynamic load at supersonic speeds (panel flutter) and high acoustic load (acoustic or sonic fatigue) due to fluctuating boundary layer or jet engine acoustic pressure. The thermal deflection of the skin panels will also alter the vehicle's configuration, thus it may affect the aerodynamic characteristics of the vehicle and lead to poor performance. Shape memory alloys (SMA) have an unique ability to recover large strains completely when the alloy is heated above the characteristic transformation (austenite finish T(sub f)) temperature. The recovery stress and elastic modulus are both temperature dependent, and the recovery stress also depends on the initial strain. An innovative concept is to utilize the recovery stress by embedding the initially strained SMA wire in a graphite/epoxy composite laminated panel. The SMA wires are thus restrained and large inplane forces are induced in the panel at elevated temeperatures. By embedding SMA in composite panel, the panel becomes much stiffer at elevated temperatures. That is because the large tensile inplane forces induced in the panel from the SMA recovery stress. A stiffer panel would certainly yield smaller dynamic responses.

Mei, Chuh↗

Investigation of the Rocket Induced Flow Field in a Rectangular Duct

Several tests were performed on a one-sixth scale Rocket Based Combined Cycle (RBCC) engine model at the University of Alabama in Huntsville. The UAH RBCC facility consists of a rectangular duct with a vertical strut mounted in the center. The scaled strut consists of two supersonic rocket nozzles with an embedded vertical turbine between the rocket nozzles. The tests included mass flow, flow visualization and horizontal pressure traverses. The mass flow test indicated a c:hoked condition when the rocket chamber pressure is between 200 psi and 300 psi. The flow visualization tests narrowed the rocket chamber pressure range from, 250 psi to 300 psi. Also, from this t.est, an assumption of a minimum <area choke was made. Horizontal pressure traverses were performed at locations both upstream and downstream of the nozzle exit plane. These traverses aided in examining the induced airflow due to the ejector effect and the growth of the: boundary layer in the strut sidewall gap.

Landrum, D. Brian↗

Theoretical study of stress concentrations at circular holes and inclusions in strain hardening materials.

Nonlinear boundary value problems of an infinite elastic-plastic plate with a circular hole subjected to pure tension and pure shear at infinity are solved by a method involving Fourier series and finite difference. On the basis of these solutions, the validity of Neuber's relationship between the stress and strain concentration factors for the plane stress problems is examined and a generalized Stowell formula for the stress concentration factor is proposed for problems in which the applied loading may be pure shear as well as pure tension and, furthermore, other stress states. By the same method of solution, the stress distributions around a rigid circular cylindrical inclusion embedded in an infinite rigid-plastic matrix subjected to uniform transverse pure shear and tension are obtained.

Huang, W. C.↗

Using a Commercial Ethernet PHY Device in a Radiation Environment

This work involved placing a commercial Ethernet PHY on its own power boundary, with limited current supply, and providing detection methods to determine when the device is not operating and when it needs either a reset or power-cycle. The device must be radiation-tested and free of destructive latchup errors. The commercial Ethernet PHY's own power boundary must be supplied by a current-limited power regulator that must have an enable (for power cycling), and its maximum power output must not exceed the PHY's input requirements, thus preventing damage to the device. A regulator with configurable output limits and short-circuit protection (such as the RHFL4913, rad hard positive voltage regulator family) is ideal. This will prevent a catastrophic failure due to radiation (such as a short between the commercial device's power and ground) from taking down the board's main power. Logic provided on the board will detect errors in the PHY. An FPGA (field-programmable gate array) with embedded Ethernet MAC (Media Access Control) will work well. The error detection includes monitoring the PHY's interrupt line, and the status of the Ethernet's switched power. When the PHY is determined to be non-functional, the logic device resets the PHY, which will often clear radiation induced errors. If this doesn't work, the logic device power-cycles the FPGA by toggling the regulator's enable input. This should clear almost all radiation induced errors provided the device is not latched up.

Parks, Jeremy↗

Plasma Circulation in the Magnetosphere

We investigate the global structure and dynamics of plasma circulation produced by prototypical solar wind disturbances of the interplanetary magnetic field and dynamic pressure. We track the global circulation and energization of solar wind, polar wind, and auroral wind plasmas throughout the magnetosphere, until they precipitate or escape into the downstream solar wind. We use the full equations of motion of the plasma ions within fields produced by a global MHD simulation of the dynamic solar wind interaction. We use the dynamic hot plasma density and Poynting energy flux specified at the inner boundary of the MHD simulation as drivers of conjugate ion outflow fluxes using local empirical relations obtained from the FAST and Polar missions. Birkeland currents computed by the MHD code are used to derive a field-parallel potential drop from a Knight-like relation [as modified by Lyons and Evans, 1980]. This potential drop is applied to each ion as an initial bulk energy, added to a thermal heating driven by the locally incident Poynting flux. The solar wind pressure increase case (B(sub Y) = 5; B(sub z) = 0 nT) produces an immediate substorm owing to compression of pre-existing plasmas. The SB(sub z), interval (embedded in NB(sub z)) produces a substorm after about one hour of development. Both disturbances enhance the auroral wind flux and heavy ion pressure of the magnetosphere substantially, with complex dynamic structuring by auroral acceleration vortexes and dynamic reconnection. Comparisons are made with observations during disturbed periods including the Halloween 2003 super-storm and other periods.

Moore, T. E.↗

Crater Flux Transfer Events: Highroad to the X Line?

We examine Cluster observations of a so-called magnetosphere crater FTE, employing data from five instruments (FGM, CIS, EDI, EFW, and WHISPER), some at the highest resolution. The aim of doing this is to deepen our understanding of the reconnection nature of these events by applying recent advances in the theory of collisionless reconnection and in detailed observational work. Our data support the hypothesis of a stratified structure with regions which we show to be spatial structures. We support the bulge-like topology of the core region (R3) made up of plasma jetting transverse to reconnected field lines. We document encounters with a magnetic separatrix as a thin layer embedded in the region (R2) just outside the bulge, where the speed of the protons flowing approximately parallel to the field maximizes: (1) short (fraction of a sec) bursts of enhanced electric field strengths (up to approximately 30 mV/m) and (2) electrons flowing against the field toward the X line at approximately the same time as the bursts of intense electric fields. R2 also contains a density decrease concomitant with an enhanced magnetic field strength. At its interface with the core region, R3, electric field activity ceases abruptly. The accelerated plasma flow profile has a catenary shape consisting of beams parallel to the field in R2 close to the R2/R3 boundary and slower jets moving across the magnetic field within the bulge region. We detail commonalities our observations of crater FTEs have with reconnection structures in other scenarios. We suggest that in view of these properties and their frequency of occurrence, crater FTEs are ideal places to study processes at the separatrices, key regions in magnetic reconnection. This is a good preparation for the MMS mission.

Farrugia, C. J.↗

Aerosol Meteorology of Maritime Continent for the 2012 7SEAS Southwest Monsoon Intensive Study - Part 2: Philippine Receptor Observations of Fine-Scale Aerosol Behavior

The largest 7 Southeast Asian Studies (7SEAS) operations period within the Maritime Continent (MC) occurred in the August-September 2012 biomass burning season. Data included were observations aboard the MY Vasco, dispatched to the Palawan Archipelago and Sulu Sea of the Philippines for September 2012. At these locations, the Vasco observed MC smoke and pollution entering the southwest monsoon (SWM) monsoonal trough. Here we describe the research cruise findings and the finer-scale aerosol meteorology of this convectively active region. This 2012 cruise complemented a 2-week cruise in 2011 and was generally consistent with previous findings in terms of how smoke emission and transport related to monsoonal flows, tropical cyclones (TC), and the covariance between smoke transport events and the atmosphere's thermodynamic structure. Biomass burning plumes were usually mixed with significant amounts of anthropogenic pollution. Also key to aerosol behavior were squall lines and cold pools propagating across the South China Sea (SCS) and scavenging aerosol particles in their path. However, the 2012 cruise showed much higher modulation in aerosol frequency than its 2011 counterpart. Whereas in 2011 large synoptic-scale aerosol events transported high concentrations of smoke into the Philippines over days, in 2012 measured aerosol events exhibited a much shorter-term variation, sometimes only 312h. Strong monsoonal flow reversals were also experienced in 2012. Nucleation events in cleaner and polluted conditions, as well as in urban plumes, were observed. Perhaps most interestingly, several cases of squall lines preceding major aerosol events were observed, as opposed to 2011 observations where these lines largely scavenged aerosol particles from the marine boundary layer. Combined, these observations indicate pockets of high and low particle counts that are not uncommon in the region. These perturbations are difficult to observe by satellite and very difficult to model. Indeed, the Navy Aerosol Analysis and Prediction System (NAAPS) simulations captured longer period aerosol events quite well but largely failed to capture the timing of high-frequency phenomena. Ultimately, the research findings of these cruises demonstrate the real world challenges of satellite-based missions, significant aerosol life cycle questions such as those the future Aerosol/Clouds/Ecosystems (ACE) will investigate, and the importance of small-scale phenomena such as sea breezes, squall lines, and nucleation events embedded within SWM patterns in dominating aerosol life cycle and potential relationships to clouds.

Reid, Jeffrey S.↗

The Sm-Nd systematics of silicate inclusions in iron meteorites: Results from Caddo (IAB)

The timing of events leading to the formation of silicate-rich and metal-rich regions in planetesimals remains an important problem in the study of planetary formation and differentiation in the early solar system. The IAB irons are especially important as they are considered to represent a magmatic differentiation series. Iron meteorites present a particular challenge for chronological studies, due to the relative paucity of phases serving as hosts for radioactive parent-daughter nuclides. Recent work using the Re-Os system, following on the pioneering work by Herr et al. and Luck and Allegre, appears promising, but investigators up to now have concentrated on whole rock isochrons. Silicate clasts enclosed within iron meteorites can provide information about the chronology and thermal history of irons. Extensive work on Rb-Sr, K-Ar, and I-Xe has been reported on silicate inclusions in iron meteorites. We report the initial results from our Sm-Nd study of an inclusion with the Caddo IAB iron, the first Sm-Nd isotopic study of a silicate clast embedded within an iron meteorite. Our results include measurements of the standard long-lived Sm-147/Nd-143 (tau = 152 AE) system, as well as the shorter-lived SM-146/Nd-142 (tau = 0.149 AE) system, which has been shown to be very useful in deciphering the history of the early solar system. The Caddo silicate clast was described by Palme et al., who kindly provided us with a major part of the inclusion. The inclusion is coarse-grained consisting predominantly of olivine, clinopyroxene, and plagioclase, with lesser amounts of orthopyroxene, Fe-Ni metal, sulfide, and phosphate. The relatively large grain size (up to 3 mm) and 120 degree grain boundaries suggest extensive metamorphism at high temperatures. Based on study of a thin section, there is evidence for metal invading along grain boundaries in some regions of the inclusion, suggesting that the Fe-Ni metal was molten when the silicate clast was incorporated. Metamorphic recrystallization may have occurred during this event.

Stewart, Brian W.↗

Application of the boundary element method to transient heat conduction

An advanced boundary element method (BEM) is presented for the transient heat conduction analysis of engineering components. The numerical implementation necessarily includes higher-order conforming elements, self-adaptive integration and a multiregion capability. Planar, three-dimensional and axisymmetric analyses are all addressed with a consistent time-domain convolution approach, which completely eliminates the need for volume discretization for most practical analyses. The resulting general purpose algorithm establishes BEM as an attractive alternative to the more familiar finite difference and finite element methods for this class of problems. Several detailed numerical examples are included to emphasize the accuracy, stability and generality of the present BEM. Furthermore, a new efficient treatment is introduced for bodies with embedded holes. This development provides a powerful analytical tool for transient solutions of components, such as casting moulds and turbine blades, which are cumbersome to model when employing the conventional domain-based methods.

Dargush, G. F.↗

Lunar regolith densification

Core tube samples of the lunar regolith obtained during the Apollo missions showed a rapid increase in the density of the regolith with depth. Various hypotheses have been proposed for the possible cause of this phenomenon, including the densification of the loose regolith material by repeated shaking from the seismic tremors which have been found to occur at regular monthly intervals when the moon and earth are closest to one another. A test bed was designed to study regolith densification. This test bed uses Minnesota Lunar Simulant (MLS) to conduct shaking experiments in the geotechnical centrifuge with an inflight shake table system. By reproducing realistic in-situ regolith properties, the experiment also serves to test penetrator concepts. The shake table system was designed and used for simulation experiments to study effects of earthquakes on terrestrial soil structures. It is mounted on a 15 g-ton geotechnical centrifuge in which the self-weight induced stresses are replicated by testing an n-th scale model in a gravity field which is n times larger than Earth's gravity. A similar concept applies when dealing with lunar prototypes, where the gravity ratio required for proper simulation of lunar gravity effects is that between the centrifugal acceleration and the lunar gravity. Records of lunar seismic tremors, or moonquakes, were obtained. While these records are being prepared for use as the input data to drive the shake table system, records from the El Centro earthquake of 1940 are being used to perform preliminary tests, using a soil container which was previously used for earthquake studies. This container has a laminar construction, with the layers free to slide on each other, so that the soil motion during the simulated earthquake will not be constrained by the otherwise rigid boundaries. The soil model is prepared by pluviating the MLS from a hopper into the laminar container to a depth of 6 in. The container is mounted on the shake table and the centrifuge is operated to generate an acceleration of 10 times Earth's gravity or 60 times the lunar gravity, thus simulating a lunar regolith thickness of 30 ft. The shake table is then operated using the scaled 'moonquake' as the input motion. One or more model moonquakes are used in each experiment, after which the soil is analyzed for its density profile with depth. This is accomplished by removing from the soil bed a column of soil contained within a thin rubber sleeve which has been previously embedded vertically in the soil during pluviation. This column of soil is transferred to a gamma ray device, in which the gamma ray transmission transversely through the soil is measured and compared with standard calibration samples. In this manner, the density profile can be determined. Preliminary results to date are encouraging, and the Center plans to study the effects of duration of shaking, intensity of the shaking motion, and the frequency of the motion.

Ko, Hon-Yim↗

Flow Trajectories in Friction Stir Welding

In the Friction Stir Welding (FSW) process, a rotating, shouldered tool with a threaded pin is inserted under pressure along the seam of two pieces of metal which are tightly clamped together, and secured against a rigid anvil underneath. The rotating pin travels along the seam and through a combination of pressure and friction heating produces a zone of plastic deformation around the pin within the workpiece on either side of the seam. As the pin is moved in the direction of welding, the plasticised material moves around the tool and bonds together behind it. The elements of the material flow behavior are a combination of three elements. There is a rotational transport of material being carried around the tool, extrusion of material being forced around the pin on both sides into the cavity created behind it, and a lifting and dropping of material as it is stirred and mixed by the rotating action of the pin. It was assumed that rotational motion of the plastic zone is the primary mechanism for transport of material around the welding tool. A kinematic mathematical model was used to compute trajectories of material movement for various distributions of rotational slip within the plastic zone. These trajectories were then compared with the results of an experiment that produced radiographs of markers embedded in a workpiece that was welded with the FSW process. It was assumed that the copper wire markers retained their original length as the aluminum material flowed around them. The kinematic model included a constraint so that the displacements were such that the total length of the wire markers did not increase after deformation. There was good agreement between the calculated trajectories for the case of localized slip at the outer surface of the plastic deformation zone and the radiographs of the copper wire markers. The trajectories differed markedly from the radiograph traces when a distributed slip zone was assumed. It was concluded that the flow field could be described generally by a rotating plug model with slip occurring on the surface at the outer boundary of the plastic zone.

Bernstein, Edward L.↗

Results of the first statistical study of Pioneer Venus Orbiter plasma observations in the distant Venus tail - Evidence for a hemispheric asymmetry in the pickup of ionospheric ions

Pioneer Venus Orbiter plasma and magnetometer observations from the first nine tail seasons of crossings of the Venus wake are used to study ion pickup in the far wake of an unmagnetized object embedded in the solar wind. This statistical study treats all of the plasma spectra containing pickup ions in the vicinity of the Venus tail. A hemispheric asymmetry is found in the pickup of ionospheric ions, with approximately four times more O(+) events observed in the 'northern' magnetic hemisphere (where the induced electric field points outward) than in the 'southern' magnetic hemisphere. Out of a total of 167 large O(+) events, 125, or 75 percent, occurred in the northern hemisphere when position is calculated in terms of Venus radii, and 129 or 77 percent occurred in the 'northern' hemisphere when position is expressed in gyroradii. This hemispheric asymmetry in ion pickup is consistent with the prediction of the Cloutier et al. (1974) mass loading model for Venusian ions above the ionopause boundary.

Intriligator, Devrie S.↗