There are multiple examples of gravitational losses in neutron stars and black holes. Protons and neutrons have been found to have enormous repulsive pressures that highly squeezed collapsing matter cannot overcome. The case against singularities follows. Galactic black hole gravitational losses can supply the missing dark energy. With highly squeezed nucleons, the big bang could begin as a hot core and a cold dark matter shell. The 21 cm. radiation data has identified baryon sized particles as cold dark matter. Highly squeezed nucleons will not decompose to produce antimatter. The flatness of the Universe is due to a baryonic bounce. The highly correlated galaxies originated from primordial black holes capturing hot core gasses.There is evidence that galaxies have not grown nor merged significantly since formation.

Using a code that employs a self-consistent method for computing the effects of photoionization on circumstellar gas dynamics, we model the formation of wind-driven nebulae around massive Wolf-Rayet (W-R) stars. Our algorithm incorporates a simplified model of the photo-ionization source, computes the fractional ionization of hydrogen due to the photoionizing flux and recombination, and determines self-consistently the energy balance due to ionization, photo-heating and radiative cooling. We take into account changes in stellar properties and mass-loss over the star's evolution. Our multi-dimensional simulations clearly reveal the presence of strong ionization front instabilities. Using various X-ray emission models, and abundances consistent with those derived for W-R nebulae, we compute the X-ray flux and spectra from our wind bubble models. We show the evolution of the X-ray spectral features with time over the evolution of the star, taking the absorption of the X-rays by the ionized bubble into account. Our simulated X-ray spectra compare reasonably well with observed spectra of Wolf-Rayet bubbles. They suggest that X-ray nebulae around massive stars may not be easily detectable, consistent with observations.

Using an ionization gasdynamics code, we simulate a model of the wind-blown bubble around a 40 solar mass star. We use this to compute the X-ray spectra from the bubble, which can be directly compared to observations. We outline our methods and techniques for these computations, and contrast them with previous calculations. Our simulated X-ray spectra compare reasonably well with observed spectra of Wolf-Rayet bubbles. They suggest that X-ray nebulae around massive stars may not be easily detectable, consistent with observations.

Eliptical and bulge galaxies share a tight correlation of velocity distribution to both luminosity and black hole mass. There are similar orbital speeds for all galaxies of a given luminosity including dark matter (DM) at large radii. The halo surface density of DM is constant for almost all types of galaxies and ranges 14 mag. down to dwarf spherical galaxies. There are supermassive black holes or giant, pure disk galaxies at high redshift inexplicable with hierarchical clustering or collapse dynamics. These and a myriad of other galaxy formation problems are explainable by an initial shell which caused the Planck cosmic microwave background radiation. A reduction in the energy-density of primordial galactic black holes is necessary to explain dark energy.

Minutes into the big bang, nucleosynthesis finds 96 percent of matter nonreactive. Massive pure disk protogalaxies must be formed in the early universe as their galaxies are found at high redshift. Surrounded by hot gas, they made little imprint on the CBR. The LHC has not found the Higgs or dark matter candidates. These phenomena can be explained by a cold baryonic shell.

A recently developed spectral-element adaptive refinement incompressible magnetohydrodynamic (MHD) code [Rosenberg, Fournier, Fischer, Pouquet, J. Comp. Phys. 215, 59-80 (2006)] is applied to simulate the problem of MHD island coalescence instability (MICI) in two dimensions. MICI is a fundamental MHD process that can produce sharp current layers and subsequent reconnection and heating in a high-Lundquist number plasma such as the solar corona [Ng and Bhattacharjee, Phys. Plasmas, 5, 4028 (1998)]. Due to the formation of thin current layers, it is highly desirable to use adaptively or statically refined grids to resolve them, and to maintain accuracy at the same time. The output of the spectral-element static adaptive refinement simulations are compared with simulations using a finite difference method on the same refinement grids, and both methods are compared to pseudo-spectral simulations with uniform grids as baselines. It is shown that with the statically refined grids roughly scaling linearly with effective resolution, spectral element runs can maintain accuracy significantly higher than that of the finite difference runs, in some cases achieving close to full spectral accuracy.

We are investigating the SN 1987A remnant by modeling the circumstellar environment of the progenitor star. Interacting stellar winds models have been reasonably successful at reproducing the gross features of the nebula, but some details, such as the early rise of the radio and x-ray emission from the supernova blast and the presence of the outer rings, are not explained in a pure wind model. In this paper we describe preliminary results from 2-D models that incorporate the effects of photoionization due to radiation from the central star. These models have successfully produced a thick HII region, as needed to produce the early radio and X-ray emission. The HII region is thickest away from the equatorial plane. The models have also produced a feature resembling the outer rings, but we suspect this to be an artifact of the 2-D calculations that would not persist in 3-D.

The tight correlation of galactic velocity distribution to both luminosity and its black hole mass and the relation of halo parameters to luminous mass distribution, can not be due to collapse dynamics. A big bang shell can solve galaxy formation problems by forming the supermassive black holes necessary to capture the initial blast wave in a coordinated pattern.

To propel the farthest galaxies at 0.5c, the big bang required a core energy of 160MeV. Its baryon-photon ratio becomes $\eta=87.6\Omega_Bh^2$ which is comparable to a supernova. Accretion photons emitted during the prior contraction supplied the CMBR planck spectrum. The finding of primordial 2^H levels varying by an order of magnitude is evidence of large scale structure during nuclesynthesis. Rotating black holes from the primordial shell were the density fluctuations for galaxy formation. Sufficient baryonic matter was recaptured during the expansion to increase the universe scale factor and allow supernovae Ia to occur in a closed universe. A low energy bounce model with zeroing of the stress-energy tensor can solve the hot big bang problems.