Contaminant and sediment transport by advection and diffusion 229
parameters from actual field data. The ability of a good qualitative
description but weak quantitative one is characteristic of sediment
transport models.
Computer code 8.6
% Example 8.6 Sediment Transport with Engelund-Hansen Method
% d = Particle diameter [m];
% fb = Bed friction coefficient;
% dr = Relative density difference;
% dwtr = Water density [kg/L];
% dsed = Sediment density [kg/L];
% Dd = Spatial step (same in both directions)[m];
% Dt = Time step [s];
% nx = Number of spatial steps in the x-direction;
% ny = Number of spatial steps in the y-direction;
% nt = Time of simulation steps [s];
clc; clear all; close all;
% Input data;
g=9.81;
d=0.001;
fb=0.001;
dwtr=1.0;
dsed=2.65;
dr=(dsed-dwtr)/dwtr;
Dd=5;
Dt=1;
nx=70;
ny=40;
nt=10368000;
up2=dr*g*d;
qss=d*sqrt(up2);
% Importing data for depths, and velocities;
load Depths.txt
load CoastU.txt
load CoastV.txt
fDep=Depths;
fCoU=CoastU;
fCoV=CoastV;
% Subroutine for depths h=fm(i,j,fDep);
% Subroutine for velocity u=fm(i,j,fCoU);
% Subroutine for velocity v=fm(i,j,fCoV);
% Initialization of variables;
for i=1:nx
for j=1:ny
qsx(i,j)=0;
qsy(i,j)=0;
dh(i,j)=0;
end