More analysis on PN

Jet in x-direction

This is mostly to address a comment from Eric, if we put in an equatorial outflow, will it still be collimated?

Instead a "disk" outflow, I put a jet in the +- x-directions, with 15-degree open angle, otherwise the same as model A.

http://www.pas.rochester.edu/~yzou5/PN_xjet_figures.pdf

On momentum (L/c)

In our simulation, the "momentum rate" (L/c) is set by the outflow's density and velocity

for high-momentum flux (models A and C), (L/c)_high = 3.835e+34 (g cm s-1 year-1), and

for low-momentum flux (models B and D), (L/c)_high = 1.917e+33 (g cm s-1 year-1)

At the last frame presented in PN paper, the total momentum ejected is then:

Observationally, the upper limit on t comes from the size and jet speed: t<R/v, and maximum momentum is estimated as:

Here I use the extend of the top lobe for R, and wind speed v=300 km/s, the momentum calculated in those two ways and their ratio are:

model time at measurement (days) P_sim (g cm s-1) P_measure (g cm s-1) P_sim / P_measure
A 920 9.665e+34 6.874e+34 1.41
B 1780 9.350e+33 3.476e+33 2.69
C 1600 1.681e+35 6.946e+34 2.42
D 4000 2.101e+34 3.473e+33 6.05

Initial Mass Flow

This seems problematic. I will check my Visit code and come back with another update. Ignore the following for now.

Outer boundary flows inward, floor density raise remove some mass, but both are small.

total change: -0.2 Msun

outer BC: 1.1e-7 Msun

floor raise: -1.4e-3 Msun

Left: including frame 0; Right: excluding frame 0.

http://www.pas.rochester.edu/~yzou5/mass_change.png http://www.pas.rochester.edu/~yzou5/mass_no_frame0.png

Comments

No comments.