Этот пример показывает, как решить уравнение Пуассона, -, u = f на 2-D геометрии, созданной как комбинация двух прямоугольников и двух кругов.
Чтобы решить эту проблему в приложении PDE Modeler, выполните следующие шаги:
Откройте приложение PDE Modeler при помощи pdeModeler
команда.
Отображение линий сетки. Для этого выберите Options > Grid Spacing и снимите флажок Auto для линейного интервала по оси X. Введите X-axis linear spacing следующим -1.5:0.25:1.5
. Затем выберите Options > Grid.
Выровнять новые формы по линиям сетки можно путем выбора Options > Snap.
Нарисуйте две окружности: одну с радиусом 0,4 и центр в (-0,5,0) и другую с радиусом 0,2 и центр в (0.5,0,2). Чтобы нарисовать круг, сначала нажмите кнопку. Затем щелкните правой кнопкой мыши источник и перетащите мышью, чтобы нарисовать круг. Щелчок правой кнопкой мыши ограничивает рисуемую форму так, чтобы она была окружностью, а не эллипсом.
Потяните два прямоугольника: один с углами (-1 0.2), (1 0.2), (1,-0.2), и (-1,-0.2) и другой с углами (0.5,1), (1,1), (1,-0.6), и (0.5,-0.6). Чтобы нарисовать прямоугольник, сначала нажмите кнопку. Затем щелкните любой угол и перетащите мышью, чтобы нарисовать прямоугольник.
Моделируйте геометрию путем ввода (R1+C1+R2)-C2
в поле Set formula.
Сохраните модель в файл, выбрав File Save As.
Удалите границы поддомена. Для этого перейдите в граничный режим путем выбора Boundary > Boundary Mode. Затем выберите Boundary > Remove All Subdomain Borders.
Задайте граничные условия для всех дуг окружности. Используя Shift + нажатие кнопки, выберите эти границы. Затем выберите Boundary > Specify Boundary Conditions и задайте граничное условие Неймана с g = -5 и q = 0. Это граничное условие означает, что решение имеет уклон -5 в нормальном направлении для этих граничных сегментов.
Для всех других контуров сохраните условие Дирихле контура умолчанию: h = 1
, r = 0
.
Задайте коэффициенты, выбрав PDE > PDE Specification или нажав кнопку PDE на панели инструментов. Задайте c = 1
, a = 0
, и f = 10
.
Инициализируйте mesh путем выбора Mesh > Initialize Mesh. Уточнить mesh можно путем выбора Mesh > Refine Mesh.
Решить УЧП можно путем выбора Solve > Solve PDE или нажатия кнопки = на панели инструментов. Тулбокс собирает задачу УЧП, решает её и строит график решения.
Постройте график решения как 3-D график:
Выберите Plot > Parameters.
В получившемся диалоговом окне выберите Height (3-D plot).
Нажмите Plot.