% Representation of the experimental dataset for 
% "Wall Repulsion during Electrophoresis: Testing the theory of
% Concentration-Polarization Electroosmosis"

close all;clear;clc
%% Figure 4(a)
% Frequencies used
f = [10e3 5e3 3e3 1.7e3 950 527 292 167 90 50]; % Hz
% Experimental separation measured
y_Exp = [9.7780   10.6237   11.0465   11.4693   12.3150   14.4292...
    15.6977   16.9662   16.5433   16.9662]*1e-6;

% Theoretical predictions from CPEO theory
y_T2 = 1.0e-04*[0.112095428775240   0.122788955932148   0.131082688124535...
    0.140405064266081   0.149558699139459   0.157759450353773...
    0.164326149920006   0.168806198005342   0.171983243812994...
    0.173706248799021];

% Representation
figure
hold on
box on
ax=gca;
ax.FontSize = 20;
ax.FontName = 'Times New Roman';
ax.XAxis.Exponent = 4;
ax.LineWidth = 1;
colorR = [212 56 37]/255;
colorO = [239 135 51]/255;
colorY = [248 206 70]/255;
colorG = [70 156 118]/255;
colorB = [56 117 176]/255;

h3=plot(f,y_T2*1e6,'-k','Linewidth',1.25,'color','k');
h4 = plot(f, y_Exp*1e6,'^r','MarkerSize',10,'Linewidth',1.25,'Color',colorR,...
    'MarkerFaceColor',colorR);
legend([h3 h4],'Predictions','3 um Carboxy','Location','SE')

xlabel('Frequency (Hz)')
ylabel('Wall Separation (\mum)')

set(gca,'XScale','log')
axis([40,1.1e4,0,20])

%% Figure 4(b)
% Frequencies used
f = [10e3 5e3 3e3 1.7e3 950 527 292 167 90 50];
% Experimental separation measured
y_Exp = [15.2748   16.1205   15.6977   14.8520   15.2748   16.9662...
    19.0803   19.9260   20.3488   20.3488]*1e-6;

% Theoretical predictions from CPEO theory
y_T2 = 1.0e-04*[0.145405550781064   0.158476118876378   0.168371410252333...
    0.179199524202033   0.189471418702061   0.198309724207252...
    0.205095189524324   0.209557502164483   0.212635213530889...
    0.214274868522463];

% Representation
figure
hold on
box on
ax=gca;
ax.FontSize = 20;
ax.FontName = 'Times New Roman';
ax.XAxis.Exponent = 4;
ax.LineWidth = 1;

h3=plot(f,y_T2*1e6,'-k','Linewidth',1.5,'color','k');
h4 = plot(f, y_Exp*1e6,'^r','Linewidth',1.25,'Color',colorR,...
    'MarkerFaceColor',colorR,'MarkerSize',10);
legend([h3 h4],'Predictions','3 um Carboxy','Location','SE')

xlabel('Frequency (Hz)')
ylabel('Wall Separation (\mum)')

set(gca,'XScale','log')
axis([40,1.1e4,0,25])

%% Figure 5(a)
% Amplitudes used
E0 = [12.571 25.143 37.714 50.286 62.857 75.429 88.000 100.57...
    113.14]*1e3; % V/m

% Experimental Data 90Hz
% 1 um
y_1um = [3.4355    5.1268    6.3953    7.2410    8.0867    8.9323...
    9.3552   10.2008   11.0465]*1e-6;
% 2 um
y_2um = [6.1839    8.9323   10.2008   11.8922   13.1607   14.4292...
    15.6977   17.3890   18.6575]*1e-6;
% 3 um
y_3um = [ 7.2410    9.7780   11.8922   13.5835   14.4292   15.6977...
    16.5433   17.8118   17.8118]*1e-6;
% 3 um PLAIN
y_3umP = [6.8182    8.9323   10.6237   11.8922   13.1607   14.0063...
    14.8520   15.2748   15.6977]*1e-6;

% Experimental Data 527Hz
% 1 um
y527_1um = [3.8584    4.2812    5.1268    6.3953    7.2410    8.0867...
    8.9323    9.7780   10.6237]*1e-6;
% 2 um
y527_2um = [4.7040    7.2410    8.9323   10.2008   11.8922   12.4695...
    13.5835   14.0063   14.4292]*1e-6;
% 3 um
y527_3um = [7.2410    8.5095   10.2008   11.0465   13.5835   14.4292...
    14.0063   13.5835   14.4292]*1e-6;
% 3 um PLAIN
y527_3umP = [5.1268    6.8182    8.5095   10.2008   11.4693   12.3150...
    12.7378   13.5835   13.5835]*1e-6;

% Figures
% 90Hz
figure
hold on
box on
ax=gca;
ax.FontSize = 20;
ax.FontName = 'Times New Roman';
ax.XAxis.Exponent = 4;

% 1umC
h1 = plot(E0, y_1um*1e6,'or','Linewidth',1.25,'Color',colorB,...
    'MarkerFaceColor',colorB,'MarkerSize',10);
% 2umC
h2 = plot(E0, y_2um*1e6,'sr','Linewidth',1.25,'Color',colorG,...
    'MarkerFaceColor',colorG,'MarkerSize',10);
% 3umC
h3 = plot(E0, y_3um*1e6,'^r','Linewidth',1.25,'Color',colorR,...
    'MarkerFaceColor',colorR,'MarkerSize',10);
% 3umP
h3P = plot(E0, y_3umP*1e6,'^k','Linewidth',1.25,'Color','k',...
    'MarkerFaceColor','k','MarkerSize',10);

xlabel('Electric Field Magnitude (V/m)')
ylabel('Wall Separation (\mum)')
legend([h1 h2 h3 h3P],'1 \mum Carboxylate', '2 \mum Carboxylate',...
    '3 \mum Carboxylate','3 \mum Plain')
set(gca,'XScale','log','YScale','log')

% 527Hz
figure
hold on
box on
ax=gca;
ax.FontSize = 20;
ax.FontName = 'Times New Roman';
ax.XAxis.Exponent = 4;

% 1umC
h1 = plot(E0, y527_1um*1e6,'or','Linewidth',1.25,'Color',colorB,...
    'MarkerFaceColor',colorB,'MarkerSize',10);
% 2umC
h2 = plot(E0, y527_2um*1e6,'sr','Linewidth',1.25,'Color',colorG,...
    'MarkerFaceColor',colorG,'MarkerSize',10);
% 3umC
h3 = plot(E0, y527_3um*1e6,'^r','Linewidth',1.25,'Color',colorR,...
    'MarkerFaceColor',colorR,'MarkerSize',10);
% 3umP
h3P = plot(E0, y527_3umP*1e6,'^k','Linewidth',1.25,'Color','k',...
    'MarkerFaceColor','k','MarkerSize',10);

xlabel('Electric Field Magnitude (V/m)')
ylabel('Wall Separation (\mum)')
legend([h1 h2 h3 h3P],'1 \mum Carboxylate', '2 \mum Carboxylate',...
    '3 \mum Carboxylate','3 \mum Plain')
set(gca,'XScale','log','YScale','log')