Mar 29, 2012 at 2:24pm UTC
this is my code:
for (i=0;i<=n;i++){
nodx[i] = Lx/((n-1)*2)+Lx/(n-1)*(i-1);}
nodx[0]=0;
nodx[n]=Lx;
for (j=0;j<=m;j++){
nody[j] = Ly/((m-1)*2)+Ly/(m-1)*(j-1);}
nody[0]=0;
nody[m]=Ly;
//condicioens de contorno ks
for (j=0;j<=m;j++)
for (i=0;i<=n;i++){
if (i==0){ // convección con fluido
if (nody[j] <= 0.4){
ac[i][j]= k3*(vcly[j]-vcly[j-1])/(nodx[i+1]-nodx[i]);}
else {
ac[i][j]=k1*(vcly[j]-vcly[j-1])/(nodx[i+1]-nodx[i]);}}
else if (j == 0){
if (nodx[i]<= 0.3){
ac[i][j]= k3*(vclx[i]-vclx[i-1])/(nody[j+1]-nody[j]);}
else {
ac[i][j]= k4*(vclx[i]-vclx[i-1])/(nody[j+1]-nody[j]);}}
else if (j==m){
if (nodx[i] <= 0.3){
ac[i][j]= k1*(vclx[i]-vclx[i-1])/(nody[j]-nody[j-1]);}
else {
ac[i][j]= k2*(vclx[i]-vclx[i-1])/(nody[j]-nody[j-1]);}}
else if (i==n){ // convección con fluido
if (nody[j]<= 0.1){
ac[i][j]= k4*(vcly[j]-vcly[j-1])/(nodx[i]-nodx[i-1]);}
else {
ac[i][j]=k2*(vcly[j]-vcly[j-1])/(nodx[i]-nodx[i-1]);}}
else {
ac[i][j]=0;
}}
ac[n][0]=0;
ac[n][m]=0;
ac[0][m]=0;
ac[0][0]=0;
for (i=0;i<=n;i++){
cout << vclx[i]<< " " ;}
cout << "\n ";
for (j=0;j<=m;j++){
cout << vcly[j]<< " " ;}
cout << "\n";
for(i=0;i<=n;i++){
cout << nodx[i] << " " ;}
when I finish, the array nodx[i] is wrong. I gave nodx[0]=0 when I started but when the program run the final value is 0.9... which is imposible because I haven´t done any assignation to nodx[i] after... any suggestions¿¿??