[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[ff3d-users] "NOT IMPLEMENTED" problem when I used "characteristic funct
From: |
Robert Li |
Subject: |
[ff3d-users] "NOT IMPLEMENTED" problem when I used "characteristic function" |
Date: |
Sun, 1 Feb 2004 09:44:44 -0800 (PST) |
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/
- [ff3d-users] "NOT IMPLEMENTED" problem when I used "characteristic function",
Robert Li <=