Разница между двухмерной и трёхмерной динамической сеткой.

[gebbels 1]

Добрый день, уважаемые специалисты.

Недавно столкнулся с проблемой,которой не могу нийти объяснение.

Есть задача: корпус дирижабля без оперения и гондолы (на первое время) движется прямолинейно равномерно через возмущённую атмосферу со скоростью 18 м/с. Порыв ветра в атмосфере моделируется в виде инжекции из сопла. Движение прямое, сетка динамическая перестраиваемая.

<noindex>http://narod.ru/disk/33585204001/grid_23112011_2D.msh.html</noindex>

Эта сетка двухмерная, на ней решено было проводить предварительные расчёты с целью экономии времени. Решение на такой сетке сходится без проблем, параметры динамической сетки приведены на рисунках. Максимальный шаг по времени, который я пробовал 0.04 секунды. То есть даже когда за один шаг по времени тело проходит 18*0.04=0.72 метра, что примерно в 3 раза больше характерного размера ячейки, решение не разваливается и отрицательные объёмы не возникают. Считалось по 40 итераций на каждый шаг.

В этом cas файле можно посмотреть сетку, которая получается после 30 секунд расчёта:

<noindex>http://narod.ru/disk/33590247001/grid_2311...000000.cas.html</noindex>

После расчётов на двухмерной сетке я построил трёхмерную сетку.

<noindex>http://narod.ru/disk/33589331001/grid_05122011_3D.msh.html</noindex>

И вот тут я столкнулся с проблемой. При таких же параметрах динамической сетки упорно и настойчиво возникают отрицательные объёмы. При чём уменьшение шага по времени не помогает. Проходит определённое число шагов и расчёт останавливается, даже при шаге по времени 0.001 сек.

Сначала была мысль, что это из-за пирамид, которые возникают на стыке призм и тетраэдров. Но полностью тетраэдральная сетка также не хочет работать. Ссылка:

<noindex>http://narod.ru/disk/33591281001/setka_3D_dyn_1.msh.html</noindex>

То есть пока что не удалось провести ни одного расчёта с именно трёхмерной динамической сеткой. У кого есть опыт в этом вопросе подскажите пожалуйста, какая именно сетка и настройки динамической сетки нужны для именно трёхмерного расчёта?

В принципе ничего критичного в своей сетке я не заметил, может быть можно свести расчёт за счёт изменения параметров. Правильно ли я выставил фактор жёсткости пружины? Судя по туториалу при значении 0 сетка получается более равномерной. Также хотелось бы узнать о применении size function.

[gebbels 1]

Отрицательные объемы можно посмотреть после срыва на самой сетке,

можно будет увидеть разорванную сетку.

Как именно посмотреть? У меня при появлении отрицательного объёма флюент виснет намертво и убить я его могу только из другого терминала. (Я по сетке подключён к удалённой линуксовой машине)

1. А Вы боковые движения дирижабля в UDF-ке исключили

применением команд ??? например:

prop[sDOF_ZERO_TRANS_Y] = TRUE; // Запрет на перемещение по оси OY

prop[sDOF_ZERO_ROT_Z] = TRUE; // Запрет на вращение относительно оси OZ

2. Было бы правильно сначала потренироваться

на обучающем примере "Falling-box" - кубик падает в воду.

С уважением, Герман

Вот моя UDF-ка:

#include"udf.h"

DEFINE_CG_MOTION(derigable, dt, cg_vel, cg_omega, time, dtime)

{

real Velocity; // Назначение типа переменной

Velocity=-18.0; // Назначение величины х-компоненты скорости дирижабля

cg_vel[0] = Velocity; /* x-velocity*/

cg_vel[1] = 0.0; /* y-velocity*/

cg_vel[2] = 0.0; /* z-velocity*/

NV_S(cg_omega, =, 0.0); /* no angular motion */

}

На самом деле причина, видимо, в слишком агрессивном сгущении сетки. На равномерной сетке всё начало работать нормально.