Week 12 - Validation studies of Symmetry BC vs Wedge BC in OpenFOAM vs Analytical H.P equation

1. Validation of Hydro-dynamic length with the numerical result / Validate maximum velocity and pressured drop for fully developed flow / Post-processing Shear stress and validating wall shear stress for fully developed flow:

wedge* angle=10 :

symmetry* angle=10 :

wedge* angle=25 :

symmetry* angle=25 :

wedge* angle=45 :

symmetry* angle=45 :

as we notice there is no difference in pressure drop, maximum velocity,shear stress between all the cases so every one of them is validated.

There is the code :

clear all 
close all 
clc
%inputs
R=0.01/2;
Vavg=0.21;



%%%%% a=10
sym_fileID='sym_10 (1)_txt.txt';
wdg_fileID='wedge_10_txt.txt';

[nm_line1_s_10,nm_line2_s_10,nm_line1_w_10,nm_line2_w_10]=f_strnum(sym_fileID,wdg_fileID);


%%%%% a=25
sym_fileID='sym_25_txt.txt';
wdg_fileID='wdg_25_txt.txt';

[nm_line1_s_25,nm_line2_s_25,nm_line1_w_25,nm_line2_w_25]=f_strnum(sym_fileID,wdg_fileID);

%%%%% a=45
sym_fileID='sym_45_txt.txt';
wdg_fileID='wedge_45_txt.txt';

[nm_line1_s_45,nm_line2_s_45,nm_line1_w_45,nm_line2_w_45]=f_strnum(sym_fileID,wdg_fileID);


hold on
plot(nm_line1_s_10,nm_line2_s_10)
plot(nm_line1_w_10,nm_line2_w_10)

plot(nm_line1_s_25,nm_line2_s_25)
plot(nm_line1_w_25,nm_line2_w_25)

plot(nm_line1_s_45,nm_line2_s_45)
plot(nm_line1_w_45,nm_line2_w_45)

legend('sym-10','wdg-10','sym-25','wdg-25','sym-45','wdg-45')




% %%%%%%% sym
% f1=fopen('sym_10_txt.txt','r'); %read the file
% line{1}=fgetl(f1);
% line{2}=fgetl(f1);%copy the first line into a dictionarry called line
% 
%  
% for i=2:1002  
%          
% 
%     num_line=str2num(line{i});
%     nm_line1=num_line(1);
%     nm_line2(i)=num_line(2);
%     v(i)=2*Vavg*(1-nm_line1^2/R^2);
%     i=i+1; %to point into the reference of the next line 
%     line{i}=fgetl(f1);% copy the next line
%     
%     
% end
% line=line(1:end-1);%delete the last element cause it includes[-1](the stop condition for the while loop)
% fclose(f1)
% v1=v;
% nm_line2_10_s=nm_line2;
% 
% %%%%%%%%%%%%%%%Wdg
% 
% f1=fopen('wedge_10_txt.txt','r'); %read the file
% line{1}=fgetl(f1);
% line{2}=fgetl(f1);%copy the first line into a dictionarry called line
% 
%  
% for i=2:1002  
%          
% 
%     num_line=str2num(line{i});
%     nm_line1=num_line(1);
%     nm_line2(i)=num_line(2);
%     v(i)=2*Vavg*(1-nm_line1^2/R^2);
%     i=i+1; %to point into the reference of the next line 
%     line{i}=fgetl(f1);% copy the next line
%     
%     
% end
% line=line(1:end-1);%delete the last element cause it includes[-1](the stop condition for the while loop)
% fclose(f1)
% nm_line2_10_w=nm_line2;
% 
% 
% %%%%%%% sym
% f1=fopen('sym_25_txt.txt','r'); %read the file
% line{1}=fgetl(f1);
% line{2}=fgetl(f1);%copy the first line into a dictionarry called line
% 
%  
% for i=2:1002  
%          
% 
%     num_line=str2num(line{i});
%     nm_line1=num_line(1);
%     nm_line2(i)=num_line(2);
%     v(i)=2*Vavg*(1-nm_line1^2/R^2);
%     i=i+1; %to point into the reference of the next line 
%     line{i}=fgetl(f1);% copy the next line
%     
%     
% end
% line=line(1:end-1);%delete the last element cause it includes[-1](the stop condition for the while loop)
% fclose(f1)
% v1=v;
% nm_line2_25_s=nm_line2;
% 
% %%%%%%%%%%%%%%%Wdg
% 
% f1=fopen('wdg_25_txt.txt','r'); %read the file
% line{1}=fgetl(f1);
% line{2}=fgetl(f1);%copy the first line into a dictionarry called line
% 
%  
% for i=2:1002  
%          
% 
%     num_line=str2num(line{i});
%     nm_line1=num_line(1);
%     nm_line2(i)=num_line(2);
%     v(i)=2*Vavg*(1-nm_line1^2/R^2);
%     i=i+1; %to point into the reference of the next line 
%     line{i}=fgetl(f1);% copy the next line
%     
%     
% end
% line=line(1:end-1);%delete the last element cause it includes[-1](the stop condition for the while loop)
% fclose(f1)
% nm_line2_25_w=nm_line2;
% 
% %%%%%%%%%%%%%%%%%%Sym
% f1=fopen('sym_45_txt.txt','r'); %read the file
% line{1}=fgetl(f1);
% line{2}=fgetl(f1);%copy the first line into a dictionarry called line
% 
%  
% for i=2:1002  
%          
% 
%     num_line=str2num(line{i});
%     nm_line1=num_line(1);
%     nm_line2(i)=num_line(2);
%     v(i)=2*Vavg*(1-nm_line1^2/R^2);
%     i=i+1; %to point into the reference of the next line 
%     line{i}=fgetl(f1);% copy the next line
%     
%     
% end
% line=line(1:end-1);%delete the last element cause it includes[-1](the stop condition for the while loop)
% fclose(f1)
% v1=v;
% nm_line2_2=nm_line2;
% 
% %%%%%%%%%%%%%%%Wdg
% 
% f1=fopen('wedge_45_txt.txt','r'); %read the file
% line{1}=fgetl(f1);
% line{2}=fgetl(f1);%copy the first line into a dictionarry called line
% 
%  
% for i=2:1002  
%          
% 
%     num_line=str2num(line{i});
%     nm_line1=num_line(1);
%     nm_line2(i)=num_line(2);
%     v(i)=2*Vavg*(1-nm_line1^2/R^2);
%     i=i+1; %to point into the reference of the next line 
%     line{i}=fgetl(f1);% copy the next line
%     
%     
% end
% line=line(1:end-1);%delete the last element cause it includes[-1](the stop condition for the while loop)
% fclose(f1)
% 
% 
% 
% %analytical velocity along y axis
% %plot(v)
% hold on
% %plot(v1)
% %numerical velocity
% plot(nm_line2)
% plot(nm_line2_2)
% 
% legend('wedge','symmetry')
% title('alpha=45')
% 
% 
% 
% 

There is difference just due to the angle which gives a maximum radium more & more shorter

So there is no difference between thes two kinds of boundary