[ff3d-users] Memory error

[Robert Li]

Dear Stephane Del Pino,
  I met a memory error when I tried to test my
example.
I need both plane boundary and whole box boundary.
Therefore, I defined the domain, which is made of six
planes, is as same as box. 
 Can I use this method to get both plane boundary and
whole box boundary?

Robert,

>>>>output<<<<
oops: Something strange happened: incorrect mesh:
        One of the faces of your mesh is related to 4
cells (vertices: )


<<<<<<<<pov
//layer1
plane  { <-1,0,0>, 2 pigment { color rgb <1,0,0> }}
plane  { <1,0,0>, 2 pigment { color rgb <0,1,0> }}
plane  { <0,-1,0>, 2 pigment { color rgb <0,0,1> }}
plane  { <0,1,0>, 2 pigment { color rgb <1,1,0> }}
plane  { <0,0,-1>, 2 pigment { color rgb <1,0,1> }}
plane  { <0,0,1>, 2 pigment { color rgb <0,1,1> }}


box { <-2, -2 , -2>, < 2, 2, 2> pigment { color rgb
<1,1,1> } }

//layer2
plane  { <-1,0,0>, 2 pigment { color rgb <0.9,0,0> }}
plane  { <1,0,0>, 2 pigment { color rgb <0,0.9,0> }}
plane  { <0,-1,0>, 2 pigment { color rgb <0,0,0.9> }}
plane  { <0,1,0>, 2 pigment { color rgb <0.9,0.9,0> }}
plane  { <0,0,-1>, -2 pigment { color rgb <0.9,0,0.9>
}}
plane  { <0,0,1>, 5 pigment { color rgb <0,0.9,0.9> }}

box { <-2, -2 , 2>, < 2, 2, 5> pigment { color rgb
<0.9,0.9,0.9> } }


<<<<<<<<<

//This example is  dpendent for time t(transient
state)
//dt(H) !=0:
//k=1920000.0 density=3000000.0  Cap=0.85 
//g!=0:
//materials=1:

vector n = (20,20,20);
vector a = (-3,-3,-2);
vector b = (3,3,5);

scene S = pov("Heat5.pov"); // the pov-ray file for
the geometry

mesh M = structured(n,a,b);

domain O = domain(S, (inside(<1,0,0>) and
inside(<0,1,0>) and inside(<0,0,1>) and
                                inside(<1,1,0>) and inside(<1,0,1>) and
inside(<0,1,1>) ) or
                                (inside(<0.9,0,0>) and inside(<0,0.9,0>) and
inside(<0,0,0.9>) and
                                inside(<0.9,0.9,0>) and inside(<0.9,0,0.9>) and
inside(<0,0.9,0.9>)) or
                                inside(<1,1,1>) or
                                inside(<0.9,0.9,0.9>)
);

mesh Omesh = surface(O,M);
mesh s1  = extract(Omesh, 1);   //xmin of bottom box
mesh s2  = extract(Omesh, 2);   //xmax of bottom box
mesh s3  = extract(Omesh, 3);   //ymin of bottom box
mesh s4  = extract(Omesh, 4);   //ymax of bottom box
mesh s5  = extract(Omesh, 5);   //zmin of bottom box
mesh s6  = extract(Omesh, 6);   //zmax of bottom box
mesh s7  = extract(Omesh, 7);   //xmin of middle box
mesh s8  = extract(Omesh, 8);   //xmax of middle box
mesh s9  = extract(Omesh, 9);   //ymin of middle box
mesh s10 = extract(Omesh, 10);  //ymax of middle box
mesh s11 = extract(Omesh, 11);  //zmin of middle box
mesh s12 = extract(Omesh, 12);  //zmax of middle box

femfunction Tn(M) = 0;  //T0; when t=0; T
double  i = 0;                  //loop times
double  dt = 24*60*60*20;               //time step 1 day
double  Cu=4184000.0;   //unfrozen volumetric heat
capacity(1000*4184)
double  Cf=1932000.0;   //frozen volumetric heat
capacity (920*2100)
double  Tpc=0.0;                //phase change temperature

function        Kx=2.2*one(<1,1,1>) +
2.9*one(<0.9,0.9,0.9>);                 //2.2*one(Tn<Tpc) +
0.57*one(Tn>=Tpc);      //conductivity in x direction
function        Ky=2.2;                 //2.2*one(Tn<Tpc) +
0.57*one(Tn>=Tpc);      //conductivity in y direction
function        Kz=2.2;                 //2.2*one(Tn<Tpc) +
0.57*one(Tn>=Tpc);      //conductivity in z direction
function        Cp=1365000.0;   //Cf*one(Tn<Tpc)  +
Cu*one(Tn>=Tpc);                //volumetric heat capacity
double          Lf=334000000.0;                                                 
        //Latent heat of
fusion of water
double          Theat=0.8;                                                      
                //Normalized volumetric
water content(0-1)
function        Theatu1=0.04;   //*one(Tn<=-2) +
(0.5*Tn+1.05)*one(Tn<=-0.5 and Tn>-2) +
                                                        // (1.6*Tn)*one(Tn<=0 
and Tn>-0.5) + 0 *
one(Tn>0);
function        factor=dt/(Cp+(Lf)*(Theat)*(Theatu1));
do {
         solve(T) in O by M {
                 pde(T)
                        T-(
dx(((Kx)*(factor))*dx(T))+dy(((Ky)*(factor))*dy(T))+dz(((Kz)*(factor))*dz(T))
) = Tn;
//                      T = -10 on s12;
//                  T = 0 on s11;
                        T = -30 on s12;
//                  T = -10 on s11;
//                  T = -10 on s6;
                    T = 0 on s5;

         }

        //Save result 
        save (medit , "heat5.00".i , T, M, dos ) ;
        save (medit , "heat5.00".i,  M, dos ) ;

        // compute deltaT here
        Tn = T;
        i=i+1;

} while (i<=10);

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