Scribd
Importer
5 vues4 pages

DM 2

Le document présente des données sur la variation horaire de l'irradiance solaire globale dans trois villes algériennes. Il inclut également un code MATLAB pour calculer la température du vitrage et le coefficient global de pertes d'énergie, ainsi que pour déterminer la surface requise pour un capteur solaire en fonction de la puissance demandée. Les résultats montrent les performances théoriques du capteur en fonction des conditions environnementales.

Transféré par

ahmedsprtae
Copyright
© All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats DOCX, PDF, TXT ou lisez en ligne sur Scribd
5 vues4 pages

DM 2

Le document présente des données sur la variation horaire de l'irradiance solaire globale dans trois villes algériennes. Il inclut également un code MATLAB pour calculer la température du vitrage et le coefficient global de pertes d'énergie, ainsi que pour déterminer la surface requise pour un capteur solaire en fonction de la puissance demandée. Les résultats montrent les performances théoriques du capteur en fonction des conditions environnementales.

Transféré par

ahmedsprtae
Copyright
© All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats DOCX, PDF, TXT ou lisez en ligne sur Scribd

Ministère du l’enseignement supérieure et de l’éducation

Ecole Nationale Polytechnique


Departement de génie mécanique

DM 02

Etudiants :
Professeur :
HANI AHMED AMINE Mr. KHALDI

Module : RENEWABLE ENERGY


Année universitaire : 2025/2026
Table des valeurs :

hourly variation of global solar irradiance w/m^2


alger batna Tamanrasset
heure 21/12 21/12 21/12
9:00 101 142 211
10:00 255 296 351
11:00 381 412 461
12:00 458 475 525
13:00 479 478 530
14:00 434 421 460
15:00 335 310 347
16:00 191 158 208
17:00 36 14 67
Total solar 2670 2706 3160
radiation: wh/m2 wh/m2 wh/m2

Le code matlab :

clc; clear; close all

g = 9.81;
sigma = 5.67e-8;

% Températures
Tp = 373;
Ta = 283;
Tsky = Ta;

% Propriétés
alpha_p = 0.97;
eps_p = 0.20;
tau_c = 0.75;
eps_c = 0.70;

% Géométrie
b = 0.025; % Espacement plaque-vitre (m)
L = 0.05; % Isolation arrière (m)
k_ins = 0.045; % Conductivité isolation (W/m.K)

% Air
k_air = 0.029; % W/m.K
nu = 1.96e-5; % m2/s
Pr = 0.7;

% Convection externe
hw = 8; % W/m2.K (vent faible)

%% ==================================================
% QUESTION 1

% Initialisation température du vitrage


Tc = 308; % 35°C
tol = 1e-3;
err = 1;

while err > tol

% --- Rayonnement plaque -> vitrage


hr_pc = sigma*(Tp^2 + Tc^2)*(Tp + Tc) / ...
(1/eps_p + 1/eps_c - 1);

% --- Rayonnement vitrage -> ciel


hr_ca = eps_c*sigma*(Tc^2 + Tsky^2)*(Tc + Tsky);

% --- Convection plaque -> vitrage


Tm = (Tp + Tc)/2;
Ra = g*(Tp-Tc)*b^3/(nu^2*Tm)*Pr;

if Ra < 1e5
Nu = 1 + 1.44*(1 - 1708/Ra);
else
Nu = 0.54*Ra^(1/4);
end

h_pc = Nu*k_air/b;

% --- Coefficient de pertes par le haut


Ut = 1 / ( 1/(h_pc+hr_pc) + 1/(hw+hr_ca) );

% --- Mise à jour température vitrage


Tc_new = Tp - Ut*(Tp-Ta)/(h_pc+hr_pc);

err = abs(Tc_new - Tc);


Tc = Tc_new;
end

% --- Pertes arrière


U_back = k_ins/L;

% --- Coefficient global


U = Ut + U_back;

fprintf('QUESTION 1\n');
fprintf('Température du vitrage Tc = %.1f °C\n', Tc-273);
fprintf('Coefficient global de pertes U = %.2f W/m2.K\n\n', U);

%% ==================================================
% QUESTION 2

% Rayonnement solaire (W/m2)


I = [300 400 500 600 700 800];

% Variable réduite
x = (Tp - Ta)./I;

% Rendement théorique
eta = tau_c*alpha_p - U*x;
% Ajustement polynomial (η = a0 - a1 x - a2 x^2)
p = polyfit(x, eta, 2);

a0 = p(3);
a1 = -p(2);
a2 = -p(1);

fprintf('QUESTION 2\n');
fprintf('a0 = %.3f\n', a0);
fprintf('a1 = %.3f\n', a1);
fprintf('a2 = %.3f\n\n', a2);

% Tracé
xx = linspace(min(x), max(x), 100);
yy = polyval(p, xx);

figure
plot(x, eta, 'o', xx, yy, 'LineWidth',2)
grid on
xlabel('(T_p - T_a)/I')
ylabel('\eta')
title('Courbe de performance du capteur')

%% ==================================================
% QUESTION 3

Q_req = 10e3; % Puissance demandée (W)

% Rayonnement moyen le 21 décembre


sites = {'Alger','Batna','Tamanrasset'};
I_sites = [450 380 520];

fprintf('QUESTION 3\n');

for i = 1:3

eta_moy = tau_c*alpha_p - U*(Tp-Ta)/I_sites(i);

if eta_moy <= 0
fprintf('%s : fonctionnement impossible à 100°C\n', sites{i});
else
Ac = Q_req / (eta_moy * I_sites(i));
fprintf('%s : Surface requise = %.1f m2\n', sites{i}, Ac);
end
end

Vous aimerez peut-être aussi