UDF для трехмерной задачи

[The_Immortal 1]

Здравствуйте!

Помогите, пожалуйста, разобраться с формированием UDF-файла для трехмерной задачи: имеется формула (во вложении). Необходимо получить профиль распределения скорости Z, по направлению движения Y.
При попытке сделать расчет получаем сообщение "Error: Divergence detected in AMG solver: x-momentum".

Подскажите, пожалуйста, как сиё поправить?

Спасибо!

[The_Immortal 1]

P.S. UDF:

#include "udf.h"
#include "math.h"
DEFINE_PROFILE(inlet_y_velocity, thread, index)

#define v_mean 4
#define z_max 750.53
#define z_min -249.47
#define B 1./7.

{
    real x[ND_ND]; /* this will hold the position vector */
    real y;
    real z;
    real u;
    face_t f;

    begin_f_loop(f, thread) /*loops over all faces in the thread passed in the DEFINE macro argument*/
    {
        F_CENTROID(x,f,thread);
        z = x[1];
        y = x[1];
        u = v_mean*(B+1);;
        F_PROFILE(f, thread, index) = u*pow((z-249.47)/10,B);
    }
    end_f_loop(f, thread)
}

 

 

[The_Immortal 1]

Господа, подскажите, пожалуйста... Задача конкретно встала :(