Contaminant and sediment transport by advection and diffusion 241
Computer code 8.8
% Example 8.8 Cross Shore Sediment Transport and Beach
Profile Formation
% H = Maximum water height at the left-hand boundary [m];
% ao = Wave height at the left-hand boundary [m];
% T = Wave period [s];
% S = Bed slope;
% dsed = Sediment density [kg/m^3];
% dwtr = Water density [kg/m^3];
% d50 = Mean particle diameter [m];
% ks = Absolute bed roughness [m];
% h() = Water depth [m];
% Dx = Spatial step [m];
% Dt = Time step [s];
% nt = Number of time steps;
clc; clear all; close all;
% Input data;
g=9.81;
H=5;
ao=2;
T=6;
S=0.02;
ks=0.001;
dsed=2.5;
dwtr=1.0;
d50=0.001;
Dx=3;
Dt=12;
nt=1500;
lo=g*T^2/(2*pi);
sed=g^2.5*((dsed/dwtr)-1)^2*d50;
% Definition of initial depths and solution domain size;
kx=1000;
for k=2:kx
h(k)=H-S*(k-1)*Dx;
if h(k)>0
nxx=k;
else
nx=nxx-1;
break
end
end
% Breaking water depth;
hb=ao/0.8;
% Estimation of water setup and modified depths;
xsi=S/sqrt(ao/lo);
G=xsi^0.17+0.08;
set=hb/(1+8/3/G^2);
% Estimation of the breaking line location;