Hi Niraj, You need to be more specific with your question. Do you get the wrong result? Or do you get an error message? In what way does the solution differ? Can you direct us to a plot of your solution compared to the one you are expecting? Etc.
-- Anders On Fri, May 24, 2013 at 11:38:59AM +0200, Niraj Jha wrote: > 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 _______________________________________________ fenics mailing list [email protected] http://fenicsproject.org/mailman/listinfo/fenics
