UDF во Fluent

[Nez 0]

Люди, а есть у когонибудь описание как писать udf'ы во Fluent?

Конкретно - нужно вставить в уравнения навье-стокса дополнительную силу, зависящую от градиента температуры. Если кто может прямо так подсказать - буду очень благодарен!)

[chupal 0]

Люди, а есть у когонибудь описание как писать udf'ы во Fluent?

Конкретно - нужно вставить в уравнения навье-стокса дополнительную силу, зависящую от градиента температуры. Если кто может прямо так подсказать - буду очень благодарен!)

Смотрите UDF Manual. Например <noindex>здесь</noindex> можно посмотреть пример задания источника, а <noindex>здесь</noindex> посмотреть информацию по функции, возвращающей градиент температуры.

[Nez 0]

Прчитал, написал - не работает...

при попытке итерировать выдает

Error:

FLUENT received fatal signal (ACCESS_VIOLATION)

1. Note exact events leading to error.

2. Save case/data under new name.

3. Exit program and restart to continue.

4. Report error to your distributor.

Error Object: ()

Вот мой udf:

#include "udf.h"

#define E 1e6

#define CONSTN 0.146e-3



DEFINE_SOURCE(electric_force_x, c, t, dS, eqn)

{

real x[ND_ND];

real eforce_x;

real deps_x;

real drho_x = C_R_G(c,t)[0];



deps_x = -2*CONSTN*(drho_x/(CONSTN*C_R(c,t)-1))-((-2*CONSTN*C_R(c,t)-1)/((CONSTN*C_R(c,t)-1)*(CONSTN*C_R(c,t)-1)))*CONSTN*drho_x;



eforce_x = 0.5*E*E*deps_x;



return eforce_x;

}

Не понимаю, в чем причина.. в мануале написано про то что нужно использовать coupled солвер при C_R_G - не помогло (...

[chupal 0]

Прчитал, написал - не работает...

при попытке итерировать выдает

Error:

FLUENT received fatal signal (ACCESS_VIOLATION)

1. Note exact events leading to error.

2. Save case/data under new name.

3. Exit program and restart to continue.

4. Report error to your distributor.

Error Object: ()

Вот мой udf:

#include "udf.h"

#define E 1e6

#define CONSTN 0.146e-3



DEFINE_SOURCE(electric_force_x, c, t, dS, eqn)

{

real x[ND_ND];

real eforce_x;

real deps_x;

real drho_x = C_R_G(c,t)[0];



deps_x = -2*CONSTN*(drho_x/(CONSTN*C_R(c,t)-1))-((-2*CONSTN*C_R(c,t)-1)/((CONSTN*C_R(c,t)-1)*(CONSTN*C_R(c,t)-1)))*CONSTN*drho_x;



eforce_x = 0.5*E*E*deps_x;



return eforce_x;

}

Не понимаю, в чем причина.. в мануале написано про то что нужно использовать coupled солвер при C_R_G - не помогло (...

Во первых, по-видимому, не стоит использовать C_R_G для pressure based solver. А можно попробовать выразить градиент плотности через градиенты давления и температуры (если, например у Вас идеальный газ).

Во вторых, сначала загрузите решение, затем дайте команду solve/set/expert и ответьте yes на вопрос Keep temporary solver memory from being freed? (на остальные вопросы жмите enter, оставляя по-умолчанию). После этого сделайте хотя бы одну итерацию не подключая udf. После того, как итерация будет сделана, подключите udf и запустите расчет снова. Думаю, что это поможет.

[Nez 0]

chupal, спасибо, за совет, помогло выражение плотности через т-ру. Ошибки не выдает, но теперь пишет что т-ра превысила 5000 К.. задача - свободная конвекция жидкого азота =( Где я напортачил... И это при том, что никаких источников тепловыделения нет, есть только объемные силы...

chupal, спасибо, за совет, помогло выражение плотности через т-ру. Ошибки не выдает, но теперь пишет что т-ра превысила 5000 К.. задача - свободная конвекция жидкого азота =( Где я напортачил... И это при том, что никаких источников тепловыделения нет, есть только объемные силы...

Сори, все нормально!!! Спасибо!