wave friction using Soulsby page 80
0001 % wave friction using Soulsby page 80 0002 0003 function [tauw,fw] = wave_fric(u_orb,t_orb,z0,nu,rho) 0004 0005 0006 ks = z0*30.; %Nikuradse roughness 0007 0008 A = u_orb*t_orb/(2.*pi); %semi-orbital excursion 0009 r = A/ks; %relative roughness 0010 Rw = u_orb*A/nu; %wave Reynolds 0011 0012 fwr = 1.39*(A/z0)^(-.52); %rough turbulent flow friction fac 0013 if(Rw > 5e5) 0014 B = .0521; N = .187; %smooth turbulent 0015 else 0016 B = 2; N = 0.5; % laminar 0017 end; 0018 fws = B*Rw^(-N);%smooth bed friction factor 0019 0020 fw = max(fwr,fws); %wave friction factor 0021 0022 tauw = 0.5*rho*fw*u_orb*u_orb; %wave shear stress