Search NASA⌕ Search

SEARCH · Search NASA

Results for “hybrid simulations”

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 307 records · Page 17

Wave and particle evolution downstream of quasi-perpendicular shocks

Distributions of ions heated in quasi-perpendicular bow shocks have large perpendicular temperature anisotropies that provide free energy for the growth of Alfven ion cyclotron (AIC) and mirror waves. These modes are often obsreved in the Earth's magnetosheath. Using two-dimensional hybrid simulations, we show that these waves are produced near the shock front and convected downstream rather than being produced locally downstream. The wave activity reduces the proton anisotropy to magnetosheath levels within a few tens of gyroradii of the shock but takes significantly longer to reduce the anisotropy of He(++) ions. The waves are primarily driven by proton anisotropy and the dynamics of the helium ions is controlled by the proton waves. Downstream of high Mach number shocks, mirror waves compete effectively with AIC waves. Downstream of low Mach number shocks, AIC waves dominate.

Mckean, M. E.↗

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.↗

Kinetic Kelvin-Helmholtz instability at a finite sized object

Two-dimensional hybrid simulations with particle ions and fluid electrons are used to calculate the kinetic evolution of the self-consistent flow around a two-dimensional obstacle with zero intrinsic magnetic field. Plasma outlfow from the obstacle is used to establish a boundary layer between the incoming solar wind and the outgoing plasma. Because the self-consistent flow solution, a velocity shear is naturally set up at this interface, and since the magnetic field for these simulations is transverse to this flow, the Kelvin-Helmholtz (K-H) instability can be excited at low-velocity shear. Simulations demonstrate the existence of the instability even near the subsolar location, which normally is thought to be stable to this instability. The apparent reason for this result is the overall time dependence at the boundary layer, which gives rise to a Rayleigh-Taylor like instability which provides seed perturbations for the K-H instability. These results are directly applicable to Venus, comets, artificial plasma releases, and laser target experiments. This result has potentially important ramifications for the interpretation of observational results as well as for an estimation of the cross-field transport. The results suggest that the K-H instability may play a role in dayside processes and the Venus ionopause, and may exist within the context of more general situations, for example, the Earth's magnetopause.

Thomas, V. A.↗

Studies with Cluster upstream and downstream of the bow shock: An experimenter's perspective

Some open questions in the physics of bow shock formation, the evolution of the particle distributions from solar wind into the magnetosheath, and the acceleration of ions at the moment of the shock are summarized. A layout of the current situation is presented in view of recent theoretical developments and the new diagnostic tools provided by the Cluster mission. The transition of ions across the quasi-perpendicular bow shock and their downstream thermalization are discussed. The processes and spatial scales are found to be species dependent and are discussed for H(+), He(2+), and He(+). The theory of particle acceleration at quasi-parallel shocks are reviewed. It is shown how Cluster can study the time variable structures of the shock as predicted by hybrid simulation. It is emphasized that high time resolution measurement with simultaneous species separation is necessary for the study of the ion acceleration. Suggestions for the spacecraft separations at the bow shock are suggested.

Moebius, E.↗

The making of an Alfvenic fluctuation: The resolution of a second-order analysis

Ulysses observations of the high speed polar streams show that they are largely occupied by very large amplitude Alfvenic fluctuations accompanied by many rotational discontinuities. These fluctuations have a nearly constant magnetic intensity or amplitude, and the magnetic field direction per wave cycle sweeps only through a limited arc, much as a car wiperblade would do. Barnes and Hollweg (JGR, 79, 2302, 1974) suggested that this unusual waveform could arise from an obliquely propagating and linearly polarized Alfven wave of finite amplitude. From a second-order analysis, they showed that the existence of a particular solution with a constant amplitude but could not resolve the outcome of the homogeneous solution which consisted of fast waves. They suggested that Landau damping of these fast waves may be needed to get the observed waveform. We present a 1 1/2 D hybrid simulation which is fully nonlinear and correctly describes the ion kinetics for an initially monochromatic and linearly polarized Alfven wave propagating obliquely to the background magnetic field. The wave has a large amplitude and a wavelength so long that it can be considered dispersionless for simulation times. At early times, the second harmonic in density and in magnetic field transverse to the initial wave magnetic field are generated and have more power than other harmonics. Steepening is observed with a weak fast shock emerging, but no rotational discontinuity is left behind, and instead a constant amplitude and an arc-shaped waveform is made. The compressional component which develops after the shocks have dissipated is to zeroth order better described as a pure acoustic wave than as a fast wave. This might be explained by the relaxing of the Alfven wave to a state where its ponderomotive force vanishes so that the compressional component can travel almost independently of it.

Vasquez, Bernard J.↗

Hybrid modeling of the formation and structure of thin current sheets in the magnetotail

Hybrid simulations are used to investigate the formation of a thin current sheet inside the plasma sheet of a magnetotail-like configuration. The initial equilibrium is subjected to a driving electric field which is qualitatively similar to what would be expected from solar wind driving. As a result, a new current sheet with the thickness of approximately the ion inertial length is formed. The current density inside the current sheet region is supplied largely by the electrons. Ion acceleration in the cross-tail direction is absent as the driving electric field fails to penetrate into the equatorial region.

Hesse, Michael↗

Application of Active Controls Technology to Aircraft Ride Smoothing Systems

A critical review of past efforts in the design and testing of ride smoothing and gust alleviation systems is presented. Design trade-offs involving sensor types, choice of feedback loops, human comfort and aircraft handling-qualities criteria are discussed. Synthesis of a system designed to employ direct-lift and side-force producing surfaces is reported. Two STOL-class aircraft and an executive transport are considered. Theoretically-predicted system performance is compared with hybrid simulation and flight test data. Pilot opinion rating, pilot workload, and passenger comfort rating data for the basic and augmented aircraft are included.

Lapins, Maris↗

Investigation of the Large Scale Evolution and Topology of Coronal Mass Ejections in the Solar Wind

This investigation is concerned with the large-scale evolution and topology of coronal mass ejections (CMEs) in the solar wind. During the course of this three-year investigation, we have undertaken a number of studies that are discussed in more detail in this report. For example, we conducted an analysis of all CMEs observed by the Ulysses spacecraft during its in-ecliptic phase between 1 and 5 AU. In addition to studying the properties of the ejecta, we also analyzed the shocks that could be unambiguously associated with the fast CMEs. We also analyzed a series of 'density holes' observed in the solar wind that bear many similarities with CMEs. To complement this analysis, we conducted a series of 1-D and 2 1/2-D fluid, MHD, and hybrid simulations to address a number of specific issues related to CME evolution in the solar wind. For example, we used fluid simulations to address the interpretation of negative electron temperature-density relationships often observed within CME/cloud intervals. As part of this investigation, a number of fruitful international collaborations were forged. Finally, the results of this work were presented at nine scientific meetings and communicated in eight scientific, refereed papers.

Riley, Pete↗

Microphysics of Waves and Instabilities in the Solar Wind and their Macro Manifestations in the Corona and Interplanetary Space

Investigations of the physical processes responsible for the acceleration of the solar wind were pursued with the development of two new solar wind codes: a hybrid code and a 2-D MHD code. Hybrid simulations were performed to investigate the interaction between ions and parallel propagating low frequency ion cyclotron waves in a homogeneous plasma. In a low-beta plasma such as the solar wind plasma in the inner corona, the proton thermal speed is much smaller than the Alfven speed. Vlasov linear theory predicts that protons are not in resonance with low frequency ion cyclotron waves. However, non-linear effect makes it possible that these waves can strongly heat and accelerate protons. This study has important implications for study of the corona and the solar wind. Low frequency ion cyclotron waves or Alfven waves are commonly observed in the solar wind. Until now, it is believed that these waves are not able to heat the solar wind plasma unless some cascading processes transfer the energy of these waves to high frequency part. However, this study shows that these waves may directly heat and accelerate protons non-linearly. This process may play an important role in the coronal heating and the solar wind acceleration, at least in some parameter space.

Habbal, Shadia R.↗

Hybrid MD-Nernst Planck Model of Alpha-hemolysin Conductance Properties

Motivated by experiments in which an applied electric field translocates polynucleotides through an alpha-hemolysin protein channel causing ionic current transient blockade, a hybrid simulation model is proposed to predict the conductance properties of the open channel. Time scales corresponding to ion permeation processes are reached using the Poisson-Nemst-Planck (PNP) electro-diffusion model in which both solvent and local ion concentrations are represented as a continuum. The diffusion coefficients of the ions (K(+) and Cl(-)) input in the PNP model are, however, calculated from all-atom molecular dynamics (MD). In the MD simulations, a reduced representation of the channel is used. The channel is solvated in a 1 M KCI solution, and an external electric field is applied. The pore specific diffusion coefficients for both ionic species are reduced 5-7 times in comparison to bulk values. Significant statistical variations (17-45%) of the pore-ions diffusivities are observed. Within the statistics, the ionic diffusivities remain invariable for a range of external applied voltages between 30 and 240mV. In the 2D-PNP calculations, the pore stem is approximated by a smooth cylinder of radius approx. 9A with two constriction blocks where the radius is reduced to approx. 6A. The electrostatic potential includes the contribution from the atomistic charges. The MD-PNP model shows that the atomic charges are responsible for the rectifying behaviour and for the slight anion selectivity of the a-hemolysin pore. Independent of the hierarchy between the anion and cation diffusivities, the anionic contribution to the total ionic current will dominate. The predictions of the MD-PNP model are in good agreement with experimental data and give confidence in the present approach of bridging time scales by combining a microscopic and macroscopic model.

Cozmuta, Ioana↗

Energy Deposition Processes in Titan's Upper Atmosphere

Most of Titan's atmospheric organic and nitrogen chemistry, aerosol formation, and atmospheric loss are driven from external energy sources such as Solar UV, Saturn's magnetosphere, solar wind and galactic cosmic rays. The Solar UV tends to dominate the energy input at lower altitudes of approximately 1100 km but which can extend down to approximately 400 km, while the plasma interaction from Saturn's magnetosphere, Saturn's magnetosheath or solar wind are more important at higher altitudes of approximately 1400 km, but the heavy ion plasma [O(+)] of approximately 2 keV and energetic ions [H(+)] of approximately 30 keV or higher from Saturn's magnetosphere can penetrate below 950km. Cosmic rays with energies of greater than 1 GeV can penetrate much deeper into Titan's atmosphere with most of its energy deposited at approximately 100 km altitude. The haze layer tends to dominate between 100 km and 300 km. The induced magnetic field from Titan's interaction with the external plasma can be very complex and will tend to channel the flow of energy into Titan's upper atmosphere. Cassini observations combined with advanced hybrid simulations of the plasma interaction with Titan's upper atmosphere show significant changes in the character of the interaction with Saturn local time at Titan's orbit where the magnetosphere displays large and systematic changes with local time. The external solar wind can also drive sub-storms within the magnetosphere which can then modify the magnetospheric interaction with Titan. Another important parameter is solar zenith angle (SZA) with respect to the co-rotation direction of the magnetospheric flow. Titan's interaction can contribute to atmospheric loss via pickup ion loss, scavenging of Titan's ionospheric plasma, loss of ionospheric plasma down its induced magnetotail via an ionospheric wind, and non-thermal loss of the atmosphere via heating and sputtering induced by the bombardment of magnetospheric keV ions and electrons. This energy input evidently drives the large positive and negative ions observed below approximately 1100 km altitude with ion masses exceeding 10,000 daltons. We refer to these ions as seed particles for the aerosols observed below 300 km altitude. These seed particles can be formed, for example, from the polymerization of acetylene (C2H2) and benzene (C6H6) molecules in Titan's upper atmosphere to form polycyclic aromatic hydrocarbons (PAH) and/or fullerenes (C60). In the case of fullerenes, which are hollow spherical carbon shells, magnetospheric keV [O(+)] ions can become trapped inside the fullerenes and eventually find themselves inside the aerosols as free oxygen. The aerosols are then expected to fall to Titan's surface as polymerized hydrocarbons with trapped free oxygen where unknown surface chemistry can take place.

Sittler, Edward C., Jr.↗

Titan's Interaction with Saturn's Magnetosphere as Viewed by T9 and T18 Flybys

New results are presented of Cassini's T9 flyby with complementary observations from T18. Based on Cassini Plasma Spectrometer (CAPS) and Cassini Magnetometer (MAG) data, compositional evidence shows the upstream flow for both T9 and T18 is composed of light ions (H(+) and H2(+), with external pressures 30 times lower than that for the earlier TA flyby where heavy ions dominated the magnetospheric plasma. When describing the plasma heating and sputtering of Titan's atmosphere, T9 and T18 can be considered interactions of low magnetospheric energy input. Perpendicular temperature over parallel temperature anisotropy varied between 1 and 3 for upstream flow. The upstream flow is perpendicular to B, indicative of local picked up ions from Titan"s H and H2 coronae extending to Titan's Hill sphere radius. Beyond this distance the corona form a neutral torus that surrounds Saturn. The T9 flyby unexpectedly resulted in observation of two "wake" crossings referred to as Events 1 and 2. Event 2 was evidently caused by draped magnetosphere field lines which are scavenging pickup ions from Titan's induced magnetopause boundary with outward flux 3.0e6 ions/cm2/s. The composition of this out flow is dominated by H2(+) and H(+) ions. Ionospheric flow away from Titan with ion flux 1.2e7 ion/cm2 /s is observed for Event 1. In between Events 1 and 2 are high energy field aligned flows of magnetosphere protons. T18 observations are much closer to Titan than T9, allowing one to probe this type of interaction down to altitudes 950 km. Comparisons with previously reported hybrid simulations are made.

Sittler, Edward↗

Mercury's Dynamic Magnetic Tail

The Mariner 10 and MESSENGER flybys of Mercury have revealed a magnetosphere that is likely the most responsive to upstream interplanetary conditions of any in the solar system. The source of the great dynamic variability observed during these brief passages is due to Mercury's proximity to the Sun and the inverse proportionality between reconnection rate and solar wind Alfven Mach number. However, this planet's lack of an ionosphere and its small physical dimensions also contribute to Mercury's very brief Dungey cycle, approx. 2 min, which governs the time scale for internal plasma circulation. Current observations and understanding of the structure and dynamics of Mercury's magnetotail are summarized and discussed. Special emphasis will be placed upon such questions as: 1) How much access does the solar wind have to this small magnetosphere as a function of upstream conditions? 2) What roles do heavy planetary ions play? 3) Do Earth-like substorms take place at Mercury? 4) How does Mercury's tail respond to extreme solar wind events such coronal mass ejections? Prospects for progress due to advances in the global magnetohydrodynamic and hybrid simulation modeling and the measurements to be taken by MESSENGER after it enters Mercury orbit on March 18, 2011 will be discussed.

Slavin, James A.↗

Saturn's Magnetospheric Interaction with Titan as Defined by Cassini Encounters T9 and T18: New Results

We present new results of Cassini s T9 flyby with complementary observations from T18. Based on Cassini plasma spectrometer (CAPS) and Cassini magnetometer (MAG), compositional evidence shows the upstream flow for both T9 and T18 appears composed of light ions (H+ and H2+), with external pressures approx.30 times lower than that for the earlier TA flyby where heavy ions dominated the magnetospheric plasma. When describing the plasma heating and sputtering of Titan s atmosphere, T9 and T18 can be considered interactions of low magnetospheric energy input. On the other hand, T5, when heavy ion fluxes are observed to be higher than typical (i.e., TA), represents the limiting case of high magnetospheric energy input to Titan s upper atmosphere. Beyond this distance the corona forms a neutral torus that surrounds Saturn. The T9 flyby unexpectedly resulted in observation of two wake crossings referred to as Events 1 and 2. Event 2 was evidently caused by draped magnetosphere field lines, which are scavenging pickup ions from Titan s induced magnetopause boundary with outward flux approx.2 x 10(exp 6) ions/sq cm/s. The composition of this out flow is dominated by H2+ and H+ ions. Ionospheric flow away from Titan with ion flux approx7 x 10(exp 6) ion/sq cm/s is observed for Event 1. In between Events 1 and 2 are high energy field aligned flows of magnetosphere protons that may have been accelerated by the convective electric field across Titan s topside ionosphere. T18 observations are much closer to Titan than T9, allowing one to probe this type of interaction down to altitudes approx.950 km. Comparisons with previously reported hybrid simulations are made.

Sittler, E. C., Jr.↗

Whistler Waves Associated with Weak Interplanetary Shocks

We analyze the properties of 98 weak interplanetary shocks measured by the dual STEREO spacecraft over approximately 3 years during the past solar minimum. We study the occurrence of whistler waves associated with these shocks, which on average are high beta shocks (0.2 < Beta < 10). We have compared the waves properties upstream and downstream of the shocks. In the upstream region the waves are mainly circularly polarized, and in most of the cases (approx. 75%) they propagate almost parallel to the ambient magnetic field (<30 deg.). In contrast, the propagation angle with respect to the shock normal varies in a broad range of values (20 deg. to 90 deg.), suggesting that they are not phase standing. We find that the whistler waves can extend up to 100,000 km in the upstream region but in most cases (88%) are contained in a distance within 30,000 km from the shock. This corresponds to a larger region with upstream whistlers associated with IP shocks than previously reported in the literature. The maximum amplitudes of the waves are observed next to the shock interface, and they decrease as the distance to the shock increases. In most cases the wave propagation direction becomes more aligned with the magnetic field as the distance to the shock increases. These two facts suggest that most of the waves in the upstream region are Landau damping as they move away from the shock. From the analysis we also conclude that it is likely that the generation mechanism of the upstream whistler waves is taking place at the shock interface. In the downstream region, the waves are irregularly polarized, and the fluctuations are very compressive; that is, the compressive component of the wave clearly dominates over the transverse one. The majority of waves in the downstream region (95%) propagate at oblique angles with respect to the ambient magnetic field (>60 deg.). The wave propagation with respect to the shock-normal direction has no preferred direction and varies similarly to the upstream case. It is possible that downstream fluctuations are generated by ion relaxation as suggested in previous hybrid simulation shocks.

Whistler waves↗

Magnetosheath Filamentary Structures Formed by Ion Acceleration at the Quasi-Parallel Bow Shock

Results from 2.5-D electromagnetic hybrid simulations show the formation of field-aligned, filamentary plasma structures in the magnetosheath. They begin at the quasi-parallel bow shock and extend far into the magnetosheath. These structures exhibit anticorrelated, spatial oscillations in plasma density and ion temperature. Closer to the bow shock, magnetic field variations associated with density and temperature oscillations may also be present. Magnetosheath filamentary structures (MFS) form primarily in the quasi-parallel sheath; however, they may extend to the quasi-perpendicular magnetosheath. They occur over a wide range of solar wind Alfvénic Mach numbers and interplanetary magnetic field directions. At lower Mach numbers with lower levels of magnetosheath turbulence, MFS remain highly coherent over large distances. At higher Mach numbers, magnetosheath turbulence decreases the level of coherence. Magnetosheath filamentary structures result from localized ion acceleration at the quasi-parallel bow shock and the injection of energetic ions into the magnetosheath. The localized nature of ion acceleration is tied to the generation of fast magnetosonic waves at and upstream of the quasi-parallel shock. The increased pressure in flux tubes containing the shock accelerated ions results in the depletion of the thermal plasma in these flux tubes and the enhancement of density in flux tubes void of energetic ions. This results in the observed anticorrelation between ion temperature and plasma density.

magnetosheath↗

On Lunar Exospheric Column Densities and Solar Wind Access Beyond the Terminator from ROSAT Soft X-Ray Observations of Solar Wind Charge Exchange

We analyze the Rontgen satellite (ROSAT) position sensitive proportional counter soft X-ray image of the Moon taken on 29 June 1990 by examining the radial profile of the surface brightness in three wedges: two 19 deg wedges (one north and one south) 13-32 deg off the terminator toward the dark side and one wedge 38 deg wide centered on the antisolar direction. The radial profiles of both the north and the south wedges show significant limb brightening that is absent in the 38 deg wide antisolar wedge. An analysis of the soft X-ray intensity increase associated with the limb brightening shows that its magnitude is consistent with that expected due to solar wind charge exchange (SWCX) with the tenuous lunar atmosphere based on lunar exospheric models and hybrid simulation results of solar wind access beyond the terminator. Soft X-ray imaging thus can independently infer the total lunar limb column density including all species, a property that before now has not been measured, and provide a large-scale picture of the solar wind-lunar interaction. Because the SWCX signal appears to be dominated by exospheric species arising from solar wind implantation, this technique can also determine how the exosphere varies with solar wind conditions. Now, along with Mars, Venus, and Earth, the Moon represents another solar system body at which SWCX has been observed.

Exospheric↗

Control of Future Air Traffic Systems via Complexity Bound Management

The complexity of the present system for managing air traffic has led to "discreteness" in approaches to creating new concepts: new concepts are created as point designs, based on experience, expertise, and creativity of the proposer. Discrete point designs may be highly successful but they are difficult to substantiate in the face of equally strong substantiation of competing concepts, as well as the state of the art in concept evaluation via simulations. Hybrid concepts may present a compromise - the golden middle. Yet a hybrid of sometimes in principle incompatible concepts forms another point design that faces the challenge of substantiation and validation. We are faced with the need to re-design the air transportation system ab initio. This is a daunting task, especially considering the problem of transitioning from the present system to any fundamentally new system. However, design from scratch is also an opportunity to reconsider approaches to new concept development. In this position paper we propose an approach, Optimized Parametric Functional Design, for systematic development of concepts for management and control of airspace systems, based on optimization formulations in terms of required system functions and states. This reasoning framework, realizable in the context of ab initio system design, offers an approach to deriving substantiated airspace management and control concepts. With growing computational power, we hope that the approach will also yield a methodology for actual dynamic control of airspace

Alexandrov, Natalia↗