LinearFree sheath BC extrapolation doesn't work from scratch

Open mikekryjak opened this issue 5 months ago • 2 comments

In sheath_boundary_simple, we support three types of extrapolating the evolved variables to the target:

https://github.com/boutproject/hermes-3/blob/0b564f8cf1300476676a12a21bb6c13da06cc958/src/sheath_boundary_simple.cxx#L124-L134

This is implemented here: https://github.com/boutproject/hermes-3/blob/0b564f8cf1300476676a12a21bb6c13da06cc958/src/sheath_boundary_simple.cxx#L35-L52

I am unable to start a case from scratch with mode = 2, i.e. LinearFree. It complains that I am setting an invalid value for e:energy_source. I am able to proceed by using the default mode = 1 and restarting with mode = 2.

@ZedThree something to keep in mind during your refactor.

Jul 09 '25 09:07 mikekryjak

What does "invalid value" imply -- NaN?

Jul 09 '25 09:07 ZedThree

NaN definitely causes it. I don't know if there are other values that can trigger it as well.

Jul 09 '25 14:07 mikekryjak