Dear All,
This is my code and unfortunately, I am getting the same error in integral line of problem.
Kindly suggest me to fix it,
codes are given below.
/////////////////////////////////////////////////////////////////////////////////////////////////////
real dt=0.50;
real tend =10;
real phio= 0.1003, rhof=1000.0, cvf = 4200, kts=2.49, cvs=1000.0, rhos=2910.2;
real qf=0.0, qs=0.0, cpf=4200, ktf=0.6, g=9.8, miu=0.2;
real kt= (1-phio)kts+phioktf;
real RhoCv= (1-phio)rhoscvs+ phiorhofcvf;
border C1(t=0,1) { x = 20t; y= 0; }
border C2(t=0,1) { x = 20; y = 20t; }
border C3(t=0,1) { x = 20 - 20t; y =20; }
border C4(t=0,1) { x = 0; y = 20 - 20t;}
mesh Th= buildmesh(C1(20)+C2(20)+C3(20)+C4(20));
fespace Vh(Th, P1);
Vh u,v;
Vh uold,f;
/macro grad(u) (dx(u)+dy (u)) /
problem prs(u,v) = int2d(Th) (
dtkt(dx(u)dx(v)+dy(u)dy(v)) + RhoCv(uv))
+ int2d(Th) (dt*cpf*rhof*qf*v*(dx(u)+dy(v)))
- int2d(Th) (dt*qs*v)
- int2d(Th) (RhoCv*uold*v)
+ on(C2, u=f)
+ on(C4, u=473.15);
uold=0;
real t=0;
for(int i=0; i<50; i++)
{ uold =u;
t=t+dt;
f= 473.15-10*t;
prs;
plot(u, wait=0,fill=1, value=1);
}