Поиск назначения граничного условия для геометрической области
BCregion = findBoundaryConditions(BCs,RegionType,RegionID)BCregion назначается указанной области.
Создайте модель 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]);
На грани 3 задайте граничное условие Неймана для уравнения 1 и граничное условие Дирихле для уравнений 2 и 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);
Задайте граничные условия Неймана с противоположными знаками на гранях 5 и 6 для всех уравнений.
applyBoundaryCondition(model,'neumann','Face',4:5,'g',[1;1;1]); applyBoundaryCondition(model,'neumann','Face',6,'g',[-1;-1;-1]);
Проверьте спецификацию граничного условия на грани 1.
findBoundaryConditions(model.BoundaryConditions,'Face',1)ans =
BoundaryCondition with properties:
BCType: 'dirichlet'
RegionType: 'Face'
RegionID: [1 2]
r: []
h: []
g: []
q: []
u: [0 0 0]
EquationIndex: []
Vectorized: 'off'
Проверьте спецификацию граничного условия на грани 3.
findBoundaryConditions(model.BoundaryConditions,'Face',3)ans =
BoundaryCondition with properties:
BCType: 'mixed'
RegionType: 'Face'
RegionID: 3
r: [3x1 double]
h: [3x3 double]
g: [3x1 double]
q: [3x3 double]
u: []
EquationIndex: []
Vectorized: 'off'
Проверьте спецификацию граничного условия на грани 5.
findBoundaryConditions(model.BoundaryConditions,'Face',5)ans =
BoundaryCondition with properties:
BCType: 'neumann'
RegionType: 'Face'
RegionID: [4 5]
r: []
h: []
g: [3x1 double]
q: []
u: []
EquationIndex: []
Vectorized: 'off'