[ff3d-users] "NOT IMPLEMENTED" problem when I used "characteristic function"

[Robert Li]

Deat Stephane Del Pino,
  I met "NOT IMPLEMENTED" problem when I used such
"characteristic function" as 
function Cp=Cu*one(Tn>Tpc) +
Cf*one(Tn<=Tpc);//volumetric heat capacity
Could you help me?

Robert


>>>ff3d output
Solving Problem:

        unknowns: List of unknowns

        Problem:
                pde(T):
                1*T -
dz(((1.92e+006*(0.1/(((4.184e+006*boolop)+(1.932e+006*bool
op))+(2.672e+008*((((0.04*boolop)+(((0.5*0)+1.05)*./language\FunctionExpression.
hpp:1104NOT IMPLEMENTED
))+((1.6*0)*./language\FunctionExpression.hpp:1104NOT
IMPLEMENTED
))+(0*boolop))))))*1)*dz(T)) -
dy(((1.92e+006*(0.1/(((4.184e+006*boolop)+(1.932e
+006*boolop))+(2.672e+008*((((0.04*boolop)+(((0.5*0)+1.05)*./language\FunctionEx
pression.hpp:1104NOT IMPLEMENTED
))+((1.6*0)*./language\FunctionExpression.hpp:1104NOT
IMPLEMENTED
))+(0*boolop))))))*1)*dy(T)) -
dx(((1.92e+006*(0.1/(((4.184e+006*boolop)+(1.932e
+006*boolop))+(2.672e+008*((((0.04*boolop)+(((0.5*0)+1.05)*./language\FunctionEx
pression.hpp:1104NOT IMPLEMENTED
))+((1.6*0)*./language\FunctionExpression.hpp:1104NOT
IMPLEMENTED
))+(0*boolop))))))*1)*dx(T)) = 0
                T = -3 on (4)

                T = -10 on (5)



>>>>>>>input file 
vector n = (20,20,20);
vector a = (-2,-2,-2);
vector b = (2,2,2);

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

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



//femfunction Tn(M) = 0;        //T0; when t=0; T
function Tn = 0;                        //T0; when t=0; T
double  i = 0;                  //loop times
double  dt = 0.1;               //time step
double  Kx=1920000.0;   //conductivity in x direction
double  Ky=1920000.0;   //conductivity in y direction
double  Kz=1920000.0;   //conductivity in z direction
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 Cp=Cu*one(Tn>Tpc) +
Cf*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 M {
                 pde(T)
                        T-(
dx(((Kx)*(factor))*dx(T))+dy(((Ky)*(factor))*dy(T))+dz(((Kz)*(factor))*dz(T))
) = Tn;

                        T = -10 on M zmax;  
                    T = -3 on M zmin;
         }

        //Save result 
        save (medit , "heat6.00".i , T, M, dos ) ;
        save (medit , "heat6.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/