Re: [mesa-users] Rapidly Rotating A-Stars

Michel Rieutord <[email protected]>
Newsgroups gmane.comp.video.mesa3d.user
Message-ID <[email protected]>
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 :[email protected]
> 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
>> [email protected] <mailto:[email protected]>
>>
>>
>> ------------------------------------------------------------------------------
>> 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
>> [email protected]
>> https://lists.sourceforge.net/lists/listinfo/mesa-users
>

_______________________________________________
mesa-users mailing list
[email protected]
http://lists.freedesktop.org/mailman/listinfo/mesa-users
lmpx.com only provides a reader for public news (NNTP) servers. It is not affiliated with the servers or forums shown here and is not responsible for the content of articles, which is written by their respective authors.