Hello Fenics team,
I am trying to solve a problem described here
http://www.caesarsystems.co.uk/NAFEMS_benchmarks/le10.html..I don't have
the exact result what i am looking for??
I tried to check my code many times but there is still error..Any
advance user can help me out with this problem..
tol = 1e-10
ABA'B'
class Left(SubDomain):
def inside(self, x, on_boundary):
return near (x[1],2.75) and between(x[2],(0.0,0.600))
#return between(x[1],(1.0,2.75)) or between(x[2],(0.0,0.600))
DCD'C'
class Right(SubDomain):
def inside(self, x, on_boundary):
return near (x[0],3.25) and between(x[2],(0.0,0.600))
#return between(x[0],(2.0,3.25)) or between(x[2],(0.0,0.600 ))
BCB'C'
class Outerall(SubDomain):
def inside(self, x, on_boundary):
return (((x[0]/3.25)**2 +(x[1]/2.75)**2) - 1.0 > tol) and
between(x[2],(0.0,0.600))
Mid plane of BCB'C'
class Outermid(SubDomain):
def inside(self, x, on_boundary):
return (((x[0]/3.25)**2 +(x[1]/2.75)**2 )-1.0 > tol) and
near(x[2],0.300)
ABCD
class Top(SubDomain):
def inside(self, x, on_boundary):
return near(x[2], 0.600 )
V = VectorFunctionSpace(mesh, "Lagrange", 1)
# Initialize sub-domain instances
left = Left()
right = Right()
outer1 = Outerall()
outer2 = Outermid()
top = Top()
facet_domains = FacetFunction('size_t', mesh, 0)
top.mark(facet_domains, 3)
bc_right = DirichletBC(V.sub(1), Constant(0.0), right, "pointwise")
bc_left = DirichletBC(V.sub(0), Constant(0.0), left, "pointwise")
bc_outer1x = DirichletBC(V.sub(0), Constant(0.0),outer1, "pointwise")
bc_outer1y = DirichletBC(V.sub(1), Constant(0.0),outer1, "pointwise")
bc_outer2x = DirichletBC(V.sub(0), Constant(0.0),outer2, "pointwise")
bc_outer2y = DirichletBC(V.sub(1), Constant(0.0),outer2, "pointwise")
bc_outer2z = DirichletBC(V.sub(2), Constant(0.0),outer2, "pointwise")
bc =
[bc_right,bc_left,bc_outer1x,bc_outer1y,bc_outer2x,bc_outer2y,bc_outer2z]
Thanks
Niraj Kumar Jha
004917634988484
_______________________________________________
fenics mailing list
[email protected]
http://fenicsproject.org/mailman/listinfo/fenics
