Homec4science
Diffusion GYACOMO (OBSOLETE: LATEST VERSION ON gitlab.epfl.ch) (half_kx_plane) Browse

gyacomo/matlab/flux_results.m
85911ff7c9f2half_kx_plane

Last Change
Raw File

flux_results.m
View Options

default_plots_options
if 1
%% Compute time average and std of the mean flow
t0 = 1000; t1 = 3000; [~,it0] = min(abs(t0-Ts2D)); [~,it1] = min(abs(t1-Ts2D));
range = it0:it1;
avg = mean(GFlux_ri(range))
stdev = std(GFlux_ri(range))^(.5)
figure
hist(GFlux_ri(range),20); xlabel('$\Gamma$')
end
if 0
%% Handwritten results for nu = 0.01
% High definition results (256x128)
SCALING = 2*sqrt(2);
Results_256x128.Gamma = [0.02, 0.03, 0.20, 0.037, 2.7, 2.25, 4, 5e-4, 2e-3, 0.03];
Results_256x128.L = [ 66, 66, 66, 50, 66, 66, 66, 66, 66, 66];
Results_256x128.P = [ 2, 3, 4, 5, 2, 3, 4, 2, 3, 4];
Results_256x128.J = [ 1, 2, 2, 3, 1, 2, 2, 1, 2, 2];
Results_256x128.etaB = [ 0.5, 0.5, 0.5, 0.5, 0.4, 0.4, 0.4, 0.6, 0.6, 0.6];
Results_256x128.mrkx = [ 'v', '>', '^', 'o', 'v', '>', '^', 'v', '>', '^'];
Results_256x128.clr = [ 'k', 'k', 'k', 'r', 'r', 'r', 'r', 'k', 'k', 'k'];
% Low definition results (128x64)
% Results_128x64.Gamma = [0.29, 0.05, 7e-4, 0.31, 3.7, 2e-3];
% Results_128x64.L = [ 25, 25, 25, 33 33, 33];
% Results_128x64.P = [ 2, 2, 2, 2, 2, 2];
% Results_128x64.J = [ 1, 1, 1, 1, 1, 1];
% Results_128x64.NU = [0.01, 0.1, 0.01, 0.01, 0.01, 0.01];
% Results_128x64.etaB = [ 0.5, 0.5, 0.67, 0.5 0.4, 0.6];
% Results_128x64.mrkx = [ 's', 's', 's', 's', 's', 's'];
% Results_128x64.clr = [ 'b', 'b', 'b', 'r', 'r', 'r'];
% Ricci_Rogers.Gamma = [2.5 1 1e-2];
% Ricci_Rogers.etaB = [0.4 0.5 1.0];
Ricci_Rogers.Gamma = [10 1e-1];
Ricci_Rogers.etaB = [0.5 1.0];
if 1
% Fig 3 of Ricci Rogers 2006
fig = figure;
semilogy(Ricci_Rogers.etaB,Ricci_Rogers.Gamma,'--','color',[0,0,0]+0.6);
hold on;
res = Results_256x128;
for i = 1:numel(res.Gamma)
semilogy(res.etaB(i),res.Gamma(i)*SCALING,...
res.mrkx(i),'DisplayName','256x128', 'color', res.clr(i));
hold on;
end
% res = Results_128x64;
% for i = 1:numel(res.Gamma)
% if res.NU(i) == 0.01
% semilogy(res.etaB(i),res.Gamma(i),...
% res.mrkx(i),'DisplayName','128x64', 'color', res.clr(i));
% end
% hold on;
% end
xlabel('$\eta_B$'); ylabel('$\Gamma^\infty_{part}$')
end
grid on; title('$\nu = 0.01$')
legend('Mix. Length, Ricci 2006','$P=2$, $J=1$','$P=3$, $J=2$','$P=4$, $J=2$','$P=5$, $J=3$')
FIGNAME = [SIMDIR,'flux_study_nu_1e-2.png']; ylim([1e-4, 10])
saveas(fig,FIGNAME);
disp(['Figure saved @ : ',FIGNAME])
end
if 0
%% Handwritten results for nu = 0.1
Results_256x128.Gamma = [0.026,0.026, 1e-2, 1, 1, 1, 2e-2, 1, 0.15, 3e-3];
Results_256x128.P = [ 2, 3, 4, 2, 3, 4, 2, 3, 4, 4];
Results_256x128.J = [ 1, 2, 2, 1, 2, 2, 1, 2, 2, 2];
Results_256x128.etaB = [ 0.5, 0.5, 0.5, 0.4, 0.4, 0.4, 0.6, 0.6, 0.6, 0.7];
Results_256x128.mrkx = [ 'v', '>', '^', 'v', '>', '^', 'v', '>', '^', '^'];
Results_256x128.clr = [ 'k', 'k', 'k', 'b', 'b', 'b', 'r', 'b', 'r', 'r'];
% Ricci_Rogers.Gamma = [2 1e-1];
% Ricci_Rogers.etaB = [0.5 1.0];
Ricci_Rogers.Gamma = [10 1e-1];
Ricci_Rogers.etaB = [0.5 1.0];
if 1
SCALING = 2*sqrt(2);
% Fig 3 of Ricci Rogers 2006
fig = figure;
semilogy(Ricci_Rogers.etaB,Ricci_Rogers.Gamma,'--','color',[0,0,0]+0.6);
hold on;
res = Results_256x128;
for i = 1:numel(res.Gamma)
semilogy(res.etaB(i),res.Gamma(i)*SCALING,...
res.mrkx(i),'DisplayName','256x128', 'color', res.clr(i));
hold on;
end
xlabel('$\eta_B$'); ylabel('$\Gamma^\infty_{part}$')
end
grid on; title('$\nu = 0.1$')
legend('Mix. Length, Ricci 2006','$P=2$, $J=1$','$P=3$, $J=2$','$P=4$, $J=2$')
FIGNAME = [SIMDIR,'flux_study_nu_1e-1.png'];
saveas(fig,FIGNAME);
disp(['Figure saved @ : ',FIGNAME])
end
if 0
%% Handwritten results for nu = 1.0, eta = 0.5, 200x100, L=100
Results_256x128.Gamma = [0.026,0.026, 1e-2, 1, 1, 1, 2e-2, 1, 0.15, 3e-3];
Results_256x128.P = [ 2, 3, 4, 2, 3, 4, 2, 3, 4, 4];
Results_256x128.J = [ 1, 2, 2, 1, 2, 2, 1, 2, 2, 2];
Results_256x128.etaB = [ 0.5, 0.5, 0.5, 0.4, 0.4, 0.4, 0.6, 0.6, 0.6, 0.7];
Results_256x128.mrkx = [ 'v', '>', '^', 'v', '>', '^', 'v', '>', '^', '^'];
Results_256x128.clr = [ 'k', 'k', 'k', 'b', 'b', 'b', 'r', 'b', 'r', 'r'];
% Ricci_Rogers.Gamma = [2 1e-1];
% Ricci_Rogers.etaB = [0.5 1.0];
Ricci_Rogers.Gamma = [10 1e-1];
Ricci_Rogers.etaB = [0.5 1.0];
if 1
% Fig 3 of Ricci Rogers 2006
SCALING = 2*sqrt(2);
fig = figure;
semilogy(Ricci_Rogers.etaB,Ricci_Rogers.Gamma,'--','color',[0,0,0]+0.6);
hold on;
res = Results_256x128;
for i = 1:numel(res.Gamma)
semilogy(res.etaB(i),res.Gamma(i)*SCALING,...
res.mrkx(i),'DisplayName','256x128', 'color', res.clr(i));
hold on;
end
xlabel('$\eta_B$'); ylabel('$\Gamma^\infty_{part}$')
end
grid on; title('$\nu = 0.1$')
legend('Mix. Length, Ricci 2006','$P=2$, $J=1$','$P=3$, $J=2$','$P=4$, $J=2$')
FIGNAME = [SIMDIR,'flux_study_nu_1e-1.png'];
saveas(fig,FIGNAME);
disp(['Figure saved @ : ',FIGNAME])
end
if 0
%% Handwritten results for nu = 1.0, eta = 0.5, 150x75, L=100, DGGK
Results_150x75.Gamma = [0.026,0.026, 1e-2, 1, 1, 1, 2e-2, 1, 0.15, 3e-3];
Results_150x75.P = [ 2, 3, 4, 2, 3, 4, 2, 3, 4, 4];
Results_150x75.J = [ 1, 2, 2, 1, 2, 2, 1, 2, 2, 2];
Results_150x75.etaB = [ 0.5, 0.5, 0.5, 0.4, 0.4, 0.4, 0.6, 0.6, 0.6, 0.7];
Results_150x75.mrkx = [ 'v', '>', '^', 'v', '>', '^', 'v', '>', '^', '^'];
Results_150x75.clr = [ 'k', 'k', 'k', 'b', 'b', 'b', 'r', 'b', 'r', 'r'];
% Ricci_Rogers.Gamma = [2 1e-1];
% Ricci_Rogers.etaB = [0.5 1.0];
Ricci_Rogers.Gamma = [10 1e-1];
Ricci_Rogers.etaB = [0.5 1.0];
if 1
% Fig 3 of Ricci Rogers 2006
SCALING = 2*sqrt(2);
fig = figure;
semilogy(Ricci_Rogers.etaB,Ricci_Rogers.Gamma,'--','color',[0,0,0]+0.6);
hold on;
res = Results_150x75;
for i = 1:numel(res.Gamma)
semilogy(res.etaB(i),res.Gamma(i)*SCALING,...
res.mrkx(i),'DisplayName','256x128', 'color', res.clr(i));
hold on;
end
xlabel('$\eta_B$'); ylabel('$\Gamma^\infty_{part}$')
end
grid on; title('$\nu = 0.1$')
legend('Mix. Length, Ricci 2006','$P=2$, $J=1$','$P=3$, $J=2$','$P=4$, $J=2$')
FIGNAME = [SIMDIR,'flux_study_nu_1e-1.png'];
saveas(fig,FIGNAME);
disp(['Figure saved @ : ',FIGNAME])
end
c4science · Help