Добавьте граничное условие в PDEModel контейнер
applyBoundaryCondition( добавляет граничное условие Дирихле к model,'dirichlet',RegionType,RegionID,Name,Value)model. Граничное условие применяется к граничным областям типа RegionType с идентификационными номерами в RegionID, и с аргументами rHU, EquationIndex заданный в Name,Value пары. Для граничных условий Дирихле задайте любые оба аргументы r и h, или аргумент u. При определении u, можно также использовать EquationIndex.
applyBoundaryCondition( добавляет Нейманово граничное условие к model,'neumann',RegionType,RegionID,Name,Value)model. Граничное условие применяется к граничным областям типа RegionType с идентификационными номерами в RegionID, и со значениями g и q заданный в Name,Value пары.
applyBoundaryCondition( добавляет отдельное граничное условие для каждого уравнения в системе УЧП. Граничное условие применяется к граничным областям типа model,'mixed',RegionType,RegionID,Name,Value)RegionType с идентификационными номерами в RegionID, и со значениями, заданными в Name,Value пары. Для смешанных граничных условий можно использовать Name,Value пары и от Дирихле и от Неймановых граничных условий по мере необходимости.
bc = applyBoundaryCondition(___)
Создайте модель PDE и геометрию.
model = createpde(1); R1 = [3,4,-1,1,1,-1,-.4,-.4,.4,.4]'; g = decsg(R1); geometryFromEdges(model,g);
Просмотрите метки ребра.
pdegplot(model,'EdgeLabels','on') xlim([-1.2,1.2]) axis equal
Примените нуль условие Дирихле на ребро 1.
applyBoundaryCondition(model,'dirichlet','Edge',1,'u',0);
На других ребрах примените условие Дирихле h*u = r, где h = 1 и r = 1.
applyBoundaryCondition(model,'dirichlet','Edge',2:4,'r',1,'h',1);
Создайте модель PDE и геометрию.
model = createpde(2); R1 = [3,4,-1,1,1,-1,-.4,-.4,.4,.4]'; g = decsg(R1); geometryFromEdges(model,g);
Просмотрите метки ребра.
pdegplot(model,'EdgeLabels','on') xlim([-1.2,1.2]) axis equal
Примените следующие Неймановы граничные условия на ребро 4.
applyBoundaryCondition(model,'neumann','Edge',4,'g',[0;.123],'q',[0;0;0;0]);
Примените оба типа граничных условий к скалярной проблеме. Во-первых, создайте модель PDE и импортируйте простую геометрию блока.
model = createpde;
importGeometry(model,'Block.stl');Просмотрите метки поверхности.
pdegplot(model,'FaceLabels','on','FaceAlpha',0.5)
Установите нуль условия Дирихле на узких поверхностях, которые помечены 1 through 4.
applyBoundaryCondition(model,'dirichlet','Face',1:4,'u',0);
Установленные Неймановы граничные условия с противоположным входят в систему, стоит 5 и 6.
applyBoundaryCondition(model,'neumann','Face',5,'g',1); applyBoundaryCondition(model,'neumann','Face',6,'g',-1);
Решите эллиптический УЧП с этими граничными условиями и постройте результат.
specifyCoefficients(model,'m',0,'d',0,'c',1,'a',0,'f',0); generateMesh(model); results = solvepde(model); u = results.NodalSolution; pdeplot3D(model,'ColorMapData',u)
Создайте модель PDE и импортируйте простую геометрию блока.
model = createpde(3);
importGeometry(model,'Block.stl');Просмотрите метки поверхности.
pdegplot(model,'FaceLabels','on','FaceAlpha',0.5)
Установите нуль условия Дирихле на поверхностях 1 и 2.
applyBoundaryCondition(model,'dirichlet','Face',1:2,'u',[0,0,0]);
Установленные Неймановы граничные условия с противоположным входят в систему, стоит 4, 5, и 6.
applyBoundaryCondition(model,'neumann','Face',4:5,'g',[1;1;1]); applyBoundaryCondition(model,'neumann','Face',6,'g',[-1;-1;-1]);
Для поверхности 3, примените обобщенное Нейманово граничное условие для первого уравнения и граничные условия Дирихле для вторых и третьих уравнений.
h = [0 0 0;0 1 0;0 0 1]; r = [0;3;3]; q = [1 0 0;0 0 0;0 0 0]; g = [10;0;0]; applyBoundaryCondition(model,'mixed','Face',3,'h',h,'r',r,'g',g,'q',q);
Решите эллиптический УЧП с этими граничными условиями и постройте результат.
specifyCoefficients(model,'m',0,'d',0,'c',1,'a',0,'f',[0;0;0]); generateMesh(model); results = solvepde(model); u = results.NodalSolution; pdeplot3D(model,'ColorMapData',u(:,1))
Когда существует несколько присвоений граничного условия на ту же геометрическую область, тулбокс использует последнюю прикладную установку.
Чтобы постараться не присваивать граничные условия неправильной области, гарантируйте, что вы используете правильные геометрические идентификаторы области путем графического вывода и визуально осмотра геометрии.
Если вы не задаете граничное условие для ребра или поверхности, значением по умолчанию является Нейманово граничное условие с нулевыми значениями для 'g' и 'q'.
Свойства BoundaryCondition | findBoundaryConditions | PDEModel