Updates on rad feedback
Run 1.
2D, periodic boundary conditions, 1282 resolution. There is a sink particle in the middle of the box that is injecting a constant 'accretion energy' into the center zones of the grid.
Plot is of Etherm (left) and Erad (right) over time,
Thermal energy is derived from subtracting the specific kinetic energy from the total energy (as there are no magnetic fields or gravity in this simulation). The total radiative energy in the box at the start of the sim is, Erad=90. By the end, Erad=10,038.
Both the planck opacity and the roseland are set to =1d-6. The right plot shows that the radiative energy diffuses through the grid much faster than it is exchanging energy with the gas (i.e. changes in thermal energy happen slower). Hardly any change in thermal energy is happening at all, in fact. Unless I made an error with my expression, not sure how the radiative energy can be increasing, yet the thermal and total energy are equal..
There is a slight 'blip' noticeable between the 1st and 3rd frame. Not immediately sure what this is being caused by either..
I would like to check that the total energy is being conserved. To do this, should I make an expression that sums E and Erad and check that it is constant over the course of the sim?
Overall though, the run seems to proceed nicely for these parameters (need to check my thermal energy expression still). I am trying to set up a simulation where I control the rate of energy injection, such that the total Erad doubles over the course of the simulation. In the meantime, experimented with differing opacities, as will show next.
Run 2.
High roseland opacity (radiation trapped, k=.001), low planck opacity (no coupling of radiation and gas, k=.1d-6).
Note, no blip in this run. I see that the thermal energy is not changing much here either. The max Erad is much greater than compared to above. The radiation being confined to a smaller space explains this.
Run 3.
High roseland opacity (radiation trapped, k=.001), High planck opacity (coupling of radiation and gas, k=10).
Now am starting to see an increase in thermal energy. Need to look at the equations and see how the opacities are coming in and how they are related to each other.
Attachments (17)
- e_rad0041.png (24.3 KB) - added by 9 years ago.
- e_rad.gif (5.8 MB) - added by 9 years ago.
- erad_etherm0027.png (28.1 KB) - added by 9 years ago.
- erad_etherm.gif (4.8 MB) - added by 9 years ago.
- total_e.gif (7.6 MB) - added by 9 years ago.
- slow2.gif (146.5 KB) - added by 9 years ago.
- 2erad_etherm.gif (1.9 MB) - added by 9 years ago.
- 2erad_etherm0027.png (25.7 KB) - added by 9 years ago.
- test2.png (24.8 KB) - added by 9 years ago.
- test.png (25.2 KB) - added by 9 years ago.
- highhigherad_etherm0027.png (22.2 KB) - added by 9 years ago.
- highhigherad_etherm.gif (1.5 MB) - added by 9 years ago.
- lowhigherad_etherm.gif (1.5 MB) - added by 9 years ago.
- lowhigherad_etherm0027.png (21.6 KB) - added by 9 years ago.
- newerad_etherm.gif (2.3 MB) - added by 9 years ago.
- newerad_etherm0100.png (30.7 KB) - added by 9 years ago.
- newslow2.gif (70.5 KB) - added by 9 years ago.
Comments
No comments.