SEARCH · Engineering Papers

Results for “Simulation Results”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS1995

Effect of Stacking Faults on the X-Ray Diffraction Profiles of Beta-SiC Powders

X-ray diffraction patterns or beta-SiC (3C or the cubic polytype or sic) powders often exhibit an additional peak at d = 0.266 nm, high background intensity around the (111) peak, and relative intensities for peaks which differ from those predicted from the crystal structure. Computer simulations were used to show that all these features are due to stacking faults in the powders and not due to the presence of other polytypes in the powders. Such simulations allow diffraction patterns to be generated for different types, frequencies, and spatial distribution or faults. Comparison of the simulation results to the XRD data indicates that the B-SiC particles consist either of heavily faulted clusters distributed irregularly between regions that have only occasional faults or twins, or the powders consist of two types of particles with different populations of faults: those with a high density of faults and those with only twins or occasional faults. Additional information is necessary to determine which description is correct. However, the simulation results can be used to rule out certain fault configurations.

Pujar, Vijay V.↗

NASA NTRS2003

Three-Dimensional MHD Modeling of The Solar Corona and Solar Wind: Comparison with The Wang-Sheeley Model

We present simulation results from a tilted-dipole steady-state MHD model of the solar corona and solar wind and compare the output from our model with the Wang-Sheeley model which relates the divergence rate of magnetic flux tubes near the Sun (inferred from solar magnetograms) to the solar wind speed observed near Earth and at Ulysses. The boundary conditions in our model specified at the coronal base and our simulation region extends out to 10 AU. We assumed that a flux of Alfven waves with amplitude of 35 km per second emanates from the Sun and provides additional heating and acceleration for the coronal outflow in the open field regions. The waves are treated in the WKB approximation. The incorporation of wave acceleration allows us to reproduce the fast wind measurements obtained by Ulysses, while preserving reasonable agreement with plasma densities typically found at the coronal base. We find that our simulation results agree well with Wang and Sheeley's empirical model.

Usmanov, A. V.↗

NASA NTRS1990

A description and evaluation of TIMER: A time-based terminal flow-control concept

A description of a time-based ATC concept called TIMER (Traffic Intelligence for the Management of Efficient Runway-scheduling) and the results of a fast time and real time computer evaluation are presented. The concept was designed to improve the efficiency of extended terminal area operations (en route approach, transition, and terminal flight to the runway). TIMER integrates en route metering, fuel efficient cruise and profile descents, terminal sequencing and spacing together with computer-generated controller aids, in order to fully use runway capacity and improve efficiency of delay absorption. The concept, by using simplified aircraft models, accommodates both 4-D and non 4-D equipped aircraft and is designed for integration into the manual, voice linked ATC system in an evolutionary manner and still be able to accommodate proposed system upgrade features such as data link and further ground automation. Fast time and real time computer simulation results identify and show the effects and interactions of such key variables as horizon of control, metering fix and final approach delivery time errors, aircraft separation requirements, delay discounting, wind, flight technical error, and knowledge of aircraft final approach speed. The current ATC system has a runway interarrival-error standard deviation of approx. 26 seconds. Simulation results indicate that, with computer aiding, the runway interarrival-error standard deviation for non 4-D equipped traffic can be reduced to the region of 8 to 12 seconds if expected final approach speed is known; however, the reduction is only in the region of 16 to 20 seconds if expected final approach speed is unknown. Another major finding is that en route metering fix delivery-error standard deviation should be kept to less than a number somewhere between 35 to 45 seconds to achieve full runway capacity. This requirement implies the need for either airborne automation or assistance to the controller since the current manual performance in today's en route metering environment is in the order of 1.5 minutes.

Credeur, Leonard↗

NASA NTRS1992

Investigation of forced unsteady separated flows using velocity-vorticity form of Navier-Stokes equations

The phenomenon of forced unsteady separation and eruption of boundary-layer vorticity is a highly-complex, high-Reynolds number flow phenomenon, which abruptly leads to the formation of a dynamic stall vortex as demonstrated earlier by the authors for a NACA 0015 airfoil undergoing constant rate pitch-up motion. This, as well as the results of other researchers, have convincingly demonstrated a complex vortical structure within the state of unsteady separation prior to the evolution of dynamic stall. This phenomenon of vortex eruption, although observed in studying dynamic stall phenomena, is also associated with transition from laminar to turbulence flow and its generic nature has been stressed by many researchers including the present investigators. An unsteady Navier-Stokes (NS) analysis is developed for arbitrarily maneuvering bodies using velocity-vorticity variables; this formulation is nearly form-invariant under a generalized non-inertial coordinate transformation. A fully-implicit uniformly second-order accurate method is used, with the nonlinear convective terms approximated using a biased third-order upwind differencing scheme to be able to simulate higher-Re flows. No explicit artificial dissipation is added. The numerical method is fully vectorized and currently achieves a computational index of 7 micro-seconds per time step per mesh point, using a single processor on a CRAY Y-MP. The simulation results show that the energetic free shear from the leading edge is responsible for the wall viscous layer to abruptly erupt near the center of the counterclockwise rotating eddy in the unsteady boundary layer. Primary, secondary, tertiary and quaternary vortices have been observed before the dynamic stall vortex evolves and gathers its maximum strength. This study will discuss the simulation results of Reynolds number up to Re = 45,000 and will also discuss the efforts of initial acceleration in a specific maneuver, on the evolution of the stall vortex.

Ghia, K. N.↗

NASA NTRS2001

Thermal Analysis of Compressible CO2 Flow for Major Equipment of Fire Detection System

A thermal analysis of the compressible CO2 flow for the Portable Fire Extinguisher (PFE) system has been performed. The purpose of this analysis is to determine the discharged CO2 mass from the PFE tank through the Temporary Sleep Station (TeSS) nozzle in reflecting to the latest design of the extended nozzle, and to evaluate the thermal issues associated to the latest nozzle configuration. A SINDA/FLUINT model has been developed for this analysis. The model includes the PFE tank and the TeSS nozzle, and both have initial temperature of 72 of. In order to investigate the thermal effect on the nozzle due to discharging C02, the PFE TeSS nozzle pipe has been divided into three segments. This model also includes heat transfer predictions for PFE tank inner and outer wall surfaces. The simulation results show that the CO2 discharge rates have fulfilled the minimum flow requirements that the PFE system discharges 3.0 Ibm CO2 in 10 seconds and 5.5 Ibm of CO2 in 45 seconds during its operation. At 45 seconds, the PFE tank wall temperature is 63 OF, and the TeSS nozzle cover wall temperatures for the three segments are 47 OF, 53 OF and 37 OF, respectively. Thermal insulation for personal protection is used for the first two segments of the TeSS nozzle. The simulation results also indicate that at 50 seconds, the remaining CO2 in the tank may be near the triple point (gas, liquid and solid) state and, therefore, restricts the flow.

Zhang, Michael Y.↗

NASA NTRS2002

Thermal Analysis of Compressible CO2 Flow for PFE TeSS Nozzle of Fire Detection System

A thermal analysis of the compressible CO2 flow for the Portable Fire Extinguisher (PFE) system has been performed. A SINDA/FLUINT model has been developed for this analysis. The model includes the PFE tank and the TeSS nozzle, and both have initial temperature of 72 F. In order to investigate the thermal effect on the nozzle due to discharging CO2, the PFE TeSS nozzle pipe has been divided into three segments. This model also includes heat transfer predictions for PFE tank inner and outer wall surfaces. The simulation results show that the CO2 discharge rates have fulfilled the minimum flow requirements that the PFE system discharges 3.0 Ibm CO2 in 10 seconds and 5.5 Ibm of CO2 in 45 seconds during its operation. At 45 seconds, the PFE tank wall temperature is 63 F, and the TeSS nozzle cover wall temperatures for the three segments are 47 F, 53 F and 37 F, respectively. Thermal insulation for personal protection is used for the first two segments of the TeSS nozzle. The simulation results also indicate that at 50 seconds, the remaining CO2 in the tank may be near the triple point (gas, liquid and solid) state and, therefore, restricts the flow.

Zhang, Y. Michael↗

NASA NTRS1995

Kinetic structure of rotational discontinuities: Implications for the magnetopause

Magnetic field rotations in the high ion beta magnetosheath that are part of the magnetopause structure are expected to have only a small normal component. We have studied the properties of rotational discontinuities (RDs) under these conditions, viewed as the limit of weak intermediate shocks (ISs), by performing hybrid simulations with a reflecting wall boundary condition (piston method). With this dynamic formation, the sense and size of rotation are not arbitrarily predetermined, but rather evolve from the given upstream (magnetosheath) and downstream (magnetospheric) boundary conditions, similar to what takes place at the magnetopause. This work focuses on several aspects: the observed minimum shear of RDs, their width, their internal signature, and their relation to ISs in isotropic plasmas. Our simulation results are in agreement with the minimum shear observations, that is, the RDs choose the sense of rotation that corresponds to the minimum angle between the upstream and downstream field vector. The RDs are stable, with a unique scale size. Typical gradient scale half widths are one to four ion inertial lengths with a total width up to ten times of that, in agreement with magnetopause observations. We develop a generalized fluid theory of RDs and discuss the characteristic internal signatures of the rotational layer, comparing the kinetic simulation results to predictions from the generalized fluid theory. The results show that ion inertia, anisotropic pressure, finite Larmor radius effects, nonzero ion heat flux, and reflected ions all contribute to the signatures of RDs on kinetic scales. The RDs may have upstream or downstream wave trains, which become weak for high ion beta and small normal components of the magnetic field. We explain the presence and direction of wave trains in terms of the kinetic properties of the Alfven/ion-cyclotron mode. Away from the RD limit there is a smooth transition to weak intermediate shocks, which have small jumps close to expected Rankine-Hugoniot values. Apart from that, there are few kinetic plasma signatures that distinguish RDs from their neighboring ISs. However, noncoplanar ISs evolve in time into thin RDs. Using the properties of RDs and ISs, we make specific suggestions how these discontinuities can be distinguished observationally in the case of an isotropic plasma.

Krauss-Varban, D.↗

NASA NTRS1997

SEDS Tether M/OD Damage Analyses

The Small Expendable Deployer System (SEDS) was designed to deploy an endmass at the end of a 20-km-long tether which acts as an upper stage rocket, and the threats from the meteoroid and orbital debris (M/OD) particle environments on SEDS components are important issues for the safety and success of any SEDS mission. However, the possibility of severing the tether due to M/OD particle impacts is an even more serious concern, since the SEDS tether has a relatively large exposed area to the M/OD environments although its diameter is quite small. The threats from the M/OD environments became a very important issue for the third SEDS mission, since the project office proposed using the shuttle orbiter as a launch platform instead of the second stage of a Delta II expendable rocket, which was used for the first two SEDS missions. A series of hyper-velocity impact tests were performed at the Johnson Space Center and Arnold Engineering Development Center to help determine the critical particle sizes required to sever the tether. The computer hydrodynamic code or hydrocode called CTH, developed by the Sandia National Laboratories, was also used to simulate the damage on the SEDS tether caused by both the orbital debris and test particle impacts. The CTH hydrocode simulation results provided the much needed information to help determine the critical particle sizes required to sever the tether. The M/OD particle sizes required to sever the tether were estimated to be less than 0.1 cm in diameter from these studies, and these size particles are more abundant in low-Earth orbit than larger size particles. Finally, the authors performed the M/OD damage analyses for the three SEDS missions; i.e., SEDS-1, -2, and -3 missions, by using the information obtained from the hypervelocity impact test and hydrocode simulations results.

Hayashida, K. B.↗

NASA NTRSDate not supplied

Achieving Equilibrium for Dense, Integrated, Vehicle Navigation

Drone usage has proliferated in recent years with many applications that have market-changing potential. Applications that include parcel delivery, wildlife protection, precision farming, law enforcement, and industrial inspection, just to name a few. In this paper, the term drone is used to mean both Unmanned Aerial Vehicle and small Unmanned Aircraft System vehicles. Flight infrastructure can currently only support a few thousand aircraft flying over the United States National Airspace System (NAS) at any given time. A delay atone airport sends ripple effects through the system, causing more delays and missed connections. Once regulations and safety policies are put in place to allow for the widespread use of unmanned drones, the number of aircraft in the NAS is expected to skyrocket to millions, potentially congesting the airspace resulting in possible separation violations. In air traffic control, separation is the concept of keeping an “ownship” aircraft outside a minimum distance from “intruder” aircraft to reduce the risk of the aircraft colliding, as well as preventing accidents due to secondary factors, such as wake turbulence. Maintaining proper separation is a safety critical property for drones in the airspace. This paper addresses separation in time and in distance for high volume corridors (en-route) and lanes(on ground). The requirements and necessary conditions for maintaining proper separation and reaching maximum throughput for a given corridor/lane are addressed assuming unidirectional corridors/lanes where an aircraft arrives from one side and departs from the opposite side. Simulation results are presented that show the presented solution guarantees a set of drones to reach the equilibrium state by adjusting their speed based on their distance to the aircraft in front of them. The equilibrium state is defined as a state when a set of n aircraft moving at a relatively constant speed and uniform spacing from each other in a congested system. A congested system is defined as a state when the aircraft cannot move at their maximum allowed speed. Unlike existing centralized and pre-planned approaches, the proposed solution is fully distributed and enables autonomous aircraft to decide to adjust their speed and distance with respect to the preceding aircraft, dynamically. Simulation results are presented that assess the feasibility of the approach using a large number of drones and evaluate the scalability of the proposed solution.

Mahyar R. Malekpour↗

NASA NTRS1975

Experimental and simulation study results of an Adaptive Video Guidance System /AVGS/

Studies relating to stellar-body exploration programs have pointed out the need for an adaptive guidance scheme capable of providing automatic real-time guidance and site selection capability. For the case of a planetary lander, without such guidance, targeting is limited to what are believed to be generally benign areas in order to ensure a reasonable landing-success probability. Typically, the Mars Viking Lander will be jeopardized by obstacles exceeding 22 centimers in diameter. The benefits of on-board navigation and real-time selection of a landing site and obstacle avoidance have been demonstrated by the Apollo lunar landings, in which man performed the surface sensing and steering functions. Therefore, an Adaptive Video Guidance System (AVGS) has been developed, bread-boarded, and flown on a six-degree-of-freedom simulator.

Schappell, R. T.↗

NASA NTRS2018

Determination of Ballistic Limit for IM7/8552 Using LS-DYNA MAT 261

The goal of the NASA ACC High Energy Dynamic Impact Project is to determine the state of the art of dynamic fracture simulations for high velocity impact for composite fuselage shielding applications. Using a building block approach, several computational models considered under NASA ACC are being validated against test data, starting at unconfigured panels and progressing to configured panels under combined out-of-plane and in-plane loading due to ballistic impact. The computational models being evaluated in this project include MAT 162, MAT 213, MAT 261, SPG, and Peridynamics. In this paper, the simulation results using LS-Dyna Material MAT 261 are presented. In particular, a series of blind predictions for unconfigured panels were performed to determine the ballistic limit or V50 velocity. MAT 261 employs failure approach that is generally physically-based using fracture toughness criteria. The overall material model relies on typical ply-level stiffness properties, similar to MAT 162 and other composite continuum damage material models. The fracture toughness values are based on standard tests, and thus are not subject to extensive calibration. This approach is more efficient than performing extensive optimization studies for calibration of parameters. Also, this approach of relying on physical properties reduces the uncertainty of results, as questions concerning the quality and extent of the calibration studies is no longer relevant. However, it was found that carefully controlled coupon-level tests are needed to accurately obtain the required fracture toughness values. Additionally, it should be noted that there is one significant parameter in MAT 261 that does appear to require calibration, and that is the overall failure strain. This is the strain at which the element is deleted, and is not the same as the strain at which damage begins to accumulate. This failure strain is termed EFS (Effective Failure Strain), and is the maximum effective stain for element failure. Simulations have shown that this value will significantly affect impact response and failure. The paper presents the effect of this element failure strain parameter, along with possible uncertainties in fracture toughness values. With an adjusted appropriate value for EFS, it is seen that simulation results compare well with impact test data for predicted penetration velocity.

Byar, Alan↗

NASA NTRS2003

[Results of Simulation Studies]

Lattice Monte Carlo and off-lattice molecular dynamics simulations of h(sub 1)t(sub 4) and h(sub 4)t(sub l) (head/tail) amphiphile solutions have been performed as a function of surfactant concentration and temperature. The lattice and off-lattice systems exhibit quite different self-assembly behavior at equivalent thermodynamic conditions. We found that in the weakly aggregating regime (no preferred-size micelles), all models yield similar micelle size distributions at the same average aggregation number, albeit at different thermodynamic conditions (temperatures). In the strongly aggregating regime, this mapping between models (through temperature adjustment) fails, and the models exhibit qualitatively different micellization behavior. Incipient micellization in a model self-associating telechelic polymer solution results in a network with a transient elastic response that decays by a two-step relaxation: the first is due to a heterogeneous jump-diffusion process involving entrapment of end-groups within well-defined clusters and this is followed by rapid diffusion to neighboring clusters and a decay (terminal relaxation) due to cluster disintegration. The viscoelastic response of the solution manifests characteristics of a glass transition and entangled polymer network.

Source record↗

NASA NTRS1986

Early time interaction of lithium ions with the solar wind in the AMPTE mission

The early time interaction of an artifically injected lithium cloud with the solar wind is simulated with a one-dimensional hybrid code. Simulation results indicate that the lithium cloud presents an obstacle to the solar wind flow, forming a shock-like interaction region. Several notable features are found: (1) The magnetic field is enhanced up to a factor of about 6, followed by a magnetic cavity downstream. (2) Solar wind ions are slowed down inside the lithium cloud, with substantial upstream reflection. (3) Most of the lithium ions gradually pick up the velocity of the solar wind and move downstream. (4) Intense and short-wavelength electric fields exist ahead of the interaction region. (5) Strong electron heating occurs within the lithium cloud. (6) The convection electric field in the solar wind is modulated in the interaction region. The simulation results are in remarkable agreement with in situ spacecraft measurements made during lithium releases in the solar wind by the AMPTE (Active Magnetospheric Particle Tracer Explorers) Program.

Lui, A. T. Y.↗

NASA NTRSDate not supplied

A Strategic Approach for Dense, Integrated, Vehicle Navigation

Drone usage has been on the rise in recent years with applications that include parcel delivery, wildlife protection, precision farming, law enforcement, and industrial inspection, just to name a few. Once regulations and safety policies are put in place to allow for the widespread use of unmanned drones, the number of aircraft in the National Airspace System (NAS) is expected to skyrocket to millions, potentially congesting the airspace which increases the likelihood of separation violations and possibly incidents. Currently, flight infrastructure can only support a few thousand aircraft flying over the United States National Airspace System (NAS) at any given time. A delay at one airport can send ripple effects throughout the system, causing more delays and missed connections. In air traffic control, separation is the concept of keeping an “ownship” aircraft outside a minimum distance from “intruder” aircraft to reduce the risk of the aircraft colliding, as well as preventing accidents due to secondary factors, such as wake turbulence. Maintaining proper separation is often a safety critical property for fixed-wing drones in the airspace. This paper addresses drone separation in time and distance for high volume corridors (en route) and lanes (on ground), merging as well as crossing intersections of multiple corridors/lanes. In this paper, the term drone is applied to both Unmanned Aerial Vehicle (UAV) and small Unmanned Aircraft System (UAS) vehicles operating autonomously. There exists a gamut of approaches to the merging and crossing problems. At one end of the extreme are the conservative yet low cost and verifiable solutions of today that deal with two drones at a time. At the other end are complex Machine Learning-based solutions with high computing requirements for fully autonomous drones of the future that are expected to handle all contentions. This paper presents a feasible and verifiable strategic approach to these problems that is based on distributed cooperation between the drones and the infrastructure. Three phases of the strategic approach (Prepare, Adjust, Commit) are presented. Simulation results are presented that show the proposed approach is stable and resilient to induced perturbations and guarantees a set of fixed-wing drones to merge and cross intersections by adjusting their speed based on their distance to the aircraft in front of them while remaining in the equilibrium state. The equilibrium state is defined as the state when a set of n aircraft move at a relatively constant speed and uniform spacing from each other in a congested system. A congested system is defined as the state when at least one aircraft cannot move at its maximum allowed speed. Unlike existing centralized and pre-planned approaches, the proposed solution is fully distributed and enables autonomous aircraft to decide to adjust their speed and distance with respect to the preceding aircraft, dynamically. Simulation results are presented that assess the feasibility of the approach.

Distributed↗

NASA NTRSDate not supplied

A Strategic Approach for Dense, Integrated, Vehicle Navigation

Distributed↗

NASA NTRS2021

Present and future aerosol impacts on Arctic climate change in the GISS-E2.1 Earth system model

The Arctic is warming 2 to 3 times faster than the global average, partly due to changes in short-lived climate forcers (SLCFs) including aerosols. In order to study the effects of atmospheric aerosols in this warming, recent past (1990–2014) and future (2015–2050) simulations have been carried out using the GISS-E2.1 Earth system model to study the aerosol burdens and their radiative and climate impacts over the Arctic (>60°N), using anthropogenic emissions from the Eclipse V6b and the Coupled Model Intercomparison Project Phase 6 (CMIP6) databases, while global annual mean greenhouse gas concentrations were prescribed and kept fixed in all simulations. Results showed that the simulations have underestimated observed surface aerosol levels, in particular black carbon (BC) and sulfate (SO2−4), by more than 50 %, with the smallest biases calculated for the atmosphere-only simulations, where winds are nudged to reanalysis data. CMIP6 simulations performed slightly better in reproducing the observed surface aerosol concentrations and climate parameters, compared to the Eclipse simulations. In addition, simulations where atmosphere and ocean are fully coupled had slightly smaller biases in aerosol levels compared to atmosphere-only simulations without nudging. Arctic BC, organic aerosol (OA), and SO2−4 burdens decrease significantly in all simulations by 10 %–60 % following the reductions of 7 %–78 % in emission projections, with the Eclipse ensemble showing larger reductions in Arctic aerosol burdens compared to the CMIP6 ensemble. For the 2030–2050 period, the Eclipse ensemble simulated a radiative forcing due to aerosol–radiation interactions (RFARI) of −0.39 ± 0.01 W/sq. m, which is −0.08 W/sq. m larger than the 1990–2010 mean forcing (−0.32 W/sq. m), of which −0.24 ± 0.01 W/sq. m was attributed to the anthropogenic aerosols. The CMIP6 ensemble simulated a RFARI of −0.35 to −0.40 W/sq. m for the same period, which is −0.01 to −0.06 W/sq. m larger than the 1990–2010 mean forcing of −0.35 W/sq. m. The scenarios with little to no mitigation (worst-case scenarios) led to very small changes in the RFARI, while scenarios with medium to large emission mitigations led to increases in the negative RFARI, mainly due to the decrease in the positive BC forcing and the decrease in the negative SO2−4 forcing. The anthropogenic aerosols accounted for −0.24 to −0.26 W/sq. m of the net RFARI in 2030–2050 period, in Eclipse and CMIP6 ensembles, respectively. Finally, all simulations showed an increase in the Arctic surface air temperatures throughout the simulation period. By 2050, surface air temperatures are projected to increase by 2.4 to 2.6 °C in the Eclipse ensemble and 1.9 to 2.6 °C in the CMIP6 ensemble, compared to the 1990–2010 mean. Overall, results show that even the scenarios with largest emission reductions leads to similar impact on the future Arctic surface air temperatures and sea-ice extent compared to scenarios with smaller emission reductions, implying reductions of greenhouse emissions are still necessary to mitigate climate change.

Aerosols↗

NASA NTRS2017

Temporal Evolution of the Plasma Sheath Surrounding Solar Cells in Low Earth Orbit

Initial results from the PIC simulation and the LEM simulation have been presented. The PIC simulation results show that more detailed study is required to refine the ISS solar array current collection model and to understand the development of the current collection in time. The initial results from the LEM demonstrate that is it possible the transients are caused by solar array interaction with the environment, but there are presently too many assumptions in the model to be certain. Continued work on the PIC simulation will provide valuable information on the development of the barrier potential, which will allow refinement the LEM simulation and a better understanding of the causes and effects of the transients.

Willis, Emily M.↗

NASA NTRS1983

Numerical simulations of counterstreaming plasmas and their relevance to interhemispheric flows

The collisionless expansion of counterstreaming plasmas has been studied in order to elucidate the basic physical processes that may be operating during the initial refilling of depleted flux tubes after a magnetic storm. The numerical technique applied is briefly described, and simulation results are presented. The simulation geometry consisted of two high-density H(+)-O(+) electron plasmas separated by a low-density H(+)-electron plasma. The temporal evolution of the expanding plasmas and the electostatic potential in the region between the two sources is described. The main interacting streams are found to be stable with respect to both the ion acoustic and ion cyclotron modes and only the suprathermal forerunner ions are unstable with respect to the ion cyclotron mode. The results also suggest that a localized potential hill can form at the equator and that this potential hill can play an important role in the subsequent trapping and thermalization of the ion streams.

Singh, N.↗