Engineering Papers⌕ Search

Engineering topics

Gombosi, T. I.

Publications and source records attributed to Gombosi, T. I..

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Magnetotail Current Sheet Thinning and Magnetic Reconnection Dynamics in Global Modeling of Substorms

Magnetotail current sheet thinning and magnetic reconnection are key elements of magnetospheric substorms. We utilized the global MHD model BATS-R-US with Adaptive Mesh Refinement developed at the University of Michigan to investigate the formation and dynamic evolution of the magnetotail thin current sheet. The BATSRUS adaptive grid structure allows resolving magnetotail regions with increased current density up to ion kinetic scales. We investigated dynamics of magnetotail current sheet thinning in response to southwards IMF turning. Gradual slow current sheet thinning during the early growth phase become exponentially fast during the last few minutes prior to nightside reconnection onset. The later stage of current sheet thinning is accompanied by earthward flows and rapid suppression of normal magnetic field component $B-z$. Current sheet thinning set the stage for near-earth magnetic reconnection. In collisionless magnetospheric plasma, the primary mechanism controlling the dissipation in the vicinity of the reconnection site is non-gyrotropic effects with spatial scales comparable with the particle Larmor radius. One of the major challenges in global MHD modeling of the magnetotail magnetic reconnection is to reproduce fast reconnection rates typically observed in smallscale kinetic simulations. Bursts of fast reconnection cause fast magnetic field reconfiguration typical for magnetospheric substorms. To incorporate nongyritropic effects in diffusion regions we developed an algorithm to search for magnetotail reconnection sites, specifically where the magnetic field components perpendicular to the local current direction approaches zero and form an X-type configuration. Spatial scales of the diffusion region and magnitude of the reconnection electric field are calculated self-consistently using MHD plasma and field parameters in the vicinity of the reconnection site. The location of the reconnection sites and spatial scales of the diffusion region are updated during the simulations. Such an approach allows quantifying the interaction between large-scale global magnetospheric dynamics and microphysical processes in diffusion regions localized near reconnection sites. To clarify the role of smallscale non-MHD effects in diffusion region on the global magnetospheric dynamic and to test different models of dissipation we perform simulations with steady southward IMF driving of the magnetosphere.

Kuznetsova, M. M.↗
NASA NTRS2006

Magnetic effects change our view of the Heliosheath

In this paper we review these recent results, and present an additional simulation with constant neutral atom background. In this case the jet is still present but with reduced intensity. Further study, e.g., including neutrals and the tilt of the solar rotation from the magnetic axis, is required before we can definitively address how the heliosheath behaves. Already we can say that this region presents remarkable dynamics, with turbulent flows, indicating that the heliosheath might be very different from what we previously thought.

turbulent flows↗
NASA NTRS2004

On The Effect of Electron Collisions in the Excitation of Cometary HCN

The electron-HCN collision rate for the excitation of rotational transitions of the HCN molecule is evaluated in comets C/1995 01 (Hale-Bopp) and C/1996 B2 (Hyakutake). Based on theoretical models of the cometary atmosphere, we show that collisions with electrons can provide a significant excitation mechanism for rotational transitions in the HCN molecule. Computed values of the cross section sigma(sub e-HCN) can be as high as 1.3 x cm2, more than 2 orders of magnitude greater than the commonly assumed HCN-H2O cross section. For the ground rotational transitions of HCN, the electron-HCN collision rate is found to exceed the HCN-H2O collision rate at distances greater than 3000 km from the cometary nucleus of Hale-Bopp and 1000 km from that of Hyakutake. Collisional excitation processes dominate over radiative excitation processes up to a distance of 160,000 km from the cometary nucleus of Hale-Bopp and 50,000 km from that of Hyakutake. Excitation models that neglect electron collisions can underestimate the HCN gas production rates by as much as a factor of 2.

Lovell, Amy J.↗
NASA NTRS2001

Io's Magnetospheric Interaction: An MHD Model with Day-Night Asymmetry

In this paper we present the results of all improved three-dimensional MHD model for Io's interaction with Jupiter's magnetosphere. We have included the day-night asymmetry into the spatial distribution of our mass-loading, which allowed us to reproduce several smaller features or the Galileo December 1995 data set. The calculation is performed using our newly modified description of the pick-up processes that accounts for the effects of the corotational electric field existing in the Jovian magnetosphere. This change in the formulation of the source terms for the MHD equations resulted in significant improvements in the comparison with the Galileo measurements. We briefly discuss the limitations of our model and possible future improvements.

Kabin, K.↗
NASA NTRS2000

Global MHD Simulations of Space Plasma Environments: Heliosphere, Comets, Magnetospheres of Plants and Satellites

Magnetohydrodynamics (MHD) provides an approximate description of a great variety of processes in space physics. Accurate numerical solutions of the MHD equations are still a challenge, but in the past decade a number of robust methods have appeared. Once these techniques made the direct solution of MHD equations feasible, a number of global three-dimensional models were designed and applied to many space physics objects. The range of these objects is truly astonishing, including active galactic nuclei, the heliosphere, the solar corona, and the solar wind interaction with planets, satellites, and comets. Outside the realm of space physics, MHD theory has been applied to such diverse problems as laboratory plasmas and electromagnetic casting of liquid metals. In this paper we present a broad spectrum of models of different phenomena in space science developed in the recent years at the University of Michigan. Although the physical systems addressed by these models are different, they all use the MHD equations as a unifying basis.

Kabin, K.↗
NASA NTRS1999

On Europa's Magnetospheric Interaction: A MHD Simulation of the E4 Flyby

The global three-dimensional interaction of Europa with the Jovian magnetosphere is modeled by using a complete set of ideal magnetohydrodynamic (MHD) equations. The model accounts for exospheric mass loading, ion-neutral charge exchange, recombination, and a possible intrinsic dipole magnetic field of Europa. The single-fluid MHD equations are solved by using a modem, finite volume, higher-order, Godunov-type method on an adoptively refined unstructured grid, which allows detailed modeling of the region near Europa while still resolving both the upstream region and the satellite's wake. The magnetic field and plasma density measured during Galileo's EGA flyby of December 19, 1996, are reproduced reasonably well in the simulation. We find the agreement between the data and our model particularly convincing if we assume that the plasma velocity during the EGA flyby deviated from the nominal corotation direction by approximately 20 deg. Evidence from the Galileo energetic particle detector also supports this assumption. In this case, we can fit the data using a dipole with orientation close to that of an induced dipole arising from the interaction of a hypothetical conducting subsurface layer on Europa with the periodically changing magnetic field of Jupiter. However, the magnitude of the dipole in our model is somewhat smaller (70%) than that suggested by Khurana et al. The total mass loading and ion-neutral charge exchange rates are consistent with the estimates of Europa's atmosphere and ionosphere.

Kabin, K.↗
NASA NTRS1998

Io's Plasma Environment During the Galileo Flyby: Global Three-Dimensional MHD Modeling with Adaptive Mesh Refinement

The first results for applying a three-dimensional multimedia ideal MHD model for the mass-loaded flow of Jupiter's corotating magnetospheric plasma past Io are presented. The model is able to consider simultaneously physically realistic conditions for ion mass loading, ion-neutral drag, and intrinsic magnetic field in a full global calculation without imposing artificial dissipation. Io is modeled with an extended neutral atmosphere which loads the corotating plasma torus flow with mass, momentum, and energy. The governing equations are solved using adaptive mesh refinement on an unstructured Cartesian grid using an upwind scheme for AHMED. For the work described in this paper we explored a range of models without an intrinsic magnetic field for Io. We compare our results with particle and field measurements made during the December 7, 1995, flyby of to, as published by the Galileo Orbiter experiment teams. For two extreme cases of lower boundary conditions at Io, our model can quantitatively explain the variation of density along the spacecraft trajectory and can reproduce the general appearance of the variations of magnetic field and ion pressure and temperature. The net fresh ion mass-loading rates are in the range of approximately 300-650 kg/s, and equivalent charge exchange mass-loading rates are in the range approximately 540-1150 kg/s in the vicinity of Io.

Combi, M. R.↗
NASA NTRS1995

Two-stream modeling of plasmaspheric refilling

Plasmaspheric refilling on an L = 4 flux tube was studied by using a time-dependent, hydrodynamic plasmaspheric flow model in which the ion streams from the two hemispheres are treated as distinct fluids. In the model the continuity, momentum, and energy equations of a two-ion (O(+) and H(+)), quasi-neutral, currentless plasma are solved along a closed geomagnetic field line; diffusive equilibrium is not assumed. collisions between all stream pairs and with neutral species are included. The model includes a corotating, tilted dipole magnetic field and neutral winds. Ionospheric sources and sinks are accounted for in a self-consistent manner. Electrons are assumed to be heated by photoelectrons. The model flux tube extends from a 200-km altitude in one hemisphere to a 200-km altitude in the other hemisphere. Initially, the upwelling streams pass through each other practically unimpeded. When the streams approach the boundary in the conjugate ionosphere, a shock develops there, which moves upward and dissipates slowly; at about the same time a reverse shock develops in the hemisphere of origin, which moves upward. After about 1 hour, large shocks develop in each stream near the equator; these shocks move toward the equator and downward after crossing the equator. However, these shocks are probably artificial, because counterstreaming flows occur in each H(+) fluid, which the model can only handle by creating shocks.

Guiter, S. M.↗
NASA NTRS1995

What is the source of observed annual variations in plasmaspheric density?

Plasmaspheric densities have been observed previously to be higher in December than in June, with the ratio varying between 1.5 and 3.0 and with larger variations at lower L shells. In order to search for the cause of the observed annual variations, we have modeled plasmaspheric density, using a time-dependent hydrodynamic model. On an L = 2 field line with geomagnetic longitude equal to 300 deg, the modeled plasmaspheric densities were a factor of 1.5 times higher in December than in June. The modeled December to June density ratio was found to increase slightly with L shell, in contrast to observations; this discrepancy may be due to the fact that outer plasmaspheric flux tubes are never completely full. In addition, for an L = 2 field line with geomagnetic longitude equal to 120 deg, the modeled plasmaspheric density was higher in June than in December by a factor of about 1.2. Various numerical tests were also performed in order to examine the sensitivity of plasmaspheric density to various parameters. In particular, a large vertical neutral wind was applied in order to raise the O(+) profile which had the effect of raising plasmaspheric density by a factor of 6. This in conjunction with a theoretical analysis suggests that plasmaspheric density levels are very sensitive to O(+) levels in the upper ionosphere. We conclude that annual variations in plasmaspheric density are due to similar variations in ionospheric O(+).

Guiter, S. M.↗
NASA NTRS1994

A unifying comparison of nearly scatter free transport models

Gombosi et al. (1993) recently derived a modified telegrapher's equation for charged particle transport under the influence of isotropic scattering. This equation obeys causality and disallows upstream diffusion for particles with random velocities smaller than the bulk flow velocity. The acausal diffusion equation was obtained to lowest order in the expansion of smallness prameters. The paper by Gombosi et al. (1993) prompted responses from Pauls et al. (1993) and Earl (1993). This paper is written to explain the differences between the methods, assumptions, and results of Gombosi et al. (1993), Pauls et al. (1993), and Earl (1993) and presents a new method of obtaining approximate solutions. It is shown that the assumptions used by Gombosi et al. (1993) and Pauls et al. (1993) are physically equivalent. In our solution method, the solution of the modified telegrapher's equation is obtained as the casual limit of solutions accurate to second order in the smallness parameter expansion. In order to investigate the coherent velocity, we have also developed `wavenumber eigenfunctions' which account for all the pitch angle dependence in our Boltzmann equation. Using truncation, Earl (1993) obtains approximations for the wavenumber dependence of the lowest two frequency modes, which correspond to two of the wavenumber eigenmodes. We find that a consequence of including only two wavenumber eigenmodes is that one obtains solutions which disobey causality at sufficiently short times. Furthermore, the coherent velocity of the two eigenmodes is strongly dependent on wavenumber and approaches the particle velocity in the limit of large wavenumber for both isotropic and anisotropic scattering processes. We conclude that Earl's (1993) solutions and solutions obtained using the new solution method implicitly assume weak acausality and reasonable behavior in the temporal regime, t less than 4 tau. The solutions are not strictly consistent with the behavior of the lowest two frequency modes but have similar behavior in the regime of low wavenumber.

Schwadron, N. A.↗
NASA NTRS1994

Thermodynamic effect of the ion sound instability in the ionosphere

During geomagnetic disturbances when the ring current interacts intensely with the plasmasphere, the plasma of this region undergoes a strong heating due to an ion cyclotron instability. This is followed by the transfer of heat along geomagnetic field lines from the heating region to the ionosphere. One of the results of this process is the formation of non-isothermal region (in which T(sub e) greater than 3.4 T(sub i) at ionospheric heights) caused by a rapid cooling the H(+) ions due to their resonant charge exchange with neutral hydrogen. Heat transfer from the top of the flux tube to the ionosphere is investigated using a hydrodynamic model for the ionosphere-plasmasphere coupling. Field-aligned currents, present in the topside ionosphere, are often accompanied by ion sound turbulence. The turbulence scatters electrons, increasing the total electron collision frequency through wave-particle effects. The influence of wave-particle interactions introduces an anomalous component to the total collision frequency, which modifies substantially the heat conduction coefficient of the plasma. As a result, the plasma is heated more intensely above than below this region of ion sound turbulence.

Khazanov, G. V.↗
NASA NTRS1994

Effect of magnetospheric convection on thermal plasma in the inner magnetosphere

The effects of E x B convection on the distribution of plasma parameters in the inner magnetosphere have been examined. Analytical solutions describing density distributions along convective trajectories in the equatorial plane have been found. These solutions suggest the following dependence of plasma concentration n on the magnetic field B along convective trajectories: n proportional to B(exp alpha), where the parameter alpha varies between alpha is less than or = between 4/3 and 2. The alpha = 2 case corresponds to disregarding transport parallel to B. The lower bound alpha = 4/3 describes the situation when parallel transport dominates over convective motions perpendicular to B. A solution has also been obtained describing the effect of convection on diffusive equilibrium profiles in the dayside plasmasphere. In addition, ion temperature variations due to adiabatic effects associated with plasma convection have been analyzed in detail. Convective drifts lead to ion temperature anisotropies with the value and sign of the anisotropy contingent on density and temperature variations, local time, and the location of the convective trajectory with respect to the location of the plasmapause. We have also found that convective motions lead to a substantial exchange of energy between the ionosphere and the plasmasphere through electron heat fluxes.

Khazanov, G. V.↗
NASA NTRS1993

A unified transport equation for both cosmic rays and thermal particles

We present a unified transport equation that is valid for particles of all energies if the particle mean free paths are much smaller than macroscopic fluid length scales. If restricted to particles with random speeds much greater than fluid flow speeds, this equation reduces to the previously discussed extended cosmic-ray transport equation. It is significant that this allows one to describe the acceleration of particles from thermal energies to cosmic-ray energies using one transport equation. This is in contrast to previous transport equations (the Parker equation and the extended cosmic-ray transport equation), which were restricted to fast particles. The close connection to the extended cosmic-ray transport equation is demonstrated.

Williams, L. L.↗
NASA NTRS1993

Pickup protons and water ions at Comet Halley - Comparisons with Giotto observations

The cometary ion pickup process along the sun-comet line at Comet Halley is investigated using a quasi-linear diffusion model including both pitch angle and energy diffusion, adiabatic compression, and convective motion with the solar wind flow. The model results are compared with energetic ion distributions observed by instruments on board the Giotto spacecraft. The observed power spectrum index of magnetic turbulence (gamma) is 2-2.5. The present simulation shows that when gamma was 2, the calculated proton distributions were much more isotropic than the observed ones. The numerical solutions of the quasi-linear diffusion equations show that the isotropization of the pickup ion distribution, particularly at the pickup velocity, is not complete even close to the bow shock. Given the observed turbulence level, quasi-linear theory yields pickup ion energy distributions that agree with the observed ones quite well and easily produces energetic ions with energies up to hundreds of keV.

Ye, G.↗
NASA NTRS1993

The telegraph equation in charged particle transport

We present a new derivation of the telegraph equation which modifies its coefficients. First, an infinite order partial differential equation is obtained for the velocity space solid angle-averaged phase-space distribution of particles which underwent at least a few collisions. It is shown that, in the lowest order asymptotic expansion, this equation simplifies to the well-known diffusion equation. The second-order asymptotic expansion for isotropic small-angle scattering results in a modified telegraph equation with a signal propagation speed of v(5/11) exp 1/2 instead of the usual v/3 exp 1/2. Our derivation of a modified telegraph equation follows from an expansion of the Boltzmann equation in the relevant smallness parameters and not from a truncation of an eigenfunction expansion. This equation is consistent with causality. It is shown that, under steady state conditions in a convecting plasma, the telegraph equation may be regarded as a diffusion equation with a modified transport coefficient, which describes a combination of diffusion and cosmic-ray inertia.

Gombosi, T. I.↗
NASA NTRS1993

Plasma hydrodynamics in view of quasilinear effects

We have obtained expressions for the moments of the quasi-linear integral of anisotropic plasma collisions. Resonant and adiabatic interactions of particles with an arbitrary wave mode are considered. Expressions for the moments describing the interactions of ion-cyclotron waves with ring current and magnetospheric helium ions have been derived as an example, using the obtained results.

Gamaiunov, K. V.↗