Нагрейте распределение в круговом цилиндрическом стержне: приложение PDE Modeler

Решите 3-D параболическую задачу УЧП путем сокращения проблемы до 2D преобразования координаты использования. Этот пример использует приложение PDE Modeler. Для решения командной строки смотрите Распределение Тепла в Круговом Цилиндрическом Стержне.

Рассмотрите цилиндрический радиоактивный стержень. Тепло постоянно добавляется в левом конце стержня, в то время как правильный конец сохранен при постоянной температуре. Во внешней границе теплом обмениваются со средой переводом. В то же время тепло однородно производится в целом стержне из-за радиоактивных процессов. Предположение, что начальная температура является нулем, приводит к следующему уравнению:

ρCut·(ku)=q

Здесь, ρ, C и k являются плотностью, тепловой способностью и теплопроводностью материала, u является температурой, и q является теплом, выработанным в стержне.

Поскольку проблема осесимметрична, удобно написать это уравнение в системе цилиндрической координаты.

ρCut1rr(krur)1r2θ(kuθ)z(kuz)=q

Здесь r, θ и  z являются тремя координатными переменными цилиндрической системы. Поскольку проблема осесимметрична, u/θ=0.

Это - цилиндрическая проблема, и Partial Differential Equation Toolbox™ требует, чтобы уравнения были в Декартовых координатах. Чтобы преобразовать уравнение к Декартовым координатам, сначала умножьте обе стороны уравнения r:

ρrCutr(krur)z(kruz)=qr

Затем задайте r как y и z как x:

ρyCut·(kyu)=qy

В данном примере используйте эти параметры:

  • Плотность, ρ = 7 800 кг/м3

  • Тепловая способность, C = 500 Вт · s/kg · ºC

  • Теплопроводность, k = 40 W/mºC

  • Радиоактивный источник тепла, q = 20 000 Вт/м3

  • Температура в правильном конце, T_right = 100 ºC

  • Нагрейте поток в левом конце, HF_left = 5 000 Вт/м2

  • Окружающая температура во внешней границе, T_outer = 100 ºC

  • Коэффициент теплопередачи, h_outer = 50 Вт/м2· ºC

Чтобы решить эту задачу в приложении PDE Modeler, выполните эти шаги:

  1. Смоделируйте стержень как прямоугольник с углами в (-1.5,0), (1.5,0), (1.5 0.2), и (-1.5 0.2). Здесь, x - ось представляет направление z и y - ось представляет направление r.

    pderect([-1.5,1.5,0,0.2])
  2. Задайте граничные условия. Для этого дважды кликните контуры, чтобы открыть диалоговое окно Boundary Condition. Приложение PDE Modeler требует граничных условий в конкретной форме. Таким образом Неймановы граничные условия должны быть в форме n·(cu)+qu=g, и граничные условия Дирихле должны быть в форме h u = r. Кроме того, потому что обе стороны уравнения умножаются на r = y, умножают коэффициенты для граничных условий y.

    • Для левого конца используйте Нейманово условие n·(ku)=HF_left=5000. Задайте g = 5000*y и q = 0.

    • Для правильного конца используйте условие Дирихле u = T_right = 100. Задайте h = 1 и r = 100.

    • Для внешней границы используйте Нейманово условие n·(ku)=h_outer(T_outeru)=50(100u). Задайте g = 50*y*100 и q = 50*y.

    • Цилиндрический r оси  = 0 не является контуром в исходной проблеме, но в 2D обработке это стало той. Используйте искусственное Нейманово граничное условие для оси, n·(ku)=0. Задайте g = 0 и q = 0.

  3. Задайте коэффициенты путем выбора PDE> PDE Specification или нажмите кнопку PDE на панели инструментов. Нагрейтесь уравнение является параболическим уравнением, так выберите тип Parabolic УЧП. Поскольку обе стороны уравнения умножаются на r = y, умножают коэффициенты на y и вводят следующие значения: c = 40*y, a = 0, f = 20000*y, и d = 7800*500*y.

  4. Инициализируйте mesh путем выбора Mesh> Initialize Mesh.

  5. Установите начальное значение к 0, время решения к 20 000 секунд, и вычисляйте решение каждые 100 секунд. Для этого выберите Solve> Parameters. В диалоговом окне Solve Parameters, время набора к 0:100:20000, и u (t 0) к 0.

  6. Решите уравнение путем выбора Solve> Solve PDE или нажатия кнопки = на панели инструментов.

  7. Постройте решение, с помощью цвета и контурного графика. Для этого выберите Plot> Parameters и выберите цвет и контурные графики в получившемся диалоговом окне.

    Temperature distribution plot with the isothermal lines

Можно исследовать решение путем варьирования параметров модели и графического вывода результатов. Например, вы можете:

  • Покажите решение, когда u не будет зависеть вовремя, то есть, решение для устойчивого состояния. Для этого откройте диалоговое окно PDE Specification и измените тип УЧП в Elliptic. Получившееся решение для устойчивого состояния находится в близком соглашении с переходным решением в 20 000 секунд.

  • Покажите решение для устойчивого состояния, не охлаждаясь на внешней границе: коэффициент теплопередачи является нулем. Для этого установите Нейманово граничное условие во внешней границе (главная сторона прямоугольника) к g = 0 и q = 0. Получившийся график показывает, что температура повышается до больше чем 2 500 на левом конце стержня.

Для просмотра документации необходимо авторизоваться на сайте