[mesa-users] Brunt Vaisala Frequency

[Jim F. <sla...@gm...>]

Hey everybody,

I've been using mesa to calculate the Brunt-Vaisala frequency in main
sequence stars of various masses. Since the brunt freq. is proportional to
the force of gravity, it should go to zero at the center of the star,
regardless of mass. However, in the models I am generating, it approaches
some finite positive value near the center of my stars (for stars of various
masses and ages). Has anybody run into this problem or can point me to the
location in the code where the brunt frequency is calculated? Thanks,

-Jim