19 Apr
2016
19 Apr
'16
3:38 p.m.
On 19 Apr 2016, at 16:35, Onno Bokhove <O.Bokhove@leeds.ac.uk> wrote:
I missed that last week in the code:
the 2nd order equation dzz(phi) needs 2 conditions, chosen here to be the same at the top and bottom; also phi_s at later stages will be one of the unknowns! So how is it possible that the equation was solved at all with only one condition? Did it take Neumann condition at the bottom, as part of an automatic procedure in Firedrake?
A neumann condition just kills one of your surface terms when integrating by parts, no? If you just completely left it out, that's equivalent to choosing u.n = 0 as a neumann condition on that side. Lawrence