Archived OpenModelica forums. Posting is disabled.

Alternative forums include GitHub discussions or StackOverflow (make sure to read the Stack Overflow rules; you need to have well-formed questions)


Forgot password? | Forgot username? | Register

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?

Edited by: lizzi - Mar-12-17 09:40:16
There are 0 guests and 0 other users also viewing this topic
You are here: