`我正在解决这里发现的问题“ http://www.caesarsystems.co.uk/NAFEMS_benchmarks/le10.html ” ..我对边界条件感到非常困惑，并尽我所能计算它，但结果完全是不同的

这是我的bcs`

class Left(SubDomain): def inside(self, x, on_boundary):
    return (between (x[1], (1.0,2.75)) and between (x[2], (0.0,0.600000000))

class Right(SubDomain):
def inside(self, x, on_boundary):
    return (between(x[0], (2.0,3.25)) and between (x[2], (0.0,0.600000000)))  

class Outerall(SubDomain):
def inside(self, x, on_boundary):
    return (((x[0]/3.25)**2 + (x[1]/2.75)**2 ) -1.0 < tol) and abs(x[2]-0.600000000)


class Outermid(SubDomain):
def inside(self, x, on_boundary):
    return  (((x[0]/3.25)**2 + (x[1]/2.75)**2 ) -1.0 < tol) and abs(x[2]-0.300000000) < tol

class Top(SubDomain):
def inside(self, x, on_boundary):
    return abs(x[2]-0.600000000) < tol