Discussion :

[Gerris]

Bonjour à tous et à toutes ,
Je sollicite vos compétences pour m'aider à résoudre un problème que j'ai eu depuis quelques jours !! en effet je devrais calculer la température d'un disque tournant dans une boite (température T1) rempli d'un mélange (air,hélium) à l'aide d'un code Matlab(j'ai utilisé un schéma de résistance thermique), le problème c'est que la densité du mélange dépend de la température dans l'enceinte . pourriez-vous m'aider ?
Merci d'avance!
Voilà le code:

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
  clear all
  close all
  T_air=288;                                        % température de l'air à 15°
  p=101325;                                         % pression en pascal
  rho=0.169;                                        % la densité de l'air
  mu=1.9388e-5 ;                                     % la viscosité dynamique de l'air
  % M=[0.029;0.004];                                  % Masse molaire respectivement de l'air et de l'hélium
  %mu=[1.802e-5;1.9388e-5];                          % viscosité dynamique respectivement de l'air et de l'hélium
  Ma=0.029;                                                     %Masse molaire de l'air
  Mhe=0.004;                                              % Masse molaire de l'hélium
  R=0.25;                                          % rayon total du disque en m
  N=10000;                                          % nombre de tours (en rpm)
  omega=2*pi*N/60;
  e= 0.010 ;                                        % epaisseur du disque  en m
  %delta=[0.01:0.005:0.03];                             %vecteur air gap
  delta=0.01;                                       % air gap
  L=0.01+2*R+0.01;                                  % Longueur de la boîte
  l=e+2*delta;                                      % largeur de la boîte
  E=2*R+0.02;                                           % Epaisseur de l'enceinte
  G=delta/R;                                        % air gap ratio
  k=0.02476;                                        % Conductivité thermique de l'air à 15°
  P_bearing=0.5*omega*0.05*(295/2)*0.015;           %Pertes produites par les roulements
  Q=P_bearing+900;                                  % Chaleur dégagée par le moteur+roulements
 
 
 
  %Introduire la densité du mélange
Xa = input('  \n entrer la fraction de l''air : \n')
 Xhe = input('\n entrer la fraction de l''hélium : \n ')
 
if Xa+Xhe==1 
    rho_m=(Xa*Ma+Xhe*Mhe)*p/(8.314*T1)     % densité du mélange(air+hélium)
else 
    input('erreur')
 
end
X=[Xa;Xhe];
%Re_theta=rho_m*omega*R^2/mu;
g1= 0.2091;
g2=1.2139;
%%%% Construction de la matrice %%%% (Appelée C au lieu de G!!!!!)
 for i=1 : 2
    for j=1 :2
        C=zeros;
        C(1,1)=(1/g1)+(1/g2);
        C(2,2)=-C(1,1);
    end
end
G1=inv(C);
T1=Q/C(1,1)+T_air
T3= T1-Ra*Q

[le fab]

Salut

on résout ce genre de problème à l'aide d'une boucle de convergence
tu pars d'une hypothèse de température
puis dans une boucle while tu recalcules densité et température
et tu sors de ta boucles lorsque tu as suffisement convergé (en général une condition du type valeur absolue de température - température au cycle précèdent plus petite qu'une certaine valeur)