[Mesa-users] difficulty with custom mass loss

[TUCHOW, NOAH WOLFE]

The main reason why I wasn't using other_wind is that the mass loss which I wanted to model wasn't exclusively due to wind. If I were to use an other_wind hook, would it matter if the mass loss isn't actually wind? Would I still modify s% star_mdot to set the mass loss rate, or is there another variable I should change? If I wanted to set a custom wind for main sequence evolution would that be setting hot_wind_scheme = 'other' in the inlist?

Thanks,
Noah Tuchow


Noah Tuchow
Penn State University
Department of Astronomy and Astrophysics
532C Davey Lab
University Park, PA 16802


On Mon, Jul 15, 2019 at 9:47 PM Aaron Dotter <aaron.dotter at gmail.com<mailto:aaron.dotter at gmail.com>> wrote:
Hi Noah,

Just a suggestion: Have you tried using other_wind?  I've used that successfully a number of times over the years.  It might be what you're looking for rather than other_adjust_mdot.

Aaron

On Mon, Jul 15, 2019 at 9:34 PM TUCHOW, NOAH WOLFE <nxt5109 at psu.edu<mailto:nxt5109 at psu.edu>> wrote:
Hello everyone,

I've been having a bit of difficulty trying to implement a custom mass loss rate for a main sequence star. I've started off simple, just slightly modifying the default work folder's inlists to make a 1 solar mass star losing mass with a mass loss rate given as a function of time. Since I want to specify a function for a custom mass loss mechanism, I have been trying to use the "other" hooks in run_star_extras to do so. However, I can't seem to get even the simplest example, a constant mass loss rate, working via run star extras. I made a simple subroutine called const_mdot, and I direct MESA to use it in extras_controls as an other_adjust_mdot. I am trying to set  s% star_mdot to be a value that I specify in the inlist, but, when I run MESA, I find that the mass of the star doesn't change at all.  I noticed that setting mass_change in the inlist controls does what I thought this would be doing, but I would eventually like to implement more complicated mass loss rates.

Am I not implementing mass loss in the correct manner? Should I not put my mass loss subroutine as an other_adjust_mdot, or is s% star_mdot not the right quantity to be modifying? I know there is an other_wind option, but I'm dealing with mass loss from other sources such as CMEs not only wind. I've attached my inlists and run_star_extras.f  which I've been using in MESA r10398.  How does one specify a custom mass loss rate in MESA?

Thanks,
Noah Tuchow

Noah Tuchow
Penn State University
Department of Astronomy and Astrophysics
532C Davey Lab
University Park, PA 16802

_______________________________________________
mesa-users at lists.mesastar.org<mailto:mesa-users at lists.mesastar.org>
https://lists.mesastar.org/mailman/listinfo/mesa-users<https://nam01.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.mesastar.org%2Fmailman%2Flistinfo%2Fmesa-users&data=02%7C01%7Cnxt5109%40psu.edu%7Cc95c63496c234f099d5908d7098f817a%7C7cf48d453ddb4389a9c1c115526eb52e%7C0%7C0%7C636988384274156909&sdata=Va99kcbM73X8xbhqtzDTFEpK9FFEXim%2FfJ%2BC%2B%2B%2Bft8Y%3D&reserved=0>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20190716/5ea88e25/attachment-0001.html>