boutproject
WIP: FV::Div_a_Grad_perp operator
Going through this operator (again). Reviewing how metric elements are interpolated to cell edges.
To do:
- [ ] Tests
- [ ] Documentation
Thanks for this @bendudson. I've put this through our puff test. For everyone else, we are facing a very tough to debug problem where our neutrals are spontaneously increasing in temperature and decreasing in density around the X-point. We have ruled out a lot of things and our current suspect is the operator. We came up with a test where we take a simplified model with a fixed plasma background of constant temperature and density and impose a neutral source near the X-point. We then observe how the diffusive transport behaves, particularly how the density falloff evolves with distance from the puff source.
I am seeing no difference with this new operator. There is a bifurcation near the source which is identical in both cases. The bifurcation can be seen near X=0 in the RHS plot we see several points of similar distance but different density. It may be that this is due to cell size after all, i.e. because the cell value represents the average across the entire cell, meaning that a larger cell covering a greater distance from the puff source would have a lower density value (I think?).
Old results:
New results:
