[Mesa-users] Using other_momentum in Riemann solver

[Jing-Ze Ma]

Hi everyone,


I'm trying to add another kind of force into the momentum equation in the ppisn test_suite, but the other_momentum and other_momentum_implicit hooks are not implemented into the Riemann solver. I can always introduce a vector change to cgrav (setting cgrav to different values in different cells) so that an additional force is equivalently added to the momentum equation, but I think the energy equation would be physically wrong according to section 4 of MESA instrument paper IV published in 2018 (or perhaps I'm getting it wrong?). Is there any other way to add another force into Riemann solver using hooks?


I'm using MESA r15140.


Thanks a lot!


Best regards,
Jing-Ze


--

Jing-Ze Ma
Undergraduate Student,
Tsien Excellence in Engineering Program (TEEP) 2018,
School of Aerospace Engineering,
Tsinghua University

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20210219/fcb89aff/attachment.htm>