[mesa-users] trying to run RLO calculations for hot Jupiter systems
Pablo Marchant
pamarca at gmail.com
Tue Feb 28 19:38:05 EST 2017
Besides from the nuclear network discrepancy, here is an issue that Evan
ran into and reported to me a couple of weeks ago. The default
implementation for the Eddington mass accretion limit applies to BH
accretors, and it breaks down and produces NaNs when evolving both stars
simultaneously with limit_retention_by_mdot_edd = .true.. This translates
into the obscure error messages you are getting Brian.
Just commited a fix that will be up in the next MESA version. For the
moment, just set limit_retention_by_mdot_edd = .false. in your inlist and
things should work.
Cheers!
On Tue, Feb 28, 2017 at 5:47 PM, Robert Farmer <rjfarmer at asu.edu> wrote:
> Hi
> The nuclear networks (and hence the isotopes each model has defined) in
> the two model files are different, can you try in the planet inlist
>
> new_net_name = 'o18_and_ne22.net'
> change_initial_net = .true.
>
> so the planet model will have the same network as the star
>
> Rob
>
>
> On Tue, Feb 28, 2017 at 4:30 PM, Brian Jackson <bjackson at boisestate.edu>
> wrote:
>
>> Dear MESA Users,
>> In my group, we're trying to continue our analysis of planetary
>> Roche-lobe overflow (RLO).
>>
>> For that, we generate starting planetary and stellar models
>> (attached). The star model is generated using the 1M_pre_ms_to_WD test
>> suite with only "max_age = 20e6" added to inlist_1.0 file. The
>> planetary model is generated using the make_planets test suite with
>> only "max_age = 20d6" changed in inlist_evolve_1.0_MJ_10.0_ME_2.0_RJ
>> file.
>>
>> We make some modifications to inlists for the binary_both_stars test
>> suite to load in the saved model files but otherwise use the standard
>> test suite.
>>
>> The following shows the last few lines from MESA before it stops
>> running our calculation:
>>
>> get_xa_for_accretion: accretion species mass fractions do not add to 1.0
>> sum(xa(1:species))
>> 0.0000000000000000D+00
>> h1 1
>> 0.0000000000000000D+00
>> he3 2
>> 0.0000000000000000D+00
>> he4 3
>> 0.0000000000000000D+00
>> c12 4
>> 0.0000000000000000D+00
>> n14 5
>> 0.0000000000000000D+00
>> o16 6
>> 0.0000000000000000D+00
>> ne20 7
>> 0.0000000000000000D+00
>> mg24 8
>> 0.0000000000000000D+00
>> 1 1 2.000000000E+07
>> 9.244518984E-04 -9.900000000E+01 -3.245769587E-01
>> 2 1 2.000000000E+07
>> 9.999853565E-01 -9.900000000E+01 -6.799606683E-01
>> have reached retry limit so now backup
>> have reached retry limit so now backup
>>
>> dt
>> 6.3566365930376634D-07
>> min_timestep_limit
>> 9.9999999999999995D-07
>>
>> stopping because of convergence problems dt < min_timestep_limit
>>
>> terminated evolution: convergence problems
>>
>> termination code: min_timestep_limit
>>
>> ---
>>
>> The species mass fractions look like nonsense, even though they look
>> fine in the initial model files. And strangely, they change depending
>> on the initial orbital period we choose.
>>
>> I am running Mac OS X 10.9.5, MESA version 8845, and the MESA SDK
>> version from April 4 2016. I have attached the binary inlists and
>> initial models we use to run the calculation.
>>
>> Please let me know if you need additional information and if you have
>> any advice.
>>
>> Many thanks,
>> Brian Jackson
>>
>> ------------------------------------------------------------
>> ------------------
>> Check out the vibrant tech community on one of the world's most
>> engaging tech sites, SlashDot.org! http://sdm.link/slashdot
>> _______________________________________________
>> mesa-users mailing list
>> mesa-users at lists.sourceforge.net
>> https://lists.sourceforge.net/lists/listinfo/mesa-users
>>
>>
>
> ------------------------------------------------------------
> ------------------
> Check out the vibrant tech community on one of the world's most
> engaging tech sites, SlashDot.org! http://sdm.link/slashdot
> _______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users
>
>
--
Pablo Marchant Campos
M.Sc on Astrophysics, Universidad Católica de Chile
PhD student, Argelander-Institut für Astronomie
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20170228/66bcd47f/attachment.html>
More information about the Mesa-users
mailing list