ff3d-users
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[ff3d-users] How to use plane?


From: Robert Li
Subject: [ff3d-users] How to use plane?
Date: Mon, 26 Jan 2004 07:02:03 -0800 (PST)

Dear Stephane Del Pino,
  I have a box with vertices  <0,  0,  0> and < 4, 4,
6  >. I tried to get the planes(z=0,2,4,6) of the box.
Then I set the boundary conditions on them. I failed
with the following errors! 
  Could you tell me why No cell has reference that
matches (4), which is bottom plane of the box! Do you
have any example used plane?

Thanks 
Robert


 
.........
Reference translation (POV-Ray -> FEM)
                (1, 0, 0) -> 1
                (0, 1, 0) -> 2
                (0, 0, 1) -> 3
                (1, 1, 0) -> 4
                (0, 1, 1) -> 5
        Done
....
Sets 1 References
No cell has reference that matches (4)!
Press any key to continue


>>>>>>>>>POV file>>>>>>>>>
box { <0,  0,  0>, < 4, 4, 6  > pigment { color rgb
<1,0,0> } }
plane  { <0,0,1>, 0 pigment { color rgb <1,1,0> }}
plane  { <0,0,1>, 2 pigment { color rgb <0,0,1> }}
plane  { <0,0,1>, 4 pigment { color rgb <0,1,0> }}
plane  { <0,0,1>, 6 pigment { color rgb <0,1,1> }}







>>>>>>>>>INPUT file>>>>>>>>>
vector n = (20, 20, 30);
//vertex a = (-2, -2, -5);
//vertex b = (2, 2, 1);
vertex a = (-1, -1, -2);
vertex b = (5, 5, 8);

mesh M = structured(n,a,b);
 
scene S = pov("Tplane.pov");    // the pov-ray file for
the geometry
domain O = domain(S, inside(<1,0,0>) or
inside(<0,1,0>) or inside(<0,0,1>) or inside(<1,1,0>)
or inside(<0,1,1>));

solve(u) in O by M  //FDM  
{
   pde(u)
     - div(grad(u)) = 0;
 u=-1 on <0,1,1>;
 u=0 on <0,1,0>;
 u=1 on <0,0,1>;
 u=2 on <1,1,0>;

};
save(medit,"TT",u,M,dos);
save(medit,"TT",M,dos);

__________________________________
Do you Yahoo!?
Yahoo! SiteBuilder - Free web site building tool. Try it!
http://webhosting.yahoo.com/ps/sb/




reply via email to

[Prev in Thread] Current Thread [Next in Thread]