It is been a year since I last wrote here. Since two weeks ago, the flux code is finally working, and flawlessly. I spent a lot of time trying to make things to work as I would like them to do, and it took me more time than I wanted and expected, but now it is working the way I expect.
The thing that took most of my time solving was the multi-level atom. The method I was using to calculate the 3-level hydrogen atom was flawed, and I was using it just to make the code run. Now I am using the output source function value of CV to calculate a 3-level hydrogen atom. CV uses a 2-level atom with continuum approach to calculate the line source function, without needing the number densities of the required levels, but using the line source function is possible to acquire the number densities. Since the medium is considered to be optically thin for Lyman lines, the population of the fundamental level can be calculated using thermodynamical equilibrium. Knowing the first level number density, and the ratio of number densities between the line transition levels, we can solve a 3-level atom. It seems simple, but this took months of my time. I tried many different ways to calculating these number densities, before I saw the way I am using now.
After that was done, I had to solve some problems with the code that were causing problems with the results. There were problem when calculating the grid used by Flux. Problems with the size of the grid boundaries. A problem after the grid re-sampling, when points that should be out of the region of diskwind and magnetospheric-accretion have values of densities, velocities, etc, different than zero. And the last problem solved, was that one of the codes was replacing the phi-component of the velocity to zero everywhere.
Now the code is running, the grid points used are correct, and also is the velocity and density fields. No changes were made in the temperature law.
No comments:
Post a Comment