[Mesa-users] [mesa-users] Rapidly Rotating A-Stars

[Michel Rieutord]

Le 05/10/2013 14:18, Michel Rieutord a écrit :
> Hi Jeremy,
>
> I think the best answer to your question is to be found with the ESTER
> code that you can freely download at
> http://code.google.com/p/ester-project/
>
> This is a code that computes the 2D structure, including differential
> rotation and meridional circulation, of rotating early-type stars at
> any rotation rate (below critical of course). With this code you'll be
> able to compare your data with those of the model including the
> temperature dependence with latitude; by adjusting the hyrdogen content
> of the core you may get an idea of the evolution status of your star on
> the main sequence. This won't give you an age in years but this can
> obtained with a 1D code afterwards. We are working out the time
> evolution but this is not yet ready; this a tough problem all the more
> that we need to also evolve the dynamics...
>
> You may find examples of the results in
> Espinosa Lara & Rieutord 2013, A&A, 552, A35
> Rieutord & Espinosa Lara 2013, LNP 865, 49, Arxive 1208.4929
> or  2013arXiv:1310.0908   or   arXiv:1310.1017
>
> Don't hesitate if you have questions.
>
> Best wishes,
>
> Michel Rieutord
>
PS third attempt; hope this is the good one.
> **************************************************************************
> Michel Rieutord
> Institut de Recherche en Astrophysique et Planétologie
> Observatoire Midi-Pyrénées      email :Michel.Rieutord at irap.omp.eu
> 14 avenue Edouard Belin                Tel +33 -(0)5 61 33 29 49
> 31400 Toulouse, France                 Fax +33 -(0)5 61 33 28 40
> **************************************************************************
> Le 04/10/2013 22:25, Jeremy Jones a écrit :
>> Hey mesa users,
>>
>> I'm working on a project concerning interferometry of A-stars. Many 
>> of my targets are rapid rotators, and I'm looking into using mesa to 
>> determine their ages. I've played around with mesa a bit (I've 
>> included one of my inlists below). I have a question about how to 
>> compare my observables with the mesa outputs, in particular radius 
>> and temperature.
>>
>> From the modeling I do of the inerferometric data, I get the full 
>> radius profile of the star (i.e. radius as a function of latitude). 
>> Likewise, I get the full temperature profile of the star.
>>
>> My question is for a rapid rotator, which radius/temperature should I 
>> compare with the mesa outputs to get the age of the star? Polar 
>> radius/temperature, equatorial radius/temperature, something more 
>> involved?
>>
>> Thanks for any insight!
>>
>> ~Jeremy
>>
>> The aforementioned inlist:
>> ================
>> &star_job
>>
>>       mesa_dir = '../..'
>>
>>       set_initial_age = .true.
>>       initial_age = 0 ! in years
>>
>>       set_initial_model_number = .true.
>>       initial_model_number = 0
>>
>>       create_pre_main_sequence_model = .true.
>>       new_omega_div_omega_crit = 0.5
>>       near_zams_relax_omega_div_omega_crit = .true.
>>
>>       pgstar_flag = .true.
>>
>>
>> / ! end of star_job namelist
>>
>>
>> &controls
>>
>>       min_timestep_limit = 1d-12 ! (seconds)
>>
>>       ! check for retries and backups as part of test_suite
>>       ! you can delete this for use outside of test_suite
>>          max_number_backups = 400
>>          max_number_retries = 900
>>          max_model_number = 1500
>>
>>          !max_allowed_nz = 20000 ! maximum number of grid points allowed
>>
>>          initial_mass = 2
>>          initial_z = 0.016
>>
>>          photostep = 10
>>          profile_interval = 10
>>          history_interval = 1
>>          terminal_interval = 2
>>          write_header_frequency = 10
>>
>> / ! end of controls namelist
>>
>>
>> &pgstar
>>
>> History_Track1_win_flag = .true.
>> History_Track1_yname = 'log_R'
>> History_Track1_yaxis_label = 'log(R/R_sun)'
>> History_Track1_xname = 'log_Teff'
>> History_Track1_xaxis_label = 'log(T/K)'
>> History_Track1_title = 'Radius vs. Effective Temperature'
>>
>> History_Track2_win_flag = .true.
>> History_Track2_yname = 'surf_avg_v_rot'
>> History_Track2_yaxis_label = 'Equatorial Velocity (km/s)'
>> History_Track2_xname = 'star_age'
>> History_Track2_xaxis_label = 'Age (yr)'
>> History_Track2_title = 'Surface Velocity vs. Age'
>>
>> / ! end of pgstar namelist
>> ================
>>
>> -- 
>> Jeremy Jones
>> Graduate Student Assistant
>> Georgia State University
>> Department of Physics & Astronomy
>> jones at chara.gsu.edu <mailto:jones at chara.gsu.edu>
>>
>>
>> ------------------------------------------------------------------------------
>> October Webinars: Code for Performance
>> Free Intel webinars can help you accelerate application performance.
>> Explore tips for MPI, OpenMP, advanced profiling, and more. Get the most from
>> the latest Intel processors and coprocessors. See abstracts and register >
>> http://pubads.g.doubleclick.net/gampad/clk?id=60134791&iu=/4140/ostg.clktrk
>>
>>
>> _______________________________________________
>> mesa-users mailing list
>> mesa-users at lists.sourceforge.net
>> https://lists.sourceforge.net/lists/listinfo/mesa-users
>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.freedesktop.org/archives/mesa-users/attachments/20131006/a2f42d9e/attachment.html>