Engineering topics

Gombosi, Tamas I.

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

At least 19 records

NASA NTRS2016

Extended Magnetohydrodynamics with Embedded Particle-in-Cell Simulation of Ganymede's Magnetosphere

We have recently developed a new modeling capability to embed the implicit particle-in-cell (PIC) model iPIC3D into the Block-Adaptive-Tree-Solarwind-Roe-Upwind-Scheme magnetohydrodynamic (MHD) model. The MHD with embedded PIC domains (MHO-EPIC) algorithm Is a two-way coupled kinetic-fluid model. As one of the very first applications of the MHD-EPIC algorithm, we simulate the Interaction between Jupiter's magnetospherlc plasma and Ganymede's magnetosphere. We compare the MHO-EPIC simulations with pure Hall MHD simulations and compare both model results with Galileo observations to assess the Importance of kinetic effects In controlling the configuration and dynamics of Ganymede's magnetosphere. We find that the Hall MHD and MHO-EPIC solutions are qualitatively similar, but there are significant quantitative differences. In particular. the density and pressure inside the magnetosphere show different distributions. For our baseline grid resolution the PIC solution is more dynamic than the Hall MHD simulation and it compares significantly better with the Galileo magnetic measurements than the Hall MHD solution. The power spectra of the observed and simulated magnetic field fluctuations agree extremely well for the MHD-EPIC model. The MHO-EPIC simulation also produced a few flux transfer events (FTEs) that have magnetic signatures very similar to an observed event. The simulation shows that the FTEs often exhibit complex 3-0 structures with their orientations changing substantially between the equatorial plane and the Galileo trajectory, which explains the magnetic signatures observed during the magnetopause crossings. The computational cost of the MHO-EPIC simulation was only about 4 times more than that of the Hall MHD simulation.

Toth, Gabor↗

NASA NTRS2011

Spatial and Temporal Signatures of Flux Transfer Events in Global Simulations of Magnetopause Dynamics

Flux transfer events (FTEs) were originally identified by transient bipolar variations of the magnetic field component normal to the nominal magnetopause centered on enhancements in the total magnetic field strength. Recent Cluster and THEMIS multi-point measurements provided a wide range of signatures that are interpreted as evidence for FTE passage (e.g., crater FTE's, traveling magnetic erosion regions). We use the global magnetohydrodynamic (MHD) code BATS-R-US developed at the University of Michigan to model the global three-dimensional structure and temporal evolution of FTEs during multi-spacecraft magnetopause crossing events. Comparison of observed and simulated signatures and sensitivity analysis of the results to the probe location will be presented. We will demonstrate a variety of observable signatures in magnetic field profile that depend on space probe location with respect to the FTE passage. The global structure of FTEs will be illustrated using advanced visualization tools developed at the Community Coordinated Modeling Center

Kuznetsova, Maria M.↗

NASA NTRS2010

Adaptive Numerical Algorithms in Space Weather Modeling

Space weather describes the various processes in the Sun-Earth system that present danger to human health and technology. The goal of space weather forecasting is to provide an opportunity to mitigate these negative effects. Physics-based space weather modeling is characterized by disparate temporal and spatial scales as well as by different physics in different domains. A multi-physics system can be modeled by a software framework comprising of several components. Each component corresponds to a physics domain, and each component is represented by one or more numerical models. The publicly available Space Weather Modeling Framework (SWMF) can execute and couple together several components distributed over a parallel machine in a flexible and efficient manner. The framework also allows resolving disparate spatial and temporal scales with independent spatial and temporal discretizations in the various models. Several of the computationally most expensive domains of the framework are modeled by the Block-Adaptive Tree Solar wind Roe Upwind Scheme (BATS-R-US) code that can solve various forms of the magnetohydrodynamics (MHD) equations, including Hall, semi-relativistic, multi-species and multi-fluid MHD, anisotropic pressure, radiative transport and heat conduction. Modeling disparate scales within BATS-R-US is achieved by a block-adaptive mesh both in Cartesian and generalized coordinates. Most recently we have created a new core for BATS-R-US: the Block-Adaptive Tree Library (BATL) that provides a general toolkit for creating, load balancing and message passing in a 1, 2 or 3 dimensional block-adaptive grid. We describe the algorithms of BATL and demonstrate its efficiency and scaling properties for various problems. BATS-R-US uses several time-integration schemes to address multiple time-scales: explicit time stepping with fixed or local time steps, partially steady-state evolution, point-implicit, semi-implicit, explicit/implicit, and fully implicit numerical schemes. Depending on the application, we find that different time stepping methods are optimal. Several of the time integration schemes exploit the block-based granularity of the grid structure. The framework and the adaptive algorithms enable physics based space weather modeling and even forecasting.

Toth, Gabor↗

NASA NTRS2003

High Performance Parallel Methods for Space Weather Simulations

This is the final report of our NASA AISRP grant entitled 'High Performance Parallel Methods for Space Weather Simulations'. The main thrust of the proposal was to achieve significant progress towards new high-performance methods which would greatly accelerate global MHD simulations and eventually make it possible to develop first-principles based space weather simulations which run much faster than real time. We are pleased to report that with the help of this award we made major progress in this direction and developed the first parallel implicit global MHD code with adaptive mesh refinement. The main limitation of all earlier global space physics MHD codes was the explicit time stepping algorithm. Explicit time steps are limited by the Courant-Friedrichs-Lewy (CFL) condition, which essentially ensures that no information travels more than a cell size during a time step. This condition represents a non-linear penalty for highly resolved calculations, since finer grid resolution (and consequently smaller computational cells) not only results in more computational cells, but also in smaller time steps.

Hunter, Paul↗

NASA NTRS2000

Two-Species, 3D, MHD Simulation of Europa's Interaction with Jupiter's Magnetosphere

The interaction of Europa with the Jovian a magnetosphere has been studied by using a two species in ideal magnetohydrodynamic (MHD) numerical model. This model considers the upstream plasma in the Jovian magnetosphere and the molecular oxygen ions in the ionosphere of Europa, separately. We present results a from simulation studies, which take into account impact ionization, recombination, and the effect of a possible induced dipole magnetic field of Europa. The total mass loading of the magnetospheric flow and the ionization frequency used in the model are consistent with the estimates of Europa's ionosphere and atmosphere. The multi-species MHD equations are solved by using a finite volume, high-order, Godunov-type method on an adoptively refined unstructured grid, which allows a detailed modeling of the region near Europa's surface, while still resolving both the upstream region and the satellite's wake. We have paid special attention to the wake of Europa, in order to be able to make comparisons with the Galileo's E4 flyby observations, as well as other model calculations. The calculated escape flux of a O2+ down the tail was found to be about 5.6 x 10(exp 25) s(sup -1).

Liu, Yifan↗

NASA NTRS2000

Plasma Flow Past Cometary and Planetary Satellite Atmospheres

The tenuous atmospheres and ionospheres of comets and outer planet satellites share many common properties and features. Such similarities include a strong interaction with their outer radiation, fields and particles environs. For comets the interaction is with the magnetized solar wind plasma, whereas for satellites the interaction is with the strongly magnetized and corotating planetary magnetospheric plasma. For this reason there are many common or analogous physical regimes, and many of the same modeling techniques are used to interpret remote sensing and in situ measurements in order to study the important underlying physical phenomena responsible for their appearances. We present here a review of various modeling approaches which are used to elucidate the basic properties and processes shaping the energetics and dynamics of these systems which are similar in many respects.

Combi, Michael R.↗

NASA NTRS1996

Three-Dimensional Multiscale MHD Model of Cometary Plasma Environments

First results of a three-dimensional multiscale MHD model of the interaction of an expanding cometary atmosphere with the magnetized solar wind are presented. The model starts with a supersonic and super-Alfvenic solar wind far upstream of the comet (25 Gm upstream of the nucleus) with arbitrary interplanetary magnetic field orientation. The solar wind is continuously mass loaded with cometary ions originating from a 10-km size nucleus. The effects of photoionization, electron impact ionization, recombination, and ion-neutral frictional drag are taken into account in the model. The governing equations are solved on an adaptively refined unstructured Cartesian grid using our new multiscale upwind scalar conservation laws-type numerical technique (MUSCL). We have named this the multiscale adaptive upwind scheme for MHD (MAUS-MHD). The combination of the adaptive refinement with the MUSCL-scheme allows the entire cometary atmosphere to be modeled, while still resolving both the shock and the diamagnetic cavity of the comet. The main findings are the following: (1) Mass loading decelerates the solar wind flow upstream of the weak cometary shock wave (M approximately equals 2, M(sub A) approximately equals 2), which forms at a subsolar standoff distance of about 0.35 Gm. (2) A cometary plasma cavity is formed at around 3 x 10(exp 3) km from the nucleus. Inside this cavity the plasma expands outward due to the frictional interaction between ions and neutrals. On the nightside this plasma cavity considerably narrows and a relatively fast and dense cometary plasma beam is ejected into the tail. (3) Inside the plasma cavity a teardrop-shaped inner shock is formed, which is terminated by a Mach disk on the nightside. Only the region inside the inner shock is the 'true' diamagnetic cavity. (4) The model predicts four distinct current systems in the inner coma: the density peak current, the cavity boundary current, the inner shock current, and finally the cross-tail current. (5) The calculated plasma parameters (magnetic field, plasma density, speed, and temperature) are in very good agreement with published Giotto observations.

Gombosi, Tamas I.↗

NASA NTRS1996

A New Axisymmetric MHD Model of the Interaction of the Solar Wind with Venus

A new two-dimensional axisymmetric MHD model is used to study the interaction of the solar wind with Venus under conditions where the interplanetary field is approximately aligned with the solar wind velocity. This numerical model solves the MHD transport equations for density, velocity, pressure, and magnetic field on an adaptively refined, unstructured grid system. This use of an adaptive grid allows high spatial resolution in regions of large density/velocity gradients and yet can be run on a workstation. The actual grid sizes vary from about 0.06 R(sub v) near the bowshock to 2 R(sub v) in the unperturbed solar wind. The results of the calculations are compared with observed magnetic field values obtained from the magnetometer on the Pioneer Venus Orbiter, at a time when the angle between the solar wind velocity vector and the interplanetary magnetic field (IMF) was only 7.6 deg. Good qualitative agreement between the observed and calculated field behavior is found. The overall results suggest that the induced magnetotail disappears when the IMF is radial for an extended time period and implies that it weakens when the field rotated through a near-radial orientation.

DeZeeuw, Darren L.↗

NASA NTRS1994

Axisymmetric modeling of cometary mass loading on an adaptively refined grid: MHD results

The first results of an axisymmetric magnetohydrodynamic (MHD) model of the interaction of an expanding cometary atmosphere with the solar wind are presented. The model assumes that far upstream the plasma flow lines are parallel to the magnetic field vector. The effects of mass loading and ion-neutral friction are taken into account by the governing equations, whcih are solved on an adaptively refined unstructured grid using a Monotone Upstream Centered Schemes for Conservative Laws (MUSCL)-type numerical technique. The combination of the adaptive refinement with the MUSCL-scheme allows the entire cometary atmosphere to be modeled, while still resolving both the shock and the near nucleus of the comet. The main findingsare the following: (1) A shock is formed approximately = 0.45 Mkm upstream of the comet (its location is controlled by the sonic and Alfvenic Mach numbers of the ambient solar wind flow and by the cometary mass addition rate). (2) A contact surface is formed approximately = 5,600 km upstream of the nucleus separating an outward expanding cometary ionosphere from the nearly stagnating solar wind flow. The location of the contact surface is controlled by the upstream flow conditions, the mass loading rate and the ion-neutral drag. The contact surface is also the boundary of the diamagnetic cavity. (3) A closed inner shock terminates the supersonic expansion of the cometary ionosphere. This inner shock is closer to the nucleus on dayside than on the nightside.

Gombosi, Tamas I.↗

NASA NTRS1993

Non-steady-state transport of superthermal electrons in the plasmasphere

Numerical solutions to the time-dependent kinetic equation, which describes the transport of superthermal electrons in the splasmasphere between the two conjugate ionospheres, are presented. The model calculates the distribution function as a function of time, field-aligned distance, energy, and pitch-angle. The processes of refilling, depleting, and establishing steady-state conditions of superthermal electrons in the plasmasphere are discussed.

Khazanov, George V.↗

NASA NTRS1993

Kinetic simulation of plasma flows in the inner magnetosphere

A one-dimensional hybrid particle code is used to study the interactions between upflowing thermal ions from conjugate ionospheres. The simulation model allows for multiple species, convection of plasmaspheric flux tubes, and Coulomb self-collisions which conserve momentum and energy locally. The model incorporates a variable-flux boundary condition where the flux, at the boundaries, approaches zero as the plasmasphere fills and equilibrium conditions are reached. The effects of two important processes on plasmaspheric refilling have been considered. The first includes convection of the plasmaspheric flux tube. The second is the interaction of ionospheric thermal plasma and particle injection from an external source. Particle injection seems to play an important role in the evolution of the total particle distribution on the early timescales (t less than 1 hour); however, for late timescales (t larger than 8 days) the thermal plasma from the ionosphere dominates the particle distribution.

Miller, Ronald H.↗

NASA NTRS1993

Multiple scattering of light in a spherical cometary atmosphere with an axisymmetric dust jet. II - Image simulation

A numerical solution for the multiple light scattering in spherical axisymmetric geometry is applied to the simulation of images of a coma as it would appear to a near-flying satellite such as Giotto. The appearance of symmetric comas and dust jets is examined in detail; the nucleus visibility is studied; the effect of forward scattering is considered; and single and multiple scattering effects are quantified. Attention is given to simulated images of a coma with a hollow cone of dust, as predicted by dust-gas hydrodynamic modeling. The cone's appearance is very similar to the northern area of activity on Comet Halley, observed by the Giotto HMC.

Chick, Kenneth M.↗

NASA NTRS1992

Analysis of the ionosphere-plasmasphere transport of superthermal electrons. I - Transport in the plasmasphere

Analytical solutions are developed for the kinetic equation which describes the transport of superthermal electrons in the terrestrial plasmasphere, together with a relationship which makes it possible to calculate the transparency of the plasmasphere to these electrons. In addition, analytic expressions are presented for the heating rate of the thermal plasma due to the passage of these superthermal electrons through the plasmasphere.

Khazanov, George V.↗

NASA NTRS1992

Multiple scattering of light in a spherical cometary atmosphere with an axisymmetric dust jet

A numerical solution has been developed for the anisotropic multiple scattering of light in a spherical shell comet atmosphere. The code has been run for a spherically symmetric coma distribution, benchmarked against past studies, and then run for the conditions of an axisymmetric dust jet at the subsolar point of the comet. The radiant flux impinging on the nucleus surface and the mean intensity of light throughout the coma were investigated.

Chick, Kenneth M.↗

NASA NTRS1992

Mass loading at Titan and comets

Present understanding of the cometary mass loading process (the addition of newborn ions to high speed magnetized plasma flows) is explored and the cometary lesson is applied to the environment of Titan. Mass loading is the dominant physical process controlling the interaction of the solar wind with comets, and it might play an important role in the interaction of Titan's neutral torus and exosphere with Saturn's magnetosphere. Freshly born ions follow cycloidal trajectories, resulting from the superposition of gyration and E/B drift. In the plasma frame, the resulting velocity space distribution is a ring beam distribution at comets and a ring distribution at Titan. This difference in free energy is expected to generate significantly different levels of low frequency broad band electromagnetic turbulence and, thus, result in important differences in the implanted ion distributions. It is predicted that mass loading creates a hot, equatorially trapped plasma population at Titan.

Gombosi, Tamas I.↗

NASA NTRS1992

Particle acceleration at comets

This paper compares calculated and measured energy spectra of implanted H(+) and O(+) ions on the assumption that the pickup geometry is quasi-parallel and about 1 percent of the waves generated by the cometary pickup process propagates backward (towards the comet). The model provides a good description of the implanted O(+) and H(+) energy distribution near the pickup energies.

Gombosi, Tamas I.↗

NASA NTRS1991

An analytic solution to the steady-state double adiabatic equations

A set of 20 generalized moment-transport equations applicable to low-beta (low random energy density/magnetic energy density) plasmas of thermal origin in rotating planetary magnetospheres. An analytic solution is obtained for a set of simplified equations characterizing the steady-state transport of collisionless plasma; the novel element of this analytic solution is a closed-form solution for the parallel-flow velocity variation along magnetic field lines. It is shown that only T(perpendicular)/B remains constant along the field line, while density, parallel temperature, and parallel Mach number vary.

Gombosi, Tamas I.↗

NASA NTRS1991

The directional dependence of magnetic fluctuations generated by cometary ion pickup

The properties of low frequency magnetic fluctuations generated by cometary ion pickup are examined by means of one-dimensional hybrid simulations, in which newborn ions are created at a constant rate. The helicity and direction of propagation of magnetic fluctuations are investigated for various cometary ion injection angles, alpha, relative to the solar wind magnetic field. The parameter eta represents the relative contribution of wave energy density propagating in the direction away from the comet, parallel to the beam. For small (quasi-parallel) injection angles, alpha is about 0 deg and eta is of order unity, while for larger (quasi-perpendicular) angles, alpha is about 90 deg and eta is about 0.5. At intermediate angles, alpha is about 60 deg, eta can vary between 0 and 1, depending on the wave number. The wave properties are consistent with the instabilities expected from linear theory.

Miller, Ronald H.↗