Уравнения, которые вы можете решить используя PDE Toolbox

Partial Differential Equation Toolbox™ решает скалярные уравнения вида

m2ut2+dut·(cu)+au=f

и собственные уравнения вида

·(cu)+au=λduили·(cu)+au=λ2mu

Для скалярных PDE существует два выбора граничных условий для каждого ребра или грани:

  • Дирихлет - На краю или лице решение u удовлетворяет уравнению

    hu = r,

    где h и r могут быть функциями пространства (x, y, и, в 3-D случае, z), u решения и времени. Часто вы принимаете h = 1 и устанавливаете r на соответствующее значение.

  • Обобщенные граничные условия Неймана - На краю или грани решение, u удовлетворяет уравнению

    n·(cu)+qu=g

    n - внешний модуль нормали. q и g - функции, определенные на ∂ Ω, и могут быть функциями x, y, и, в 3D случае, z, решение u, и, для уравнений с временной зависимостью, время.

Тулбокс также решает системы уравнений вида

m2ut2+dut·(cu)+au=f

и собственные системы собственных значений вида

·(cu)+au=λduили·(cu)+au=λ2mu

Система PDE с N компонентами N связана PDE со связанными граничными условиями. Скалярные PDE являются такими, которые имеют N = 1, что означает только один УЧП. Системы PDE обычно означают N > 1. Документация иногда относится к системам как к многомерным PDE или к PDE с вектором u решения. Во всех случаях системы PDE имеют одну геометрию и mesh. Это только N, количество уравнений, которые могут варьироваться.

Коэффициенты m, d, c, a и f могут быть функциями местоположения (x, y, и, в 3-D, z), и, кроме собственного значения задач, они также могут быть функциями решения u или его градиента. Для задач с собственным значением коэффициенты не могут зависеть от решения u или его градиент.

Для скалярных уравнений все коэффициенты, кроме c, скаляра. c коэффициентов представляет матрицу 2 на 2 в 2-D геометрии или матрицу 3 на 3 в 3-D геометрии. Для систем N уравнений коэффициенты m, d и a являются N -by- N матрицами, f является вектором N -by-1, и c является 2 N -by-2 N тензором (2-D геометрия) или 3 N -by-3 N тензором (3-D геометрия). Для значенияcu, см. c Коэффициент для sefectCofficients.

Когда и m, и d 0, УЧП является стационарным. Когда m или d ненулевые, задача зависит от времени. Когда любой коэффициент зависит от u решения или его градиента, задача называется нелинейной.

Для систем PDE существуют обобщённые версии граничных условий Дирихле и Неймана:

  • hu = r представляет матрицу h, умножающую вектор решения u и равную вектору r.

  • n·(cu)+qu=g. Для 2-D систем обозначение n·(cu) означает N -by-1 матрицу с (i, 1) - компонент

    j=1N(cos(α)ci,j,1,1x+cos(α)ci,j,1,2y+sin(α)ci,j,2,1x+sin(α)ci,j,2,2y)uj

    где вектор внешней нормали контура n=(cos(α),sin(α)).

    Для 3-D систем обозначение n·(cu) означает вектор N -by-1 с (i, 1) - компонент

    j=1N(sin(φ)cos(θ)ci,j,1,1x+sin(φ)cos(θ)ci,j,1,2y+sin(φ)cos(θ)ci,j,1,3z)uj+j=1N(sin(φ)sin(θ)ci,j,2,1x+sin(φ)sin(θ)ci,j,2,2y+sin(φ)sin(θ)ci,j,2,3z)uj+j=1N(cos(θ)ci,j,3,1x+cos(θ)ci,j,3,2y+cos(θ)ci,j,3,3z)uj

    где вектор внешней нормали контура n=(sin(φ)cos(θ),sin(φ)sin(θ),cos(φ)).

    Для каждого ребра или грани существует в общей сложности N граничных условий.

Похожие темы