[mame]

Hi to everyone.
My Ansys get me this error
chip-exec: calcolo: wrong return type: void udf function expected
Error: received a fatal signal (Segmentation fault).
Error: received a fatal signal (Segmentation fault).
Error Object: #f
I try this program in Matlab and in C, and it is ok...Ansys in my problem
Can you help me??

#include "udf.h"
/*COSTANTI FISICHE*/
#define Kb 1.3806488e-23 /*costante di Boltzmann [J/K]*/
#define q 1.60217657e-19 /*carica elettrone [C]*/
#define sigma 5.67044e-8 /*costante di Stefan-Boltzmann [W/m2K4]*/

/*DATI di TARGA MODULO*/
#define Voc_stc 22.8
#define Isc 4.73
#define Im_stc 4.45
#define Vm_stc 17.8
#define Ncelle 40
#define alfa_Isc 0.038

#define T_STC 25
#define size_Tc 36 // rendo la dimensione di tc una COSTANTE che posso modificare a piacimento
/*DATI NECESSARI PER IL MODULO*/
#define aa 4.73e-4
#define bb 636
#define Eg0 1.21
#define PASSI 2000

/*CONDIZIONI AMBIENTALI*/
#define G 1000
#define Ta 32

/*PARAMETRI TEORICI*/
#define mCella 1.092024
#define Kph 3.284889e-1
#define Rs 1.022002e-02 /*resistenza serie*/
#define Rsh 200
/*TRASMITTANZE*/
#define etaPC 0.93
#define etaTPU 0.90
#define etaARC 0.85
/* PROTOTIPI FUNZIONI*/
real CalcoloElettrico (real Tc,int m, real I0_Stc,real Iph0, real A);
////////////////////////////////////////////////////////////////////////////////////////////////
// SVILUPPO FUNZIONI//
///////////////////////////*Funzione CALCOLO ELETTRICO*//////////////////////////////////////

real CalcoloElettrico(real Tc,int m, real I0_Stc, real Iph0,real A){

int colonna, dim_Tc,ii, ii1, ii2,i,iv;
real I[PASSI]=0, V[PASSI]=0, Step=0, csi[PASSI]=0, delta=0,Eg=0, Eg1=0, Appoggio=0, Iph=0, Potenza[PASSI]=0, PMax=0, potenza=0, esp_corrente=0, moltip=0, Voc=0, Vm=0,I0=0;
real Pmaxcella=0, etaN=0;

Voc=Voc_stc/Ncelle;
Vm=Vm_stc/Ncelle;
//Inizializzazione dei vettori corrente I e tensione V
for (colonna=0;colonna<PASSI; colonna++)
I[colonna]=0;

for (colonna=0;colonna<PASSI; colonna++)
V[PASSI]=0;

potenza=pow(Tc,2);
Eg=Eg0-((potenza*aa)/(Tc+bb));
Eg1=q*Eg;

I0=I0_Stc*(pow(Tc/298,3))*exp(-Eg1/Kb/Tc)/exp(-Eg1/Kb/298);
Iph=Iph0*(1+alfa_Isc/100*(Tc-298.15))*G/1000;/*1000 irradianza standard*/

/*linsapce come si fa*/
Step=((0.98*Voc)-(Rs*Iph))/(PASSI-1);
csi[0]=Rs*Iph;
csi[PASSI-1]=0.98*Voc;
for(i=1;i<PASSI-1;i++){
csi[i]=csi[i-1]+Step;
}

delta=q/(mCella*Kb*Tc);
for(iv=0;iv<PASSI;iv++){
esp_corrente=delta*csi[iv];
moltip=I0*(exp(esp_corrente)-1);
I[iv]=Iph-moltip-(csi[iv]/Rsh);
V[iv]= csi[iv]-(Rs*I[iv]);
}
for(ii1=0;ii1<i;ii1++){
Potenza[ii1]=I[ii1]*V[ii1];
}

for(ii2=1;ii2<i;ii2++){
if (Potenza[ii2]>=Potenza[ii2-1])
PMax=Potenza[ii2];
}
Pmaxcella=PMax;
etaN=Pmaxcella/(G*A);

return etaN;
}

/////////////
////////MAIN/////////
////////////
DEFINE_SOURCE(calcolo,c, t, dS, eqn)
{
int m=0,iniz=0;
real eta_STC=0, Vm=0, A=0, P_STC=0, I0_Stc=0, Iph0=0, Voc=0, T_STC_K=0, Ta_K=0, etaN=0, denominatore_esp=0, arg_esp=0, generazione=0;
real temp=0;
cell_t c;

P_STC=Im_stc*Vm_stc;/* lo posso anche metter come DEFINE*/
/*Dati geometrici celle*/
A=12*12;
A=A*pow(10,-4);
Iph0=A*Kph*1000;/*1000 è l' irradianza GStc*/
Voc=Voc_stc/Ncelle;
Vm=Vm_stc/Ncelle;

/*rendimento condizioni standard*/
eta_STC= Vm_stc*Im_stc/(1000*A);
T_STC_K=T_STC+273.15;/*conversione in K*/
Ta_K=Ta+273.15;/*conversione in K, CONTROLLARE IN ANSYS*/
denominatore_esp=mCella*Kb*(25+273.15);
arg_esp=q*Voc/denominatore_esp;
I0_Stc=(Iph0-Voc/Rsh)/exp(arg_esp);
/* calcolo elettrico*/
temp=C_T(c,t);
etaN=CalcoloElettrico(temp, m, I0_Stc, Iph0, A);
generazione=((etaPC*etaTPU*etaARC)-etaN)*G;
return generazione;
}

[Bruno Machado]

Hi mame,

The code you wrote is not in the same language that is used in Fluent. Have a look at the UDF Manual so you can make the changes accordingly.

[mame]

i used thr udf manuel and i used the DEFINE_ Source and a function. When i compiled or interperted the udf i have no problem.
i also change the udf :

#include "udf.h"

/*COSTANTI FISICHE*/
#define Kb 1.3806488e-23 /*costante di Boltzmann [J/K]*/
#define q 1.60217657e-19 /*carica elettrone [C]*/
#define sigma 5.67044e-8 /*costante di Stefan-Boltzmann [W/m2K4]*/

/*DATI di TARGA MODULO*/
#define Voc_stc 22.8
#define Isc 4.73
#define Im_stc 4.45
#define Vm_stc 17.8
#define Ncelle 40
#define alfa_Isc 0.038

#define T_STC 25
#define size_Tc 36 // rendo la dimensione di tc una COSTANTE che posso modificare a piacimento
/*DATI NECESSARI PER IL MODULO*/
#define aa 4.73e-4
#define bb 636
#define Eg0 1.21
#define PASSI 2000

/*CONDIZIONI AMBIENTALI*/
#define G 1000
#define Ta 32

/*PARAMETRI TEORICI*/
#define mCella 1.092024
#define Kph 3.284889e-1
#define Rs 1.022002e-02 /*resistenza serie*/
#define Rsh 200
/*TRASMITTANZE*/
#define etaPC 0.93
#define etaTPU 0.90
#define etaARC 0.85

/////////////
////////MAIN/////////
////////////
DEFINE_SOURCE(calcolo,c, t, dS, eqn)
{
int m=0,iniz=0;
real eta_STC=0, Vm=0, A=0, P_STC=0, I0_Stc=0, Iph0=0, Voc=0, T_STC_K=0, Ta_K=0, etaN=0, denominatore_esp=0, arg_esp=0, generazione=0;
real temp=0;
int colonna, dim_Tc,ii, ii1, ii2,i,iv;
real I[PASSI], V[PASSI], Step=0, csi[PASSI], delta=0,Eg=0, Eg1=0, Appoggio=0, Iph=0, Potenza[PASSI], PMax=0, potenza=0, esp_corrente=0, moltip=0,I0=0;
real Pmaxcella=0;
//cell_t c;

P_STC=Im_stc*Vm_stc;/* lo posso anche metter come DEFINE*/
/*Dati geometrici celle*/
A=12*12;
A=A*pow(10,-4);
Iph0=A*Kph*1000;/*1000 è l' irradianza GStc*/
Voc=Voc_stc/Ncelle;
Vm=Vm_stc/Ncelle;

/*rendimento condizioni standard*/
eta_STC= Vm_stc*Im_stc/(1000*A);
T_STC_K=T_STC+273.15;/*conversione in K*/
Ta_K=Ta+273.15;/*conversione in K, CONTROLLARE IN ANSYS*/
denominatore_esp=mCella*Kb*(25+273.15);
arg_esp=q*Voc/denominatore_esp;
I0_Stc=(Iph0-Voc/Rsh)/exp(arg_esp);

/* calcolo elettrico*/
temp=C_T(c,t);

//FUNZIONE
Voc=Voc_stc/Ncelle;
Vm=Vm_stc/Ncelle;
//Inizializzazione dei vettori corrente I e tensione V
for (colonna=0;colonna<PASSI; colonna++)
I[colonna]=0;

for (colonna=0;colonna<PASSI; colonna++)
V[PASSI]=0;

potenza=pow(temp,2);
Eg=Eg0-((potenza*aa)/(temp+bb));
Eg1=q*Eg;

I0=I0_Stc*(pow(temp/298,3))*exp(-Eg1/Kb/temp)/exp(-Eg1/Kb/298);
Iph=Iph0*(1+alfa_Isc/100*(temp-298.15))*G/1000;/*1000 irradianza standard*/

/*linsapce come si fa*/
Step=((0.98*Voc)-(Rs*Iph))/(PASSI-1);
csi[0]=Rs*Iph;
csi[PASSI-1]=0.98*Voc;
for(i=1;i<PASSI-1;i++){
csi[i]=csi[i-1]+Step;
}

delta=q/(mCella*Kb*temp);
for(iv=0;iv<PASSI;iv++){
esp_corrente=delta*csi[iv];
moltip=I0*(exp(esp_corrente)-1);
I[iv]=Iph-moltip-(csi[iv]/Rsh);
V[iv]= csi[iv]-(Rs*I[iv]);
}
for(ii1=0;ii1<i;ii1++){
Potenza[ii1]=I[ii1]*V[ii1];
}

for(ii2=1;ii2<i;ii2++){
if (Potenza[ii2]>=Potenza[ii2-1])
PMax=Potenza[ii2];
}
Pmaxcella=PMax;
etaN=Pmaxcella/(G*A);

generazione=((etaPC*etaTPU*etaARC)-etaN)*G;
return generazione;
}

I hope you understand that even if the variable names are in Italian.
maybe i must use tha define_profile or define specific heat??

[pakk]

Does it work without errors if you replace

temp=C_T(c,t);

by

temp=300;

?

The reason I ask is because your error message indicates that you are trying to access information that does not exist, and the only place I quickly see in your code where you try to access information is there.

If the error disappears, your problem is that the temperature is not defined everywhere. If the problem stays, it is elsewhere.

[mame]

thanks for your replay.
I try this....maybe it want a "begin_c_loop" although DEFINE SOURCE however access to the cells
My problem is that the ansys calculate the temperature of every cell, use the udf for calculate the efficienty and then calculate the heat generation for every cell