% Relation HQ
%
% Auteurs : Thibault Prévost.
% NOMA(s)   :0525 1500.
% Date    : 06/05/2021

clc;
clear all;
%close all;

DebitOpenfoam=[1 1 1 1 1 1];
HauteurOpenFoam=[2.1633 2.1786 2.1921 2.2071 2.21 2.4087];
HauteurAval= [1.614 1.8196 1.9126 2.001 2.0491 2.4087];

%Dimensionnement
H0=0.5;
p=3*H0;

%valeurs théoriques
mu0=0.496;
%creager y/H0 = 0.47(x/H0)^(1.8)
L=1;
g=9.81;


%% Calcul valeurs expérimentales de mu

for i=1:length(DebitOpenfoam)

muExpe(i)= DebitOpenfoam(i)/(L * sqrt(2*g) * (HauteurOpenFoam(i)-p)^(3/2));

end

%calcul du mu0 expériementale
mu0Expe=0;

for i=2:length(DebitOpenfoam)
    if HauteurOpenFoam(i)>2 && mu0Expe==0
        delta=(2-HauteurOpenFoam(i-1))/(HauteurOpenFoam(i)-HauteurOpenFoam(i-1));
        DebitDimensionnement=DebitOpenfoam(i-1)+delta*(DebitOpenfoam(i)-DebitOpenfoam(i-1));
        mu0Expe = DebitDimensionnement/(L * sqrt(2*g) * (2-p)^(3/2))
    else
    end    
end

%mu calculé pour le cas de base (sans influence aval)

muCreager = 0.432;

Abscisses = (HauteurOpenFoam-HauteurAval)./(HauteurOpenFoam-p)
Ordonnees = muExpe ./muCreager


figure

plot(Abscisses,Ordonnees,'ro-')
xlabel(' h_{d} /H')
ylabel('\xi = \mu_{Ressault}/ \mu_0')
legend('évolution du coefficient de débis en focntion de la hauteur aval','Location','southeast')
hold off