Ut+fx=0 в 3d

[Noctum 0]

Господа!

Никогда не писал солверы для неструктурированных сеток - всегда делал преобразование координат. Но вот тут понадобилось.

Помогите дураку, д.ф-м.н. Есть ДУЧП конвекции вида:

dU/dT+CxdU/dx+CydU/dy+CzdU/dz=0, или, в короткой форме dU/dt+dFi/dxi=0, с корректными НУ и ГУ.

Надо его решит на сетке тетраедров.

Геометрия тетраеров, объем, площади сторон, внешние нормали и их составляющие вычислены.

Пишу следующую функцию для расчета Fi с помощью простого метода Годунова (просто проверить работоспособность):

double compressible_flow::calculate_advection()

{

double ret=0.0;

// for each side

for(int i=0;i<4;i++){

double Ux=0.0;

double Uy=0.0;

double Uz=0.0;

double Cx=-0.5;

double Cy=0.0;

double Cz=0.0;

if(cosX*Cx>0)

Ux=cv[number].U; // Cos... - unit normals, that cosX^2+cosY^2+cosZ^2=1

else

Ux=U_t;

if(cosY*Cy>0)

Uy=cv[number].U;

else

Uy=U_t;

if(cosZ*Cz>0)

Uz=cv[number].U;

else

Uz=U_t;

ret=ret+S*( Cx*Ux*cosX+Cy*Uy*cosY+Cz*Uz*cosZ ); // S - face squre, i - side No.

}

return ret;

}

Получаю ерунду! Что посоветуете. Так-же пробовал схему Roe (ARP) - все равно та же фигня.

Спасибо!