[Mesa-users] wd_diffusion doesn't work?

[Evan Bauer]

Hi Michael,

Short answer: yes, this is expected behavior. You’re right to call attention to this sort of thing, but I wouldn’t worry about this particular test right now.

Long answer: this test case is out of date. It isn’t meant to test diffusion timescales, but conditions for diffusive equilibrium, so the length of timesteps shouldn’t be relevant for the physics being tested. It starts from a model that is already near diffusive equilibrium, and then runs for long enough (a few Gyr) that many diffusion timescales should have occurred in the atmosphere and H/He should have achieved diffusive equilibrium. It then tests against an analytic expression for H/He abundance profiles in diffusive equilibrium. See figure 23 in the second instrument paper for more details.

All of this is out of date because the wd_diffusion routines assume ideal gas for diffusion, as do the analytics that they test against. They agree very nicely, but the ideal gas EOS is not the correct thing to use in a WD, so this test isn’t very useful.

If you want to do diffusion in a WD, I strongly recommend going through the release notes for 9793: https://sourceforge.net/p/mesa/mailman/message/35869775/ <https://sourceforge.net/p/mesa/mailman/message/35869775/>

If you still have convergence issues after that, let me know and I can send more suggestions.

Cheers,
Evan


> On Oct 2, 2017, at 11:42 AM, Zhang, Zhaoxi M. (Michael) <mz at astro.caltech.edu> wrote:
> 
> Hello everyone,
> 
> In MESA 10000, I ran wd_diffusion in the test suites.  It worked, but its timesteps were almost 100 million years, which I thought was far higher than the diffusion timescale.    So I decreased the timesteps to see how the elements diffuse initially:
> 
> max_years_for_timestep = 1d3
> 
> and the test stopped working.  It struggles to converge, and fails with:
> 
> stopping because of problems dt < min_timestep_limit
> terminated evolution: cannot find acceptable model
> 
> Is it expected behaviour for this test to fail due to smaller timesteps?  Why would having a smaller timestep create convergence problems?
> 
> 
> Cheers,
> 
> Michael

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20171002/4b9128fd/attachment.html>