
We studied the role of electron physics in 3D twofluid 10moment simulation
of the Ganymede's magnetosphere. The model captures nonideal physics like the
Hall effect, the electron inertia, and anisotropic, nongyrotropic pressure
effects. A series of analyses were carried out: 1) The resulting magnetic field
topology and electron and ion convection patterns were investigated. The
magnetic fields were shown to agree reasonably well with insitu measurements
by the Galileo satellite. 2) The physics of collisionless magnetic reconnection
were carefully examined in terms of the current sheet formation and
decomposition of generalized Ohm's law. The importance of pressure anisotropy
and nongyrotropy in supporting the reconnection electric field is confirmed.
3) We compared surface "brightness" morphology, represented by surface electron
and ion pressure contours, with oxygen emission observed by the Hubble Space
Telescope (HST). The correlation between the observed emission morphology and
spatial variability in electron/ion pressure was demonstrated. Potential
extension to multiion species in the context of Ganymede and other
magnetospheric systems is also discussed.

The growth and saturation of magnetic fields due to the Weibel instability
(WI) have important implications for laboratory and astrophysical plasmas, and
this has drawn significant interest recently. Since the WI can generate a large
magnetic field from no initial field, the maximum magnitudes achieved can have
significant consequences for a number of applications. Hence, an understanding
of the detailed dynamics driving the nonlinear saturation of the WI is
important. This work considers the nonlinear saturation of the WI when
counterstreaming populations of initially unmagnetized electrons are perturbed
by a magnetic field oriented perpendicular to the direction of streaming.
Previous works have found magnetic trapping to be important and connected
electron skin depth spatial scales to the nonlinear saturation of the WI. 2
Results presented in this work are consistent with these findings for a
hightemperature case. However, using a highorder continuum kinetic simulation
tool, this work demonstrates that, when the electron populations are colder, a
significant electrostatic potential develops that works with the magnetic field
to create potential wells. The electrostatic field develops due to transverse
flows induced by the WI, and in some cases is strengthened by a secondary
instability. This field plays a key role in saturation of the WI for colder
populations. The role of the electrostatic potential in Weibel instability
saturation has not been studied in detail previously.

The existence and properties of low Machnumber electrostatic collisionless
shocks are investigated with a semianalytical solution for the shock
structure. We show that the properties of the shock obtained in the
semianalytical model can be well reproduced in fully kinetic Eulerian
VlasovPoisson simulations, where the shock is generated by the decay of an
initial density discontinuity. Using this semianalytical model, we study the
effect of electrontoion temperature ratio and presence of impurities on both
the maximum shock potential and Mach number. We find that even a small amount
of impurities can influence the shock properties significantly, including the
reflected light ion fraction, which can change several orders of magnitude.
Electrostatic shocks in heavy ion plasmas reflect most of the hydrogen impurity
ions.

The kinetic study of plasma sheaths is critical, among other things, to
understand the deposition of heat on walls, the effect of sputtering, and
contamination of the plasma with detrimental impurities. The plasma sheath also
provides a boundary condition and can often have a significant global impact on
the bulk plasma. In this paper, kinetic studies of classical sheaths are
performed with the continuum code, Gkeyll, that directly solves the
VlasovPoisson/Maxwell equations. The code uses a novel version of the
finiteelement discontinuous Galerkin (DG) scheme that conserves energy in the
continuoustime limit. The electrostatic field is computed using the Poisson
equation. Ionization and scattering collisions are included, however, surface
effects are neglected. The aim of this work is to introduce the
continuumkinetic method and compare its results to those obtained from an
already established finitevolume multifluid model also implemented in Gkeyll.
Novel boundary conditions on the fluids allow the sheath to form without
specifying wall fluxes, so the fluids and fields adjust selfconsistently at
the wall. The work presented here demonstrates that the kinetic and fluid
results are in agreement for the momentum flux, showing that in certain
regimes, a multifluid model can be a useful approximation for simulating the
plasma boundary. There are differences in the electrostatic potential between
the fluid and kinetic results. Further, the direct solutions of the
distribution function presented here highlight the nonMaxwellian distribution
of electrons in the sheath, emphasizing the need for a kinetic model.

As modeling of collisionless magnetic reconnection in most space plasmas with
realistic parameters is beyond the capability of today's simulations, due to
the separation between global and kinetic length scales, it is important to
establish scaling relations in model problems so as to extrapolate to realistic
scales. Recently, large scale particleincell (PIC) simulations of island
coalescence have shown that the time averaged reconnection rate decreases with
system size, while fluid systems at such large scales in the Hall regime have
not been studied. Here we perform the complementary resistive MHD, Hall MHD and
two fluid simulations using a tenmoment model with the same geometry. In
contrast to the standard Harris sheet reconnection problem, Hall MHD is
insufficient to capture the physics of the reconnection region. Additionally,
motivated by the results of a recent set of hybrid simulations which show the
importance of ion kinetics in this geometry, we evaluate the efficacy of the
tenmoment model in reproducing such results.

Idealized merging argon plasma jets are simulated in 2D using both gas
dynamic and MHD models. Results indicate that peak pressures of several hundred
kilobar can be achieved for high Mach number jets. Including a simple optically
thin Brehmstrahlung radiation model and plasma targets shows that extremely
high densities and magnetic fields can be achieved during jet merging on the
order of ~1000 times the initial density/field. Further investigations should
include detailed ionization processes and more accurate radiation modeling to
properly capture the radiation transport and subsequent target compression.

Coupling separately developed codes offers an attractive method for
increasing the accuracy and fidelity of the computational models. Examples
include the earth sciences and fusion integrated modeling. This paper describes
the Framework Application for CoreEdge Transport Simulations (FACETS).

A discontinuous Galerkin method for the ideal 5 moment twofluid plasma
system is presented. The method uses a second or third order discontinuous
Galerkin spatial discretization and a third order TVD RungeKutta time stepping
scheme. The method is benchmarked against an analytic solution of a dispersive
electron acoustic square pulse as well as the twofluid electromagnetic shock
and existing numerical solutions to the GEM challenge magnetic reconnection
problem. The algorithm can be generalized to arbitrary geometries and three
dimensions. An approach to maintaining small gauge errors based on error
propagation is suggested.