We have decided to analyze the accretion onto a neutron star just up until the hydrogen burning begins. We'll first attempt this as an adaptation to the Bondi accretion problem with a solid sphere near the origin (also relevant to the problem of infall onto the proto-neutron star during core collapse). If we're successful, we'll try an accretion disk in 2D with polar geometry.
Determined to be outside scope of project.
* [ ] Determine the input abundances
* [ ] Build the nuclear network for the problem
* [ ] Take output from the codes
* [ ]determine when hydrogen burning begins
* [ ]analyze data values to optimize scale(linear vs log) and ranges(0-?) for visualization
1D on a solid surface was successful so the 2D case was the bulk of the project work. The solid case was unstable and resulted in an non-physical explosion. Both the infall and reflective cases were attempted, we found the  to be optimal. We were able to model a non-viscous and viscous case using hydronamic force and energy equation. Terms for the extremely large B field (10^7-8 T) were not included. Both accreted around the central mass with a large high pressure/empty barrier for the non-viscous case, as the angular momentum is not dissipated. The viscous case correctly accreted on to the surface.
(further resources in presentation)