Re: [ff3d-users] Using complex geometry....error using "cylinder"??

[Arun Viswanathan]

Hey Stephane, thank you for your reply. I have
attached the 2 files -- p2.ff and p2.pov -- at the end
of this mail. Besides getting a syntax error (for
"cylinder") in the povray file, another issue I am
unsure of is using the "intersection" function to
define a container consisting of 12 drums in it.

Also, incase any of the drums leaks, the chemical
would diffuse into different parts of the container,
moving around the other drums....I am not sure what
sort of boundary conditions I should give for this
sort of motion?!

Please give me any suggestions, the code is given
below. Thanks, best regards -- Arun.

******************p2.ff (below)********************
// Aim: To solve Fick's law and determine the 
// concentration of the gas escaping from the drum(s)
// in various parts of the container.

function D = 10;
function J = 20;
function C0 = x*y*z;

double dt = 0.1;
double i;

vector n = (10,10,10);
vector a = (0,0,0);
vector b = (2,2,6);
mesh M = structured(n,a,b);
scene S = pov("p2.pov");
domain O = domain(S);

for (i=0; i<5; i++)
{       solve(C) in O by M
        {
                pde(C)
                        C - div(D*dt*grad(C)) = C0;
/*** Equation: dC/dt = D * d/dx(dC/dx) ***/
        };
//      cout << "C = " << C << "\n";
        C0 = C;
}

save(opendx, "p1.dat", C, M);

******************p2.pov (below)*******************

intersection {
        box {
                <0,0,0>, <2,2,6>
        }
        pigment { color rgb <1,0,0> }

        cylinder {
                <0.5,0,0.5>, <0.5,1.4,0.5>, 0.4
        }
        pigment { color rgb <0,1,0> }

        cylinder {
                <1.5,0,0.5>, <1.5,1.4,0.5>, 0.4
        }
        pigment { color rgb <0,1,0> }

        cylinder {
                <0.5,0,1.5>, <0.5,1.4,1.5>, 0.4
        }
        pigment { color rgb <0,1,0>

        cylinder {
                <1.5,0,1.5>, <1.5,1.4,1.5>, 0.4
        }
        pigment { color rgb <0,1,0> }

        cylinder {
                <0.5,0,2.5>, <0.5,1.4,2.5>, 0.4
        }
        pigment { color rgb <0,1,0> }

        cylinder {
                <1.5,0,2.5>, <1.5,1.4,2.5>, 0.4
        }
        pigment { color rgb <0,1,0> }

        cylinder {
                <0.5,0,3.5>, <0.5,1.4,3.5>, 0.4
        }
        pigment { color rgb <0,1,0> }

        cylinder {
                <1.5,0,3.5>, <1.5,1.4,3.5>, 0.4
        }
        pigment { color rgb <0,1,0> }

        cylinder {
                <0.5,0,4.5>, <0.5,1.4,4.5>, 0.4
        }
        pigment { color rgb <0,1,0> }

        cylinder {
                <1.5,0,4.5>, <1.5,1.4,4.5>, 0.4
        }
        pigment { color rgb <0,1,0> }

        cylinder {
                <0.5,0,5.5>, <0.5,1.4,5.5>, 0.4
        }
        pigment { color rgb <0,1,0> }

        cylinder {
                <1.5,0,5.5>, <1.5,1.4,5.5>, 0.4
        }
        pigment { color rgb <0,1,0> }
}

****************************************************

--- Stephane Del Pino <address@hidden> wrote:
> Dear Arun.
> 
> Could you send your .pov and .ff files so that we
> can help you to solve
> your problem?
> 
> Regards,
> Stephane.
> 
> PS: send them to the list, I will not read this
> address tomorrow.
> 
> 
> _______________________________________________
> ff3d-users mailing list
> address@hidden
> http://mail.nongnu.org/mailman/listinfo/ff3d-users



        
                
__________________________________
Do you Yahoo!?
Friends.  Fun.  Try the all-new Yahoo! Messenger.
http://messenger.yahoo.com/