In the Galactic fountain scenario, supernovae and/or stellar winds propel material into the Galactic halo. As the material cools, it condenses into clouds. By using FLASH three-dimensional magnetohydrodynamic simulations, we model and study the dynamical evolution of these gas clouds after they form and begin to fall toward the Galactic plane. In our simulations, we assume that the gas clouds form at a height of z=5 kpc above the Galactic midplane, then begin to fall from rest. We investigate how the cloud's evolution, dynamics, and interaction with the interstellar medium (ISM) are affected by the initial mass of the cloud. We find that clouds with sufficiently large initial densities (> 0.1 hydrogen atoms per cc) accelerate sufficiently and maintain sufficiently large column densities as to be observed and identified as high-velocity clouds (HVCs) even if the ISM is weakly magnetized (1.3 micro Gauss). We also investigate the effects of various possible magnetic field configurations. As expected, the ISM's resistance is greatest when the magnetic field is strong and perpendicular to the motion of the cloud. The trajectory of the cloud is guided by the magnetic field lines in cases where the magnetic field is oriented diagonal to the Galactic plane. The model cloud simulations show that the interactions between the cloud and the ISM can be understood via analogy to the shock tube problem which involves shock and rarefaction waves. We also discuss accelerated ambient gas, streamers of material ablated from the clouds, and the cloud's evolution from a sphere-shaped to a disk- or cigar-shaped object.

In order to determine the circumstances under which isolated SNRs are capable of rising into and enriching the thick disk and galactic halo, simulations of supernova remnants are performed with the FLASH magnetohydrodynamic code. We performed simulations in which the interstellar magnetic field is parallel to or perpendicular to the galactic plane as well as a simulation without a magnetic field. The ambient gas density distribution and gravitational potential are based on observations of our galaxy. We evolve the remnants to ages of roughly 10,000,000 years. For our simulation without a magnetic field, we compare the evolution of the hot bubble's velocity with the velocity evolution calculated from the buoyant and drag accelerations. We found surprisingly small vertical velocities of the hot gas, from which we estimated the drag coefficient to be ten for the non-magnetic simulation. Although we found little buoyant motion of the hot gas during the remnant's lifetime, we found rapid vertical motion of the associated cool dense gas near the end of the remnants life. This motion deformed the remnant into a mushroom cloud structure similar to those found in previous simulations. The simulation in which we have a 4 micro-Gauss magnetic field parallel to the galactic mid-plane shows a dramatically elongated bubble parallel to the magnetic field. The magnetic field pins the supernova remnant preventing it from rising. In the simulation with the 4 micro-Gauss magnetic field perpendicular to the midplane the hot bubble rises more, indicating that having the magnetic field in the same direction as the gravitational force enhances the rise of the bubble.

We have developed a formalism for analyzing the stability of a fluid in the presence of neutrinos of all flavors and in the presence of a gravitational field. When applied to an extensive two-dimensional grid of core radii and fluid element sizes for each of several time slices of a number of proto-supernovae, we find no evidence for the neutron finger instability as described by the Livermore group. We find, instead, that the rate of lepton equilibration always exceeds that of thermal equilibration. Furthermore, we find that thermal diffusion driven by a lepton fraction difference and lepton diffusion driven by an entropy difference are both nonzero and that the first of these tends to be large in magnitude. An important consequence of this is the presence of a doubly diffusive instability, which we refer to as ``lepto-entropy fingers,'' in an extensive region below the neutrinosphere where the lepton number, $Y_{\ell}$, is small. This instability is driven by a mechanism very different from that giving rise to neutron fingers, and may play an important role in enhancing the neutrino emission. Deep in the core where the entropy is low and the lepton number higher, our analysis indicates a region unstable to another instability which we refer to as ``lepto-entropy semiconvection.'' These instabilities, particularly lepto-entropy fingers, may have already been seen in some multi-dimensional core collapse simulations described in the literature.