mesa-users

[mesa-users] irradiation

[張世昕 <ss...@as...>]

Hi,

Two parameters are needed in the &controls namelist
when the irradiation is involved.
I could figure out the irradiation_flux, but don't know
how to determine the column_depth_for_irradiation (i.e. m_irradiation).

ex.
a planet (0.001Msun & atm=simple_photosphere) with a 3 day period
around a solar-type star. What's the column_depth_for_irradiation?

Thanks!!
Shih-Hsin

[mesa-users] irradiation

[ztt <ztt...@12...>]

  Dear Mesa users,

Do you have experience with the irradiated_planet test_suite? I tried with 2 different inlist:
      !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

and the other one:
      !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

both use which_atm_option = 'grey_and_kap'

then do ./star_irrad, and the results are the same.
Could you give me some hints about how to make irradiation work??
Thank you!

Best wishes, Tingtao

Re: [mesa-users] irradiation

[Jean-Claude P. <jc...@gm...>]

Hi Tingtao,

you have to remove the exclamation mark before 'set_irradiation = .true.'

Best,

JC

On 21/07/11 21:54, ztt wrote:
> Dear Mesa users,
>
> Do you have experience with the irradiated_planet test_suite? I tried
> with 2 different inlist:
> !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
>
> and the other one:
> !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
>
> both use which_atm_option = 'grey_and_kap'
>
> then do ./star_irrad, and the results are the same.
> Could you give me some hints about how to make irradiation work??
> Thank you!
>
> Best wishes, Tingtao
>
>

[mesa-users] irradiation

[ztt <ztt...@12...>]

  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
            atm_grey_irradiated_max_tries = 50 ! limit on iterations
            trace_atm_grey_irradiated = .false.

Thank you very much!

Best wishes, Tingtao

Re: [mesa-users] irradiation

[Jean-Claude P. <jc...@gm...>]

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
>
>

Re: [mesa-users] irradiation

[Bill P. <pa...@ki...>]

Hi Shih-Hsin,

I don't know if there is a "first prinicples" way to decide
on the value for column_depth_for_irradiation.
As a rule of thumb, you're looking for a value that is
small enough so that the energy is deposited in a region
of "not too large" optical depth (something < 10 for sure),
but at the same time large enough so that mesa can
deal with it (too tiny means it will try to deposit
all of the energy in one or two cells and that might
be tough for convergence.)   The test suite example
uses a column depth of 1 g cm^2 -- you might try that
and check to see the optical depth corresponding to that.

BTW: just a reminder to use
      kappa_file_prefix = 'gs98f'
so that you get the new Freedman low T opacities.

Also, you should know that there is another way
to deal with irradiated planets that uses a constant
pressure surface boundary condition following
Guillot & Havel, A&A 527, A20 (2011).
I haven't added a test_suite case for it yet (I'll do that soon!),
but let me know if you are interested in trying it
and I'll provide some more info for you.

Cheers,
Bill






On Jun 9, 2011, at 12:48 AM, 張世昕 wrote:

> Hi,
> 
> Two parameters are needed in the &controls namelist
> when the irradiation is involved.
> I could figure out the irradiation_flux, but don't know 
> how to determine the column_depth_for_irradiation (i.e. m_irradiation).
> 
> ex. 
> a planet (0.001Msun & atm=simple_photosphere) with a 3 day period 
> around a solar-type star. What's the column_depth_for_irradiation?
> 
> Thanks!!
> Shih-Hsin
> 
> 

Re: [mesa-users] irradiation

[Bill P. <pa...@ki...>]

Hi,

Oops -- I lied when I said we didn't have a test suite case showing the new way to do irradiated planets.
What a surprise to find that the test called 'irradiated_planet' is actually demonstrating the new way!
I answered your earlier email before I'd had a morning cup of coffee -- bad idea.  ;-)

The irradiated_planet example works in 2 steps.  
The 1st (inlist1) prepares a model for changing to the constant pressure surface boundary condition, 
and the 2nd (inlist2) switches to the new BC and evolves the planet with it.  

-B



On Jun 9, 2011, at 8:24 AM, Bill Paxton wrote:

> Also, you should know that there is another way
> to deal with irradiated planets that uses a constant
> pressure surface boundary condition following
> Guillot & Havel, A&A 527, A20 (2011).
> I haven't added a test_suite case for it yet (I'll do that soon!),
> but let me know if you are interested in trying it
> and I'll provide some more info for you.

[mesa-users] Error compiling EOS when installing MESA

[Ehsan M. <mor...@ia...>]

<html><div>
Dear MESA users.
Hi.

I tried to install MESA on a rather
old machine, but with a recent Ubuntu distribution, hence the last
up-to-dated gfortran compiler, but I received the following error:

gfortran -fopenmp -o ../sample sample_eos.o -L../../make -leos
-L../../../lib -lchem -lnum -lmtx -lmesalapack
-lmesablas -lalert -linterp_2d -linterp_1d -lutils -lconst
../utils/build_and_test: line 33: 21696 Segmentation
fault ./test_quietly

/home/emoravveji/mesa/eos/test
FAILED

I would appreciate
if someone help me with this.
With kind regards.
Ehsan
Moravveji.

-- 
This message has been scanned for viruses and
dangerous content by MailScanner, and is
believed to be clean.