[Mesa-users] Query regarding implementing hydrodynamics

[Abha Vishwakarma]

Dear MESA Users,

I have been trying to switch on hydrodynamics in a 12 solar mass red giant
star. I first ran the hse_reimann test suite which implements
hydrodynamics in a 12 solar mass star with the core excised. But, I wish to
implement hydrodynamics with the core intact.

I switched on hydrodynamics by adding the following lines to my inlist:

        change_initial_u_flag = .true.
    change_initial_v_flag = .true.
    new_v_flag = .false.
    new_u_flag = .true.

My program takes as input a 12 solar mass red giant. But, the program
crashes because it hits the timestep limit or because of crossing some
tolerances. But, if I instead take the 12 solar mass star with the core
removed, as is done in the hse_reimann module, the code works fine.

I am not sure what I am missing here, but I feel it has something to do
with the core of the red giant, or some boundary conditions or physics that
I might have missed. I have attached my inlist to this email, and any help
in this regard would be greatly appreciated!

Sincerely,
Abha Vishwakarma
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20240308/60161d1c/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist_project
Type: application/octet-stream
Size: 1871 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20240308/60161d1c/attachment.obj>