Создайте модель PDE, импортируйте геометрию и установите начальное условие к 50 на целой геометрии.

model = createpde();
importGeometry(model,'BracketWithHole.stl');
setInitialConditions(model,50);

Установите различные начальные условия для каждого компонента системы УЧП.

Создайте модель PDE для системы с пятью компонентами. Импортируйте геометрию Block.stl.

model = createpde(5);
importGeometry(model,'Block.stl');

Установите начальные условия для каждого компонента к дважды номеру компонента.

u0 = [2:2:10]';
setInitialConditions(model,u0)

ans = 
  GeometricInitialConditions with properties:

           RegionType: 'cell'
             RegionID: 1
         InitialValue: [5x1 double]
    InitialDerivative: []

Установите различные начальные условия на каждом фрагменте L-образной мембранной геометрии.

Создайте модель, установите функцию геометрии и просмотрите метки субдомена.

model = createpde();
geometryFromEdges(model,@lshapeg);
pdegplot(model,'FaceLabels','on')
axis equal
ylim([-1.1,1.1])

Установите субдомен 1 на начальное значение-1, субдомен 2 к начальному значению 1 и субдомен 3 к начальному значению 5.

setInitialConditions(model,-1);
setInitialConditions(model,1,'Face',2);
setInitialConditions(model,5,'Face',3);

Начальная установка применяется к целой геометрии. Последующие настройки заменяют начальные установки для областей 2 и 3.

Установите начальные условия для L-образной мембранной геометрии, чтобы быть $$x^2+y^2$$, кроме нижнего левого квадрата, где это $$x^2-y^4$$.

model = createpde();
geometryFromEdges(model,@lshapeg);
pdegplot(model,'FaceLabels','on')
axis equal
ylim([-1.1,1.1])

Установите начальные условия к $$x^2+y^2$$.

initfun = @(location)location.x.^2 + location.y.^2;
setInitialConditions(model,initfun);

Установите начальные условия на области 2 к $$x^2-y^4$$. Эта установка заменяет первую установку, потому что вы применяете его после первой установки.

initfun2 = @(location)location.x.^2 - location.y.^4;
setInitialConditions(model,initfun2,'Face',2);

Гиперболические уравнения имеют ненулевой коэффициент m, таким образом, необходимо установить обоих аргументы u0 и ut0.

Импортируйте Block.stl к модели PDE с N = 3 компонента.

model = createpde(3);
importGeometry(model,'Block.stl');

Установите начальное значение условия, чтобы быть 0 для всех компонентов. Установите начальную производную.

Чтобы создать этот начальный градиент, запишите файл функции и гарантируйте, что функция находится на вашем пути MATLAB.

function ut0 = ut0fun(location)

M = length(location.x);

ut0 = zeros(3,M);

denom = location.x.^2+location.y.^2+location.z.^2;

ut0(1,:) = 4 + location.x./denom;

ut0(2,:) = 5 - tanh(location.z);

ut0(3,:) = 10*location.y./denom;

end

Установите начальные условия.

setInitialConditions(model,0,@ut0fun)

ans = 
  GeometricInitialConditions with properties:

           RegionType: 'cell'
             RegionID: 1
         InitialValue: 0
    InitialDerivative: @ut0fun

Установите начальные условия с помощью решения от предыдущего анализа той же геометрии и mesh.

Создайте и просмотрите геометрию: квадрат с круговым субдоменом.

% Square centered at (1,1), circle centered at (1.5,0.5).
rect1 = [3;4;0;2;2;0;0;0;2;2];
circ1 = [1;1.5;.75;0.25];
% Append extra zeros to the circle;
circ1 = [circ1;zeros(length(rect1)-length(circ1),1)];
gd = [rect1,circ1];
ns = char('rect1','circ1');
ns = ns';
sf = 'rect1+circ1';
[dl,bt] = decsg(gd,sf,ns);
pdegplot(dl,'EdgeLabels','on','FaceLabels','on')
axis equal
ylim([-0.1,2.1])

Включайте геометрию в модель PDE, границу множества и начальные условия, и задайте коэффициенты.

model = createpde();
geometryFromEdges(model,dl);

% Set boundary conditions that the upper and left edges are at temperature 10.
applyBoundaryCondition(model,'dirichlet','Edge',[2,3],'u',10);

% Set initial conditions that the square region is at temperature 0,
% and the circle is at temperature 100.
setInitialConditions(model,0);
setInitialConditions(model,100,'Face',2);

specifyCoefficients(model,'m',0,...
                          'd',1,...
                          'c',1,...
                          'a',0,...
                          'f',0);

Решите проблему в течение многих времен 0 через 1/2 на шагах 0,01.

generateMesh(model,'Hmax',0.05);
tlist = 0:0.01:0.5;
results = solvepde(model,tlist);

Постройте решение в течение многих времен 0.02, 0.04, 0.1, и 0.5.

sol = results.NodalSolution;

subplot(2,2,1)
pdeplot(model,'XYData',sol(:,3))
title('Time 0.02')
subplot(2,2,2)
pdeplot(model,'XYData',sol(:,5))
title('Time 0.04')
subplot(2,2,3)
pdeplot(model,'XYData',sol(:,11))
title('Time 0.1')
subplot(2,2,4)
pdeplot(model,'XYData',sol(:,51))
title('Time 0.5')

Теперь, возобновите анализ и решите проблему в течение многих времен от 1/2 до 1. Используйте ранее полученное решение в течение времени 1/2 как начальное условие. Поскольку 1/2 является последним элементом в tlist, вы не должны задавать индекс времени решения. По умолчанию setInitialConditions использует последний индекс решения.

setInitialConditions(model,results)

ans = 
  NodalInitialConditions with properties:

         InitialValue: [7289x1 double]
    InitialDerivative: []

Решите проблему в течение многих времен 1/2 до 1 на шагах 0,01.

tlist1 = 0.5:0.01:1.0;
results1 = solvepde(model,tlist1);

Постройте решение в течение многих времен 0.5, 0.7, 0.9, и 1.

sol1 = results1.NodalSolution;

figure

subplot(2,2,1)
pdeplot(model,'XYData',sol1(:,1))
title('Time 0.5')
subplot(2,2,2)
pdeplot(model,'XYData',sol1(:,21))
title('Time 0.7')
subplot(2,2,3)
pdeplot(model,'XYData',sol1(:,41))
title('Time 0.9')
subplot(2,2,4)
pdeplot(model,'XYData',sol1(:,51))
title('Time 1.0')

Чтобы использовать ранее полученное решение в течение конкретного времени решения вместо последнего, задайте индекс времени решения как третий параметр setInitialConditions. Например, используйте решение во время 0.2, который является 21-м элементом в tlist.

setInitialConditions(model,results,21)

ans = 
  NodalInitialConditions with properties:

         InitialValue: [7289x1 double]
    InitialDerivative: []

Решите проблему в течение многих времен 0.2 до 1 на шагах 0,01.

tlist2 = 0.2:0.01:1.0;
results2 = solvepde(model,tlist2);