- Index
- » Programming
- » Modelica Language
- » error in thr output of openmodelica
error in thr output of openmodelica
error in thr output of openmodelica
model betondechanvrehome
parameter Real rho=1000;//material density
parameter Real rhosurface=413;//material density
parameter Real Rhoambiant=0.007;//material density
parameter Real Cp=1000;//specific heat capacity
parameter Real L=0.1;//length of the wall
parameter Real lambda=0.11;//heat conductivity
parameter Real Tinit=293.15;//the wall is at 20 degree initialy
parameter Real Text=293.15;//the wall is at 20 degree initialy
parameter Real Tint=293.15;//the wall is at 20 degree initialy
parameter Real Thi=293.15;//the wall is at 20 degree initialy
parameter Real tetainit=0.4;//the wall is at 40% final volume specifique d'humidite
parameter Real tetaext=0.4;//the wall is at 40% final volume specifique d'humidite
parameter Real Tfinal=303.15;//the wall is at 30 degree final
parameter Real tetafinal=0.7;//the wall is at 70% final volume specifique d'humidite
parameter Integer N=10;
parameter Real deltax=L/N;//Discretisation
parameter Real Dteta=1.16e-9;
parameter Real Dtemp=1.02e-12;
parameter Real Dtetav=1.07e-9;
parameter Real Dtempv=1.02e-12;
parameter Real Lv=1000;//specific heat capacity
parameter Real hMassique=0.003;//coefficient d'echange massique
parameter Real coefficientThermique=3;//coefficient d'echange thermique
Real T[N];
Real teta[N];
initial equation//initialise the state of variable of the module
for i in 1:N-1 loop
T[i]=Tinit;
teta[i]=tetainit;
end for;
equation
-rho*(T[1]*Dtemp+Dteta*teta[1])=0;
-rho*(T[N]*Dtemp+Dteta*teta[N])=hMassique*(rhosurface-Rhoambiant);
-lambda*T[1]-Lv*rho*(Dtempv*T[1]+Dtetav*teta[1])=0;
-lambda*T[N]-Lv*rho*(Dtempv*T[N]+Dtetav*teta[N])=coefficientThermique*(Tinit-Tfinal)+Lv*hMassique*(rhosurface-Rhoambiant);
for i in 2:N-1 loop
der(teta[i])=Dtemp*(T[i+1]-2*T[i]+T[i-1])/(L/N)^2 +Dteta*(teta[i+1]-2*teta[i]+teta[i-1])/(L/N)^2;
rhosurface*Cp*der(T[i])=lambda*(T[i+1]-2*T[i]+T[i-1])/(L/N)^2 +Lv*Dtempv*(T[i+1]-2*T[i]+T[i-1])/(L/N)^2 + Dteta*(teta[i+1]-2*teta[i]+teta[i-1])/(L/N)^2;
end for;
end betondechanvrehome;
the output error is
[1] 11:28:24 Translation Notification
It was not possible to analyze the given system symbolically, because the relevant equations are part of an algebraic loop. This is not supported yet.
[2] 11:28:24 Translation Warning
Assuming redundant initial conditions for the following 2 initial equations:
T[1] = Tinit
teta[1] = tetainit
what is the problem? can anyone help?
- Index
- » Programming
- » Modelica Language
- » error in thr output of openmodelica