PPN

I used radmc3d to produce a radiative transfer simulation based on the result of my hydro simulation. The radiative transfer code use Monte Carlo method. Below are the wavelength picture correspondingly from observation and my simulation.

This is a synthetic radiative transfer simulation. It corresponds to the dense shell model​dense shell hydro simulation. I also analyzed the density and velocity properties of constant wind and dense shell model here. ​What defines a disk or a torus?

The viewing angle is the same as L2 pup, that is 83 degree. In this radiative transfer simulation, I assumed a 1 AU radius, 3000 K black body in the center of the simulation zone. I also assumed the dust/(dust+gas) mass ratio is 0.05. This is a bit high dust gas ratio but since the actual binary separation is 2 AU (in the simulation it is 4 AU), the optical depth may go higher if the separation shrink to 2 AU.

I do not know the particular species of the dust (probably silicates). I adopted the default dust species setting in radmc3d, I make the scattering to be isotropic even though the dust scattering is likely to be forward.

Further research can be done on:

  1. Redo 2 AU separation simulation.
  1. Figure out the dust species and their scattering, absorption properties.
  1. If there is enough computational resources, I want to do a larger simulation zone. Currently, for 2 AU separation, I can only go 40AU*40AU*40AU with 100*100*100 box plus 3 level AMR and it takes 120CPU*7days.

Attachments (5)

Download all attachments as: .zip

Comments

No comments.