[mesa-users] trying to run RLO calculations for hot Jupiter systems

[Robert Farmer]

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
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20170228/e71f4e17/attachment.html>