Contaminant and sediment transport by advection and diffusion 249
Computer code 8.9
% Example 8.9 Air Bubbles Induced Flow
% qg = Air discharge [m^3/s];
% wb = Bubbles/particle rising velocity [m/s];
% cp = Pseudo compressibility coefficient [m/s];
% xo, yo = Particles origin coordinates [m];
% tau = Drag force divided by the density per particle [m^3/
s^2];
% ed = Eddy diffusion coefficient [m^2/s];
% ip = Total number of particles;
% Dd = Spatial step same in both directions [m];
% Dt = Time step [s];
% nx = Number of steps in the x-direction;
% ny = Number of steps in the y-direction;
% nt = Number of time steps;
clc; clear all; close all;
% Input data;
qg=0.001;
wb=0.35;
cp=4;
xo=14.5;
yo=1.5;
tau=0.000755;
ed=0.1;
ip=8000;
Dd=1;
30
25
20
15
10
5
0 0
5
10
15
20
25
30
Horizontal distance (m)
Vertical distance (m)
Figure 8.15 Quasi-steady-state velocities field.