On 24 May 2013 11:38:59 +0200
"Niraj Jha" <[email protected]> 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))
As I said on stackexchange, you're missing here x[0]<tol
> #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))
Missing x[1] < tol
> #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))
That's obviously wrong, you need
(((x[0]/3.25)**2 +(x[1]/2.75)**2 ) > 1.0-tol)
instead of
(((x[0]/3.25)**2 +(x[1]/2.75)**2 ) > 1.0+tol)
>
> 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)
Here this should be: (((x[0]/3.25)**2 +(x[1]/2.75)**2 ) < 1.0+tol)
> 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
