Changes between Version 10 and Version 11 of u/erica/CoreCollapseBlog
- Timestamp:
- 03/15/18 18:37:04 (7 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
u/erica/CoreCollapseBlog
v10 v11 12 12 I’m not too sure why the astrobear’s output does not perfectly line up along the density progenitor input profile towards the center of the core. I suspect it has to do with the interpolation scheme and the calculated mass enclosed vs. average density for the cell center. 13 13 14 Figures 3 & 4 are zoom-ins of the density and pressure profiles between the innermost zone (located at dx=5.6e+7 cm) and the radius at which the profiles begin to noticeably deviate. Note the progenitor profile has data sampling down to radius of r=7.84e+5 cm, so about 2 orders of magnitude finer resolution than in this example simulation. 15 14 16 [[Image(2D_rho_zoomin.png, 50%)]][[br]] 15 17 '''Figure 3.''' Zoom-in of density profiles described in text. … … 18 20 '''Figure 4.''' Zoom-in of pressure profiles described in text. 19 21 20 Figure 3 is a zoom-in of the profiles between the inner most zone (located at dx=5.6e+7 cm) and the radius at which the profiles begin to noticeably deviate. Note the progenitor profile has data sampling down to radius of r=7.84e+5 cm, so about 2 orders of magnitude finer resolution than in this example simulation.22 In addition to the given input progenitor profile (red curve) and the lineout from the simulation (blue curve), I’ve plotted two profile objects in these figures. The first is the initial profile object (populated with the progenitor input values, green curve). This is just the input progenitor values after they have been read into a readable format for astrobear’s HSE solver (i.e. a “profile object”). 21 23 22 In addition to the input progenitor profile (red curve) and the lineout from the simulation (blue curve), I’ve plotted two profile objects in Fig. 3. The first is the initial profile object (populated with the progenitor input values, green curve). This is just the input progenitor values after they have been read into a readable format for astrobear’s HSE solver (i.e. a “profile object”). 23 24 The second is the profile object after it has been sent through the HSE solver (yellow curve). Note, although it can’t be seen here, the initial and final density profile objects are exactly identical (green and yellow curves) and they perfectly line up with the progenitor density profile (red), but the green and yellow pressure profiles are marginally different (~<1%) near the center of the core with the green curve being identical to the initial progenitor pressure profile. This is because the HSE solver calculates what the pressure profile should be in HSE, given the input density profile. Since the input and output pressure profiles deviate slightly indicates that the initial progenitor profile wasn’t exactly in HSE (to within numerical uncertainty). 24 The second is the profile object after it has been sent through the HSE solver (yellow curve). Note, although it can’t be seen here, the initial and final ''density'' profile objects are exactly identical (green and yellow curves) and they perfectly line up with the progenitor density profile, but the green and yellow ''pressure'' profiles are marginally different (~<1%) near the center of the core with the green curve being identical to the initial progenitor pressure profile. This is because the HSE solver calculates what the pressure profile should be in HSE, given the input density profile. Since the input and output pressure profiles deviate slightly indicates that the initial progenitor profile wasn’t exactly in HSE (to within numerical uncertainty). 25 25 26 26 When we are talking overall differences of ~$10^8$ in the density profile, perhaps this relatively small discrepancy shouldn’t be worried about too much… Although, it is curious. I am going to run a simulation with a finer resolution to compare to these results here.