Discussion :

[Pausedejeuner]

Bonjour à tous,

Je suis en train de modeliser à l'aide de PDE tool la répartition de chaleur dans une pièce rectangulaire avec une source chaude et une source froide modélisées toutes les deux par des conditions aux bords du type T=0°c et T=20°C. J'obtiens alors la répartition de la chaleur en régime stationnaire sur tout mon espace comme sur l'image jointe. Je me demandais alors s'il était possible de poser une condition au bord du type T=constante inconnue pour la source chaude et ainsi de visualiser la température dans mon tout mon espace en fonction de cette constante? Car Matlab m'envoie un message d'erreur quand je rentre en condition au bord pour la source chaude une constante k.

Merci d'avance !

[Gooby]

Bonjour,

Peux-tu nous montrer le code associé?

[Pausedejeuner]

Bonjour Gooby,

Bien sur, le voici:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
 
function pdemodel
[pde_fig,ax]=pdeinit;
pdetool('appl_cb',1);
set(ax,'DataAspectRatio',[1 1 1]);
set(ax,'PlotBoxAspectRatio',[1.5 1 1]);
set(ax,'XLim',[-1.5 1.5]);
set(ax,'YLim',[-1 1]);
set(ax,'XTickMode','auto');
set(ax,'YTickMode','auto');
 
% Geometry description:
pderect([-1.0292792792792791 0.92567567567567544 0.4346846846846848 -0.36711711711711725],'R1');
pderect([-0.98423423423423428 -0.87162162162162149 0.13738738738738743 -0.11036036036036023],'R2');
pderect([0.72747747747747793 0.83558558558558627 0.14189189189189211 -0.10585585585585577],'R3');
set(findobj(get(pde_fig,'Children'),'Tag','PDEEval'),'String','R1-R2-R3')
 
% Boundary conditions:
pdetool('changemode',0)
pdesetbd(12,...
'dir',...
1,...
'1',...
'0')
pdesetbd(11,...
'dir',...
1,...
'1',...
'20')
pdesetbd(10,...
'dir',...
1,...
'1',...
'20')
pdesetbd(9,...
'dir',...
1,...
'1',...
'0')
pdesetbd(8,...
'neu',...
1,...
'0',...
'0')
pdesetbd(7,...
'neu',...
1,...
'0',...
'0')
pdesetbd(6,...
'dir',...
1,...
'1',...
'0')
pdesetbd(5,...
'dir',...
1,...
'1',...
'0')
pdesetbd(4,...
'dir',...
1,...
'1',...
'20')
pdesetbd(3,...
'dir',...
1,...
'1',...
'20')
pdesetbd(2,...
'neu',...
1,...
'0',...
'0')
pdesetbd(1,...
'neu',...
1,...
'0',...
'0')
 
% Mesh generation:
setappdata(pde_fig,'Hgrad',1.3);
setappdata(pde_fig,'refinemethod','regular');
setappdata(pde_fig,'jiggle',char('on','mean',''));
setappdata(pde_fig,'MesherVersion','preR2013a');
pdetool('initmesh')
pdetool('refine')
 
% PDE coefficients:
pdeseteq(1,...
'1.0',...
'0.0',...
'10.0',...
'1.0',...
'0:10',...
'0.0',...
'0.0',...
'[0 100]')
setappdata(pde_fig,'currparam',...
['1.0 ';...
'0.0 ';...
'10.0';...
'1.0 '])
 
% Solve parameters:
setappdata(pde_fig,'solveparam',...
char('0','1926','10','pdeadworst',...
'0.5','longest','0','1E-4','','fixed','Inf'))
 
% Plotflags and user data strings:
setappdata(pde_fig,'plotflags',[1 1 1 1 1 1 1 1 0 0 0 1 1 0 0 0 0 1]);
setappdata(pde_fig,'colstring','');
setappdata(pde_fig,'arrowstring','');
setappdata(pde_fig,'deformstring','');
setappdata(pde_fig,'heightstring','');
 
% Solve PDE:
pdetool('solve')