[mesa-users] irradiation

[ztt]

Hi Jean,

I tried with your suggestion, below is related part of my inlist:

      load_saved_model = .true.
      saved_model_name = 'results_star_irrad/irrad/irrad1.mod'
      set_irradiation = .true.
      relax_irradiation = .true.
      relax_irradiation_min_steps = 200
      relax_irradiation_max_yrs_dt = 1d-1

      new_irrad_flux = 1d7 ! erg s^-1 cm^-2
      new_irrad_col_depth = 1 ! g cm^-2

      !relax_tau_factor = .true. ! gradually change tau_factor, reconverging at each step.
      !relax_to_this_tau_factor = 150 ! tau_factor of 1 puts outer cell at photosphere.
      !dlogtau_factor = 0.1 ! change log10(tau_factor) by at most this amount per step
      set_tau_factor = .true. ! change tau_factor without reconverging.
      set_to_this_tau_factor = 100
However, it doesn't seem to be working, since I compare the effective temperature evolution with another plot from the run without irradiation and they are not quite different, while the irradiation heat flux is certainly not a negligible value here. Please see the attachment plots(1st is with irradiation). Million Thanks.

Best wishes, Tingtao

At 2011-07-23 12:48:09,"Jean-Claude Passy" <jcpassy at gmail.com> wrote:
Hi,

no you need both. The 'relax_irradiation' helps you getting started by relaxing the star as you turn on irradiation progressively. The section in your inlist should look like:

----------------------------------------------------
set_irradiation = .true.

relax_irradiation = .true.
relax_irradiation_min_steps = 200
relax_irradiation_max_yrs_dt = 1d-1

new_irrad_flux = 1d7 ! erg s^-1 cm^-2
new_irrad_col_depth = 1 ! g cm^-2
----------------------------------------------------

Hope that helps.
Good luck,

JC


On 22/07/11 21:59, ztt wrote:
  Hi Jean,

Thanks for your reply. My mail box didn't get your email but I saw it on the mailing list.
So is this the correct procedure example for turning on irradiation??
1. do relax with:
      relax_irradiation = .true.
      relax_irradiation_min_steps = 200
      relax_irradiation_max_yrs_dt = 1d-1

      new_irrad_flux = 1d7 ! erg s^-1 cm^-2
      new_irrad_col_depth = 1 ! g cm^-2
with 'set_irradiation = .true.' commented;
2. turn on irradiation with:
      set_irradiation = .true.
      new_irrad_flux = 1d7 ! erg s^-1 cm^-2
      new_irrad_col_depth = 1 ! g cm^-2
but with all the relax_* commented now;

while keep the setup for tau & atmosphere all the time:
      set_tau_factor = .true. ! change tau_factor without reconverging.
      set_to_this_tau_factor = 100

      which_atm_option = 'grey_and_kap'
      !which_atm_option = 'grey_irradiated'

            atm_grey_irradiated_atol = 1d-4
            atm_grey_irradiated_rtol = 1d-4
            atm_grey_irradiated_T_eq = 1000 ! equilibrium temperature based on irradiation
               ! e.g., might use T_eq = Teff_star*(R_star/(2*Distance_to_star))^(1/2)
            atm_grey_irradiated_kap_nu = 4d-3 ! opacity for irradiation
            atm_grey_irradiated_P_surf = 1d6 ! 1 bar in cgs units
    &nb! sp;  ;      atm_grey_irradiated_max_tries = 50 ! limit on iterations
            trace_atm_grey_irradiated = .false.

Thank you very much!

Best wishes, Tingtao



-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20110723/44adaecb/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: log_Teff.pdf
Type: application/pdf
Size: 21981 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20110723/44adaecb/attachment.pdf>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: log_Teff.pdf
Type: application/pdf
Size: 12560 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20110723/44adaecb/attachment-0001.pdf>